gurobi_fractional_scheduler.h

#ifndef __GUROBIFRACTIONALSCHEDULER_H__
#define __GUROBIFRACTIONALSCHEDULER_H__

#include <string>
#include <sstream>
#include <map>
#include <typeinfo>
#include "fractional_scheduler.h"
#include "../../../problems/scheduling/filtering.h"
#include "gurobi_model.h"
#include <iostream>
#include <algorithm>

//using std::cerr;
//using std::endl;

namespace delphos{ 

//#ifndef GRB_MAXIMIZE
//#define GRB_MAXIMIZE -1
//#define GRB_MINIMIZE 1
//#endif

class GurobiFractionalScheduler : public FractionalScheduler {
private:
    GurobiModel * _gmodel;
    
    
    bool _created_vars;
    bool _updated_ctrs;
    bool _stocks;
    bool _bmodelcreated;
    
    Model::OptimizationSense _optimization_sense;
    string _logfile;

    bool CreateVariables();
    bool CreateGoalFunction();
    bool CreateCoreConstraints();
    bool CreateUserConstraints();
    
    const bool CreatedVars() const;
    const bool UpdatedConstrs() const;
    
    void ApplyPrecedenceConstraint(const PrecedenceConstraint& pc);
    void ApplyAttributeSumConstraint(const AttributeSumConstraint& asc);
    void ApplyAttributeAvgConstraint(const AttributeAvgConstraint& asc);
    void ApplyDoNotMineBeforeConstraint(const DoNotMineBeforeConstraint& dc);
    void ApplyForceMiningByConstraint(const ForceMiningByConstraint& fc);
    void ApplyDoNotSendToConstraint(const DoNotSendToConstraint& dnsc);
    void ApplyGeometricConstraint(const GeometricConstraint& gc);
    void ApplyMinExtractPeriodConstraint(const MinExtractPeriodConstraint& mepc);
    void ApplyBlockSelectionPrecedenceConstraint(const BlockSelectionPrecedenceConstraint& bspc);
    void ApplyConstraint(const Constraint& ctr);
    
    vector<BlockIndexType> _index_in_bs;
    Variable ProcVar(
        const BlockIndexType& block_index, 
        const unsigned int& period, 
        const unsigned int& destination);
        
    Variable StockVar(
        const unsigned int& stock_index, 
        const unsigned int& period, 
        const unsigned int& destination);

    LinearExpr RealProcVar(
        const BlockIndexType& block_index, 
        const unsigned int& period, 
        const unsigned int& destination);
        
    LinearExpr RealProcVar2P(
        const BlockIndexType & block_index, 
        const unsigned int & period, 
        const unsigned int & destination);
    
    bool SetFractionalSchedule();
    unsigned long _nb_prec; 
        
public:
    GurobiFractionalScheduler(delphos::SchedulingInstance & si);

    virtual ~GurobiFractionalScheduler();
    
    virtual bool Run();

    bool CreateModel();

    #ifndef DOXYGEN

    void Optimize();
    #endif
    
    void WriteLP(const string outfile);
    
    void WriteSol(const string outfile);

    double GetSol(const BlockIndexType& block_index, 
        const unsigned int& period, 
        const unsigned int& destination);
        
    GurobiModel& GetMathematicalModel();

    double GetObjValue();
    
    void SetOptimizationSense(const Model::OptimizationSense& sense);

    void SetLogFile(const string& name);
    
    const bool IsStock();

    #ifndef DOXYGEN
    static string ProcVarName(
        const BlockIndexType& block_index, 
        const unsigned int& period, 
        const unsigned int& destination);
    static string ProcVarName2P(
        const BlockIndexType & block_index, 
        const unsigned int & period, 
        const unsigned int & destination);  
    static string StockVarName(
        const unsigned long& stock_index, 
        const unsigned int& period, 
        const unsigned int& destination);   
    #endif
    
    const double GetDual(const string & cname);

    void ComputeExpectedTime(const string column);

    const delphos::FractionalSchedule& GetFractionalSchedule();
    
};
}
#endif