MATH_Integration.h

Go to the documentation of this file.

#ifndef MATH_Integration_H
#define MATH_Integration_H

#include "CORE_Object.h"

#include "CORE_Array.h"

DEFINE_SPTR(MATH_IntegrandFunction);

DEFINE_SPTR(MATH_Integration);

class MATH_Integration: public CORE_Object {
    SP_OBJECT(MATH_Integration);

    // ATTRIBUTES

public:
private:

 
private:
    //bounds of intergration
    CORE_RealArray mBounds;

    // ASSOCIATIONS

    //function to integrate
    WP::MATH_IntegrandFunction mFunction;
    
protected: 
    // METHODS
    
    // CONSTRUCTORS 
  
    MATH_Integration(void) {
        mBounds.setSize(2);
        mBounds[0]=-1;
        mBounds[1]=1;
    }
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MATH_Integration(void) {
    }
  

public:

    // NEW methods
 
    
    // COPY methods
    
    virtual void copy(const MATH_Integration& integral) {
        mBounds.copy(integral.getBounds());
       
    }
 
  
    
    
public:
  
    
    inline void setFunction(SP::MATH_IntegrandFunction f) {
        mFunction=f;
    }; 
   

    void setBounds(const tUSInt& n,const tReal I[]);
    
    inline void setBounds(const tReal& a,const tReal& b) {
        mBounds.setSize(2);
        mBounds[0]=(a<b)?a:b;
        mBounds[1]=(a>b)?a:b;
        setBoundsCalled();
    }
    inline void setBounds(const tReal& a1,const tReal& b1,const tReal& a2,const tReal& b2) {
        mBounds.setSize(4);
        mBounds[0]=(a1<b1)?a1:b1;
        mBounds[1]=(a1>b1)?a1:b1;
        mBounds[2]=(a2<b2)?a2:b2;
        mBounds[3]=(a2>b2)?a2:b2;
        setBoundsCalled();
    }
    inline void setBounds(const tReal& a1,const tReal& b1,const tReal& a2,const tReal& b2,const tReal& a3,const tReal& b3) {
        mBounds.setSize(6);
        mBounds[0]=(a1<b1)?a1:b1;
        mBounds[1]=(a1>b1)?a1:b1;
        mBounds[2]=(a2<b2)?a2:b2;
        mBounds[3]=(a2>b2)?a2:b2;
        mBounds[4]=(a3<b3)?a3:b3;
        mBounds[5]=(a3>b3)?a3:b3;
        setBoundsCalled();
    }
protected:
    virtual void setBoundsCalled()=0;
    
    // GET methods 

public:



    inline tBoolean hasFunction() const {
        return (mFunction.lock().get()!=null);
    }
  
    inline const MATH_IntegrandFunction& getFunction() const {
        SPC::MATH_IntegrandFunction p=mFunction.lock();
        return *p.get();
    };
    inline MATH_IntegrandFunction& getFunction() {
        SP::MATH_IntegrandFunction p=mFunction.lock();
        return *p.get();
    };

   
    inline const CORE_RealArray& getBounds() const {
        return mBounds;
    } 
    inline const tReal& getBound(const tUSInt& i) const {
        return mBounds[i];
    }

    inline tUSInt getDimension() const {
        return mBounds.getSize()/2;
        
    }
public:
    // ------------------
    // INTEGRATE function
    // ------------------

    virtual tReal integrate(SP::MATH_IntegrandFunction f) {
        setFunction(f);
        return integrate();
    }
    
    virtual tReal integrate()=0;
    
    
    
};


#endif