# single H atom, test of fftalg and istwfk. Use a non-cubic unit cell.

#ndtset 2   jdtset 1 13   getwfk 1
 ndtset 15  getwfk 1

fftalg1  112    istwfk1  2   getwfk1 0
fftalg2  112    istwfk2  1
fftalg3  111    istwfk3  1
fftalg4  111    istwfk4  2
fftalg5  110    istwfk5  2
fftalg6  110    istwfk6  1
fftalg7  100    istwfk7  1
fftalg8  112    istwfk8  2   nstep8  0
fftalg9  112    istwfk9  1   nstep9  0
fftalg10 400    istwfk10 1
fftalg11 401    istwfk11 1
fftalg12 402    istwfk12 1
fftalg13 410    istwfk14 2
fftalg14 411    istwfk15 2
fftalg15 412    istwfk16 2

acell 3*10
diemac 1.0d0
diemix 1.0d0
ecut 5
intxc 1

 kptopt 0

kpt 3*0
natom 1
nband 1
nkpt 1
nline 3
nstep 20
nsym 1
ntypat 1
rprim 1 1 0
      0 1 0
      0 1 1
tolwfr 1.0d-14
typat 1
xred 3*0
wtk 1 znucl 1

 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosTM_pwteter/1h.pspnc"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t88.out, tolnlines = 0, tolabs = 0.000e+00, tolrel = 0.000e+00
#%% [paral_info]
#%% max_nprocs = 1
#%% [extra_info]
#%% keywords = NC
#%% authors =  Unknown
#%% description = 
#%%   H atom, similar to test 87, except that the box is not cubic,
#%%   but only parallelipipedic (not even a rectangle parallelipiped).
#%%   However, the primitive vectors describe the same lattice as
#%%   in case 87, so that the results must be identical
#%%   to those of test 87, even if a different FFT grid is used.
#%%   That is what is observed.
#%%<END TEST_INFO>