rams 0.1.2 → 0.1.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 8340e74253257a2e32d343f0c1c9a17e50ddeb47
-  data.tar.gz: 62eb080153e76cfa116a34c0afcb326e023484f6
+  metadata.gz: 6ef086c84ca5a0a24442527f9cfea81b147d521d
+  data.tar.gz: 36dc61c090667e648edc3c17cd1339093899cf92
 SHA512:
-  metadata.gz: c4cf8db5cc61366fcf80cdabde0ba1f87d93a837f2cff11ca62de37088f4afb02d12df01d4fa2899d3cf35fbe18f9747da7b1ba0bdfea44e116ec319abc494b9
-  data.tar.gz: d7ea6f0faf629e241cb205fc6e28057b6639d51b1ce2bff4fc654594e0d84bce5fd87c7901c2ac57a6eb1dcc0c8275f239f1df61c03bec0a9d9c009db4d4c68d
+  metadata.gz: fe7659717c3e04f74b0ab2e624e0b23cb83f8592d6cabfbf09162c87dafe57ebab179c1af2a1107497c78e6a3ed9ec5ea3b55132c26c33cc1151f7a46032b907
+  data.tar.gz: 008681815a50ec96d8d85120a2863977b535e338871b624d49b9c8fb800e9aedcd444dbc764d24de3a4047ae90ab153ea274fd55fe6fd019b37de37f95c4a63e

data/lib/rams/model.rb

@@ -4,6 +4,7 @@ require_relative 'solvers/cbc'
 require_relative 'solvers/clp'
 require_relative 'solvers/cplex'
 require_relative 'solvers/glpk'
+require_relative 'solvers/scip'
 
 module RAMS
   # A Model is a collection of:
@@ -35,7 +36,8 @@ module RAMS
       cbc: RAMS::Solvers::CBC.new,
       clp: RAMS::Solvers::CLP.new,
       cplex: RAMS::Solvers::CPLEX.new,
-      glpk: RAMS::Solvers::GLPK.new
+      glpk: RAMS::Solvers::GLPK.new,
+      scip: RAMS::Solvers::SCIP.new
     }.freeze
 
     def initialize

data/lib/rams/solvers/cbc.rb

@@ -13,10 +13,10 @@ module RAMS
       def parse_status(_model, lines)
         return :undefined if lines.count < 1
         status = lines.first
-        return :optimal if status =~ /Optimal/
-        return :feasible if status =~ /Stopped/
-        return :infeasible if status =~ /Infeasible/
-        return :unbounded if status =~ /Unbounded/
+        return :optimal if status =~ /optimal/i
+        return :feasible if status =~ /stopped/i
+        return :infeasible if status =~ /infeasible/i
+        return :unbounded if status =~ /unbounded/i
         :undefined
       end
 

data/lib/rams/solvers/clp.rb

@@ -13,10 +13,10 @@ module RAMS
       def parse_status(_model, lines)
         return :undefined if lines.count < 1
         status = lines.first
-        return :optimal if status =~ /optimal/
-        return :feasible if status =~ /stopped/
-        return :infeasible if status =~ /infeasible/
-        return :unbounded if status =~ /unbounded/
+        return :optimal if status =~ /optimal/i
+        return :feasible if status =~ /stopped/i
+        return :infeasible if status =~ /infeasible/i
+        return :unbounded if status =~ /unbounded/i
         :undefined
       end
 

data/lib/rams/solvers/cplex.rb

@@ -29,9 +29,9 @@ module RAMS
 
       def parse_status(_model, xml_doc)
         status = xml_doc.css('CPLEXSolution').css('header').first['solutionStatusString']
-        return :optimal if status =~ /optimal/
-        return :feasible if status =~ /feasible/
-        return :unbounded if status =~ /unbounded/
+        return :optimal if status =~ /optimal/i
+        return :feasible if status =~ /feasible/i
+        return :unbounded if status =~ /unbounded/i
         :unknown
       end
 

data/lib/rams/solvers/glpk.rb

@@ -12,10 +12,10 @@ module RAMS
 
       def parse_status(_model, lines)
         status = lines.select { |l| l =~ /^Status/ }.first
-        return :optimal if status =~ /OPTIMAL/
-        return :feasible if status =~ /FEASIBLE/
-        return :infeasible if status =~ /EMPTY/
-        return :unbounded if status =~ /UNBOUNDED/
+        return :optimal if status =~ /optimal/i
+        return :feasible if status =~ /feasibe/i
+        return :infeasible if status =~ /empty/i
+        return :unbounded if status =~ /unbounded/i
         :undefined
       end
 

data/lib/rams/solvers/scip.rb

@@ -0,0 +1,63 @@
+require 'open3'
+require_relative 'solver'
+
+module RAMS
+  module Solvers
+    # Interface to SCIP
+    class SCIP < Solver
+      def solve_and_parse(model, model_path, solution_path)
+        output = call_solver model, model_path, solution_path
+        parse_solution model, output
+      end
+
+      def solver_command(model_path, _solution_path, args)
+        ['scip', '-c', "read #{model_path}"] + args +
+          ['-c', 'optimize', '-c', 'display solution', '-c', 'display dualsolution', '-c', 'quit']
+      end
+
+      private
+
+      def parse_status(_model, lines)
+        return :undefined if lines.count < 1
+        status = lines.select { |l| l =~ /^SCIP Status/ }.first
+        return :optimal if status =~ /optimal/i
+        return :unknown if status =~ /(unknown|infeasible or unbounded)/i
+        return :infeasible if status =~ /infeasible/i
+        return :unbounded if status =~ /unbounded/i
+        :feasible
+      end
+
+      def parse_objective(_model, lines)
+        status = (lines.select { |l| l =~ /^objective value:/ }.first || '').split
+        return nil if status.size < 3
+        status[2].to_f
+      end
+
+      def parse_primal(model, lines)
+        primal = model.variables.values.map { |v| [v, 0.0] }.to_h
+        primal.merge(solution_lines(lines).map do |l|
+          comps = l.split
+          model.variables[comps[0]].nil? ? [] : [model.variables[comps.first], comps[1].to_f]
+        end.reject(&:empty?).to_h)
+      end
+
+      def parse_dual(model, lines)
+        dual = model.constraints.values.map { |c| [c, 0.0] }.to_h
+        dual.merge(solution_lines(lines).map do |l|
+          comps = l.split
+          model.constraints[comps[0]].nil? ? [] : [model.constraints[comps[0]], dual_value(model, comps[1].to_f)]
+        end.reject(&:empty?).to_h)
+      end
+
+      def solution_lines(lines)
+        start_idx = lines.index { |l| l =~ /^objective value:/ }
+        return [] unless start_idx
+        lines[start_idx, lines.count]
+      end
+
+      def dual_value(model, dual)
+        model.sense == :max ? -dual : dual
+      end
+    end
+  end
+end

data/lib/rams/solvers/solver.rb

@@ -48,6 +48,7 @@ module RAMS
           error = output + (stderr.gets(nil) || '')
           puts output if model.verbose && output != ''
           raise error unless exit_code.value == 0
+          return output
         ensure
           stdout.close
           stderr.close

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: rams
 version: !ruby/object:Gem::Version
-  version: 0.1.2
+  version: 0.1.3
 platform: ruby
 authors:
 - Ryan J. O'Neil
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2017-01-10 00:00:00.000000000 Z
+date: 2017-01-16 00:00:00.000000000 Z
 dependencies: []
 description: A library for solving MILPs in Ruby.
 email:
@@ -27,6 +27,7 @@ files:
 - lib/rams/solvers/clp.rb
 - lib/rams/solvers/cplex.rb
 - lib/rams/solvers/glpk.rb
+- lib/rams/solvers/scip.rb
 - lib/rams/solvers/solver.rb
 - lib/rams/variable.rb
 homepage: https://github.com/ryanjoneil/rams