1) Run DMFT
2) Obtain matrix elements:

   change TOT to FERMI in case.in2c
   change TETRA with value=101 in case.in2c for metallic systems
   > run[sp]_lapw [-so] -s lapw1 -e lcore
   
   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 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
    gamma-Ce   # case
    0.01       # gamma  -- broadening of all bands
    0.01       # gammac -- broadening of the correlated bands
    10.        # ommax  -- maximum frequency
    1e-3       # delta  -- minimum domega for logarithmic mesh in frequency
    5          # Nd     -- number of points in each linear mesh, 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.33333 0.0 0.0    # alphaV(1,:)
    0.0 0.33333 0.0    # alphaV(2,:)
    0.0 0.0 0.33333    # alphaV(3,:)
    0          # Project to correlated only

    

   run DMFT optics
   > 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".
   
