Function MprodcutA

Applications/HHL_Algorithm/BasicFunctionByOrigin.cpp:246–265

}

SparseQMatrix<double> MprodcutA(SparseQMatrix<double>& M, SparseQMatrix<double>& A, int SparseIndex) {
	int size = A.size;
 
	vector<vector<pair<int, double>>> sparse_dataM(size);

#pragma omp for
	for (int i = 0; i < size; i++) {
		auto rowV = M.getSparseEigenRowVector(i);
		for (int j = 0; j < size; j++) {
			auto colV = A.getSparseEigenColVector(j);
			auto val = rowV.dot(colV);
			if (val != 0) {
				sparse_dataM[i].emplace_back(j, val);
			}
		}
	}
	SparseQMatrix<double> MA;
	FastTransfer2SparseQMat(sparse_dataM.data(), size, MA);
	return MA;
}

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