text: //==============================================================================
  //         Copyright 2003 - 2012   LASMEA UMR 6602 CNRS/Univ. Clermont II
  //         Copyright 2009 - 2012   LRI    UMR 8623 CNRS/Univ Paris Sud XI
  //
  //          Distributed under the Boost Software License, Version 1.0.
  //                 See accompanying file LICENSE.txt or copy at
  //                     http://www.boost.org/LICENSE_1_0.txt
  //==============================================================================
  #define NT2_UNIT_MODULE "nt2 trigonometric toolbox - secant/simd Mode"
  
  //////////////////////////////////////////////////////////////////////////////
  // unit test behavior of trigonometric components in simd mode
  //////////////////////////////////////////////////////////////////////////////
  /// created  by jt the 11/02/2011
  ///
  #include <nt2/trigonometric/include/functions/secant.hpp>
  #include <boost/simd/sdk/simd/native.hpp>
  #include <nt2/trigonometric/constants.hpp>
  extern "C" {extern long double cephes_cosl(long double);}
  
  #include <boost/type_traits/is_same.hpp>
  #include <nt2/sdk/functor/meta/call.hpp>
  #include <nt2/sdk/meta/as_integer.hpp>
  #include <nt2/sdk/meta/as_floating.hpp>
  #include <nt2/sdk/meta/as_signed.hpp>
  #include <nt2/sdk/meta/upgrade.hpp>
  #include <nt2/sdk/meta/downgrade.hpp>
  #include <nt2/sdk/meta/scalar_of.hpp>
  #include <boost/dispatch/meta/as_floating.hpp>
  #include <boost/type_traits/common_type.hpp>
  #include <nt2/sdk/unit/tests.hpp>
  #include <nt2/sdk/unit/module.hpp>
  #include <nt2/sdk/memory/buffer.hpp>
  #include <nt2/constant/constant.hpp>
  #include <nt2/sdk/meta/cardinal_of.hpp>
  #include <nt2/include/functions/splat.hpp>
  #include <nt2/include/functions/load.hpp>
  
  
  NT2_TEST_CASE_TPL ( secant_real__1_0,  NT2_SIMD_REAL_TYPES)
  {
    using nt2::secant;
    using nt2::tag::secant_;
    using nt2::load;
    using boost::simd::native;
    using nt2::meta::cardinal_of;
    typedef NT2_SIMD_DEFAULT_EXTENSION  ext_t;
    typedef typename nt2::meta::upgrade<T>::type   u_t;
    typedef native<T,ext_t>                        n_t;
    typedef n_t                                     vT;
    typedef typename nt2::meta::as_integer<T>::type iT;
    typedef native<iT,ext_t>                       ivT;
    typedef typename nt2::meta::call<secant_<nt2::medium_>(vT)>::type r_t;
    typedef typename nt2::meta::call<secant_<nt2::medium_>(T)>::type sr_t;
    typedef typename nt2::meta::scalar_of<r_t>::type ssr_t;
  
    // specific values tests
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(-nt2::Pi<vT>())[0], nt2::Mone<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(-nt2::Pi<vT>()/nt2::splat<vT>(4))[0], nt2::Sqrt_2<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Inf<vT>())[0], nt2::Nan<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Minf<vT>())[0], nt2::Nan<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Nan<vT>())[0], nt2::Nan<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Pi<vT>())[0], nt2::Mone<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Pi<vT>()/nt2::splat<vT>(4))[0], nt2::Sqrt_2<sr_t>(), 0.5);
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Zero<vT>())[0], nt2::One<sr_t>(), 0.5);
  } // end of test for floating_
  
  NT2_TEST_CASE_TPL ( secant_int_convert__1_0,  NT2_SIMD_INT_CONVERT_TYPES)
  {
    using nt2::secant;
    using nt2::tag::secant_;
    using nt2::load;
    using boost::simd::native;
    using nt2::meta::cardinal_of;
    typedef NT2_SIMD_DEFAULT_EXTENSION  ext_t;
    typedef typename nt2::meta::upgrade<T>::type   u_t;
    typedef native<T,ext_t>                        n_t;
    typedef n_t                                     vT;
    typedef typename nt2::meta::as_integer<T>::type iT;
    typedef native<iT,ext_t>                       ivT;
    typedef typename nt2::meta::call<secant_<nt2::medium_>(vT)>::type r_t;
    typedef typename nt2::meta::call<secant_<nt2::medium_>(T)>::type sr_t;
    typedef typename nt2::meta::scalar_of<r_t>::type ssr_t;
  
    // specific values tests
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Zero<vT>())[0], nt2::One<sr_t>(), 0.5);
  } // end of test for int_convert_
  
  NT2_TEST_CASE_TPL ( secant_uint_convert__1_0,  NT2_SIMD_UINT_CONVERT_TYPES)
  {
    using nt2::secant;
    using nt2::tag::secant_;
    using nt2::load;
    using boost::simd::native;
    using nt2::meta::cardinal_of;
    typedef NT2_SIMD_DEFAULT_EXTENSION  ext_t;
    typedef typename nt2::meta::upgrade<T>::type   u_t;
    typedef native<T,ext_t>                        n_t;
    typedef n_t                                     vT;
    typedef typename nt2::meta::as_integer<T>::type iT;
    typedef native<iT,ext_t>                       ivT;
    typedef typename nt2::meta::call<secant_<nt2::medium_>(vT)>::type r_t;
    typedef typename nt2::meta::call<secant_<nt2::medium_>(T)>::type sr_t;
    typedef typename nt2::meta::scalar_of<r_t>::type ssr_t;
  
    // specific values tests
    NT2_TEST_ULP_EQUAL(secant<nt2::medium_>(nt2::Zero<vT>())[0], nt2::One<sr_t>(), 0.5);
  } // end of test for uint_convert_
  