fresnel_unittest.cpp

Go to the documentation of this file.

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

#include <gtest/gtest.h>
#include <mrpt/math/fresnel.h>
#include <mrpt/system/os.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}};

    for (auto test_value : test_values)
    {
        const double x = test_value[0], val_good = test_value[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}};

    for (auto test_value : test_values)
    {
        const double x = test_value[0], val_good = test_value[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";
    }
}