rglpk 0.1.1 → 0.2.0

data/swig/glpk.i

@@ -23,7 +23,7 @@
 *Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */
 
-%module rglpk
+%module glpk_wrapper
 
 %{
 #include "glpk.h"
@@ -34,6 +34,7 @@
 %array_functions(int, intArray)
 %array_functions(double, doubleArray)
 
+%ignore glp_vprintf;
 %include "glpk.h"
 
 

data/test/helper.rb

@@ -0,0 +1,3 @@
+require 'test/unit'
+$LOAD_PATH << './ext'
+require 'lib/rglpk'

data/test/test_all.rb

@@ -1,188 +1,6 @@
-require 'test/unit'
-require 'glpk'
+here = File.expand_path(File.dirname(__FILE__))
 
-class TestGLPK < Test::Unit::TestCase
-
-  def test_create
-    assert_instance_of GLPK::Problem, GLPK::Problem.new
-  end
-
-  def test_name
-    p = GLPK::Problem.new
-    p.name = 'test'
-    assert_equal 'test', p.name
-  end
-
-  def test_obj_fun_name
-    p = GLPK::Problem.new
-    p.obj.name = 'test'
-    assert_equal 'test', p.obj.name
-  end
-
-  def test_obj_fun_dir
-    p = GLPK::Problem.new
-    p.obj.dir = GLPK::GLP_MIN
-    assert_equal GLPK::GLP_MIN, p.obj.dir
-    p.obj.dir = GLPK::GLP_MAX
-    assert_equal GLPK::GLP_MAX, p.obj.dir
-    assert_raise(ArgumentError){p.obj.dir = 3}    
-  end
-
-  def test_add_rows
-    p = GLPK::Problem.new
-    p.add_rows(2)
-    assert_equal 2, p.rows.length
-    p.add_rows(2)
-    assert_equal 4, p.rows.length
-  end
-
-  def test_add_cols
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    assert_equal 2, p.cols.length
-    p.add_cols(2)
-    assert_equal 4, p.cols.length
-  end
-
-  def test_set_row_name
-    p = GLPK::Problem.new
-    p.add_rows(10)
-    p.rows[1].name = 'test'
-    assert_equal 'test', p.rows[1].name
-    assert_nil p.rows[2].name
-  end
-
-  def test_set_col_name
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.cols[1].name = 'test'
-    assert_equal 'test', p.cols[1].name
-    assert_nil p.cols[2].name
-  end
-
-  def test_set_row_bounds
-    p = GLPK::Problem.new
-    p.add_rows(2)
-    p.rows[1].set_bounds(GLPK::GLP_FR,nil,nil)
-    assert_equal [GLPK::GLP_FR, nil, nil], p.rows[1].bounds
-  end
-
-  def test_set_col_bounds
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.cols[1].set_bounds(GLPK::GLP_FR,nil,nil)
-    assert_equal [GLPK::GLP_FR, nil, nil], p.cols[1].bounds
-  end
-
-  def test_obj_coef
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.obj.set_coef(1,2)
-    assert_equal [2,0], p.obj.coefs
-    p.obj.coefs = [1,2]
-    assert_equal [1,2], p.obj.coefs    
-  end
-
-  def test_set_row
-    p = GLPK::Problem.new
-    p.add_rows(2)
-    assert_raise(RuntimeError){p.rows[1].set([1,2])}    
-    p.add_cols(2)
-    p.rows[1].set([1,2])
-    assert_equal [1,2], p.rows[1].get
-  end
-
-  def test_set_col
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    assert_raise(RuntimeError){p.cols[1].set([1,2])}    
-    p.add_rows(2)
-    p.cols[1].set([1,2])
-    assert_equal [1,2], p.cols[1].get
-  end
-
-  def test_set_mat
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    assert_equal [1,2], p.rows[1].get
-    assert_equal [3,4], p.rows[2].get
-    assert_equal [1,3], p.cols[1].get
-    assert_equal [2,4], p.cols[2].get
-  end
-
-  def test_del_row
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    assert_equal [1,2], p.rows[1].get
-    p.del_rows([1])
-    assert_equal [3,4], p.rows[1].get
-    assert_equal [3],   p.cols[1].get
-  end
-
-  def test_del_col
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    assert_equal [1,3], p.cols[1].get
-    p.del_cols([1])
-    assert_equal [2,4], p.cols[1].get
-    assert_equal [2],   p.rows[1].get
-  end
-
-  def test_nz
-    p = GLPK::Problem.new
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    assert_equal 4, p.nz
-  end
-
-  def test_row_get_by_name
-    p = GLPK::Problem.new
-    assert_raises(RuntimeError){ p.rows['test'] }
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    assert_raises(ArgumentError){ p.rows['test'] }
-    p.rows[1].name = 'test'
-    assert_equal [1,2], p.rows['test'].get
-  end
-
-  def test_col_get_by_name
-    p = GLPK::Problem.new
-    assert_raises(RuntimeError){ p.cols['test'] }
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    assert_raises(ArgumentError){ p.cols['test'] }
-    p.cols[1].name = 'test'
-    assert_equal [1,3], p.cols['test'].get
-  end
-
-  def test_solve
-    p = GLPK::Problem.new
-    assert_raises(RuntimeError){ p.cols['test'] }
-    p.add_cols(2)
-    p.add_rows(2)
-    p.set_matrix([1,2,3,4])
-    p.simplex({:msg_lev => 1})
-  end
-
-  class D < GLPK::Problem
-    attr_accessor :species
-    def initialize
-      @species = []
-      super
-    end
-  end  
-
-  def test_derived
-    D.new.add_rows(10)
-  end
-
-end
+require here+'/helper'
+require here+'/test_basic'
+require here+'/test_problem_kind'
+require here+'/test_brief_example'

data/test/test_basic.rb

@@ -0,0 +1,186 @@
+require File.expand_path('helper', File.dirname(__FILE__))
+
+class TestRglpk < Test::Unit::TestCase
+
+  def test_create
+    assert_instance_of Rglpk::Problem, Rglpk::Problem.new
+  end
+
+  def test_name
+    p = Rglpk::Problem.new
+    p.name = 'test'
+    assert_equal 'test', p.name
+  end
+
+  def test_obj_fun_name
+    p = Rglpk::Problem.new
+    p.obj.name = 'test'
+    assert_equal 'test', p.obj.name
+  end
+
+  def test_obj_fun_dir
+    p = Rglpk::Problem.new
+    p.obj.dir = Rglpk::GLP_MIN
+    assert_equal Rglpk::GLP_MIN, p.obj.dir
+    p.obj.dir = Rglpk::GLP_MAX
+    assert_equal Rglpk::GLP_MAX, p.obj.dir
+    assert_raise(ArgumentError){p.obj.dir = 3}
+  end
+
+  def test_add_rows
+    p = Rglpk::Problem.new
+    p.add_rows(2)
+    assert_equal 2, p.rows.size
+    p.add_rows(2)
+    assert_equal 4, p.rows.size
+  end
+
+  def test_add_cols
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    assert_equal 2, p.cols.size
+    p.add_cols(2)
+    assert_equal 4, p.cols.size
+  end
+
+  def test_set_row_name
+    p = Rglpk::Problem.new
+    p.add_rows(10)
+    p.rows[1].name = 'test'
+    assert_equal 'test', p.rows[1].name
+    assert_nil p.rows[2].name
+  end
+
+  def test_set_col_name
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.cols[0].name = 'test'
+    assert_equal 'test', p.cols[0].name
+    assert_nil p.cols[1].name
+  end
+
+  def test_set_row_bounds
+    p = Rglpk::Problem.new
+    p.add_rows(2)
+    p.rows[1].set_bounds(Rglpk::GLP_FR, nil, nil)
+    assert_equal [Rglpk::GLP_FR, nil, nil], p.rows[1].bounds
+  end
+
+  def test_set_col_bounds
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.cols[1].set_bounds(Rglpk::GLP_FR, nil, nil)
+    assert_equal [Rglpk::GLP_FR, nil, nil], p.cols[1].bounds
+  end
+
+  def test_obj_coef
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.obj.set_coef(1, 2)
+    assert_equal [2, 0], p.obj.coefs
+    p.obj.coefs = [1, 2]
+    assert_equal [1, 2], p.obj.coefs
+  end
+
+  def test_set_row
+    p = Rglpk::Problem.new
+    p.add_rows(2)
+    assert_raise(RuntimeError){p.rows[1].set([1, 2])}
+    p.add_cols(2)
+    p.rows[1].set([1, 2])
+    assert_equal [1, 2], p.rows[1].get
+  end
+
+  def test_set_col
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    assert_raise(RuntimeError){p.cols[1].set([1, 2])}
+    p.add_rows(2)
+    p.cols[1].set([1, 2])
+    assert_equal [1, 2], p.cols[1].get
+  end
+
+  def test_set_mat
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    assert_equal [1, 2], p.rows[0].get
+    assert_equal [3, 4], p.rows[1].get
+    assert_equal [1, 3], p.cols[0].get
+    assert_equal [2, 4], p.cols[1].get
+  end
+
+  def test_del_row
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    assert_equal [1, 2], p.rows[0].get
+    p.del_rows([1])
+    assert_equal [3, 4], p.rows[0].get
+    assert_equal [3], p.cols[0].get
+  end
+
+  def test_del_col
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    assert_equal [1, 3], p.cols[0].get
+    p.del_cols([1])
+    assert_equal [2, 4], p.cols[0].get
+    assert_equal [2], p.rows[0].get
+  end
+
+  def test_nz
+    p = Rglpk::Problem.new
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    assert_equal 4, p.nz
+  end
+
+  def test_row_get_by_name
+    p = Rglpk::Problem.new
+    assert_raises(RuntimeError){ p.rows['test'] }
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    assert_raises(ArgumentError){ p.rows['test'] }
+    p.rows[0].name = 'test'
+    assert_equal [1, 2], p.rows['test'].get
+  end
+
+  def test_col_get_by_name
+    p = Rglpk::Problem.new
+    assert_raises(RuntimeError){ p.cols['test'] }
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    assert_raises(ArgumentError){ p.cols['test'] }
+    p.cols[0].name = 'test'
+    assert_equal [1, 3], p.cols['test'].get
+  end
+
+  def test_solve
+    p = Rglpk::Problem.new
+    assert_raises(RuntimeError){ p.cols['test'] }
+    p.add_cols(2)
+    p.add_rows(2)
+    p.set_matrix([1, 2, 3, 4])
+    p.simplex({:msg_lev => 1})
+  end
+
+  class D < Rglpk::Problem
+    attr_accessor :species
+    def initialize
+      @species = []
+      super
+    end
+  end
+
+  def test_derived
+    D.new.add_rows(10)
+  end
+end

data/test/test_brief_example.rb

@@ -0,0 +1,56 @@
+require File.expand_path('helper', File.dirname(__FILE__))
+
+class TestExample < Test::Unit::TestCase
+  def test_example
+    # The same Brief Example as found in section 1.3 of 
+    # glpk-4.44/doc/glpk.pdf.
+    #
+    # maximize
+    #   z = 10 * x1 + 6 * x2 + 4 * x3
+    #
+    # subject to
+    #   p:      x1 +     x2 +     x3 <= 100
+    #   q: 10 * x1 + 4 * x2 + 5 * x3 <= 600
+    #   r:  2 * x1 + 2 * x2 + 6 * x3 <= 300
+    #
+    # where all variables are non-negative
+    #   x1 >= 0, x2 >= 0, x3 >= 0
+    #    
+    p = Rglpk::Problem.new
+    p.name = "sample"
+    p.obj.dir = Rglpk::GLP_MAX
+    
+    rows = p.add_rows(3)
+    rows[0].name = "p"
+    rows[0].set_bounds(Rglpk::GLP_UP, 0, 100)
+    rows[1].name = "q"
+    rows[1].set_bounds(Rglpk::GLP_UP, 0, 600)
+    rows[2].name = "r"
+    rows[2].set_bounds(Rglpk::GLP_UP, 0, 300)
+    
+    cols = p.add_cols(3)
+    cols[0].name = "x1"
+    cols[0].set_bounds(Rglpk::GLP_LO, 0.0, 0.0)
+    cols[1].name = "x2"
+    cols[1].set_bounds(Rglpk::GLP_LO, 0.0, 0.0)
+    cols[2].name = "x3"
+    cols[2].set_bounds(Rglpk::GLP_LO, 0.0, 0.0)
+    
+    p.obj.coefs = [10, 6, 4]
+    
+    p.set_matrix([
+       1, 1, 1,
+      10, 4, 5,
+       2, 2, 6
+    ])
+    
+    p.simplex
+    z = p.obj.get
+    x1 = cols[0].get_prim
+    x2 = cols[1].get_prim
+    x3 = cols[2].get_prim
+    
+    result = "z = %g; x1 = %g; x2 = %g; x3 = %g" % [z, x1, x2, x3]
+    assert_equal "z = 733.333; x1 = 33.3333; x2 = 66.6667; x3 = 0", result
+  end
+end