15 long int lower, iupper;
29 lower =
MIN2( n1 , n2 );
30 iupper =
MAX2( n1, n2 );
31 if( lower < 1 || lower == iupper )
33 fprintf(
ioQQQ,
" HydroEinstA called with impossible ns, =%li %li\n", lower, iupper);
50 double z4 =
POW4((
double)nelem+1.);
53 if( ipHi >=
iso_sp[ipISO][nelem].numLevels_max-
iso_sp[ipISO][nelem].nCollapsed_max )
55 if( ipLo >=
iso_sp[ipISO][nelem].numLevels_max-
iso_sp[ipISO][nelem].nCollapsed_max )
69 iso_sp[ipISO][nelem].CachedAs[
N_(ipHi)-
iso_sp[ipISO][nelem].n_HighestResolved_max-1 ][ ipLo ][0] = (
realnum)Aul;
71 Aul *= (2.*
L_(ipLo)+3.) * 2. / (2.*(
double)
N_(ipHi)*(
double)
N_(ipHi));
78 iso_sp[ipISO][nelem].CachedAs[
N_(ipHi)-
iso_sp[ipISO][nelem].n_HighestResolved_max-1 ][ ipLo ][1] = (
realnum)Aul1;
81 Aul += Aul1*(2.*
L_(ipLo)-1.) * 2. / (2.*(
double)
N_(ipHi)*(
double)
N_(ipHi));
84 iso_sp[ipISO][nelem].CachedAs[
N_(ipHi)-
iso_sp[ipISO][nelem].n_HighestResolved_max-1 ][ ipLo ][1] = 0.f;
92 if(
N_(ipHi) ==
N_(ipLo) )
99 else if( ipLo == 0 && ipHi == 1 )
102 Aul = 2.46e-6*pow((
double)(nelem+1.),10.);
105 else if( ipLo == 0 && ipHi == 2 )
110 else if( abs(
L_(ipLo) -
L_(ipHi) )== 1 )
119 ( abs(
L_(ipHi)-
L_(ipLo)) > 1) );
#define DEBUG_ENTRY(funcname)
double H_Einstein_A(long int n, long int l, long int np, long int lp, long int iz)
double HydroEinstA(long int n1, long int n2)
realnum hydro_transprob(long nelem, long ipHi, long ipLo)
double HydroOscilStr(double xLower, double Upper)
t_iso_sp iso_sp[NISO][LIMELM]
void iso_put_error(long ipISO, long nelem, long ipHi, long ipLo, long whichData, realnum errorOpt, realnum errorPess)
UNUSED const double TRANS_PROB_CONST
UNUSED const double RYD_INF