Changeset 1967

Show
Ignore:
Timestamp:
04/24/08 06:19:37 (3 weeks ago)
Author:
gary
Message:

add sum of absolute values of residuals from case B assert - can see at a glance how good or bad the result is.

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • trunk/source/assertresults.h

    r1732 r1967  
    2727EXTERN bool lgAssertsOK , lgBigBotch, lgPrtSciNot; 
    2828 
     29EXTERN struct t_assertresults { 
     30        double SumErrorCaseAssert; 
     31        long int nSumErrorCaseAssert; 
     32 
     33}       assertresults; 
     34 
    2935#endif /* _ASSERTRESULTS_H_ */ 
  • trunk/source/assert_results.cpp

    r1962 r1967  
    19251925                                                AssertQuantity[i] = 0.; 
    19261926                                                RelError[i] = MAX2( RelError[i] , fabs(error) ); 
     1927 
     1928                                                /* save sum which we will report later */ 
     1929                                                assertresults.SumErrorCaseAssert += RelError[i]; 
     1930                                                ++assertresults.nSumErrorCaseAssert; 
     1931 
    19271932                                        } 
    19281933                                } 
     
    19641969                                        AssertQuantity[i] = 0.; 
    19651970                                        RelError[i] = MAX2( RelError[i] , fabs(error) ); 
     1971 
     1972                                        /* save sum which we will report later */ 
     1973                                        assertresults.SumErrorCaseAssert += RelError[i]; 
     1974                                        ++assertresults.nSumErrorCaseAssert; 
    19661975                                } 
    19671976                                fprintf(ioQQQ,"\n"); 
     
    26042613                } 
    26052614 
     2615                if( assertresults.nSumErrorCaseAssert>0 ) 
     2616                { 
     2617                        fprintf(ioQQQ,"\n The sum of the %li Case A and B assert relative " 
     2618                                "residuals add up to %.2f\n\n" ,  
     2619                                assertresults.nSumErrorCaseAssert, 
     2620                                assertresults.SumErrorCaseAssert /assertresults.nSumErrorCaseAssert ); 
     2621                } 
     2622 
    26062623                /* explain how we were compiled, but only if printing time */ 
    26072624                if( prt.lgPrintTime ) 
  • trunk/source/init_coreload_postparse.cpp

    r1776 r1967  
    44* after parser, will be called one time per core load */ 
    55#include "cddefines.h"  
     6#include "assertresults.h" 
    67#include "init.h"  
    78 
     
    2425        ++nCalled; 
    2526 
     27        assertresults.SumErrorCaseAssert = 0.; 
     28        assertresults.nSumErrorCaseAssert = 0; 
     29 
    2630        return; 
    2731}