simplex.h

Go to the documentation of this file.

#ifndef SIMPLEXH
#define SIMPLEXH

#include "box.h"
#include "interval.h"
#include "linear_constraint.h"
#include "equation.h"

#if (_SIMPLEX_ENGINE==_GLPK)
  #include <glpk.h>
  #define SimplexType LPX
#endif

#if (_SIMPLEX_ENGINE==_CLP)
  #include <coin/Clp_C_Interface.h>  
  #define SimplexType __attribute__ ((aligned(8))) Clp_Simplex
#endif

#if (_SIMPLEX_ENGINE==_LP_SOLVE)
  #include <lp_solve/lp_lib.h>
  #define SimplexType lprec
#endif

/************************************************************************************/
#define SIMPLEX_TIMEOUT  10
/************************************************************************************/

typedef struct {

  SimplexType *lp; 
  #if (_SIMPLEX_ENGINE==_CLP)
    unsigned int fakeRows; 
    double *lower; 
    double *upper; 
    double *obj; 
  #endif

  #if (_SIMPLEX_ENGINE==_LP_SOLVE)
    double inf; 
  #endif

} TSimplex;

/* Functions implemented in particular for each simplex engine
   (GLPK,CLP,...) in simplex_glpk.c simplex_clp.c,.... */

void SimplexCreate(double epsilon,unsigned int ncols,TSimplex *s);


void ResetSimplex(TSimplex *s);

unsigned int SimplexNColumns(TSimplex *s);


unsigned int SimplexNRows(TSimplex *s);

void SimplexSetColBounds(unsigned int ncol,Tinterval *i,TSimplex *s);

void SimplexGetColBounds(unsigned int ncol,Tinterval *i,TSimplex *s);

void SimplexGetColConstraint(unsigned int ncol,TLinearConstraint *lc,TSimplex *s);

boolean SimplexColEmpty(unsigned int ncol,TSimplex *s);

double SimplexGetColPrimal(unsigned int ncol,TSimplex *s);


double SimplexGetColDual(unsigned int ncol,TSimplex *s);

void SimplexSetRowBounds(unsigned int nrow,Tinterval *i,TSimplex *s);

void SimplexGetRowBounds(unsigned int nrow,Tinterval *i,TSimplex *s);

void SimplexGetRowConstraint(unsigned int nrow,TLinearConstraint *lc,TSimplex *s);

double SimplexGetRowPrimal(unsigned int nrow,TSimplex *s);

double SimplexGetRowDual(unsigned int nrow,TSimplex *s);

void SimplexAddNewConstraintRaw(TLinearConstraint *lc,
                                TSimplex *s);

void SimplexSetOptimizationFunction(TLinearConstraint *obj,TSimplex *s);

void SimplexGetOptimizationFunction(TLinearConstraint *obj,TSimplex *s);

double SimplexGetOptimalValueRaw(TSimplex *s);

unsigned int SimplexSolve(TSimplex *s);


void SimplexDelete(TSimplex *s);


void SimplexExpandBounds(unsigned int eq_type,Tinterval *b);


void SetSimplexBounds(Tbox *b,TSimplex *lp);

boolean SimplexAddNewConstraint(double epsilon,
                                unsigned int safeSimplex,
                                TLinearConstraint *lc,
                                unsigned int eq_type,
                                Tinterval *is,
                                TSimplex *s);

double SimplexGetOptimalValue(unsigned int safeSimplex,double m,Tbox *b,TSimplex *s);

unsigned int ReduceRange(double epsilon,unsigned int safeSimplex,unsigned int nv,Tbox *b,TSimplex *lp);
#endif