# Crystalline silicon
# Same physical test as in tutorial t31.in, but here
# test the DOS using the tetrahedron method
#

 ndtset 2

#Dataset 1 : SCF calculation
 ngkpt1   2 2 2   # Use a rather loose grid for SCF calculation
 toldfe1  1.0d-10
 prtden1  1

#Dataset 2 : DOS calculation
 iscf2   -3
 getden2  1
 getwfk2  1
 prtdos2  2
 dosdeltae 0.00005  # This is also the default value for prtdos2=2,
                    # but the presence of dosdeltae is there for checking purposes
 ngkpt2   4 4 4     # This is not sufficient for a nice DOS.
                    # One should use ngkpt2 8 8 8
 tolwfr2  1.0d-16
 shiftk2 0.0 0.0 0.0 # Use a non-shifted k point grid, giving a much
         0.0 0.5 0.5 # better-looking DOS graph than the shifted one,
         0.5 0.0 0.5 # because tetrahedra with all four same energies
         0.5 0.5 0.0 # are avoided


#Definition of the unit cell
acell 3*10.18
rprim  0.0  0.5  0.5
       0.5  0.0  0.5
       0.5  0.5  0.0

#Definition of the atom types
ntypat 1
znucl 14

#Definition of the atoms
natom 2
typat 1 1
xred   0.0  0.0  0.0
       1/4  1/4  1/4

#Definition of the planewave basis set
ecut  8.0

#Definition of the k-point grid
kptopt 1
nshiftk 4
shiftk 0.5 0.5 0.5
       0.5 0.0 0.0
       0.0 0.5 0.0
       0.0 0.0 0.5

#Definition of the SCF procedure
nstep 10
diemac 12.0


 pp_dirpath "$ABI_PSPDIR"
 pseudos "PseudosTM_pwteter/14si.pspnc"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t46.out, tolnlines = 0, tolabs = 0.000e+00, tolrel = 0.000e+00;
#%%   t46o_DS2_DOS, tolnlines = 0, tolabs = 0.000e+00, tolrel = 0.000e+00
#%% [paral_info]
#%% max_nprocs = 2
#%% [extra_info]
#%% authors = Unknown
#%% keywords = 
#%% description = 
#%%   2 Silicon atoms, in a diamond structure.
#%%   Test the computation of the DOS, using the linear tetrahedron method.
#%% topics = ElecDOS
#%%<END TEST_INFO>