Changeset 1985

Show
Ignore:
Timestamp:
04/29/08 01:03:48 (2 weeks ago)
Author:
gary
Message:

only reevaluate H branching ratios when static opacities are reevaluated. this speeds up code by 1/3 with no impact on asserted results.

with this change the trunk is now about 2x slower than C07 was.

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/source/iso_solve.cpp

    r1969 r1985  
    88#include "conv.h" 
    99#include "dense.h" 
     10#include "opacity.h" 
    1011#include "elementnames.h" 
    1112#include "h2.h" 
    …… 
    114115                                iso_radiative_recomb( ipISO , nelem ); 
    115116 
    116                                 if( ipISO==ipH_LIKE && nelem<=ipHELIUM ) 
    117                                 { 
    118                                         iso_collapsed_bnl_set( ipISO, nelem ); 
    119  
    120                                         //iso_collapsed_bnl_print( ipISO, nelem ); 
    121  
    122                                         iso_collapsed_Aul_update( ipISO, nelem ); 
    123  
    124                                         iso_collapsed_lifetimes_update( ipISO, nelem ); 
    125  
    126                                         iso_cascade( ipISO, nelem ); 
    127  
    128                                         iso_radiative_recomb_effective( ipISO, nelem ); 
    129                                 } 
    130                                 else 
    131                                 { 
    132                                         iso_cascade( ipISO, nelem ); 
    133  
    134                                         iso_radiative_recomb_effective( ipISO, nelem ); 
     117                                if( opac.lgRedoStatic ) 
     118                                { 
     119                                        if( ipISO==ipH_LIKE && nelem<=ipHELIUM ) 
     120                                        { 
     121                                                iso_collapsed_bnl_set( ipISO, nelem ); 
     122 
     123                                                //iso_collapsed_bnl_print( ipISO, nelem ); 
     124 
     125                                                iso_collapsed_Aul_update( ipISO, nelem ); 
     126 
     127                                                iso_collapsed_lifetimes_update( ipISO, nelem ); 
     128 
     129                                                iso_cascade( ipISO, nelem ); 
     130 
     131                                                iso_radiative_recomb_effective( ipISO, nelem ); 
     132                                        } 
     133                                        else 
     134                                        { 
     135                                                iso_cascade( ipISO, nelem ); 
     136 
     137                                                iso_radiative_recomb_effective( ipISO, nelem ); 
     138                                        } 
    135139                                } 
    136140