SMLBlockIP.h

Go to the documentation of this file.

#ifndef SMLBlockIP_H
#define SMLBlockIP_H

#include <string>
#include "ExpandedModelInterface.h"
#include "BlockIP.h"
#include "MatrixBlockIP.h"

using namespace std;

class SMLBlockIP {
  private:
    struct objBlock {
      bool nonlinear;
      int numVars;
      int numQuadVars;
      double *linearCoefs;
      double *quadraticCoefs;
      ExpandedModelInterface *model;
    };

    struct objFunction {
      int numBlocks;
      bool slacks;
      int *slacksOrder;
      objBlock *blocks;
    };
    
    struct Order_vector {
      string* *v;
      bool operator() (int i,int j) {
        return (*v[i] < *v[j]);
      }
    };
  
    int numBlocks; 
    int numVars; 
    int numCons; 
    int numLinCons; 
    bool slacks; 
    bool nonLinear; 
    bool linear; 
    ExpandedModelInterface *root; 
    ExpandedModelInterface* slacksBlock; 
    objFunction *objF; 
    BlockIP *bip; 
    string *blockNames; 
    string *varNames; 
    int *varsOrder; 
    string *consNames; 
    
    MatrixBlockIP *blocks; 
    MatrixBlockIP *linking_constraints; 
    
    double inf; 
    
    static void fobjnonlin(int n, double x[], double &fx, double Gx[], double Hx[], void *params); 
    
  public:
    // Constructor
    SMLBlockIP();
    
    // Constructor
    SMLBlockIP(ExpandedModelInterface *root);
    
    // Destructor
    ~SMLBlockIP();
    
    // Load a problem from ampl files
    void readAmpl(const char *modelFilename, const char *dataFilename);
    
    // Converts a problem in ampl format to mps format
    void amplToMps(const char *modelFilename, const char *dataFilename, const char *mpsFilename, bool convertToStd=false);
    
    // Converts a problem in ampl format to BlockIP format
    void amplToBlockIPFormat(const char *modelFilename, const char *dataFilename, const char *BlockIPFormatFilename, bool convertToStd=false);
    
    // Get BlockIP interface
    BlockIP* getBlockIPInterface();
    
    // Print the whole problem through standard output
    void printProblem();
        
    // Return true if the objective function is non-linear
    bool isNonLinear();
    
    // Return true if the objective function is quadratic
    bool isQuadratic();
    
    // Return the correct order of the variables, ampl change the order
    const int* getVarsOrder();
    
  private:
    void initialize();
    
    // Load a problem from SML structure to BlockIP
    void loadProblemInBlockIP();
    
    void freeMemory();
};

inline bool SMLBlockIP::isNonLinear()

{
  return nonLinear;
}

inline bool SMLBlockIP::isQuadratic()

{
  return !linear && !nonLinear;
}

inline const int* SMLBlockIP::getVarsOrder()

{
  return varsOrder;
}

#endif //SMLBlockIP_H