3 subroutine mattest(nx
,ny
,nz
,band1
,band2
,band3
,stiffness
,maxiter
,targrms
,method
)
4 integer,parameter :: dpkind
=kind(1.0D0
)
6 real(dpkind
),allocatable
,dimension(:) :: ad
,au1
,au2
,au3
,x
,b
7 allocate(ad(nxyz
),au1(nxyz
),au2(nxyz
),au3(nxyz
),x(nxyz
),b(nxyz
))
9 if ( method
=='NFCG' ) then
11 end subroutine mattest