---

Ma27Solver.h

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

#ifndef MA27LINSYS_H
#define MA27LINSYS_H

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



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

  void ma27ad_( int * n,        int * nz,       
                int irn[],      int icn[],      
                int iw[],       int * liw,
                int ikeep[],    int iw1[],      
                int * nsteps,   int * iflag,
                int icntl[],    double cntl[],
                int info[],     double * ops );

  void ma27bd_( int * n,        int * nz,
                int irn[],      int icn[],
                double a[],     int * la,
                int iw[],       int * liw,
                int ikeep[],    int * nsteps,
                int * maxfrt,   int iw1[],
                int icntl[],    double cntl[], 
                int  info[] );

  void ma27cd_( int * n,        
                double a[],     int * la,
                int iw[],       int * liw,
                double w[],     int * maxfrt,
                double rhs[],
                int iw2[],      int * nsteps,
                int icntl[],    int info[] );
}

class Ma27SolverBase : public DoubleLinearSolver {
protected:
  Ma27SolverBase() {};
protected:
  int     icntl[30];
  int     info[20];
  double  cntl[5];


  int ierror() { return info[1]; }
  int minimumRealWorkspace() { return info[4]; }
  int minimumIntWorkspace() { return info[5]; }
  int ma27ErrFlg() { return info[0]; }

  double  precision;

  int     *irowM,    *jcolM;

  double  *fact;

  int      n;

  int      nnz;

  int la;

  int     *ikeep, *iw, liw, *iw1, *iw2, nsteps, maxfrt;

  double  *w;

  double  ipessimism, rpessimism;

  virtual void firstCall();
public:

  Ma27SolverBase( int n, int nnz );
  
  double thresholdPivoting() { return cntl[0]; }
  void   setThresholdPivoting( double piv  ) { cntl[0] = piv; }

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

  virtual void copyMatrixElements( double fact[], int lfact ) = 0;

  virtual void getIndices( int irowM[], int jcolM[] ) = 0;

  virtual void diagonalChanged( int idiag, int extent );
  virtual void matrixChanged();

  virtual void basicSolve( double drhs[], int n );

  virtual void solve( OoqpVector& rhs ) = 0;

  virtual ~Ma27SolverBase();
};

class Ma27Solver : public Ma27SolverBase {
protected:
  Ma27Solver() {};

  SparseSymMatrixHandle mMat;

public:
  Ma27Solver( SparseSymMatrix * sgm );

  virtual void copyMatrixElements( double fact[], int lfact );

  virtual void getIndices( int irowM[], int jcolM[] );

  virtual void solve( OoqpVector& rhs );
};

#endif

---