---

SparseSymSolver.hh

Go to the documentation of this file.

#ifndef SPARSESYMSOLVER
#define SPARSESYMSOLVER

#include "LinearSymSolver.hh"
#include "SimpleArrayData.hh"
#include "SparseSymMatrix.hh"



class SparseSymSolver : public LinearSymSolver {
protected:
  int mFirstCall;
  int ictnl[30], info[20];
  double ctnl[5];
  SimpleArrayDataHandle mFactorization;
  SparseSymMatrixHandle  mMat;
  int mSize;
  int mN, mNz;
  int *irow, *jcol;
  int *iw, liw, *ikeep, *iw1;
  int nsteps, iflag, ops;
  int nrlnc, la, nirnec, maxfrt;
  double *Afact;
  int * ipiv;
  int mSingular;
public:
  SparseSymSolver( SparseSymMatrix * Mat );
  virtual void firstCall();
  ~SparseSymSolver();
  // Specializes LinearSolver::solve
  virtual void solve( IotrVector & b );
  // Specializes LinearSolver::matrixChanged
  virtual void matrixChanged();
  virtual void getIndices( int irow[], int jcol[] );
  virtual void copyMatrixElements( double afact[], int lafact ); 
  virtual void inspect();
  virtual void inertia( int & npos, int & nneg, int & nzero ) const;

  virtual int isSingular() const;

  double thresholdPivoting() { return ctnl[0]; }
  void   setThresholdPivoting( double piv  ) { ctnl[0] = piv; }

  double   treatAsZero() { return ctnl[2]; }
  void     setTreatAsZero( double tol ) { ctnl[2] = tol; }

};

typedef Handle<SparseSymSolver> SparseSymSolverHandle;

#endif

---