#Distorted rocksalt LaAs to test proper handling of nonlocal f psp

  ndtset  2

#Dataset 1 : self-consistent ground state calculation

#Dataset 2 : response function strain perturbation calculation
  rfstrs2  3
   rfdir2  1 0 0
    nqpt2  1
     qpt2  0 0 0
  kptopt2  2
  getwfk2  1

# common input data

   acell  3*11.81
  diemac  10.0
    ecut  8.0
  ecutsm  0.0
     ixc  2
  kptopt  1
  mqgrid  1501
   natom  2
   nband  4
   ngkpt  2 2 2
nloc_alg  2
 nshiftk  1
   nstep  50
    nsym  1
  ntypat  2
  prtvol  10
   rprim -0.007  0.548  0.473
          0.468  0.038  0.512
          0.481  0.514 -0.032
  shiftk  0.5 0.5 0.5
  tolvrs  1.0d-12
   typat  1 2
    xred  0.000  0.000  0.000
          0.543  0.493  0.508
   znucl  57  33

 pp_dirpath "$ABI_PSPDIR"
 pseudos "57la.drh, 33as.drh"

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test = 
#%%   t62.out, tolnlines = 0, tolabs = 0.000e+00, tolrel = 0.000e+00, fld_options = -easy
#%% [paral_info]
#%% max_nprocs = 4
#%% [extra_info]
#%% authors = D. R. Hamann
#%% keywords = NC, DFPT
#%% description = 
#%%   Test of the strain perturbation.
#%%   LaAs in a randomly distorted version of its normal rocksalt structure
#%%   using a non-standard norm conserving potential produced by DRH
#%%   for this test.  Distorted structures were used for most development
#%%   work so that no components of the strain 2nd derivatives would be
#%%   zero by symmetry.  A La compound was chosen because its f nonlocal
#%%   potential makes a significant contribution to the result (10%
#%%   for a well-converged result although only 1% for the low value
#%%   of ecut used here), so it tests that the very complex routines
#%%   for the strain derivatives of the nonlocal energy are behaving
#%%   properly.  (Choices in constructing this La potential were
#%%   restricted by current limitations of the treatment of fhi-
#%%   format potentials, and it should not be used for physics results.)
#%%<END TEST_INFO>