# Silicon diatomic molecule for test calculations
# Exp bond=4.244 bohr; freq=511 cm^-1 (?)
# Binggeli reports that large box (18 on side) needed to get physical answers.
# Allow relaxation.
 acell 12 8 8
 diemac 1.0d0
 diemix 1.0d0
 ecut 7.0
 enunit 2
 intxc 1
 irdwfk 1
 iscf 1
 kptopt 0

 kpt   3*0
 natom  2
 nband 5 3
 nkpt 1
 nline 3
 nsppol 2
 nstep 15
 nsym 1
 ntime 10
 ntypat  1
 occ 8*1
 occopt 2
 prtvol 10
 rprim 1 0 0  0 1 0  0 0 1
 toldfe 1.0d-12
 typat  2*1
 wtk  1
 xred  -0.1846858603     -1.0293890254E-07 -1.3165386139E-08
        0.1846858603      1.0293890254E-07  1.3165386139E-08
 znucl  14.0

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

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% input_prefix = t46o 
#%% test_chain = t45.in, t46.in, t47.in
#%% [files]
#%% files_to_test = t47.out, tolnlines = 0, tolabs = 1.001e-05, tolrel = 2.000e-05, fld_options = -easy
#%% [paral_info]
#%% max_nprocs = 1
#%% [extra_info]
#%% keywords = NC
#%% authors = Unknown
#%% description = 
#%%   Si2 molecule, spin-polarized. Test the option iscf=1 (eigenvalues of the SCF cycle). 
#%%   Start from the wavefunctions of test 46.
#%%   Cannot be executed in parallel:
#%%         --- !ERROR
#%%         message: |
#%%         the number of bands in the spin up case must be equal to
#%%         the number of bands in the spin down case.
#%%         This is not the case for the k point number : 1 
#%%         The number of bands spin up and down are : 5 3 Action : change nband, or use the sequential
#%%         version of ABINIT.
#%%         src_file: chkinp.F90
#%%         src_line: 1181
#%%         ...
#%%<END TEST_INFO>