Changeset 2034 for branches/c08_branch/source/hydroeinsta.cpp
- Timestamp:
- 05/10/08 09:03:02 (2 months ago)
- Files:
-
- branches/c08_branch/source/hydroeinsta.cpp (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
- Modified
- Copied
- Moved
branches/c08_branch/source/hydroeinsta.cpp
r1732 r2034 48 48 long ipISO = ipH_LIKE; 49 49 /* charge to 4th power, needed for scaling laws for As*/ 50 double z4 = POW 2((double)nelem+1);50 double z4 = POW4((double)nelem+1.); 51 51 DEBUG_ENTRY( "hydro_transprob()" ); 52 53 z4 *= z4;54 52 55 53 if( ipHi >= iso.numLevels_max[ipISO][nelem]-iso.nCollapsed_max[ipISO][nelem] ) … … 69 67 Aul = H_Einstein_A( N_(ipHi), L_(ipLo)+1, N_(ipLo), L_(ipLo), nelem+1 ); 70 68 69 iso.CachedAs[ipISO][nelem][ N_(ipHi)-iso.n_HighestResolved_max[ipISO][nelem]-1 ][ ipLo ][0] = (realnum)Aul; 70 71 71 Aul *= (2.*L_(ipLo)+3.) * 2. / (2.*(double)N_(ipHi)*(double)N_(ipHi)); 72 72 … … 76 76 Aul1 = H_Einstein_A( N_(ipHi), L_(ipLo)-1, N_(ipLo), L_(ipLo), nelem+1 ); 77 77 78 iso.CachedAs[ipISO][nelem][ N_(ipHi)-iso.n_HighestResolved_max[ipISO][nelem]-1 ][ ipLo ][1] = (realnum)Aul1; 79 78 80 /* now add in this part after multiplying by stat weight for lHi = lLo-1. */ 79 81 Aul += Aul1*(2.*L_(ipLo)-1.) * 2. / (2.*(double)N_(ipHi)*(double)N_(ipHi)); 80 82 } 83 else 84 iso.CachedAs[ipISO][nelem][ N_(ipHi)-iso.n_HighestResolved_max[ipISO][nelem]-1 ][ ipLo ][1] = 0.f; 85 81 86 iso_put_error(ipISO,nelem,ipHi,ipLo,IPRAD,0.01f); 82 87 ASSERT( Aul > 0.);
