---

MpsReader.h

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

#ifndef MPSREADER
#define MPSREADER

#include "OoqpVectorHandle.h"
#include "DoubleMatrixHandle.h"
#include "SimpleVectorHandle.h"
#include "stdio.h"
#include "hash.h"

struct MpsRowInfo;
struct MpsColInfo;

#ifdef TESTING
class MpsReaderTester;
#endif

class MpsReader {
#ifdef TESTING
  friend MpsReaderTester;
#endif
private:
  void insertElt( int irow[], int len, int jcol[],
                  double dval[], int& ne,
                  int row, int col, double val, int& ier );
  void stuffMatrix( GenMatrix& A,
                    int irow[], int nnz, int jcol[], double dA[] );
  void stuffMatrix( SymMatrix& A,
                    int irow[], int nnz, int jcol[], double dA[] );
protected:
  char * infilename;
  int iline;

  int nnzA, nnzC, nnzQ; 
  int my, mz;

  FILE * file;

  char    *boundType;

  MpsRowInfo * rowInfo;
  int * rowRemap;
  int totalRows;
  
  MpsColInfo * colInfo;
  int totalCols;
  int firstColumnLine;
  int columnFilePosition;

  char problemName[16];
  char objectiveName[16];
  char RHSName[16];
  //  char rangeName[10];
  char boundName[16];

  HashTable *rowTable;

  HashTable *colTable;

  double objminus;
  MpsReader() {};
  
  MpsReader( FILE * file );

  virtual void scanFile( int& iErr );

  virtual int GetLine_old(char * line );
  virtual int GetLine(char * line );

  virtual int ParseHeaderLine( char line[], char entry1[] );
  virtual int ParseHeaderLine2( char line[], char entry1[] );
  virtual int ParseBoundsLine( char line[], int& code, char name1[],
                               char name2[], double * val );
  virtual int ParseBoundsLine2( char line[], int& code, char name1[],
                               char name2[], double * val );
  virtual int ParseRowsLine( char line[], char code[], char name1[] );
  virtual int ParseRowsLine2( char line[], char code[], char name1[] );
  virtual void expectHeader( int kindOfLine, const char expectName[],
                            char line[], int& ierr );
  virtual void expectHeader2( int kindOfLine, const char expectName[],
                            char line[], int& ierr );
  virtual void remapRows();
  virtual int acceptHeader( int kindOfLine, const char expectName[],
                            char line[], int& ierr );
  virtual int acceptHeader2( int kindOfLine, const char expectName[],
                            char line[], int& ierr );
  virtual int ParseDataLine2( char line[], char code[], char name1[], 
                             char name2[], double * val1, int& hasSecondValue,
                             char name3[], double * val2);

  virtual int ParseDataLine( char line[], char code[], char name1[], 
                             char name2[], double * val1, int& hasSecondValue,
                             char name3[], double * val2);

  virtual int string_copy( char dest[], char string[], int max);
  
  virtual void readProblemName( char line[], int& iErr, int kindOfLine );
  virtual void readProblemName2( char line[], int& iErr, int kindOfLine );

  virtual void readObjectiveSense( char line[], int& iErr, int kindOfLine );
  virtual void readRowsSection( char line[62], 
                int& iErr, int& return_getline );
  virtual void scanColsSection( char line[62], 
                int& iErr, int& return_getline );
  virtual void scanRangesSection( char line[62], 
                  int& iErr, int& return_getline );
  virtual void rowHasRange( int rownum, double val, int& iErr );
  virtual void scanHessSection( char line[62], 
                int& iErr, int& return_getline );

  virtual void readColsSection( OoqpVector& c,
                                GenMatrix& A, GenMatrix& C,
                                char line[62], 
                                int& iErr, int& return_getline );
  virtual void readColsSection( double  c[],
                                int irowA[], int jcolA[], double dA[],
                                int irowC[], int jcolC[], double dC[],
                                char line[62], 
                                int& iErr, int& return_getline );
  virtual void readRHSSection( OoqpVector&  b,
                               SimpleVector& clow, OoqpVector& iclow,
                               SimpleVector& cupp, OoqpVector& icupp,
                               char line[], int& ierr, int& kindOfLine );
  virtual void readRHSSection( double b[],
                               double clow[], char iclow[],
                               double cupp[], char icupp[],
                               char line[], int& ierr, int& kindOfLine );
  virtual void readRangesSection( SimpleVector& clow, SimpleVector& cupp,
                                  char line[], int& ierr, int& kindOfLine );
  virtual void readRangesSection( double clow[], double cupp[],
                                  char line[], int& ierr, int& kindOfLine );
  virtual void readBoundsSection( OoqpVector& xlow, OoqpVector& ixlow,
                                  OoqpVector& xupp, OoqpVector& ixupp,
                                  char line[], int& ierr, int& kindOfLine );
  virtual void defaultBounds( double xlow[], char ixlow[],
                              double xupp[], char ixupp[] );
  virtual void defaultBounds( OoqpVector& xlow, OoqpVector& ixlow,
                              OoqpVector& xupp, OoqpVector& ixupp );
  virtual void readBoundsSection( double xlow[], char ixlow[],
                                  double xupp[], char ixupp[],
                                  char line[], int& ierr, int& kindOfLine );
  virtual void readHessSection( SymMatrix& Q,
                                char line[], int& ierr, int& kindOfLine );
  virtual void readHessSection( int irowQ[], int jcolQ[], double dQ[],
                                char line[], int& ierr, int& kindOfLine );

public:
  int scalingOption; 

  char objectiveSense[3]; /* MAX or MIN */

   
  static MpsReader * newReadingFile( char filename[], int& iErr );

  static  void findFile( FILE*& file, char *& resolvedName,
                         char filename[] );
  void numberOfNonZeros( int& nnzQ, int& nnzA, int& nnzC );
  void numbersOfNonZeros( int nnzQ[], int nnzA[], int nnzC[] );
  virtual void readQpBound( OoqpVector& c, SymMatrix& Q,
                            OoqpVector& xlow, OoqpVector& ixlow,
                            OoqpVector& xupp, OoqpVector& ixupp,
                            int &ierr );
  virtual void readQpGen( OoqpVector& c, SymMatrix& Q,
                          OoqpVector& xlow, OoqpVector& ixlow,
                          OoqpVector& xupp, OoqpVector& ixupp,
                          GenMatrix& A, OoqpVector& b,
                          GenMatrix& C, 
                          OoqpVector& clow, OoqpVector& iclow,
                          OoqpVector& cupp, OoqpVector& icupp,
                          int& ierr );

  virtual void readQpGen( double     c[],  
                          int    irowQ[],  int  jcolQ[],  double dQ[],
                          double  xlow[],  char ixlow[],
                          double  xupp[],  char ixupp[],
                          int    irowA[],  int  jcolA[],  double dA[],
                          double     b[],
                          int    irowC[],  int  jcolC[],  double dC[],
                          double  clow[],  char  iclow[],
                          double  cupp[],  char  icupp[],
                          int&    ierr );
  virtual void getSizes( int& nx, int& my, int& mz );

  virtual void releaseFile( int& ierr );
  
  virtual ~MpsReader();

  // Output section

  double objconst() { return -objminus; }

  void printSolution( double x[],      int     nx,
                      double xlow[],   char    ixlow[],
                      double xupp[],   char    ixupp[],
                      double gamma[],  double  phi[],
                      double y[],      int     my,
                      double s[],      int     mz,
                      double clow[],   char    iclow[],
                      double cupp[],   char    icupp[],
                      double lambda[], double  pi[],
                      double z[],      double  objectiveValue,
                      int& iErr );
  char * defaultOutputFilename( int& iErr );

};

enum { mpsok = 0,
       mpsunknownerr    = -1, // Someone didn't set the error code properly.
       mpsfileopenerr   = -2, // Couldn't open the given file
       mpsioerr         = -3, // There was an i/o error.
       mpssyntaxerr     = -4, // There was a mps syntax error.
       mpsmemoryerr     = -5  // We couldn't allocate the necessary memory
};

#endif

---