Parent Directory | Revision Log | Revision Graph | 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 | } | } |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |