Artifact 04698e671746a950997b9abf8c1d41acdf8e1a0870ed4924dfb4f4336b689042:

• File mttroot/mtt/lib/cc/mtt_AlgebraicSolver.cc — part of check-in [b94fab9d42] at 2001-08-08 02:15:00 on branch origin/master — Rationalisation of solver code, beginning with algebraic solvers. (user: geraint@users.sourceforge.net, size: 566) [annotate] [blame] [check-ins using] [more...]

#include "mtt_AlgebraicSolver.hh"

ColumnVector
MTT::AlgebraicSolver::solve (const ColumnVector	&x,
			     const ColumnVector	&u,
			     const double	&t,
			     const ColumnVector	&par)
{
  if (_nyz > 0)
    {
      _x = x;
      _uui.insert(u,0);
      _t = t;
      _par = par;
      _ui = ColumnVector(_nyz,1.0);
      Solve();
      _uui.insert(_ui,_nu);
    }
  else
    {
      _uui = u;
    }
  return _uui;
}

ColumnVector
MTT::AlgebraicSolver::eval (const ColumnVector	&ui)
{
  if (_nyz > 0)
    _uui.insert(ui,_nu);
  return mtt_ae(_x,_uui,_t,_par);
}