OscillatoryFunctor.h

Go to the documentation of this file.

#pragma once

#include "AQ_Util.h"

#include <cmath>
#include <functional>

namespace UA_CoMP {
namespace Num {
using std::cos;
using std::exp;
using std::sin;
class OscillatoryFunctor : public std::unary_function<double, double>
{
public:
        using Arg = argument_type;
        using Res = result_type;

        class Integral : public std::binary_function<Arg, Arg, Res>
        {
        public:
                Integral(Arg a, Arg b) : fA(a), fB(b) {}

                result_type operator()(Arg a, Arg b) const
                {
                        return ((exp(fA * b) * (fA * sin(fB * b) - fB * cos(fB * b)) + fB) -
                                (exp(fA * (double)a) * (fA * sin(fB * a) - fB * cos(fB * a)) + fB)) /
                               (fA * fA + fB * fB);
                }

        private:
                Arg fA;
                Arg fB;
        };

        OscillatoryFunctor(argument_type a, argument_type b) : fA(a), fB(b), fIntegral(a, b) {}

        result_type operator()(argument_type x) const { return exp(fA * x) * sin(fB * x); }

        Integral getIntegral() const { return fIntegral; }

private:
        double fA;
        double fB;
        Integral fIntegral;
};

} // namespace Num
} // namespace UA_CoMP