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.);
71 Aul *= (2.*
L_(ipLo)+3.) * 2. / (2.*(double)
N_(ipHi)*(double)
N_(ipHi));
81 Aul += Aul1*(2.*
L_(ipLo)-1.) * 2. / (2.*(double)
N_(ipHi)*(double)
N_(ipHi));
92 if(
N_(ipHi) ==
N_(ipLo) )
99 else if( ipLo == 0 && ipHi == 1 )
102 Aul = 8.226*pow((
double)(nelem+1.),6.);
105 else if( ipLo == 0 && ipHi == 2 )
110 else if( abs(
L_(ipLo) -
L_(ipHi) )== 1 )
119 ( abs(
L_(ipHi)-
L_(ipLo)) > 1) );