rams 0.1.1 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 892ffe4b05f309fee78d179745713f1880362428
-  data.tar.gz: a659e1f65de80be3dae780ea30cb3b3ed504e5f9
+  metadata.gz: 8340e74253257a2e32d343f0c1c9a17e50ddeb47
+  data.tar.gz: 62eb080153e76cfa116a34c0afcb326e023484f6
 SHA512:
-  metadata.gz: dcd53a0bd69bf259195e74ad33162b4176a89d55388ca400516da6b43916d819d2a613ade587906d3b26bf07897b105d20f7098487e839bfbbb0f99dcc1e723b
-  data.tar.gz: 45cd3f47128cb9c4ee68ece0e9e309495b0b2049e53522aa3914aa404147164183ce67b0b0e4ceb49dbc4118f31174d192215ee7407f328a403edec651063fc3
+  metadata.gz: c4cf8db5cc61366fcf80cdabde0ba1f87d93a837f2cff11ca62de37088f4afb02d12df01d4fa2899d3cf35fbe18f9747da7b1ba0bdfea44e116ec319abc494b9
+  data.tar.gz: d7ea6f0faf629e241cb205fc6e28057b6639d51b1ce2bff4fc654594e0d84bce5fd87c7901c2ac57a6eb1dcc0c8275f239f1df61c03bec0a9d9c009db4d4c68d

data/lib/rams/model.rb

@@ -2,6 +2,7 @@ require 'tempfile'
 require_relative 'variable'
 require_relative 'solvers/cbc'
 require_relative 'solvers/clp'
+require_relative 'solvers/cplex'
 require_relative 'solvers/glpk'
 
 module RAMS
@@ -33,6 +34,7 @@ module RAMS
     SOLVERS = {
       cbc: RAMS::Solvers::CBC.new,
       clp: RAMS::Solvers::CLP.new,
+      cplex: RAMS::Solvers::CPLEX.new,
       glpk: RAMS::Solvers::GLPK.new
     }.freeze
 

data/lib/rams/solvers/cbc.rb

@@ -4,13 +4,13 @@ module RAMS
   module Solvers
     # Interface to COIN-OR Branch-and-Cut
     class CBC < Solver
-      def solver_command(model_file, solution_file, args)
-        ['cbc', model_file.path] + args + ['printingOptions', 'all', 'solve', 'solution', solution_file.path]
+      def solver_command(model_path, solution_path, args)
+        ['cbc', model_path] + args + ['printingOptions', 'all', 'solve', 'solution', solution_path]
       end
 
       private
 
-      def parse_status(model, lines)
+      def parse_status(_model, lines)
         return :undefined if lines.count < 1
         status = lines.first
         return :optimal if status =~ /Optimal/

data/lib/rams/solvers/clp.rb

@@ -4,13 +4,13 @@ module RAMS
   module Solvers
     # Interface to COIN-OR Linear Programming
     class CLP < Solver
-      def solver_command(model_file, solution_file, args)
-        ['clp', model_file.path] + args + ['printingOptions', 'all', 'solve', 'solution', solution_file.path]
+      def solver_command(model_path, solution_path, args)
+        ['clp', model_path] + args + ['printingOptions', 'all', 'solve', 'solution', solution_path]
       end
 
       private
 
-      def parse_status(model, lines)
+      def parse_status(_model, lines)
         return :undefined if lines.count < 1
         status = lines.first
         return :optimal if status =~ /optimal/

data/lib/rams/solvers/cplex.rb

@@ -0,0 +1,55 @@
+require 'nokogiri'
+require_relative 'solver'
+
+module RAMS
+  module Solvers
+    # Interface to CPLEX
+    class CPLEX < Solver
+      def solve_and_parse(model, model_path, solution_path)
+        call_solver model, model_path, solution_path
+        return RAMS::Solution.new(:infeasible, nil, {}, {}) unless File.exist? solution_path
+        parse_solution model, File.read(solution_path)
+      end
+
+      def solver_command(model_path, solution_path, args)
+        ['cplex', '-c', "read #{model_path}"] + args + ['optimize', "write #{solution_path}"]
+      end
+
+      private
+
+      def parse_solution(model, solution_text)
+        xml_doc = Nokogiri::XML solution_text
+        RAMS::Solution.new(
+          parse_status(model, xml_doc),
+          parse_objective(model, xml_doc),
+          parse_primal(model, xml_doc),
+          parse_dual(model, xml_doc)
+        )
+      end
+
+      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/
+        :unknown
+      end
+
+      def parse_objective(_model, xml_doc)
+        xml_doc.css('CPLEXSolution').css('header').first['objectiveValue'].to_f
+      end
+
+      def parse_primal(model, xml_doc)
+        xml_doc.css('CPLEXSolution').css('variables').css('variable').map do |v|
+          [model.variables[v['name']], v['value'].to_f]
+        end.to_h
+      end
+
+      def parse_dual(model, xml_doc)
+        xml_doc.css('CPLEXSolution').css('linearConstraints').css('constraint').map do |c|
+          [model.constraints[c['name']], c['dual'].to_f]
+        end.to_h
+      end
+    end
+  end
+end

data/lib/rams/solvers/glpk.rb

@@ -4,8 +4,8 @@ module RAMS
   module Solvers
     # Interface to the GNU Linear Programming Kit
     class GLPK < Solver
-      def solver_command(model_file, solution_file, args)
-        ['glpsol', '--lp', model_file.path, '--output', solution_file.path] + args
+      def solver_command(model_path, solution_path, args)
+        ['glpsol', '--lp', model_path, '--output', solution_path] + args
       end
 
       private

data/lib/rams/solvers/solver.rb

@@ -8,7 +8,7 @@ module RAMS
       def solve(model)
         model_file = write_model_file model
         begin
-          get_solution model, model_file
+          get_solution model, model_file.path
         ensure
           model_file.unlink
         end
@@ -23,23 +23,24 @@ module RAMS
         model_file
       end
 
-      def get_solution(model, model_file)
-        solution_file = Tempfile.new ['', '.sol']
+      def get_solution(model, model_path)
+        solution_path = model_path + '.sol'
         begin
-          solve_and_parse model, model_file, solution_file
+          solve_and_parse model, model_path, solution_path
         ensure
-          solution_file.unlink
+          File.delete(solution_path) if File.exist?(solution_path)
         end
       end
 
-      def solve_and_parse(model, model_file, solution_file)
-        call_solver model, model_file, solution_file
-        parse_solution model, File.read(solution_file)
+      def solve_and_parse(model, model_path, solution_path)
+        call_solver model, model_path, solution_path
+        return RAMS::Solution.new(:unknown, nil, {}, {}) unless File.exist? solution_path
+        parse_solution model, File.read(solution_path)
       end
 
       # rubocop:disable MethodLength
-      def call_solver(model, model_file, solution_file)
-        command = solver_command(model_file, solution_file, model.args)
+      def call_solver(model, model_path, solution_path)
+        command = solver_command(model_path, solution_path, model.args)
         _, stdout, stderr, exit_code = Open3.popen3(*command)
 
         begin
@@ -54,7 +55,7 @@ module RAMS
       end
       # rubocop:enable MethodLength
 
-      def solver_command(_model_file, _solution_file, _args)
+      def solver_command(_model_file, _solution_path, _args)
         raise NotImplementedError
       end
 

metadata

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