# ================================================================
# Fe bcc structure (ferromagnetic metal)
#
# Test for Abinit: PAW spin-polarized + GGA
#                  Shape function for compensation charge Bessel
#                                 analytical and numerical
# Low cut-off for test
# ================================================================

# Datasets
   ndtset 2
   typat1 1  ! PAW dataset with analytical shape function
   typat2 2  ! PAW dataset with numerical shape function
  ntypat  2

# Unit cell
  acell   3*5.42
  rprim   -0.5  0.5  0.5
           0.5 -0.5  0.5
           0.5  0.5 -0.5
  natom   1
   xred   0.0  0.0  0.0
  znucl   2*26
  nband   14
  nsppol  2
  spinat  0. 0. 4.

# K-points and symmmetries
 kptopt   1
  ngkpt   6 6 6
nshiftk   1
 shiftk   0.5 0.5 0.5
 occopt   7
 tsmear   0.008

# Self-consistent cycle
 toldfe   1.0d-10
  nstep   30

# Basis sets
   ecut   8.
pawecutdg 10.

# Miscelaneous
   prtwf1 1 getwfk2 1
   prtwf2 0
   prtden 0
   prteig 0


 pp_dirpath "$ABI_PSPDIR"
 pseudos "26fe.paw, 26fe.num.paw"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t08.out, tolnlines = 40, tolabs = 1.010e-05, tolrel = 7.322e-02, fld_options = -medium
#%% [paral_info]
#%% max_nprocs = 14
#%% [extra_info]
#%% authors = Unknown
#%% keywords = PAW
#%% description = 
#%%   Fe in BCC structure. Low cut-offs and few k-points.
#%%   This test is done within GGA and (ferro-)magnetism in PAW formalism.
#%%   The same job is done twice with two differents PAW datasets.
#%%   First PAW dataset is Fe with a "Bessel" compensation charge
#%%   shape function expressed analytically.
#%%   SecondPAW dataset is Fe with a "Bessel" compensation charge
#%%   shape function expressed numerically.
#%% topics = PAW
#%%<END TEST_INFO>