1
2
3
4
5
6
7
8
9
10
11
|
program NF
end program NF
subroutine mattest(nx,ny,nz,band1,band2,band3,stiffness,maxiter,targrms,method)
integer,parameter :: dpkind=kind(1.0D0)
character(*) :: method
real(dpkind),allocatable,dimension(:) :: ad,au1,au2,au3,x,b
allocate(ad(nxyz),au1(nxyz),au2(nxyz),au3(nxyz),x(nxyz),b(nxyz))
au1(nx:nxyz:nx) = 0.0
if ( method=='NFCG' ) then
endif
end subroutine mattest
|