MATH_Solver.h

Go to the documentation of this file.

#ifndef MATH_Solver_H
#define MATH_Solver_H

#include "CORE_Object.h"


#include "MATH_Matrix.h"


DEFINE_SPTR(MATH_Solver);

class MATH_Solver : public CORE_Object  {
    SP_OBJECT(MATH_Solver);


private:
    tReal mTolerance;
    tReal mError;
    tUInt mMaxIterationsNumber;
    tUInt mIterationsNumber;
    
    SP::MATH_Vector mW;
    
protected:
    MATH_Solver(void) {
        mMaxIterationsNumber=100;
        mTolerance=1.e-12;
        mIterationsNumber=0;
        mError=0;
    }
    
    
    
    // DESTRUCTORS 
    
  
    virtual ~MATH_Solver(void) {
    }

public:
    // copy
    inline void copy(const MATH_Solver& solver) {
        mMaxIterationsNumber=solver.getMaximumIterationsNumber();
        mTolerance=solver.getTolerance();
        mIterationsNumber=solver.getIterationsNumber();
        mError=solver.getError();
    }
      
    //SET Methods
    //==========

    inline void setMaximumIterationsNumber(const tUInt& nIterations) {
        mMaxIterationsNumber=nIterations;
    }
    inline void setTolerance(const tReal& eps) {
        mTolerance=eps;
    }
    //GET Methods
    //==========

    inline const tUInt& getMaximumIterationsNumber() const {
        return mMaxIterationsNumber;
    }
    inline const tReal&  getTolerance() const {
        return mTolerance;
    }

    inline const tUInt& getIterationsNumber() const {
        return mIterationsNumber;
    }

    inline tUInt& getIterationsNumber() {
        return mIterationsNumber;
    }
    inline tReal& getError() {
        return mError;
    }
    inline const tReal& getError() const {
        return mError;
    }
    
    
    //SOLVER
public:
    virtual tBoolean solve(const MATH_Matrix& A,MATH_Vector& b,MATH_Vector& x,MATH_Vector& w)=0;

    inline tBoolean solve(const MATH_Matrix& A,MATH_Vector& b,MATH_Vector& x) {
        if (mW.get()==null) mW=b.NewInstance();
        return solve(A,b,x,*mW.get());
    }
    
    
  

    
};
#endif