1) Run DMFT
2) Obtain matrix elements:

   Prepare "case.inop" file. Check Wien2K manual.
   
   in case.inop change MME to ON
   > x optic[c] [-so] [-up]

   You should get "case.mommat" and "case.symop"
   
3) run DMFT in mode "u":
   a) execute
   > x_dmft.py dmftu
   You should get two files: Udmft[updn].0 and BasicArrays.dat

4) prepare input file named "dmftopt.in", which contains:

    0          # updn: [0|up|dn] -- 0 for non-magnetic calculation
    sr2ruo4    # case
    0.01       # gamma  -- broadening of all bands / we recommend to keep nonzero value
    0.0        # gammac -- broadening of the correlated bands (in addition to the self-energy)
    4          # ommax  -- maximum frequency for the optics calculation
    1e-2       # delta  -- minimum separation of frequency for logarithmic mesh in frequency
    5          # Nd     -- number of points in each linear mesh, a subset of logarithmic mesh
    0          # Qsym: [0|1]. Do we need to symmetrize? (over all irreducible only or over all.)
    0          # InterbandOnly [0|1] (0 -- all, 1--interband)
    0.5        # dwindow -- Not all bands are used in computation of optics, but only those from [-omeg-dwindow, omega+dwindow].
    2          # Ndirection -- How many optics type do you want to compute: xx, yy, zz,...
    0.5 0.0 0.0    # alphaV(1,:)
    0.0 0.5 0.0    # alphaV(2,:)
    0.0 0.0 0.0    # alphaV(3,:)
    0.0 0.0 0.0    # alphaV(1,:)
    0.0 0.0 0.0    # alphaV(2,:)
    0.0 0.0 1.0    # alphaV(3,:)


   prepare self-energy sig.inp1 by
   > ssplit.py

   and finally run DMFT optics by
   > dmftopt
   
   You will need "case.energy" , "case.mommat", "case.symop" ,
   "Udmft.0" and "BasicArrays.dat".

   Result is stored in "optics.dat". The corresponding total density
   of states is in "optdos.dat".
