ParametrizedFunctionTest.cpp File Reference

#include "DerivativeTest.hpp"
#include "TruncatedLorentzian.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"
#include <cstring>

Go to the source code of this file.

Classes
class	ParametrizedCosine

Functions
void	testDifferenceQuotient ()
void	testDerivatives ()
void	testErrorFunction ()
void	testErrorLorentzian ()
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0

---

Function Documentation

void testDifferenceQuotient

(

)

Definition at line 38 of file ParametrizedFunctionTest.cpp.

References e(), epsilon, os_, unit_assert_matrices_equal, and x.

Referenced by main().

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

    using namespace DerivativeTest;

    class TestFunction : public VectorFunction<double>
    {
        public:

        // f(x,y) = (x^2, xy, y^2)

        virtual unsigned int argumentCount() const {return 2;}
        virtual unsigned int valueCount() const {return 3;}

        virtual ublas::vector<double> operator()(ublas::vector<double> x) const
        {
            if (x.size() != argumentCount())
                throw logic_error("[TestFunction::()] Wrong argument count.");

            ublas::vector<double> result(3);
            result(0) = x(0)*x(0);
            result(1) = x(0)*x(1);
            result(2) = x(1)*x(1);
            return result;
        }
    };


    TestFunction f;
    ublas::vector<double> args(2);
    args(0) = 5; args(1) = 7;
    if (os_) *os_ << "f(5,7): " << f(args) << endl;

    if (os_) f.printDifferenceQuotientSequence(args, *os_);

    // f'(x,y) = ((2x, y, 0), (0, x, 2y))
    // f'(5,7) = ((10, 7, 0), (0, 5, 14))

    ublas::matrix<double> d(2,3);
    d(0,0) = 10;
    d(0,1) = 7;
    d(0,2) = 0;
    d(1,0) = 0;
    d(1,1) = 5;
    d(1,2) = 14;

    const double delta = 1e-9;
    const double epsilon = 1e-5;
    unit_assert_matrices_equal(d, f.differenceQuotient(args,delta), epsilon);
}

void testDerivatives

(

)

Definition at line 160 of file ParametrizedFunctionTest.cpp.

References os_.

Referenced by main().

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

    ParametrizedCosine f;

    ublas::vector<double> p(2);
    p(0) = 5;
    p(1) = M_PI/4;

    for (int i=0; i<8; i++)
        DerivativeTest::testDerivatives(f, i, p, os_);
}

void testErrorFunction

(

)

Definition at line 176 of file ParametrizedFunctionTest.cpp.

References e(), and os_.

Referenced by main().

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

    ParametrizedCosine f;

    ublas::vector<double> p(2);
    p(0) = 4;
    p(1) = 30;

    ParametrizedCosine::ErrorFunction::Data data;
    typedef ParametrizedCosine::ErrorFunction::Datum Datum;
    data.push_back(Datum(0,3));
    data.push_back(Datum(M_PI/2,0));

    ParametrizedCosine::ErrorFunction e(f, data);
    if (os_) *os_ << "error: " << e(p) << endl;

    DerivativeTest::testDerivatives<double>(e, p, os_);

    if (os_) *os_ << "8*pi^2: " << 8*M_PI*M_PI << endl;
}

void testErrorLorentzian

(

)

Definition at line 201 of file ParametrizedFunctionTest.cpp.

References AlphaI, AlphaR, e(), os_, and Tau.

Referenced by main().

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

    TruncatedLorentzian f(1);

    ublas::vector<double> p(4);
    p(TruncatedLorentzian::AlphaR) = 1;
    p(TruncatedLorentzian::AlphaI) = 5;
    p(TruncatedLorentzian::Tau) = 2;
    p(TruncatedLorentzian::F0) = 0;

    TruncatedLorentzian::ErrorFunction::Data data;
    typedef TruncatedLorentzian::ErrorFunction::Datum Datum;
    data.push_back(Datum(0,3));
    data.push_back(Datum(M_PI/2,0));

    TruncatedLorentzian::ErrorFunction e(f, data);
    if (os_) *os_ << "error: " << e(p) << endl;

    DerivativeTest::testDerivatives< complex<double> >(e, p, os_);
}

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

Definition at line 226 of file ParametrizedFunctionTest.cpp.

References e(), os_, TEST_EPILOG, TEST_FAILED, TEST_PROLOG, testDerivatives(), testDifferenceQuotient(), testErrorFunction(), and testErrorLorentzian().

{
    TEST_PROLOG(argc, argv)

    try
    {
        if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
        if (os_) *os_ << "ParametrizedFunctionTest\n";
        testDifferenceQuotient();
        testDerivatives();
        testErrorFunction();
        testErrorLorentzian();
    }
    catch (exception& e)
    {
        TEST_FAILED(e.what())
    }
    catch (...)
    {
        TEST_FAILED("Caught unknown exception.")
    }

    TEST_EPILOG
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 35 of file ParametrizedFunctionTest.cpp.