---

NlpLinsysSol.hh

#ifndef NLPLINSYSSOL
#define NLPLINSYSSOL

#include "IotrRefCount.hh"
#include "IotrHandle.hh"
#include "IotrVector.hh"

class NlpIterate;
class NlpResiduals;
class NlpResidualShift;
class NlpLinsysDiag;
class IotrVector;
class NlpGrads; 

typedef Handle<IotrVector> IotrVectorHandle;

class NlpLinsysSol : public IotrRefCount {
protected:
  IotrVectorHandle mXsol;
  IotrVectorHandle mYsol;
  IotrVectorHandle mLambdaSol;
  IotrVectorHandle mPiSol;
  IotrVectorHandle mGammaSol;
  IotrVectorHandle mPhiSol;
public:
  IotrVector & x()      { return *mXsol; }
  IotrVector & y()      { return *mYsol; }
  IotrVector & lambda() { return *mLambdaSol; }
  IotrVector & pi()     { return *mPiSol; }
  IotrVector & gamma()  { return *mGammaSol; }
  IotrVector & phi()    { return *mPhiSol; }

  NlpLinsysSol( IotrVector * x,      IotrVector * y,
                IotrVector * lambda, IotrVector * pi,
                IotrVector * gamma,  IotrVector * phi );
  virtual ~NlpLinsysSol();
  virtual void calcRhs( NlpIterate & vars, 
                        NlpResiduals & resids, NlpResidualShift & shift,
                        NlpGrads & grads );
  // Scales multipliers by -1/(1+2sigma)
  virtual void getp( double sigma );
  virtual void scaleWithID( NlpLinsysDiag & diag );
  virtual void scaleWithIDinv( NlpLinsysDiag &diag );
  virtual void copyFrom( NlpLinsysSol & sol );
  virtual void copyInto( IotrVector & vec );
  virtual void copyFrom( IotrVector & vec );
  virtual double dot( NlpLinsysSol & sol );
  virtual void negateMultipliers();
  virtual void scale( double factor );
  virtual void axpy( double alpha, NlpLinsysSol & sol );

  double dualFeasibility() { return mXsol->infnorm(); };
  double penaltyNorm() { return mYsol->infnorm(); }

  virtual void inspect();
  virtual void inspectAsMatlab();
};

typedef Handle<NlpLinsysSol> NlpLinsysSolHandle;

#endif

---