资源描述：
《求一个用于稀疏矩阵求解的iccg的fortran程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、求一个用于稀疏矩阵求解的ICCG的fortran程序implicitreal(a-h,o-z)reala(200000),af(200000),b(28000),x(28000),$anrm,xnrm,resd,g(28000),p(28000),t1(28000),t2(28000)integeria(28001),ja(200000),nit(3)reald,c2,a2,cfcommon/tomm/a,af,ia,jacommon/area2/d,c2,a2,cfcommon/order/iorflgdatanit/1,5,10
2、/ccthisisamainprogramforrunningtheincompleteLUcconjugategradientforsolvinglinearsystems.cthematrixneednotbesymmetric.ccilucgthisflagif1willcauseincompletefactoizationcif0nofactorizationjuststraightcgccwrite(6,*)'startinmain'cccalltrapov(100)loca=2000001continuecwrite(6
3、,*)'inputineum5,ineum6,irot,iorflg'cwrite(6,*)'setirot<0tostop'cread(5,*)ineum5,ineum6,irot,iorflgcif(irot.lt.0)thencwrite(6,*)'endtest'cstopcendifcwrite(6,*)'ineum5,ineum6,irot,iorflg',ineum5,ineum6,irot,iorflgcwrite(6,*)'inputmeshcellsforx,y,andz'cread(5,*)imx,imy,im
4、zcwrite(6,*)'inner,middle,andouterfill-inintegers'cread(5,*)ifli,iflm,ifloccwrite(6,*)'ising(normallyzero)'ccread(5,*)isingineum5=1ineum6=1irot=1iorflg=2imx=20imy=20imz=20ccif(ineum5.ne.1.or.ineum6.ne.1)ising=0ising=0ccimx=5ccimy=5ccimz=5n=imx*imy*imzcwrite(6,*)'******
5、*****problemtype*************'cif(ineum5.eq.0)thencwrite(6,*)'*****dironbottom*****'celsecwrite(6,*)'*****neumanonbottom*****'cendifcif(ineum6.eq.0)thencwrite(6,*)'*****dirontop*****'celsecwrite(6,*)'*****neumanontop*****'cendifcif(ising.eq.0)thencwrite(6,*)'*****regul
6、ar*****'celsecwrite(6,*)'*****singularifneumanontopandbottom*****'cendifcwrite(6,*)'*****meshcellsinxdir',imx,'*****'cwrite(6,*)'*****meshcellsinydir',imy,'*****'cwrite(6,*)'*****meshcellsinzdir',imz,'*****'ccifliisthenumberoffillinstripesallowedforeachofthecinnerstrip
7、es;ifzeronofillin.ciflmisthenumberoffillinstripesallowedforeachofthecmiddlestripes;ifzeronofillin.cifloisthenumberoffillinstripesallowedforeachofthecouterstripes;ifzeronofillinhere.cinum5if1thenneumannonbottomcif0thendirchletonbottomcinum6if1thenneumannontopcif0thendir
8、chletontopcirotif0thennorotflowfieldcif1thenrotflowfieldcisingifzeroregularcif1thensignularifneumanntopandbottomcimx#
