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