Line data    Source code

       1              : #include "Test.H"
       2              : #include "Linear.H"
       3              : #include "Util/Util.H"
       4              : #include "Set/Set.H"
       5              : #include "Util/Color.H"
       6              : 
       7              : namespace Numeric
       8              : {
       9              : namespace Interpolator
      10              : {
      11              : 
      12              : template<>
      13            2 : int Test<Linear<Set::Scalar> >::Match(int verbose)
      14              : {
      15              :     // Test function:
      16              :     //
      17              :     // y(x) = { -1           x  < -1
      18              :         //        { 2*x + 1      -1 < x  < 0
      19              :     //        { 1-x          0  < x  < 1
      20              :     //        { 0            1  < x
      21              : 
      22            2 :     std::vector<Set::Scalar> xs;
      23            2 :     std::vector<Set::Scalar> ys;
      24            2 :     xs.push_back(-1.0); ys.push_back(-1.0);
      25            2 :     xs.push_back(0.0);  ys.push_back(1.0);
      26            2 :     xs.push_back(1.0);  ys.push_back(0.0);
      27              : 
      28            2 :     Linear<Set::Scalar> interp(ys,xs);
      29              : 
      30            2 :     Set::Scalar normsq = 0.0;
      31            2 :     const Set::Scalar dx = 0.01;
      32          202 :     for (Set::Scalar x = -2.0; x < -1.0; x+=dx)
      33              :     {
      34          200 :         Set::Scalar exact = -1.0;
      35          200 :         normsq += pow((interp(x) - exact)/dx,2.0);
      36          200 :         if (verbose>0) Util::Message(INFO,(normsq>1E-8 ? Color::FG::Red : Color::Reset), "x = ", x , "\texact = ", exact, "\tinterp = ", interp(x), " normsq = ", normsq,Color::Reset);
      37              :     }
      38          202 :     for (Set::Scalar x = -1.0; x < -0; x+=dx)
      39              :     {
      40          200 :         Set::Scalar exact = 2*x + 1;
      41          200 :         normsq += pow((interp(x) - exact)/dx,2.0);
      42          200 :         if (verbose>0) Util::Message(INFO,(normsq>1E-8 ? Color::FG::Red : Color::Reset), "x = ", x , "\texact = ", exact, "\tinterp = ", interp(x), " normsq = ", normsq,Color::Reset);
      43              :     }
      44          202 :     for (Set::Scalar x = 0; x < 1.0; x+=dx)
      45              :     {
      46          200 :         Set::Scalar exact = 1.0 - x;
      47          200 :         normsq += pow((interp(x) - exact)/dx,2.0);
      48          200 :         if (verbose>0) Util::Message(INFO,(normsq>1E-8 ? Color::FG::Red : Color::Reset), "x = ", x , "\texact = ", exact, "\tinterp = ", interp(x), " normsq = ", normsq,Color::Reset);
      49              :     }
      50          202 :     for (Set::Scalar x = 1.0; x < 2.0; x+=dx)
      51              :     {
      52          200 :         Set::Scalar exact = 0.0;
      53          200 :         normsq += pow((interp(x) - exact)/dx,2.0);
      54          200 :         if (verbose>0) Util::Message(INFO,(normsq>1E-8 ? Color::FG::Red : Color::Reset), "x = ", x , "\texact = ", exact, "\tinterp = ", interp(x), " normsq = ", normsq,Color::Reset);
      55              :     }
      56              : 
      57            2 :     if (normsq > 1E-8) return 1;
      58            2 :     else return 0;
      59            2 : }
      60              : }
      61              : }