erfTest.cpp File Reference

#include "erf.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"
#include <limits>
#include <cstring>
#include <cmath>

Go to the source code of this file.

Functions
void	test_real ()
void	test_series ()
void	test_real_wrapper ()
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0
const double	epsilon_ = numeric_limits<double>::epsilon()

---

Function Documentation

void test_real

(

)

Definition at line 107 of file erfTest.cpp.

References e(), pwiz::math::erf(), os_, unit_assert_equal, and x.

Referenced by main().

{
    // tests our complex erf against the not-so-standard gcc-provided double erf(double)

    if (os_) *os_ << "test_real()\n";

    double a = 10;
    for (double x=-a; x<=a; x+=a/100)
    {
        complex<double> resultComplex = erf(complex<double>(x));
        double resultReal = ((double(*)(double))::erf)(x);
        if (os_) *os_ << x << " -> " << resultComplex << " " << resultReal << endl;
        unit_assert_equal(resultComplex.real(), resultReal, 1e-12);
    }

    if (os_) *os_ << endl;
}

void test_series

(

)

Definition at line 127 of file erfTest.cpp.

References e(), pwiz::math::erf(), pwiz::math::erf_series2(), os_, unit_assert_equal, x, and y.

Referenced by main().

{
    if (os_) *os_ << "test_series()\n";

    // erf_series2 matches erf in region [-2,2]x[-2,2] within 1e-10
    double a = 2;
    for (double x=-a; x<=a; x+=a/5.)
    for (double y=-a; y<=a; y+=a/5.)
    {
        complex<double> z(x,y);
        complex<double> result1 = erf(z);
        complex<double> result2 = erf_series2(z);
        if (os_) *os_ << z << ": " << abs(result1-result2) << endl;
        unit_assert_equal(abs(result1-result2), 0, 1e-10);
    }

    if (os_) *os_ << endl;
}

void test_real_wrapper

(

)

Definition at line 147 of file erfTest.cpp.

References e(), pwiz::math::erf(), os_, unit_assert_equal, and x.

Referenced by main().

{
    if (os_) *os_ << "test_real_wrapper()\n";

    double a = 10;
    for (double x=-a; x<=a; x+=a/100)
    {
        double result_pwiz = pwiz::math::erf(x);
        double result_std = ::erf(x);
        if (os_) *os_ << x << " -> " << result_pwiz << " " << result_std << endl;
        unit_assert_equal(result_pwiz, result_std, 1e-12);
    }

    if (os_) *os_ << endl;
}

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 164 of file erfTest.cpp.

References e(), os_, TEST_EPILOG, TEST_FAILED, TEST_PROLOG, test_real(), test_real_wrapper(), and test_series().

{
    TEST_PROLOG(argc, argv)

    try
    {
        if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
        if (os_) *os_ << "erfTest\n" << setprecision(20);
        test_real();
        test_series();
        test_real_wrapper();
    }
    catch (exception& e)
    {
        TEST_FAILED(e.what())
    }
    catch (...)
    {
        TEST_FAILED("Caught unknown exception.")
    }

    TEST_EPILOG
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 38 of file erfTest.cpp.

const double epsilon_ = numeric_limits<double>::epsilon()

Definition at line 39 of file erfTest.cpp.