genmodel 0.0.20 → 0.0.21

data/ext/Genmodel/InterfaceObject.cpp

@@ -0,0 +1,68 @@
+//
+//  InterfaceObject.cpp
+//  
+//
+//  Created by Mathieu Bouchard on 2014-03-19.
+//
+//
+
+#include "InterfaceObject.h"
+
+template<class T> InterfaceVector<T>::InterfaceVector()
+{
+    val = NULL;
+    size = 0;
+    //printf("\tC++: Adress at creation is %p size = %d\n", val, size);
+}
+
+template<class T> InterfaceVector<T>::InterfaceVector(int _size)
+{
+    size = _size;
+    val = new T[size];
+    //printf("\tC++: Adress at creation is %p size = %d\n", val, size);
+    //memset(val, 0, sizeof(T));
+}
+
+template<class T> InterfaceVector<T>::~InterfaceVector()
+{
+    Delete();
+}
+
+template<class T> void InterfaceVector<T>::SetSize(int _size)
+{
+    Delete();
+    size = _size;
+    val = new T[size];
+}
+template<class T> void InterfaceVector<T>::Delete()
+{
+    size = 0;
+    if(val != NULL)
+        delete[] val;
+}
+
+template<class T> void InterfaceVector<T>::Set(int index, T inval)
+{
+    if(index < 0 || index >= size)
+        printf("Index out of bound (%d not in 0..%d)\n", index, size);
+    else
+        val[index] = inval;
+}
+
+template<class T> T InterfaceVector<T>::Get(int index)
+{
+    if(index < 0 || index >= size)
+        printf("Index out of bound (%d not in 0..%d)\n", index, size);
+    else
+        return val[index];
+    return 0.0;
+}
+
+template<class T> T* InterfaceVector<T>::Ptr()
+{
+    return val;
+}
+
+template class InterfaceVector<int>;
+template class InterfaceVector<long>;
+template class InterfaceVector<double>;

data/ext/Genmodel/InterfaceObject.h

@@ -0,0 +1,28 @@
+//
+//  InterfaceObject.h
+//  
+//
+//  Created by Mathieu Bouchard on 2014-03-19.
+//
+//
+
+#ifndef _InterfaceObject_h
+#define _InterfaceObject_h
+
+template<class T> class InterfaceVector {
+public:
+    InterfaceVector();
+    InterfaceVector(int size);
+    ~InterfaceVector();
+    void SetSize(int size);
+    void Delete();
+    void Set(int index, T val);
+    T Get(int index);
+    T* Ptr();
+    
+private:
+    int size;
+    T * val;
+};
+
+#endif

data/ext/Genmodel/ProblemReaderOsi.cpp

@@ -0,0 +1,78 @@
+/*
+ * ProblemReaderOsi.cpp
+ *
+ *  Created on: 2012-10-02
+ *      Author: mbouchard
+ */
+
+#include "ProblemReaderOsi.h"
+
+int ReadFromFile(GenModel* pgm, string fn, int type, string dn)
+{
+    OsiSolverInterface* pm = new OsiClpSolverInterface();
+    //OsiSolverInterface* pm = new OsiGlpkSolverInterface();
+
+    switch (type)
+    {
+        case 0: pm->readMps(fn.c_str()); break;
+        case 1: pm->readLp(fn.c_str()); break;
+        case 2: pm->readGMPL(fn.c_str(), (dn == "" ? NULL : dn.c_str())); break;
+        default: pm->readMps(fn.c_str()); break;
+    }
+    
+    return ReadFromObject(pgm, pm);
+}
+
+int ReadFromObject(GenModel* pgm, OsiSolverInterface* pm)
+{
+    char tmp[64];
+    int numr = pm->getNumRows();
+    int numc = pm->getNumCols();
+    pgm->nz = 0;
+    
+    const double* obj = pm->getObjCoefficients();
+    const double* lb = pm->getColLower();
+    const double* ub = pm->getColUpper();
+    const char* type = pm->getColType();
+    for(int i = 0; i < numc; i++)
+    {
+        char gmtype;
+        snprintf(tmp,64,"C%d", i);
+        switch (type[i])
+        {
+            case 0: gmtype = 'C'; break;
+            case 1: gmtype = 'B'; break;
+            case 2: gmtype = 'I'; break;
+            default: gmtype = 'C'; break;
+        }
+        pgm->AddVar(string(tmp), obj[i],
+                    (lb[i]==COIN_DBL_MAX?numeric_limits<double>::infinity():(lb[i]==-COIN_DBL_MAX?-numeric_limits<double>::infinity():lb[i])),
+                    (ub[i]==COIN_DBL_MAX?numeric_limits<double>::infinity():(ub[i]==-COIN_DBL_MAX?-numeric_limits<double>::infinity():ub[i])), gmtype);
+    }
+    const char* sense = pm->getRowSense();
+    const double* lrhs = pm->getRowLower();
+    const double* urhs = pm->getRowUpper();
+    
+    const CoinPackedMatrix* pcpm = pm->getMatrixByRow();
+    const double* elem = pcpm->getElements();
+    const int* ind = pcpm->getIndices();
+    
+    for(int i = 0; i < numr; i++)
+    {
+        snprintf(tmp,64,"R%d", i);
+        if(sense[i] == 'L')
+            pgm->AddConst(string(tmp),
+                (urhs[i]==COIN_DBL_MAX?numeric_limits<double>::infinity():(urhs[i]==-COIN_DBL_MAX?-numeric_limits<double>::infinity():urhs[i])), sense[i]);
+        else
+            pgm->AddConst(string(tmp),
+                (lrhs[i]==COIN_DBL_MAX?numeric_limits<double>::infinity():(lrhs[i]==-COIN_DBL_MAX?-numeric_limits<double>::infinity():lrhs[i])), sense[i]);
+        pgm->consts.back().urhs = (urhs[i]==COIN_DBL_MAX?numeric_limits<double>::infinity():(urhs[i]==-COIN_DBL_MAX?-numeric_limits<double>::infinity():urhs[i]));
+        //if(fabs(lrhs[i]) > 1e-4 || fabs(urhs[i]) > 1e-4) //What is that ???
+        for(int j = pcpm->getVectorFirst(i); j < pcpm->getVectorLast(i); j++)
+        {
+            pgm->AddNzToLast(ind[j], elem[j]);
+        }
+    }
+    
+    return 0;
+}

data/ext/Genmodel/ProblemReaderOsi.h

@@ -0,0 +1,31 @@
+/*
+ * ProblemReaderOsi.h
+ *
+ *  Created on: 2012-10-02
+ *      Author: mbouchard
+ */
+
+#ifndef PROBLEMREADER_H_
+#define PROBLEMREADER_H_
+
+#ifndef CBC_THREAD
+    #define CBC_THREAD
+#endif
+
+#if defined WIN64 || defined WIN32
+	#ifndef snprintf
+		#define snprintf sprintf_s
+	#endif
+#endif
+
+#include "GenModel.h"
+#include "OsiClpSolverInterface.hpp"
+//#include "OsiGlpkSolverInterface.hpp"
+#include <string>
+
+using namespace std;
+
+int ReadFromFile(GenModel* pgm, string fn, int type=0, string dn = "");
+int ReadFromObject(GenModel* pgm, OsiSolverInterface* pm);
+
+#endif /* PROBLEMREADER_H_ */

data/ext/Genmodel/extconf.rb

@@ -0,0 +1,162 @@
+require 'mkmf'
+require 'rbconfig'
+
+extension_name = 'GenModel/GenModel'
+dir_config(extension_name)
+
+#abort "missing malloc()" unless have_func "malloc"
+#abort "missing free()"   unless have_func "free"
+
+parts = RUBY_DESCRIPTION.split(' ')
+type = parts[0]
+type = type[4..-1] if type.start_with?('tcs-')
+type = 'ree' if 'ruby' == type && RUBY_DESCRIPTION.include?('Ruby Enterprise Edition')
+is_windows = RbConfig::CONFIG['host_os'] =~ /(mingw|mswin)/
+is_linux = RbConfig::CONFIG['host_os'] =~ /(linux)/
+is_darwin = RbConfig::CONFIG['host_os'] =~ /(darwin??)/
+platform = RUBY_PLATFORM
+version = RUBY_VERSION.split('.')
+
+
+if is_darwin
+    puts "DEFINE Darwin"
+    $defs.push("-DDarwin")
+    elsif is_linux
+    puts "DEFINE Linux"
+    $defs.push("-DLinux")
+    elsif is_windows
+    puts "DEFINE WIN64"
+    $defs.push("-DWIN64")
+    else
+    puts "Unknown platform"
+    $defs.push("-DUnknownPlatform")
+end
+puts ">>>>> Creating Makefile for #{type} version #{RUBY_VERSION} on #{platform} (#{RbConfig::CONFIG['host_os']}) <<<<<"
+
+
+puts "#{CONFIG["dldflags"]}"
+CONFIG["DLDFLAGS"].gsub!(/-multiply_definedsuppress/, '')
+$DLDFLAGS = CONFIG["DLDFLAGS"]
+puts "#{CONFIG["DLDFLAGS"]}, #{$DLDFLAGS}"
+
+is_cplex = true;
+is_osi = true;
+if is_darwin
+    path = CONFIG['HOME']+"/Applications/IBM/ILOG/CPLEX_Studio126/cplex/include:/opt/ibm/ILOG/CPLEX_Studio126/cplex/include"
+    puts "Looking for ilcplex/cplex.h in "+path
+    if(is_cplex && find_header("ilcplex/cplex.h",path))
+        puts "found"
+    else
+        puts "not found"
+        is_cplex = false
+    end
+    
+    path = CONFIG['HOME']+"/Applications/IBM/ILOG/CPLEX_Studio126/cplex/lib/x86-64_osx/static_pic"
+    puts "Looking for libcplex (function CPXopenCPLEX) in "+path
+    if(is_cplex && find_header("libcplex","CPXopenCPLEX",path))
+        puts "found"
+    else
+        puts "not found"
+        is_cplex = false
+    end
+elsif is_linux
+    path = CONFIG['HOME']+"/ibm/ILOG/CPLEX_Studio126/cplex/include:/opt/ibm/ILOG/CPLEX_Studio126/cplex/include"
+    puts "Looking for ilcplex/cplex.h in "+path
+    if(is_cplex && find_header("ilcplex/cplex.h",path))
+        puts "found"
+    else
+        puts "not found"
+        is_cplex = false
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for pthread (function main) in "+path
+    if(find_header("pthread",nil,path))
+        puts "found"
+    else
+        puts "not found"
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for Clp (function main) in "+path
+    if(is_osi && find_header("Clp",nil,path))
+        puts "found"
+        else
+        puts "not found"
+        is_osi = false
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for Cbc (function main) in "+path
+    if(is_osi && find_header("Cbc",nil,path))
+        puts "found"
+        else
+        puts "not found"
+        is_osi = false
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for Cgl (function main) in "+path
+    if(is_osi && find_header("Cgl",nil,path))
+        puts "found"
+    else
+        puts "not found"
+        is_osi = false
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for Osi (function main) in "+path
+    if(is_osi && find_header("Osi",nil,path))
+        puts "found"
+    else
+        puts "not found"
+        is_osi = false
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for OsiClp (function main) in "+path
+    if(is_osi && find_header("Osi",nil,path))
+        puts "found"
+    else
+        puts "not found"
+        is_osi = false
+    end
+    
+    path = "/usr/lib:/usr/local/lib/"
+    puts "Looking for CoinUtils (function main) in "+path
+    if(is_cplex && find_header("CoinUtils",nil,path))
+        puts "found"
+    else
+        puts "not found"
+        is_cplex = false
+    end
+    
+    
+else
+    is_cplex = false;
+    is_osi = false;
+end
+
+[#r"/opt/ibm/ILOG/CPLEX_Studio126/cplex/bin/x86-64_sles10_4.1/",
+#os.environ['HOME']+"/source/gurobi510/linux64/lib",
+
+"Darwin": [os.environ['HOME']+"/Applications/IBM/ILOG/CPLEX_Studio126/cplex/lib/x86-64_osx/static_pic",
+r"/Library/gurobi560/mac64/lib",
+r"/Users/mbouchard/source/scipoptsuite-3.0.2/scip-3.0.2/lib/",
+r"../lib"]
+}
+
+
+if(is_cplex)
+    puts "*** With CPLEX_MOUDLE ***"
+    $defs.push("-DCPLEX_MODULE")
+end
+if(is_osi)
+    puts "*** With OSI_MOUDLE ***"
+    $defs.push("-DOSI_MODULE")
+end
+
+CONFIG["CFLAGS"] = '-O2 -Wall -pthread -fmessage-length=0 -fPIC -std=c++0x'
+CONFIG["CXX"] = 'g++'
+
+create_makefile (extension_name)

data/lib/Genmodel.rb

@@ -0,0 +1,5 @@
+module Genmodel
+    VERSION = "0.0.21"
+end
+
+require "Genmodel/Genmodel"

metadata

@@ -1,47 +1,59 @@
 --- !ruby/object:Gem::Specification
 name: genmodel
 version: !ruby/object:Gem::Version
-  version: 0.0.20
-  prerelease: 
+  version: 0.0.21
 platform: ruby
 authors:
 - Mathieu Bouchard
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2014-03-05 00:00:00.000000000 Z
+date: 2014-03-19 00:00:00.000000000 Z
 dependencies: []
 description: A generic tool to solve optimization model with different solvers
 email: mathbouchard@gmail.com
 executables: []
-extensions: []
+extensions:
+- ext/Genmodel/extconf.rb
 extra_rdoc_files: []
 files:
-- lib/genmodel.so
-- lib/genmodel.bundle
-homepage: https://github.com/mathbouchard/GenModel
+- ext/Genmodel/GenModel.h
+- ext/Genmodel/GenModelBase.cpp
+- ext/Genmodel/GenModelCplex.cpp
+- ext/Genmodel/GenModelCplex.h
+- ext/Genmodel/GenModelOsi.cpp
+- ext/Genmodel/GenModelOsi.h
+- ext/Genmodel/GenModelOsiInterface.cpp
+- ext/Genmodel/GenModelOsiInterface.h
+- ext/Genmodel/Genmodel.cpp
+- ext/Genmodel/InterfaceObject.cpp
+- ext/Genmodel/InterfaceObject.h
+- ext/Genmodel/ProblemReaderOsi.cpp
+- ext/Genmodel/ProblemReaderOsi.h
+- ext/Genmodel/extconf.rb
+- lib/Genmodel.rb
+homepage: https://github.com/mathbouchard/GenModelGem
 licenses:
 - MIT
+metadata: {}
 post_install_message: 
 rdoc_options: []
 require_paths:
 - lib
 required_ruby_version: !ruby/object:Gem::Requirement
-  none: false
   requirements:
-  - - ! '>='
+  - - '>='
     - !ruby/object:Gem::Version
       version: '0'
 required_rubygems_version: !ruby/object:Gem::Requirement
-  none: false
   requirements:
-  - - ! '>='
+  - - '>='
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 1.8.23
+rubygems_version: 2.2.2
 signing_key: 
-specification_version: 3
+specification_version: 4
 summary: GenModel Gem
 test_files: []