30 double rel_photon_energy;
42 rel_photon_energy = EgammaRyd /
EthRyd;
43 rel_photon_energy =
MAX2( rel_photon_energy , 1. + FLT_EPSILON*2. );
45 cs =
H_photo_cs(rel_photon_energy ,
N_(ipLev),
L_(ipLev), nelem + 1 );
47 ASSERT( cs > 0. && cs < 1.E-8 );
55 double alpha,RecomIntegral=0.,b,E1,E2,step,OldRecomIntegral,TotChangeLastFive;
56 double change[5] = {0.,0.,0.,0.,0.};
83 OldRecomIntegral = RecomIntegral;
88 change[4] = change[3];
89 change[3] = change[2];
90 change[2] = change[1];
91 change[1] = change[0];
92 change[0] = (RecomIntegral - OldRecomIntegral)/RecomIntegral;
93 TotChangeLastFive = change[0] + change[1] + change[2] + change[3] + change[4];
98 }
while ( ((E2-
EthRyd) < 100.*
kTRyd) && ( TotChangeLastFive > 0.0001) );
101 alpha = b * RecomIntegral;
114 x1 = ERyd * ERyd * exp(-1.0 * ( ERyd -
EthRyd ) /
kTRyd);