I am trying to use ZGETRF+ZGETRS function to solve a square matrix equation in form of Ax=b but the results are totally different comparign with results from GMM++. I use the code below in C++ and Linux:
// MKL solution to solve "YV=I" matrix equation
zgetrf(&nNodes, &nNodes, MKL_Y, &nNodes, ipiv, &info);
zgetrs(&transa, &nNodes, &nrhs, MKL_Y, &nNodes, ipiv, &I, &nNodes, &info);
where I is defined as
vector > I(nNodes);
and filled before above functions are called. nNodes also indicates the dimension of "A" which is square.
Have you had same experience before? Any help would be appreciated.