fresnel_unittest.cpp

Go to the documentation of this file.

/* +------------------------------------------------------------------------+
   |                     Mobile Robot Programming Toolkit (MRPT)            |
   |                          http://www.mrpt.org/                          |
   |                                                                        |
   | Copyright (c) 2005-2018, Individual contributors, see AUTHORS file     |
   | See: http://www.mrpt.org/Authors - All rights reserved.                |
   | Released under BSD License. See details in http://www.mrpt.org/License |
   +------------------------------------------------------------------------+ */

#include <mrpt/math/fresnel.h>
#include <mrpt/system/os.h>
#include <gtest/gtest.h>
#include <cmath>

TEST(fresnel, fresnelc)
{
    // fresnelc()
    const double test_values[][2] = {{0.0, .0},
                                     {0.79788456080286541, 0.721705924292605},
                                     {1.0, 0.779893400376823},
                                     {0.4, 0.397480759172359},
                                     {1.5, 0.445261176039822},
                                     {2.0, 0.488253406075341},
                                     {2.4, 0.554961405856428},
                                     {3.34, 0.407099627096608},
                                     {50.0, 0.499999189430728},
                                     {-0.4, -0.397480759172359},
                                     {-1.5, -0.445261176039822},
                                     {-2.0, -0.488253406075341},
                                     {-2.4, -0.554961405856428},
                                     {-3.34, -0.407099627096608},
                                     {-50.0, -0.499999189430728}};

    const unsigned int nTests = sizeof(test_values) / sizeof(test_values[0]);

    for (unsigned int i = 0; i < nTests; i++)
    {
        const double x = test_values[i][0], val_good = test_values[i][1];
        const double val = mrpt::math::fresnel_cos_integral(x);
        EXPECT_NEAR(val, val_good, 1e-5)
            << " x: " << x << "\n val_good: " << val_good << "\n val: " << val
            << "\n";
    }
}

TEST(fresnel, fresnels)
{
    // fresnelc()
    const double test_values[][2] = {{0.0, .0},
                                     {1.0, 0.438259147390355},
                                     {1.5, 0.697504960082093},
                                     {2.4, 0.619689964945684},
                                     {50.0, 0.493633802585939},
                                     {-2.0, -0.343415678363698},
                                     {-2.4, -0.619689964945684},
                                     {-3.34, -0.479600423968308},
                                     {-50.0, -0.493633802585939}};

    const unsigned int nTests = sizeof(test_values) / sizeof(test_values[0]);

    for (unsigned int i = 0; i < nTests; i++)
    {
        const double x = test_values[i][0], val_good = test_values[i][1];
        const double val = mrpt::math::fresnel_sin_integral(x);
        EXPECT_NEAR(val, val_good, 1e-5)
            << " x: " << x << "\n val_good: " << val_good << "\n val: " << val
            << "\n";
    }
}