src/test/java/com/seaofnodes/simple/progs/newtonFloat.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

#include <stdio.h>
#include <math.h>
#include <stdlib.h>

extern double test_sqrt(double);

int main( int argc, char** argv ) {
  for( int i=0; i<10; i++ ) {
    double d = test_sqrt((double)i), expect = sqrt((double)i);
    double delta = fabs(d-expect);
    printf("%d  %f   (%g)\n",i,d,delta);
    if( delta > 1e-15 )
      return 1;
  }
  return 0;
}