Search  

---

C:/temp/src/j2k/Beta/Math/Vector/NotDone/Gauss/Eqsystem.hpp

Go to the documentation of this file.

//  EqSystem.h
//  A system of equations and unknowns.

#include "Vector.h" // Added by ClassView
#include "Matrix.h" // Added by ClassView

#include <vector>

// Represent an unknown.
// The solver gets and sets its values.
class IEqUnknown
{
public:
    virtual double GetValue() = 0;
    virtual void SetValue(double dValue) = 0;
};

// Represents an equation.
// The solver calculates equations and
// queries its dependencies.
class IEqEquation
{
public:
    virtual double CalculateValue() = 0;
    virtual bool DependsUpon(IEqUnknown *pUnknown) = 0;
};

// The solver engine.
// The application implements the above interfaces
// and adds them to the system object, then kicks
// off the solution.
class CEqSystem  
{
public:
    CEqSystem();

    void AddUnknown(IEqUnknown *pUnknown);
    void AddEquation(IEqEquation *pEquation);
    bool Solve();

private:
    void CalcNorm( CVector &vY );
    bool Converged( const CVector &vY );
    CMatrix CalcJacobian(
        const CVector &vX,
        const CVector &vY );
    bool LineSearch(
        const CVector &vX,
        CVector &vXNew,
        CVector &vY,
        const CVector &vStep,
        const CVector &vGradF);

private:
    std::vector<IEqUnknown *>   m_vUnknowns;
    std::vector<IEqEquation *>  m_vEquations;

    double                      m_dMaxLength;
    double                      m_dNorm;
};

---