mirror of
https://github.com/paboyle/Grid.git
synced 2024-11-09 23:45:36 +00:00
36 lines
823 B
Bash
36 lines
823 B
Bash
|
#!/bin/bash
|
||
|
|
fn=$1
|
||
|
|
|
||
|
|
grep "double zmobius_" $fn |
|
||
|
|
awk 'BEGIN{ m["zmobius_b_coeff"]=0; m["zmobius_c_coeff"]=1; }{ val[m[substr($2,0,15)]][substr($2,17)+0]=$4; }END{
|
||
|
|
|
||
|
|
ls=length(val[0])/2;
|
||
|
|
|
||
|
|
print "ls = " ls
|
||
|
|
|
||
|
|
bmc=-111;
|
||
|
|
|
||
|
|
for (s=0;s<ls;s++) {
|
||
|
|
br[s] = val[0][2*s + 0];
|
||
|
|
bi[s] = val[0][2*s + 1];
|
||
|
|
cr[s] = val[1][2*s + 0];
|
||
|
|
ci[s] = val[1][2*s + 1];
|
||
|
|
|
||
|
|
t=br[s] - cr[s];
|
||
|
|
if (bmc == -111)
|
||
|
|
bmc=t;
|
||
|
|
else if (bmc != t)
|
||
|
|
print "Warning: b-c is not constant!";
|
||
|
|
|
||
|
|
omegar[s] = (-1.0 + 2.0* br[s])/(4.0*bi[s]**2.0 + (1.0 - 2.0* br[s])**2);
|
||
|
|
omegai[s] = - 2.0* bi[s]/(4.0*bi[s]**2.0 + (1.0 - 2.0* br[s])**2);
|
||
|
|
}
|
||
|
|
|
||
|
|
print "b-c = " bmc
|
||
|
|
|
||
|
|
for (s=0;s<ls;s++) {
|
||
|
|
printf( "omega.push_back( std::complex<double>(%.15g,%.15g) );\n",omegar[s],omegai[s]);
|
||
|
|
}
|
||
|
|
|
||
|
|
}'
|