testsuite/supernova/approximation_test.cpp

   1 #include <iostream>
   2 #include <boost/test/unit_test.hpp>
   3 #include <boost/test/floating_point_comparison.hpp>
   4
   5 #include <cmath>
   6
   7 #include "utilities/fast_approximations.hpp"
   8
   9 using namespace nova;
  10 using namespace std;
  11
  12 namespace
  13 {
  14
  15 template <typename F>
  16 void run_fast_atan2_tests(void)
  17 {
  18     for (F y = -4; y < 4; y+=0.01)
  19         for (F x = -4; x < 4; x+=0.01)
  20         {
  21             if (abs(x) < 1e-8)
  22                 continue;
  23
  24             F accurate = atan2(y, x);
  25             F fast = fast_atan2(y, x);
  26             F faster = faster_atan2(y, x);
  27
  28             BOOST_CHECK_SMALL( accurate - fast,   F(0.011) );
  29             BOOST_CHECK_SMALL( accurate - faster, F(0.072) );
  30         }
  31 }
  32
  33 }
  34
  35 BOOST_AUTO_TEST_CASE( atan2_tests_float )
  36 {
  37     run_fast_atan2_tests<float>();
  38 }
  39
  40 BOOST_AUTO_TEST_CASE( atan2_tests_double )
  41 {
  42     run_fast_atan2_tests<double>();
  43 }