tr1_incl_test.cpp

Boost provides free peer-reviewed portable C++ source libraries. We emphasize libraries that work

//  Copyright John Maddock 2008.
//  Use, modification and distribution are subject to the
//  Boost Software License, Version 1.0. (See accompanying file
//  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
//
// Basic sanity check that header <boost/math/tr1.hpp>
// #includes all the files that it needs to.
//
#include <boost/math/tr1.hpp>
//
// Note this header includes no other headers, this is
// important if this test is to be meaningful:
//
#include "test_compile_result.hpp"

void check()
{
   float f = 0;
   double d = 0;
   long double l = 0;
   int i = 0;
   unsigned ui = 0;

   check_result<float>(boost::math::tr1::assoc_laguerre(ui, ui, f));
   check_result<float>(boost::math::tr1::assoc_laguerref(ui, ui, f));
   check_result<double>(boost::math::tr1::assoc_laguerre(ui, ui, d));
   check_result<long double>(boost::math::tr1::assoc_laguerre(ui, ui, l));
   check_result<long double>(boost::math::tr1::assoc_laguerrel(ui, ui, l));
   check_result<double>(boost::math::tr1::assoc_laguerre(ui, ui, i));
   check_result<double>(boost::math::tr1::assoc_laguerre(ui, ui, ui));

   check_result<float>(boost::math::tr1::assoc_legendre(ui, ui, f));
   check_result<float>(boost::math::tr1::assoc_legendref(ui, ui, f));
   check_result<double>(boost::math::tr1::assoc_legendre(ui, ui, d));
   check_result<long double>(boost::math::tr1::assoc_legendre(ui, ui, l));
   check_result<long double>(boost::math::tr1::assoc_legendrel(ui, ui, l));
   check_result<double>(boost::math::tr1::assoc_legendre(ui, ui, i));
   check_result<double>(boost::math::tr1::assoc_legendre(ui, ui, ui));

   check_result<float>(boost::math::tr1::beta(f, f));
   check_result<float>(boost::math::tr1::betaf(f, f));
   check_result<double>(boost::math::tr1::beta(d, d));
   check_result<long double>(boost::math::tr1::beta(l, l));
   check_result<long double>(boost::math::tr1::betal(l, l));
   check_result<double>(boost::math::tr1::beta(ui, ui));
   check_result<double>(boost::math::tr1::beta(i, ui));
   check_result<double>(boost::math::tr1::beta(f, d));
   check_result<long double>(boost::math::tr1::beta(l, d));

   check_result<float>(boost::math::tr1::comp_ellint_1(f));
   check_result<float>(boost::math::tr1::comp_ellint_1f(f));
   check_result<double>(boost::math::tr1::comp_ellint_1(d));
   check_result<long double>(boost::math::tr1::comp_ellint_1(l));
   check_result<long double>(boost::math::tr1::comp_ellint_1l(l));
   check_result<double>(boost::math::tr1::comp_ellint_1(ui));
   check_result<double>(boost::math::tr1::comp_ellint_1(i));

   check_result<float>(boost::math::tr1::comp_ellint_2(f));
   check_result<float>(boost::math::tr1::comp_ellint_2f(f));
   check_result<double>(boost::math::tr1::comp_ellint_2(d));
   check_result<long double>(boost::math::tr1::comp_ellint_2(l));
   check_result<long double>(boost::math::tr1::comp_ellint_2l(l));
   check_result<double>(boost::math::tr1::comp_ellint_2(ui));
   check_result<double>(boost::math::tr1::comp_ellint_2(i));

   check_result<float>(boost::math::tr1::comp_ellint_3(f, f));
   check_result<float>(boost::math::tr1::comp_ellint_3f(f, f));
   check_result<double>(boost::math::tr1::comp_ellint_3(d, d));
   check_result<long double>(boost::math::tr1::comp_ellint_3(l, l));
   check_result<long double>(boost::math::tr1::comp_ellint_3l(l, l));
   check_result<double>(boost::math::tr1::comp_ellint_3(ui, ui));
   check_result<double>(boost::math::tr1::comp_ellint_3(i, ui));
   check_result<double>(boost::math::tr1::comp_ellint_3(f, d));
   check_result<long double>(boost::math::tr1::comp_ellint_3(l, d));

   check_result<float>(boost::math::tr1::cyl_bessel_i(f, f));
   check_result<float>(boost::math::tr1::cyl_bessel_if(f, f));
   check_result<double>(boost::math::tr1::cyl_bessel_i(d, d));
   check_result<long double>(boost::math::tr1::cyl_bessel_i(l, l));
   check_result<long double>(boost::math::tr1::cyl_bessel_il(l, l));
   check_result<double>(boost::math::tr1::cyl_bessel_i(ui, ui));
   check_result<double>(boost::math::tr1::cyl_bessel_i(i, ui));
   check_result<double>(boost::math::tr1::cyl_bessel_i(f, d));
   check_result<long double>(boost::math::tr1::cyl_bessel_i(l, d));

   check_result<float>(boost::math::tr1::cyl_bessel_j(f, f));
   check_result<float>(boost::math::tr1::cyl_bessel_jf(f, f));
   check_result<double>(boost::math::tr1::cyl_bessel_j(d, d));
   check_result<long double>(boost::math::tr1::cyl_bessel_j(l, l));
   check_result<long double>(boost::math::tr1::cyl_bessel_jl(l, l));
   check_result<double>(boost::math::tr1::cyl_bessel_j(ui, ui));
   check_result<double>(boost::math::tr1::cyl_bessel_j(i, ui));
   check_result<double>(boost::math::tr1::cyl_bessel_j(f, d));
   check_result<long double>(boost::math::tr1::cyl_bessel_j(l, d));

   check_result<float>(boost::math::tr1::cyl_bessel_k(f, f));
   check_result<float>(boost::math::tr1::cyl_bessel_kf(f, f));
   check_result<double>(boost::math::tr1::cyl_bessel_k(d, d));
   check_result<long double>(boost::math::tr1::cyl_bessel_k(l, l));
   check_result<long double>(boost::math::tr1::cyl_bessel_kl(l, l));
   check_result<double>(boost::math::tr1::cyl_bessel_k(ui, ui));
   check_result<double>(boost::math::tr1::cyl_bessel_k(i, ui));
   check_result<double>(boost::math::tr1::cyl_bessel_k(f, d));
   check_result<long double>(boost::math::tr1::cyl_bessel_k(l, d));

   check_result<float>(boost::math::tr1::cyl_neumann(f, f));
   check_result<float>(boost::math::tr1::cyl_neumannf(f, f));
   check_result<double>(boost::math::tr1::cyl_neumann(d, d));
   check_result<long double>(boost::math::tr1::cyl_neumann(l, l));
   check_result<long double>(boost::math::tr1::cyl_neumannl(l, l));
   check_result<double>(boost::math::tr1::cyl_neumann(ui, ui));
   check_result<double>(boost::math::tr1::cyl_neumann(i, ui));
   check_result<double>(boost::math::tr1::cyl_neumann(f, d));
   check_result<long double>(boost::math::tr1::cyl_neumann(l, d));

   check_result<float>(boost::math::tr1::ellint_1(f, f));
   check_result<float>(boost::math::tr1::ellint_1f(f, f));
   check_result<double>(boost::math::tr1::ellint_1(d, d));
   check_result<long double>(boost::math::tr1::ellint_1(l, l));
   check_result<long double>(boost::math::tr1::ellint_1l(l, l));
   check_result<double>(boost::math::tr1::ellint_1(ui, ui));
   check_result<double>(boost::math::tr1::ellint_1(i, ui));
   check_result<double>(boost::math::tr1::ellint_1(f, d));
   check_result<long double>(boost::math::tr1::ellint_1(l, d));

   check_result<float>(boost::math::tr1::ellint_2(f, f));
   check_result<float>(boost::math::tr1::ellint_2f(f, f));
   check_result<double>(boost::math::tr1::ellint_2(d, d));
   check_result<long double>(boost::math::tr1::ellint_2(l, l));
   check_result<long double>(boost::math::tr1::ellint_2l(l, l));
   check_result<double>(boost::math::tr1::ellint_2(ui, ui));
   check_result<double>(boost::math::tr1::ellint_2(i, ui));
   check_result<double>(boost::math::tr1::ellint_2(f, d));
   check_result<long double>(boost::math::tr1::ellint_2(l, d));

   check_result<float>(boost::math::tr1::ellint_3(f, f, f));
   check_result<float>(boost::math::tr1::ellint_3f(f, f, f));
   check_result<double>(boost::math::tr1::ellint_3(d, d, d));
   check_result<long double>(boost::math::tr1::ellint_3(l, l, l));
   check_result<long double>(boost::math::tr1::ellint_3l(l, l, l));
   check_result<double>(boost::math::tr1::ellint_3(ui, ui, i));
   check_result<double>(boost::math::tr1::ellint_3(i, ui, f));
   check_result<double>(boost::math::tr1::ellint_3(f, d, i));
   check_result<long double>(boost::math::tr1::ellint_3(l, d, f));

   check_result<float>(boost::math::tr1::expint(f));
   check_result<float>(boost::math::tr1::expintf(f));
   check_result<double>(boost::math::tr1::expint(d));
   check_result<long double>(boost::math::tr1::expint(l));
   check_result<long double>(boost::math::tr1::expintl(l));
   check_result<double>(boost::math::tr1::expint(ui));
   check_result<double>(boost::math::tr1::expint(i));

   check_result<float>(boost::math::tr1::hermite(ui, f));
   check_result<float>(boost::math::tr1::hermitef(ui, f));
   check_result<double>(boost::math::tr1::hermite(ui, d));
   check_result<long double>(boost::math::tr1::hermite(ui, l));
   check_result<long double>(boost::math::tr1::hermitel(ui, l));
   check_result<double>(boost::math::tr1::hermite(ui, i));
   check_result<double>(boost::math::tr1::hermite(ui, ui));

   check_result<float>(boost::math::tr1::laguerre(ui, f));
   check_result<float>(boost::math::tr1::laguerref(ui, f));
   check_result<double>(boost::math::tr1::laguerre(ui, d));
   check_result<long double>(boost::math::tr1::laguerre(ui, l));
   check_result<long double>(boost::math::tr1::laguerrel(ui, l));
   check_result<double>(boost::math::tr1::laguerre(ui, i));
   check_result<double>(boost::math::tr1::laguerre(ui, ui));

   check_result<float>(boost::math::tr1::legendre(ui, f));
   check_result<float>(boost::math::tr1::legendref(ui, f));
   check_result<double>(boost::math::tr1::legendre(ui, d));
   check_result<long double>(boost::math::tr1::legendre(ui, l));
   check_result<long double>(boost::math::tr1::legendrel(ui, l));
   check_result<double>(boost::math::tr1::legendre(ui, i));
   check_result<double>(boost::math::tr1::legendre(ui, ui));

   check_result<float>(boost::math::tr1::riemann_zeta(f));
   check_result<float>(boost::math::tr1::riemann_zetaf(f));
   check_result<double>(boost::math::tr1::riemann_zeta(d));
   check_result<long double>(boost::math::tr1::riemann_zeta(l));
   check_result<long double>(boost::math::tr1::riemann_zetal(l));
   check_result<double>(boost::math::tr1::riemann_zeta(ui));
   check_result<double>(boost::math::tr1::riemann_zeta(i));

   check_result<float>(boost::math::tr1::sph_bessel(ui, f));