ossimTerm.h

Go to the documentation of this file.

#ifndef OSSIM_TERM_H
#define OSSIM_TERM_H 1

#include <ossim/base/ossimConstants.h>

#include <vector>
#include <cmath>



enum termType : ossim_uint32 {
    monomial,
    polynomial,
    variable,
    constant,
    product
};



class OSSIMDLLEXPORT ossimTerm
{
public:


   virtual ~ossimTerm() {}


   virtual ossim_float64    evaluate(std::vector<ossim_float64> /* point */ ) const = 0;
   virtual ossimTerm*       derivative(const ossim_uint32& variable) const = 0;


   virtual ossim_uint32 getType() const = 0;
   virtual ossimTerm*       copy() const = 0;


   inline ossim_float64 operator()(std::vector<ossim_float64> point) const{ return this->evaluate(point); }


   inline ossimTerm*    operator+(ossimTerm* rhs) const;
   inline ossimTerm*    operator*(ossimTerm* rhs) const;

};



class OSSIMDLLEXPORT ossimMonom : public ossimTerm{
public:


    ossimMonom(const ossim_float64& coeffient, ossimTerm* variable, const ossim_float64& exponential) : coef(coeffient), var(variable), exp(exponential) {}
    ossimMonom(const ossimMonom& src) : coef(src.coef), var(src.var->copy()), exp(src.exp) {}
    ~ossimMonom() { delete var; }


    ossim_float64   evaluate(std::vector<ossim_float64> point) const { return (coef * pow(var->evaluate(point), exp)); }
    ossimTerm*      derivative(const ossim_uint32& variable) const;


    ossim_uint32    getType() const { return termType::monomial; }
    ossimTerm*      copy() const;

private:
    ossim_float64   coef;       // The coefficent of the monomial 
    ossimTerm*      var;        // The variable number of the monomial (starts at 0 - n)
    ossim_float64   exp;        // The exponenent of the monomial

};



class OSSIMDLLEXPORT ossimPolynomial : public ossimTerm{
public:


    ossimPolynomial(ossimTerm* leftTerm, ossimTerm* rightTerm) : lTerm(leftTerm), rTerm(rightTerm) {}
    ossimPolynomial(const ossimPolynomial& src) : lTerm(src.lTerm->copy()), rTerm(src.rTerm->copy()) {}
    ~ossimPolynomial() { delete lTerm; delete rTerm; }


    ossim_float64   evaluate(std::vector<ossim_float64> point) const { return lTerm->evaluate(point) + rTerm->evaluate(point); }
    ossimTerm*      derivative(const ossim_uint32& variable) const;


    ossim_uint32    getType() const { return termType::polynomial; }
    ossimTerm*      copy() const;

private:
    ossimTerm* lTerm;       // The left term of the polynomial
    ossimTerm* rTerm;       // The right term of the polynomial

};



class OSSIMDLLEXPORT ossimProd : public ossimTerm{
public:


    ossimProd(ossimTerm* leftTerm, ossimTerm* rightTerm) : lTerm(leftTerm), rTerm(rightTerm) {}
    ossimProd(const ossimProd& src) : lTerm(src.lTerm->copy()), rTerm(src.rTerm->copy()) {}
    ~ossimProd() { delete lTerm; delete rTerm; }


    ossim_float64   evaluate(std::vector<ossim_float64> point) const { return lTerm->evaluate(point) * rTerm->evaluate(point); }
    ossimTerm*      derivative(const ossim_uint32& variable) const;


    ossim_uint32    getType() const { return termType::product; }
    ossimTerm*      copy() const;

private:
    ossimTerm* lTerm;       // The left term of the product
    ossimTerm* rTerm;       // The right term of the product

};



class OSSIMDLLEXPORT ossimVar : public ossimTerm{
public:


    ossimVar(const ossim_uint32& variable) : var(variable) {}
    ossimVar(const ossimVar& src) : var(src.var) {}
    ~ossimVar() {}


    ossim_float64   evaluate(std::vector<ossim_float64> point) const{ return point[var]; }
    ossimTerm*      derivative(const ossim_uint32& variable) const;


    ossim_uint32    getType() const { return termType::variable; }
    ossimTerm*      copy() const;

private:
    ossim_uint32 var;       // The number of the variable (0-n)

};



class OSSIMDLLEXPORT ossimConst : public ossimTerm
{
public:


ossimConst(const ossim_float64& value) : val(value) {}
ossimConst(const ossimConst& src) : val(src.val) {}
   ~ossimConst() {}


   ossim_float64    evaluate(std::vector<ossim_float64> /* point */ ) const{ return val; }
   ossimTerm*       derivative(const ossim_uint32& variable) const;


   ossim_uint32 getType() const { return termType::constant; }
   ossimTerm*       copy() const;

private:
   ossim_float64 val;   // The value of the constant

};

#endif // OSSIM_TERM_H