MATH_TestFunction.h

Go to the documentation of this file.

#ifndef MATH_TestFunction_H
#define MATH_TestFunction_H

#include "MATH_IntegrandFunction.h"


DEFINE_SPTR(MATH_TestFunction);

class MATH_TestFunction: public MATH_IntegrandFunction {
    SP_OBJECT(MATH_TestFunction);
    
    // ATTRIBUTES
    
public:

    
private:
    tFlag mType;
    
private:

   

    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    MATH_TestFunction(void) {
        mType=0;
    }
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MATH_TestFunction(void) {
    }
  

public:
    
    // NEW methods
 
    static SP::MATH_TestFunction New() {
        SP::MATH_TestFunction p(new MATH_TestFunction(),
                                MATH_TestFunction::Delete());
        p->setThis(p);
        return p;
    };
  
    
    
    // COPY methods

    // ----------------
    // FUNCTION METHOD
    // ----------------

    inline void setType(const tFlag& t) {
        mType=t;
    }
    virtual tReal computeFunction(const tUSInt& n,const tReal X[])  {
        tReal ret=0;
        const tReal *x=X;
        tUSInt i;
        
        switch(mType) {
        case 0:
            //f(x,y,z)=xyz
            ret=(*x);
            for (i=1;i<n;i++) {
                x++;
                ret*=(*x);
            }
            break;
        case 1:
            //f(x,y,z)=1/sqrt(x^2+y^2+z^2)^3-3*x^2/sqrt(x^2+y^2+z^2)^5
            ret=(*x)*(*x);
            for (i=1;i<n;i++) {
                x++;
                ret+=(*x)*(*x);
            }
            ret=pow(ret,-1.5)-3*(*X)*(*X)*pow(ret,-2.5);
            break;
        case 2:
            //f(x,y,z)=-3*x*y/sqrt(x^2+y^2+z^2)^5
            ret=(*x)*(*x);
            for (i=1;i<n;i++) {
                x++;
                ret+=(*x)*(*x);
            }
            ret=pow(ret,-2.5);
            x=X;
            ret*=(*x);
            x++;
            ret*=(*x);
            ret*=-3.;
            break;
        }   
        return ret;
    }
    
    

    
    
};


#endif