First of all I deal with plural matrix is plural square, about thirty thousand order
I have when I was in call pardiso china-africa zero in front of the plural matrix elements into CSR storage format, set the mtype=13, feel all aspects are ready, but in the end at run time always said access conflict? I also try to iparm (60) is set to 2, but no effect, and made more than a day or not as a result, since the this, every brother please see what went wrong? My other iparm setting is the default
CodePudding user response:
Or for processing complex asymmetric sparse matrix in what way? I before using conjugate gradient method, but I want to know is there any other methods? I tried Zgetrf + Zgetri, other preconditioning conjugate gradient, but is not does not conform to the required level is not enough, don't change the program, ask everybody to help,,,CodePudding user response:
Asymmetric can try double conjugate gradient or the stability of the double conjugate gradient gmres is also possible