---

Ma57Solver.h

/* OOQP                                                               *
 * Authors: E. Michael Gertz, Stephen J. Wright                       *
 * (C) 2001 University of Chicago. See Copyright Notification in OOQP */

#ifndef MA57LINSYS_H
#define MA57LINSYS_H

#include "DoubleLinearSolver.h"
#include "SparseSymMatrixHandle.h"
#include "SparseStorageHandle.h"
#include "OoqpVectorHandle.h"

extern "C" {
  void ma57id_( double cntl[],  int icntl[] );

  void ma57ad_( int * n,        int * ne,       int irn[],     
                int jcn[],      int * lkeep,    int keep[],
                int iwork[],    int icntl[],    int info[],
                double rinfo[] );

  void ma57bd_( int * n,        int * ne,       double a[],
                double fact[],  int * lfact,    int ifact[],
                int * lifact,   int * lkeep,    int keep[],
                int ppos[],     int * icntl,    double cntl[],
                int info[],     double rinfo[] );
  void ma57cd_( int * job,      int * n,        double fact[],
                int * lfact,    int ifact[],    int * lifact,
                int * nrhs,     double rhs[],   int * lrhs,
                double w[],     int * lw,       int iw1[],
                int icntl[],    int info[]);
  void ma57dd_( int * job,      int * n,        int * ne,
                double a[],     int irn[],      int jcn[],
                double fact[],  int * lfact,    int ifact[],
                int * lifact,   double rhs[],   double x[],
                double resid[], double w[],     int iw[],
                int icntl[],    double cntl[],  int info[],
                double rinfo[] );
  void ma57ed_( int * n,        int * ic,       int keep[],
                double fact[],  int * lfact,    double * newfac,
                int * lnew,     int  ifact[],   int * lifact,
                int newifc[],   int * linew,    int * info );
}

class Ma57Solver : public DoubleLinearSolver {
private:
  Ma57Solver() {};
protected:
  int     icntl[20];
  int     info[40];
  double  cntl[5];
  double  rinfo[20];

  double   kThresholdPivoting;

  double   kThresholdPivotingMax;

  double   kThresholdPivotingFactor;

  double   kTreatAsZero;

  double  kPrecision;

  int     *irowM,    *jcolM;

  double  *M;

  int      n;

  int      nnz;

  int     lkeep, *keep;

  int     lifact, *ifact, lfact;

  /* storage for the factors */
  double *fact;

  double  ipessimism, rpessimism;

  int     freshFactor;

  SparseStorageHandle mStorage;

  virtual void firstCall();
public:
  Ma57Solver( SparseSymMatrix * sgm );
  
  virtual void diagonalChanged( int idiag, int extent );
  virtual void matrixChanged();
  virtual void solve( OoqpVector& rhs );

  double setTreatAsZero() { return cntl[1] = kTreatAsZero; }

  double setThresholdPivoting() { return cntl[0] = kThresholdPivoting; }  

  virtual ~Ma57Solver();
};

#endif

---