htswanalysis/MACS/lib/gsl/gsl-1.11/doc/examples/ieeeround.c

   1 #include <stdio.h>
   2 #include <gsl/gsl_math.h>
   3 #include <gsl/gsl_ieee_utils.h>
   4
   5 int
   6 main (void)
   7 {
   8   double x = 1, oldsum = 0, sum = 0;
   9   int i = 0;
  10
  11   gsl_ieee_env_setup (); /* read GSL_IEEE_MODE */
  12
  13   do
  14     {
  15       i++;
  16
  17       oldsum = sum;
  18       sum += x;
  19       x = x / i;
  20
  21       printf ("i=%2d sum=%.18f error=%g\n",
  22               i, sum, sum - M_E);
  23
  24       if (i > 30)
  25          break;
  26     }
  27   while (sum != oldsum);
  28
  29   return 0;
  30 }