---

Solver.h

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

#ifndef ALGORITHM_H
#define ALGORITHM_H

class Data;
class Variables;
class Residuals;
class LinearSystem;
class Status;
class OoqpMonitor;
class OoqpStartStrategy;
class ProblemFormulation;

class Solver
{
protected:
  OoqpMonitor * itsMonitors;
  Status *status;
  OoqpStartStrategy * startStrategy;

  double dnorm;

  double     mutol, artol;

  double gamma_f, gamma_a;

  double phi;

  int        maxit;

  double *mu_history;

  double *rnorm_history;
  
  double *phi_history;

  double *phi_min_history;

public:

  int iter;

  Solver();
  virtual ~Solver();

  virtual void start( ProblemFormulation * formulation, 
                      Variables * iterate, Data * prob,
                      Residuals * resid, Variables * step);

  virtual void defaultStart( ProblemFormulation * formulation,
                             Variables * iterate, Data * prob,
                             Residuals * resid, Variables * step);

  virtual void stevestart ( ProblemFormulation * formulation,
                            Variables * iterate, Data * prob, 
                            Residuals * resid, Variables * step );

  virtual void dumbstart ( ProblemFormulation * formulation,
                           Variables * iterate, Data * prob, 
                           Residuals * resid, Variables * step );

  virtual int solve(Data *prob, Variables *iterate, Residuals * resids) = 0;

  virtual double finalStepLength( Variables *iterate, Variables *step );

  virtual void doMonitor( Data * data, Variables * vars,
                                                  Residuals * resids,
                                                  double alpha, double sigma,
                                                  int i, double mu,
                                                  int stop_code,
                                                  int level );
  virtual void defaultMonitor( Data * data, Variables * vars,
                                                           Residuals * resids,
                                                           double alpha, double sigma,
                                                           int i, double mu, 
                                                           int stop_code,
                                                           int level ) = 0;

  virtual int doStatus( Data * data, Variables * vars, 
                        Residuals * resids,
                        int i, double mu, 
                        int level );

  virtual int defaultStatus( Data * data, Variables * vars, 
                             Residuals * resids,
                             int i, double mu, 
                             int level );

  void addMonitor( OoqpMonitor * );

  void useStatus( Status * s ) { status = s; }

  void useStartStrategy( OoqpStartStrategy * s ) { startStrategy = s; }

  void monitorSelf();
 
  void setMuTol( double m ) { mutol = m; }
  double getMuTol() { return mutol; }

  void setArTol( double ar ) { artol = ar; }
  double getArTol() { return artol; }

  double dataNorm() { return dnorm; }
  LinearSystem * getLinearSystem() { return sys; };
protected:
  LinearSystem *sys;
};


#endif

---