#!/bin/bash

cp Delta.dat Delta.dat.prev

math <<EOF

sigma = Import["resigma.dat"];
sigim = Import["imsigma.dat"];
sigma[[All, 2]] += I * sigim[[All, 2]];
fsigma = Interpolation[sigma];

htDOS0[z_] := 2  (z - I Sign[Im[z]] Sqrt[1 - z^2]);
EPS = 10^-20;
htDOS[z_] := htDOS0[Re[z] + I (If[Im[z] > 0.0, Im[z], EPS])];

gf = Map[{#[[1]], htDOS[ #[[1]] - #[[2]] ]}&, sigma];

gfim = Map[{#[[1]], Im[#[[2]]]}&, gf];
gfre = Map[{#[[1]], Re[#[[2]]]}&, gf];

awim = Map[{#[[1]], -1/Pi #[[2]] }&, gfim];
awre = Map[{#[[1]], -1/Pi #[[2]] }&, gfre];

Export["imaw.dat", awim];
Export["reaw.dat", awre];

g0inv = Table[{gf[[i,1]], gf[[i,2]]^-1+sigma[[i,2]]}, {i, Length[gf]}];
Delta = Table[{g0inv[[i,1]], Im[g0inv[[i,2]]]}, {i, Length[gf]}];

Export["Delta.dat", Delta];

EOF

${NRGPREFIX} kk Delta.dat Delta-re.dat