MATH_GaussLegendreIntegration.h

Go to the documentation of this file.

#ifndef MATH_GaussLegendreIntegration_H
#define MATH_GaussLegendreIntegration_H

#include "MATH_Integration.h"

#include "CORE_Array.h"

DEFINE_SPTR(MATH_GaussLegendreIntegration);

class MATH_GaussLegendreIntegration: public MATH_Integration {
    SP_OBJECT(MATH_GaussLegendreIntegration);
    
    // ATTRIBUTES
    
public:

    
private:
    
    
private:

    //weights per dimension
    SV::CORE_RealArray mWeights;

    //points per dimension
    SV::CORE_RealArray mPoints;


   

  
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    MATH_GaussLegendreIntegration(void);
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MATH_GaussLegendreIntegration(void);   
  

public:
    
    // NEW methods
  
    static SP::MATH_GaussLegendreIntegration New() {
        SP::MATH_GaussLegendreIntegration p(new MATH_GaussLegendreIntegration(),
                                            MATH_GaussLegendreIntegration::Delete());
        p->setThis(p);
        return p;
    };
  
    
    
    // COPY methods
    
    virtual void copy(const MATH_Integration& integral);
    
 
  
    
    
    // SET methods 
    
public:
  
    void setOrder(const tUIndex& n);

    void computeIntegrationData(const tUSInt& n,
                                const tReal& a,const tReal& b,
                                tReal* w,
                                tReal* x) const;
    
protected:
    virtual void setBoundsCalled() {
        if (mPoints.getSize()>0) setOrder(mPoints.get(0)->getSize());
    }
    
    // GET methods
public:
    inline const CORE_RealArray& getPoints(const tUIndex& k) const {
        return *mPoints.get(k);
    }
    inline CORE_RealArray& getPoints(const tUIndex& k)  {
        return *mPoints.get(k);
    }
    inline const CORE_RealArray& getWeights(const tUIndex& k) const {
        return *mWeights.get(k);
    }
    inline CORE_RealArray& getWeights(const tUIndex& k)  {
        return *mWeights.get(k);
    }
  
    
    inline tUIndex getOrder() const {
        if (mPoints.getSize()>0) return mPoints.get(0)->getSize();
        return 0;
    }
    
public:

    // ------------------
    // INTEGRATE function
    // ------------------

    virtual tReal integrate();

    
    
};


#endif