Function MproductbDen

Applications/HHL_Algorithm/BasicFunctionByOrigin.cpp:267–287

}

VectorXd MproductbDen(const SparseQMatrix<double>& M, const VectorXd &b) {
	int size = M.size;
	VectorXd Mb(size);
	auto data = M.data;
	auto col_ptr = M.nnzcol_ptr;
	auto col_ind = M.nnzcol_ind;
	if (size == 0) { 
		cout << "wrong is here" << endl; 
	}
#pragma omp parallel for
	for (int i = 0; i < size; i++) {
		double sum = 0;
		for (int j = col_ptr[i]; j < col_ptr[i + 1]; j++) {
			auto temp = data[j] * b(col_ind[j]);
			sum += temp;
		}	
		Mb(i) = sum;
	}

	return Mb;
}

vector<double> MproductbDen(const SparseQMatrix<double>& M, const vector<double>& b) {
	int size = M.size;