Show
Ignore:
Timestamp:
05/10/08 09:03:02 (2 months ago)
Author:
peter
Message:

Merging all changes from mainline upto r2033, except r1902.

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • branches/c08_branch/source/helike_einsta.cpp

    r1875 r2034  
    2121#include "helike_einsta.h" 
    2222#include "hydroeinsta.h" 
    23  
    24 /*lint -e662 creation of  out of bound pointer */ 
    25 /*lint -e661 creation of  out of bound pointer */ 
    2623 
    2724/* the array of transitions probabilities read from data file.  */ 
    …… 
    353350                        /* This is an extrapolation to the data given above.  The fit reproduces 
    354351                         * the above values to 10% or better. */ 
    355                         A= 7.22E-9*pow(nelem, 9.33); 
     352                        A= 7.22E-9*pow((double)nelem, 9.33); 
    356353                        iso_put_error(ipHE_LIKE,nelem,ipHi,ipLo,IPRAD,0.3f); 
    357354                } 
    …… 
    512509                                else if( lHi>=2 && lLo>=2 && nHi>nLo ) 
    513510                                { 
    514                                         /*lint -e790 integral to float */ 
    515511                                        Aul = H_Einstein_A(nHi ,lHi , nLo , lLo , nelem); 
    516                                         /*lint +e790 integral to float */ 
    517512                                        ASSERT( Aul > 0.); 
    518513 
    …… 
    724719                                else if( (nHi > nLo) && ((lHi > 2) || (lLo > 2)) ) 
    725720                                { 
    726                                         /*lint -e790 integral to float */ 
    727721                                        Aul = H_Einstein_A(nHi ,lHi , nLo , lLo , nelem); 
    728                                         /*lint +e790 integral to float */ 
    729722                                        ASSERT( Aul > 0.); 
    730723                                } 
    …… 
    10221015                        Aul = he_1trans( nelem, Enerwn, n_eff_hi, L_(ipLo)+1, 
    10231016                                S_(ipLo), -1, n_eff_lo, L_(ipLo), S_(ipLo), ipLo-3); 
    1024                         Aul *= (2.*L_(ipLo)+3.) * S_(ipLo) /  
    1025                                 (4.*(double)N_(ipHi)*(double)N_(ipHi)); 
     1017 
     1018                        iso.CachedAs[ipISO][nelem][ N_(ipHi)-iso.n_HighestResolved_max[ipISO][nelem]-1 ][ ipLo ][0] = (realnum)Aul; 
     1019 
     1020                        Aul *= (2.*L_(ipLo)+3.) * S_(ipLo) / (4.*(double)N_(ipHi)*(double)N_(ipHi)); 
     1021 
    10261022                        if( L_(ipLo) != 0 ) 
    10271023                        { 
    …… 
    10291025                                Aul1 = he_1trans( nelem, Enerwn, n_eff_hi, L_(ipLo)-1, 
    10301026                                        S_(ipLo), -1, n_eff_lo, L_(ipLo), S_(ipLo), ipLo-3); 
     1027 
     1028                                iso.CachedAs[ipISO][nelem][ N_(ipHi)-iso.n_HighestResolved_max[ipISO][nelem]-1 ][ ipLo ][1] = (realnum)Aul1; 
     1029 
    10311030                                /* now add in this part after multiplying by stat weight for lHi = lLo-1.       */ 
    1032                                 Aul += Aul1*(2.*L_(ipLo)-1.) * S_(ipLo) / 
    1033                                         (4.*(double)N_(ipHi)*(double)N_(ipHi)); 
     1031                                Aul += Aul1*(2.*L_(ipLo)-1.) * S_(ipLo) / (4.*(double)N_(ipHi)*(double)N_(ipHi)); 
    10341032                        } 
     1033                        else 
     1034                                iso.CachedAs[ipISO][nelem][ N_(ipHi)-iso.n_HighestResolved_max[ipISO][nelem]-1 ][ ipLo ][1] = 0.f; 
     1035 
    10351036                        iso_put_error(ipHE_LIKE,nelem,ipHi,ipLo,IPRAD,0.01f); 
    10361037                        ASSERT( Aul > 0.); 
    …… 
    12031204        return; 
    12041205} 
    1205 /*lint +e662 creation of  out of bound pointer */ 
    1206 /*lint +e661 creation of  out of bound pointer */