GammaLib

To improve the efficiency of the computations it should be possible to apply numerical integration to vectors of functions, returning hence the a vector of integrals in a single method call. This allows for example to integrate in one shot a function and its derivatives.

For this purpose a new class

class GFunctions {
public:

    // Constructors and destructors
    GFunction(void);
    GFunction(const GFunction& function);
    virtual ~GFunction(void);

    // Operators
    GFunction& operator=(const GFunction& function);

    // Methods
    virtual GVector eval(const double& x) = 0;
}

needs to be introduced which an eval() method that returns a vector of values for a given argument x.

Methods

explicit GIntegral(GFunctions* kernel);
void              kernel(GFunctions* kernel);
const GFunctions* kernel(void) const;
GVector           trapzd(const double& a, const double& b, const int& n, GVector result);
GVector           romberg_vector(const double& a, const double& b, const int& order);

need to be added to the GIntegral class to support these vectors.