/[MITgcm]/MITgcm/verification/testreport
ViewVC logotype

Diff of /MITgcm/verification/testreport

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.63 by jmc, Wed Jan 26 00:24:06 2005 UTC revision 1.64 by edhill, Mon Jan 31 19:26:24 2005 UTC
# Line 529  int main( int argc, char** argv )  { Line 529  int main( int argc, char** argv )  {
529      abave = 0.5*(fabs(a)+fabs(b));      abave = 0.5*(fabs(a)+fabs(b));
530      if (abave > 0.0) {      if (abave > 0.0) {
531        relerr=fabs(a-b)/abave;        relerr=fabs(a-b)/abave;
532        if (relerr > 0.0) { linnum = (int)lrint(log10(relerr)); }        if (relerr > 0.0) { linnum = (int)rint(log10(relerr)); }
533        else { linnum = -16 ; }        else { linnum = -16 ; }
534        best = (best > linnum) ? best : linnum;        best = (best > linnum) ? best : linnum;
535      }      }

Legend:
Removed from v.1.63  
changed lines
  Added in v.1.64

  ViewVC Help
Powered by ViewVC 1.1.22