common.h

   1 /*
   2  * Generally useful stuff.
   3  */
   4
   5 #ifndef COMMON_H
   6 #define COMMON_H
   7
   8 #define _GNU_SOURCE
   9 #include <math.h>
  10 #include <limits.h>
  11
  12 double hypotD(const double p[D3], const double q[D3]);
  13 double hypotD2(const double p[D3], const double q[D3]);
  14 double hypotD2plus(const double p[D3], const double q[D3], double add);
  15
  16 #ifdef FP_FAST_FMA
  17 # define fma_fast fma
  18 #else
  19 # define fma_fast(f1,f2,t) ((f1)*(f2)+(t))
  20 #endif
  21 #define ffsqa(factor,term) fma_fast((factor),(factor),(term))
  22
  23 #endif /*COMMON_H*/