#test silicon crystal
 ndtset 4

#First dataset : SC run with LDA and effmass_free=1.00
 ixc1     32

#Second dataset : SC run with mGGA fake1 (alpha=1-1/effmass_free) effmass=1.01
 ixc2     33

#Third dataset : (same as dataset 1 but with spin)
 ixc3     32
 nsppol3  2
 nspden3  2     spinmagntarget3   0.0d0

#Fourth dataset : (same as dataset 2 but with spin)
 ixc4     33
 nsppol4  2
 nspden4  2     spinmagntarget4   0.0d0

#Common data
 acell   3*10.26
 amu 28.086
 ecut 25.00
 kptopt 0
 nkpt 2
 kpt 0.25000   0.25000   0.25000
     0.25000   0.50000   0.50000
 wtk 0.25 0.75
 natom  2 nband 4
 nstep 8
 ntypat  1
 occopt 1
 prtden 0
 rprim   0 .5 .5  .5 0 .5  .5 .5 0
 xred 3*0.00d0 3*0.25d0
 typat  1 1  tolwfr  1.e-28
 znucl  14
 usekden 1
 #prtkden 1
 prtwf 0

 pp_dirpath "$ABI_PSPDIR"
 pseudos "14si.Hamann_mod"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t32.abo, tolnlines = 0, tolabs = 0.0, tolrel = 0.0
#%% [paral_info]
#%% max_nprocs = 2
#%% [extra_info]
#%% authors = Unknown
#%% keywords = 
#%% description = 
#%%   Test the mGGA implementation. Silicon bulk. Two datasets are run.
#%%   One with a fake native mGGA (fake2, i.e. ixc 32) involving laplacian of the density.
#%%   One with a fake native mGGA (fake3, i.e. ixc 33) involving gradient of the density.
#%%   Both calculation should give the same total and exchange-correlation energies.
#%%   (Two other datasets are also run for spin polarized case)
#%%<END TEST_INFO>