Changeset 685 for trunk/source/rt_stark.cpp

Timestamp:

12/05/06 18:14:43 (2 years ago)

Author:

gary

Message:

dayslow_master.pl - use icl rather than gcc for slow test suite
ion_solver.cpp - undo bug introduced in rev 622. now looks at matrix elements to see if singular
rt_stark.cpp - resolved ambiguous overload
blr_n09_p18_Z20.in - add punch convergence error command
limit_hi_ion.in - renamed func_hi_ion.in to this since it tests code in high ionization limit

Files:

• trunk/source/rt_stark.cpp (modified) (1 diff)

trunk/source/rt_stark.cpp

@@ -58 +58 @@
          * with negative optical depths,
          * NB did not understand why neg optical depths started */
-        aa = pow(SDIV(EmisLines[ipH_LIKE][ipHYDROGEN][ipH2p][ipH1s].TauIn),-0.75f);
+        aa = (float)SDIV(EmisLines[ipH_LIKE][ipHYDROGEN][ipH2p][ipH1s].TauIn);
+        aa = pow( aa ,-0.75);
         hydro.pestrk[ipH1s][ipH2p] = strkla/2.*MAX2(1.,aa);