temporary. epinv is expandend

epinv is allowed form of the invariant dielectric tensor. See mkqg.F
called from qg4gw. Last outcome is dmlx, where
  !! epinv is expanded as
  !!   <ehat| epinv|ehat> = \sum_lm dmlx(iq0i,lm) *Y_lm(ehat)
See L328 mkqg.F around.
temporary. 1shot exchange part for gas_pw_gw_lmfh case works OK for newaniso=T.

It reproduces results for newansiso=F.
A little confusing part is in sxcf_fal2.F
This reproduce offset-Gamma method in diagonal Coulomb MPB. Temporary.
temporary. off-set gamma method is reproduced.
test case of newmode (Christoph method) works OK.
Only 1/epsilon is added to q=0.
test is just for GaAs.
fermi_ds is added for drude weight in tetwt5 (unused yet. Not tested yet).
gw_lmfh works for newaniso mode. (not gwsc yet).

Divergent part of hx0fp0 is only epsPP part.
because other term in chrpstophs' paper PRB81 125102 gives little effect 
(but I tested only for Si ans GaAs).
Look into cmmm if you like to recover.
(then it is necessary to comment out sxcf_fal2.F:L479
          iqend=nqibz !no sum for offset-Gamma points.)

Takao think second term in Eq.(36) gives little effect except
4pi/k**2/eLe term; (note that Eq.36 is slightliy wrong... need to
multiple sqrt(vcoul(igb)).
In addition, we skip third term in Eq.(40).
Thus L(omega) is epsilon without local field correction.
Temporary. MTO/IPW rotator is given in readeigen.F (it shold be applicable also to MPB).

gw_lmfh for gas and si works OK.
minor correction. cputid() is inserted. Q0P delta_scale in mkqg.
make check OK. newaniso2 works ok. si, gas, waln are tested; debug finished. 

* still need to implement EIBZ (or inversion symmetry) method.
* how to make this to QSGW scheme?
ctrlgen2.py (including setting of dimer) and so. TestHomoDimerAtom.

TestHomoDimerAtom constains systems and info to reproduce
data for diatomic molecule calculations.

ctrlgen2.py srtio3 (e.g. for ctrls.srtio3)
will give some reasonable ctrlgen.ctrl.srtio3.
R are not calculated but given by atomlist embedded in ctrlgen2.py.

ctrlgen2.py is not so perfect, especially above Kr.
Fix PNU and QNU part. atomic part(freeat.F and m_rdctrl). Add samples of dimer.

Takao think previous version was problematic for MMOM calculation for
atom. But, this should not affects to previous dimer result (FSMOM
method with MMOM=0 0 0 0).

TESTsamples/Fe2dimer and so are added.

*Reconsider MMOM mechanism. It might be better not to use it for
 atomic calculations for comparison...
*FRZWF is worth to consider especially for the case of small R.
ctrlgen2.py and atomsr and so are modified. temporary.
samples and ctrlgen2.py are modified.
indentation for *.F for lmv7 and NMCORE

SPEC_ATOM_NMCORE=1 means non-spin-polarized core density.
minor fix for m_rdctrl.F(just a messege). ctrlgen2.py RSMH>0.5 restriction
temporary. going to add eibz mechanism in hx0fp0.F (calling x0kf_v4h).
eibz mode works. But inefficient symmetrization in x0kf.F
make check OK. gw_lmfh OK for newaniso. But 1.samples not organized well. 2. EIBZdeveloping

Look into explanation for eibz4x0() in fpgw/gwsrc/switch.F.
Current implementation of EIBZ mode is still developing.
(only timereversal=.false. works; it is for development and very
 inefficient for symmetrization at the bottom of x0kf_v4h.)

So, TODO are
(1) EIBZ mode. probably it is better to start from q <---> -q symmetric q points.
(2) gwsc is not for newaniso mode.
(3) Tests for anisotropic case.
Temporary. Mainly eibzmode (see switch.F). It works(but slow).

This is a temporary version of EIBZ for gw_lmfh.
Symmetrization at the end of x0kf_v4h.F works OK (for test_si_gw_lmfh
case). But it is slow; I have to improve the mechanism on rotMPB.
(this can be given as a simple index changes with phase
multiplication).