solversdp.h

Go to the documentation of this file.

#ifndef SOLVERSDP
#define SOLVERSDP
#include <map>
#include <vector>
//#include <borderbasix/v3sdp/mosek.h>


struct solversdp {
public:
  
    typedef std::pair<int,int>  Pair;
    typedef std::map<Pair,double> Matrix;
    typedef std::vector<double> Vector;
    typedef std::map<int, int>::iterator   dim_iterator;
    typedef std::map<Pair, Matrix>::iterator mat_iterator;
    typedef std::map<int, double>::iterator obj_iterator;

    solversdp(): slv_prec(32) {}

    int& dim()            {return slv_dim[0];}
    int& dim(int i)       {return slv_dim[i];}
    int  nb_block() const {return slv_dim.size();}

    dim_iterator dim_begin() { return slv_dim.begin(); }
    dim_iterator dim_end()   { return slv_dim.end(); }

    int& prec()       {return slv_prec;}
    int  prec() const {return slv_prec;}

    mat_iterator begin() { return slv_input.begin(); }
    mat_iterator end()   { return slv_input.end(); }

    double& coeff(int i, int j, int k)  { return slv_input[Pair(i,0)][Pair(j,k)]; }
    double& coeff(int i, int b, int j, int k)  { return slv_input[Pair(i,b)][Pair(j,k)]; }

    int& nentries(int i)       {return slv_nentries[i];}

    virtual void run() = 0;

    int     obj_size()    { return slv_objfct.size();}
    double& obj_coeff(int i)    { return slv_objfct[i]; }
    obj_iterator obj_begin()    { return slv_objfct.begin(); }
    obj_iterator obj_end()      { return slv_objfct.end(); }

    void clear() { slv_input.clear(); slv_objfct.clear(); slv_dim.clear();}

private:
    int slv_prec;
    std::map<int, int>     slv_dim;
    std::map<Pair, Matrix> slv_input;
    std::map<int, double>  slv_objfct;
    std::map<int, int>     slv_nentries;
public:
    Vector output;
    double fminprimal;
    double fmindual;
};

#endif // SOLVERSDP