ceph/src/boost/libs/math/test/test_expint_hooks.hpp

   1 //  (C) Copyright John Maddock 2006.
   2 //  Use, modification and distribution are subject to the
   3 //  Boost Software License, Version 1.0. (See accompanying file
   4 //  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
   5
   6 #ifndef BOOST_MATH_TEST_ZETA_OTHER_HOOKS_HPP
   7 #define BOOST_MATH_TEST_ZETA_OTHER_HOOKS_HPP
   8
   9
  10 #ifdef TEST_CEPHES
  11 namespace other{
  12 extern "C" {
  13    double expn(int, double);
  14    float expnf(int, float);
  15    long double expnl(int, long double);
  16 }
  17 inline float expint(unsigned n, float a)
  18 { return expnf(n, a); }
  19 inline double expint(unsigned n, double a)
  20 { return expn(n, a); }
  21 inline long double expint(unsigned n, long double a)
  22 {
  23 #ifdef BOOST_MSVC
  24    return expn(n, (double)a);
  25 #else
  26    return expnl(n, a);
  27 #endif
  28 }
  29 // Ei is not supported:
  30 template <class T>
  31 inline T expint(T){ return 0; }
  32 }
  33 #define TEST_OTHER
  34 #endif
  35
  36 #ifdef TEST_GSL
  37 #include <gsl/gsl_sf_expint.h>
  38
  39 namespace other{
  40 inline float expint(float a)
  41 { return (float)gsl_sf_expint_Ei(a); }
  42 inline double expint(double a)
  43 { return gsl_sf_expint_Ei(a); }
  44 inline long double expint(long double a)
  45 { return gsl_sf_expint_Ei(a); }
  46 // En is not supported:
  47 template <class T>
  48 inline T expint(unsigned, T){ return 0; }
  49 }
  50 #define TEST_OTHER
  51 #endif
  52
  53 #ifdef TEST_SPECFUN
  54 namespace other{
  55 extern "C" int calcei_(double *arg, double *result, int*);
  56 inline float expint(float a)
  57 {
  58    double r, a_(a);
  59    int v = 1;
  60    calcei_(&a_, &r, &v);
  61    return (float)r;
  62 }
  63 inline double expint(double a)
  64 {
  65    double r, a_(a);
  66    int v = 1;
  67    calcei_(&a_, &r, &v);
  68    return r;
  69 }
  70 inline long double expint(long double a)
  71 {
  72    double r, a_(a);
  73    int v = 1;
  74    calcei_(&a_, &r, &v);
  75    return r;
  76 }
  77 // En is not supported:
  78 template <class T>
  79 inline T expint(unsigned, T){ return 0; }
  80 }
  81 #define TEST_OTHER
  82 #endif
  83
  84 #ifdef TEST_OTHER
  85 namespace other{
  86    boost::math::concepts::real_concept expint(unsigned, boost::math::concepts::real_concept){ return 0; }
  87    boost::math::concepts::real_concept expint(boost::math::concepts::real_concept){ return 0; }
  88 }
  89 #endif
  90
  91
  92 #endif
  93
  94