MATH_CubicPolynom.h

Go to the documentation of this file.

#ifndef MATH_CubicPolynom_H
#define MATH_CubicPolynom_H


#include "MATH_Function.h"


DEFINE_SPTR(MATH_CubicPolynom);
class MATH_CubicPolynom : public MATH_Function {
    SP_OBJECT(MATH_CubicPolynom);
  
public:
 

    // ATTRIBUTES
private:
    
    //
    // ASSOCIATIONS
    tReal mCoef[3];
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    MATH_CubicPolynom(void);
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~MATH_CubicPolynom(void);   

    // NEW methods

public:
    static inline SP::MATH_CubicPolynom New() {
        SP::MATH_CubicPolynom p(new MATH_CubicPolynom(),MATH_CubicPolynom::Delete());
        p->setThis(p);
        return p;
    };
  
  

    //SET &  GET methods
    void copy(const MATH_CubicPolynom& poly);
    
    tReal& operator[](const int& i) {
        ASSERT_IN(i<=3);
        return mCoef[i];
    }
    const tReal& operator[](const int& i) const {
        ASSERT_IN(i<=3);
        return mCoef[i];
    }

    /* ! \brief solve the eqtion a*x*x+b*x+c=0
     * the method used is explain as follow:
     * https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Cardan
     */
    template<class T>
    static void solve(const tReal& a,const tReal& b,const tReal& c,const tReal& d,
                      int& n,T solutions[3]);
    
    /* ! \brief solve this=0
     */
    inline void solve(int& n,long double solutions[3]) const  {
        MATH_CubicPolynom::solve<long double>(mCoef[3],mCoef[2],mCoef[1],mCoef[0],n,solutions);
    }
    /* ! \brief solve this=0
     */
    inline void solve(int& n,double solutions[3]) const  {
        MATH_CubicPolynom::solve<double>(mCoef[3],mCoef[2],mCoef[1],mCoef[0],n,solutions);
    }
    

    //SAVER & LOADER

    void loadFromUIClass(const UI_Class& mclass);
    void saveToUIClass(UI_Class& mclass) const;
    
    // OTHERS methods

    
    virtual tReal computeFunction(const tReal& x) {
        return ComputeFunction(mCoef[3],mCoef[2],mCoef[1],mCoef[0],x);
    }
   
    static  tReal ComputeFunction(const tReal& a,const tReal& b, const tReal& c,const tReal& d,const tReal& x) {
        return  ((a*x+b)*x+c)*x+d;
    }

    virtual tString toString() const;
 
};
#include "MATH_CubicPolynom.hpp"
#endif