1program NF
2end program NF
3subroutine mattest(nx,ny,nz,band1,band2,band3,stiffness,maxiter,targrms,method)
4  integer,parameter :: dpkind=kind(1.0D0)
5  character(*) :: method
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))
8  au1(nx:nxyz:nx) = 0.0
9  if ( method=='NFCG' ) then
10  endif
11end subroutine mattest
12