求一个用于稀疏矩阵求解的ICCG的fortran程序.docVIP

求一个用于稀疏矩阵求解的ICCG的fortran程序 implicit real(a-h,o-z) real a(200000),af(200000),b(28000),x(28000), $ anrm,xnrm,resd,g(28000),p(28000),t1(28000),t2(28000) integer ia(28001),ja(200000),nit(3) real d,c2,a2,cf common /tomm/ a,af,ia,ja common /area2/ d,c2,a2,cf common /order/ iorflg data nit/1,5,10/ c c this is a main program for running the incomplete LU c conjugate gradient for solving linear systems. c the matrix need not be symmetric. c c ilucg this flag if 1 will cause incomplete factoization c if 0 no factor