# H compressed simple cubic, metallic.
# Large smearing width

 ndtset 4

# Defaults for phonon calculations
    nqpt   1
  rfphon   1
 rfatpol   1 1
   rfdir   0 1 0


#DATASET 1 : SC run with 4 kpoints
  getwfk1  0
    nqpt1  0
  rfphon1  0
  kptopt1  1
  tolwfr1  1.0d-22

#DATASET 2 : RF at Gamma
     qpt2   0.0d0 0.0d0 0.0d0
  rfphon2   1
  toldfe2   1.0d-10

#DATASET 3 : RF at 1/2 0 0
     qpt3   0.5d0 0.0d0 0.0d0
  rfphon3   1
  toldfe3   1.0d-10

#DATASET 4 : RF at 1/4 0 0
     qpt4   0.25d0 0.0d0 0.0d0
  rfphon4   1
  toldfe4   1.0d-10

#Common data
 acell  3*2.60

 ntypat 1
 znucl 1.0
 natom 1
 typat  1
 xred 0.0 0.0 0.0

 nstep 20

 ecut 20.0

 nband  3
 nbdbuf 1
 getwfk 1

 kptopt 3
 occopt 7
 tsmear 0.2

 ngkpt  4 4 4
 shiftk
   0.5 0.5 0.5


 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosGTH_pwteter/01h.pspgth"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t90.out, tolnlines = 2, tolabs = 1.100e-09, tolrel = 9.000e-10
#%% [paral_info]
#%% max_nprocs = 4
#%% [extra_info]
#%% keywords = NC, DFPT
#%% authors = Unknown
#%% description = 
#%%   H, compressed simple cubic, to test accurately RF in the metallic case. 
#%%   No linear XC core correction.
#%%   The total energy of the 1-atom cell is -.559361014386 Ha
#%%   The 2DTE with respect to a q(0 0 0) phonon is 4.6E-10 .
#%%   The 2DTE with respect to a transverse q(1/2 0 0) phonon
#%%   is 1.75340771040435E-02. The 2DTE with respect to a transverse
#%%   q(1/4 0 0) phonon is 9.38555654208462E-03.
#%%<END TEST_INFO>