---

NlpLinsysDiag.hh

#ifndef NLPLINSYSDIAG
#define NLPLINSYSDIAG

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

class NlpIterate;
class NlpResiduals;
class NlpResidualShift;
class NlpLinsysDiag;
class NlpLinsysSol;
class IotrVector;

typedef Handle<IotrVector> IotrVectorHandle;

class NlpLinsysDiag : public IotrRefCount {
protected:
  IotrVectorHandle mDx;
  IotrVectorHandle mDy;
  IotrVectorHandle mDgamma;
  IotrVectorHandle mDphi;
  IotrVectorHandle mDlambda;
  IotrVectorHandle mDpi;
public:
  IotrVector & Dx()         { return *mDx; }
  IotrVector & Dy()         { return *mDy; }
  IotrVector & Dgamma()  { return *mDgamma; }
  IotrVector & Dphi()    { return *mDphi; }
  IotrVector & Dlambda() { return *mDlambda; }
  IotrVector & Dpi()     { return *mDpi; }

  NlpLinsysDiag( IotrVector * x,      IotrVector * y,
                 IotrVector * lambda, IotrVector * pi,
                 IotrVector * gamma,  IotrVector * phi );

  virtual ~NlpLinsysDiag();
  virtual void calcDiagonal( NlpIterate & vars, NlpResiduals & resids,
                             NlpResidualShift & rshift );
  virtual void calcTrustDiagonal( NlpLinsysDiag & diag, double sigma );

  virtual void copyFrom( NlpLinsysDiag & diag );
  virtual void invertDual();
  virtual void squareRoot();
  virtual void copyInto( IotrVector & avec );
  virtual void inspect();
};

typedef Handle<NlpLinsysDiag> NlpLinsysDiagHandle;

#endif

---