Changeset 1988

Timestamp:

04/29/08 14:19:04 (2 weeks ago)

Author:

gary

Message:

assert_results.cpp - print dashes in individual Case B line assert to see relative error, as in final set of asserts

rt_line_one.cpp - when Case B in place do not assert that damping constant is greater than 0, only that it is not negative

Files:

• trunk/source/assert_results.cpp (modified) (5 diffs)
• trunk/source/rt_line_one.cpp (modified) (1 diff)

trunk/source/assert_results.cpp

@@ -1924 +1924 @@
-                        if( nISOCaseB == ipH_LIKE )
-                        {
-
+        
-                                /* limit of 10 is because that is all we printed and saved in prt_lines_hydro */
-                                for( long int ipLo=1+iCase; ipLo< MIN2(10,nHighestPrinted-1); ++ipLo )
@@ -1946 +1946 @@
-                                                        absint = pow( 10., absint );
-                                                        error = (CBabsint - absint)/MAX2(CBabsint , absint);
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-                                                                chElemLabelCaseB , ipHi , ipLo );
-                                                        prt_wl(ioQQQ, wl );
@@ -1973 +2011 @@
-                        else if( nISOCaseB == ipHE_LIKE )
-                        {
+                                if( !dense.lgElmtOn[ipHELIUM] )
+                                {
+                                        fprintf(ioQQQ,"DISASTER assert case B for a He is requested but He is not "
+                                                "included.\n");
+                                        fprintf(ioQQQ,"Do not turn off He is you want to assert its spectrum.\n");
+                                        cdEXIT(EXIT_FAILURE);
+                                }
-#                               define N_CASEB_HEI      11
-                                realnum CaseBWlHeI[N_CASEB_HEI]=
@@ -1978 +2023 @@
-                                4026.f, 6678.f, 4922.f };
-                                /* do He I as special case */
-Species
+          
-                                for( long int ipLine=0; ipLine< N_CASEB_HEI ; ++ipLine )
-                                {
@@ -1991 +2036 @@
-                                                absint = pow( 10., absint );
-                                                error = (CBabsint - absint)/MAX2(CBabsint , absint);
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-                                                prt_wl(ioQQQ, wl );
-                                                fprintf(ioQQQ," %.2e %.2e %10.3f", 

trunk/source/rt_line_one.cpp

@@ -294 +294 @@
-                fprintf(ioQQQ,"DEBUG hit it \n");
-        }*/
-
+
+
-
-        /* static solution - which type of line will determine