integral_demo.cpp

Go to the documentation of this file.

#include "abs_diff.h"
#include "AGL_Quadrature.h"
#include "Integral.h"
#include "OscillatoryFunctor.h"
#include "PowerFunctor.h"
#include "SimpleRules.h"

#include <functional>

using std::bind;
using std::cout;
using std::endl;
using std::function;
using UA_CoMP::Num::abs_diff;
using UA_CoMP::Num::AGL_Quadrature;
using UA_CoMP::Num::OscillatoryFunctor;

int main()
{
        cout << "Verify Integral: " << endl;
        {
                cout << endl << "ComboApproximation and GaussLegendre 16 vs 32:" << endl;
                OscillatoryFunctor const ftor(-1.0, 20.0);
                AGL_Quadrature<16, 32> quad(1.0e-7, 1, 1000);

                function<double(double)> f1{bind(ftor.getIntegral(), 0.0, std::placeholders::_1)};
                function<double(double, double)> ff{make_integral(quad, ftor)};
                function<double(double)> f2{bind(ff, 0.0, std::placeholders::_1)};

                const double d{abs_diff(0.0, 3.0, 0.1, f1, f2)};
                cout << "abs diff for integral function is: " << d << endl;
        }
        return 0;
}