HPC Magazine novembre 2013 - Des outils pour calculer avec précision.

Listing 2 - Calcul de la formule de Rump avec la bibliothèque CADNA.

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

main()
{  
  cadna_init(-1);
  float_st a = 77617.;
  float_st b = 33096.;
  float_st res;

  res = (333.0 + 3.0/4.0)*(b*b*b*b*b*b) + (a*a)*(11.0*(a*a)*(b*b) - (b*b*b*b*b*b) - 
    121.0*(b*b*b*b) -2.0) + 11.0/2.0*(b*b*b*b*b*b*b*b) + a/(2.0*b);
  printf("res=%s\n",strp(res));

  res=333.75*b*b*b*b*b*b+a*a*(11.0*a*a*b*b-b*b*b*b*b*b-121.0*b*b*b*b-2.0) +
    5.5*b*b*b*b*b*b*b*b+a/(2.0*b);
  printf("res=%s\n",strp(res));

  cadna_end();
}