gsl4r 0.0.2 → 0.0.3

data/test/gsl_vector_tests_gen.c

@@ -0,0 +1,203 @@
+/* Auto generated by GSL4r::Vector::Harness */
+#include <stdio.h>
+#include "gsl/gsl_vector.h"
+int main( int argc, char **argv )
+{
+  puts("$: << File.join('..','lib')\nrequire 'test/unit'\nrequire 'test/unit/autorunner'\nrequire 'gsl4r/vector'\ninclude GSL4r::Vector\nclass VectorTests < Test::Unit::TestCase\n  EPSILON = 5.0e-15");
+ /* c_test_gsl_vector_sub */
+  gsl_vector * v76;
+  v76 = gsl_vector_alloc(3); 
+  gsl_vector_set(v76, 0, 1.0 ); gsl_vector_set(v76, 1, 2.0 ); gsl_vector_set(v76, 2, 3.0 );
+  gsl_vector * v77;
+  v77 = gsl_vector_alloc(3); 
+  gsl_vector_set(v77, 0, 1.0 ); gsl_vector_set(v77, 1, 2.0 ); gsl_vector_set(v77, 2, 3.0 );
+  gsl_vector_sub(v76,v77);
+  puts("def test_gsl_vector_sub()");
+  puts("  v76 = GSL_Vector.create(3)");
+  puts("  v76.set( [1.0,2.0,3.0] )");
+  puts("  v77 = GSL_Vector.create(3)");
+  puts("  v77.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_sub(v76,v77)");
+  printf("  assert v76.values == ");
+  printf("[%.15g,%.15g,%.15g]\n",gsl_vector_get(v76,0),gsl_vector_get(v76,1),gsl_vector_get(v76,2));
+  puts("end");
+ /* c_test_gsl_vector_scale */
+  gsl_vector * v82;
+  v82 = gsl_vector_alloc(3); 
+  gsl_vector_set(v82, 0, 1.0 ); gsl_vector_set(v82, 1, 2.0 ); gsl_vector_set(v82, 2, 3.0 );
+double v83 = (double)2.0;
+  gsl_vector_scale(v82,v83);
+  puts("def test_gsl_vector_scale()");
+  puts("  v82 = GSL_Vector.create(3)");
+  puts("  v82.set( [1.0,2.0,3.0] )");
+  puts("");
+  puts("  v83 = 2.0");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_scale(v82,v83)");
+  printf("  assert v82.values == ");
+  printf("[%.15g,%.15g,%.15g]\n",gsl_vector_get(v82,0),gsl_vector_get(v82,1),gsl_vector_get(v82,2));
+  puts("end");
+ /* c_test_gsl_vector_max_index */
+  gsl_vector * v87;
+  v87 = gsl_vector_alloc(3); 
+  gsl_vector_set(v87, 0, 1.0 ); gsl_vector_set(v87, 1, 2.0 ); gsl_vector_set(v87, 2, 3.0 );
+  size_t c_r86;
+  c_r86 = gsl_vector_max_index(v87);
+  puts("def test_gsl_vector_max_index()");
+  puts("  v87 = GSL_Vector.create(3)");
+  puts("  v87.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_max_index(v87)");
+  printf("  assert_in_delta r_r1, %lu, EPSILON\n", c_r86);
+  puts("end");
+ /* c_test_gsl_vector_isneg */
+  gsl_vector * v92;
+  v92 = gsl_vector_alloc(3); 
+  gsl_vector_set(v92, 0, 1.0 ); gsl_vector_set(v92, 1, 2.0 ); gsl_vector_set(v92, 2, 3.0 );
+  int c_r91;
+  c_r91 = gsl_vector_isneg(v92);
+  puts("def test_gsl_vector_isneg()");
+  puts("  v92 = GSL_Vector.create(3)");
+  puts("  v92.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_isneg(v92)");
+  printf("  assert_in_delta r_r1, %lu, EPSILON\n", c_r91);
+  puts("end");
+ /* c_test_gsl_vector_mul */
+  gsl_vector * v78;
+  v78 = gsl_vector_alloc(3); 
+  gsl_vector_set(v78, 0, 1.0 ); gsl_vector_set(v78, 1, 2.0 ); gsl_vector_set(v78, 2, 3.0 );
+  gsl_vector * v79;
+  v79 = gsl_vector_alloc(3); 
+  gsl_vector_set(v79, 0, 1.0 ); gsl_vector_set(v79, 1, 2.0 ); gsl_vector_set(v79, 2, 3.0 );
+  gsl_vector_mul(v78,v79);
+  puts("def test_gsl_vector_mul()");
+  puts("  v78 = GSL_Vector.create(3)");
+  puts("  v78.set( [1.0,2.0,3.0] )");
+  puts("  v79 = GSL_Vector.create(3)");
+  puts("  v79.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_mul(v78,v79)");
+  printf("  assert v78.values == ");
+  printf("[%.15g,%.15g,%.15g]\n",gsl_vector_get(v78,0),gsl_vector_get(v78,1),gsl_vector_get(v78,2));
+  puts("end");
+ /* c_test_gsl_vector_add_constant */
+  gsl_vector * v84;
+  v84 = gsl_vector_alloc(3); 
+  gsl_vector_set(v84, 0, 1.0 ); gsl_vector_set(v84, 1, 2.0 ); gsl_vector_set(v84, 2, 3.0 );
+double v85 = (double)2.0;
+  gsl_vector_add_constant(v84,v85);
+  puts("def test_gsl_vector_add_constant()");
+  puts("  v84 = GSL_Vector.create(3)");
+  puts("  v84.set( [1.0,2.0,3.0] )");
+  puts("");
+  puts("  v85 = 2.0");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_add_constant(v84,v85)");
+  printf("  assert v84.values == ");
+  printf("[%.15g,%.15g,%.15g]\n",gsl_vector_get(v84,0),gsl_vector_get(v84,1),gsl_vector_get(v84,2));
+  puts("end");
+ /* c_test_gsl_vector_isnull */
+  gsl_vector * v90;
+  v90 = gsl_vector_alloc(3); 
+  gsl_vector_set(v90, 0, 1.0 ); gsl_vector_set(v90, 1, 2.0 ); gsl_vector_set(v90, 2, 3.0 );
+  int c_r89;
+  c_r89 = gsl_vector_isnull(v90);
+  puts("def test_gsl_vector_isnull()");
+  puts("  v90 = GSL_Vector.create(3)");
+  puts("  v90.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_isnull(v90)");
+  printf("  assert_in_delta r_r1, %lu, EPSILON\n", c_r89);
+  puts("end");
+ /* c_test_gsl_vector_min */
+  gsl_vector * v88;
+  v88 = gsl_vector_alloc(3); 
+  gsl_vector_set(v88, 0, 1.0 ); gsl_vector_set(v88, 1, 2.0 ); gsl_vector_set(v88, 2, 3.0 );
+  double c_r87;
+  c_r87 = gsl_vector_min(v88);
+  puts("def test_gsl_vector_min()");
+  puts("  v88 = GSL_Vector.create(3)");
+  puts("  v88.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_min(v88)");
+  printf("  assert_in_delta r_r1, %.15g, EPSILON\n", c_r87);
+  puts("end");
+ /* c_test_gsl_vector_isnonneg */
+  gsl_vector * v93;
+  v93 = gsl_vector_alloc(3); 
+  gsl_vector_set(v93, 0, 1.0 ); gsl_vector_set(v93, 1, 2.0 ); gsl_vector_set(v93, 2, 3.0 );
+  int c_r92;
+  c_r92 = gsl_vector_isnonneg(v93);
+  puts("def test_gsl_vector_isnonneg()");
+  puts("  v93 = GSL_Vector.create(3)");
+  puts("  v93.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_isnonneg(v93)");
+  printf("  assert_in_delta r_r1, %lu, EPSILON\n", c_r92);
+  puts("end");
+ /* c_test_gsl_vector_add */
+  gsl_vector * v74;
+  v74 = gsl_vector_alloc(3); 
+  gsl_vector_set(v74, 0, 1.0 ); gsl_vector_set(v74, 1, 2.0 ); gsl_vector_set(v74, 2, 3.0 );
+  gsl_vector * v75;
+  v75 = gsl_vector_alloc(3); 
+  gsl_vector_set(v75, 0, 1.0 ); gsl_vector_set(v75, 1, 2.0 ); gsl_vector_set(v75, 2, 3.0 );
+  gsl_vector_add(v74,v75);
+  puts("def test_gsl_vector_add()");
+  puts("  v74 = GSL_Vector.create(3)");
+  puts("  v74.set( [1.0,2.0,3.0] )");
+  puts("  v75 = GSL_Vector.create(3)");
+  puts("  v75.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_add(v74,v75)");
+  printf("  assert v74.values == ");
+  printf("[%.15g,%.15g,%.15g]\n",gsl_vector_get(v74,0),gsl_vector_get(v74,1),gsl_vector_get(v74,2));
+  puts("end");
+ /* c_test_gsl_vector_div */
+  gsl_vector * v80;
+  v80 = gsl_vector_alloc(3); 
+  gsl_vector_set(v80, 0, 1.0 ); gsl_vector_set(v80, 1, 2.0 ); gsl_vector_set(v80, 2, 3.0 );
+  gsl_vector * v81;
+  v81 = gsl_vector_alloc(3); 
+  gsl_vector_set(v81, 0, 1.0 ); gsl_vector_set(v81, 1, 2.0 ); gsl_vector_set(v81, 2, 3.0 );
+  gsl_vector_div(v80,v81);
+  puts("def test_gsl_vector_div()");
+  puts("  v80 = GSL_Vector.create(3)");
+  puts("  v80.set( [1.0,2.0,3.0] )");
+  puts("  v81 = GSL_Vector.create(3)");
+  puts("  v81.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_div(v80,v81)");
+  printf("  assert v80.values == ");
+  printf("[%.15g,%.15g,%.15g]\n",gsl_vector_get(v80,0),gsl_vector_get(v80,1),gsl_vector_get(v80,2));
+  puts("end");
+ /* c_test_gsl_vector_max */
+  gsl_vector * v86;
+  v86 = gsl_vector_alloc(3); 
+  gsl_vector_set(v86, 0, 1.0 ); gsl_vector_set(v86, 1, 2.0 ); gsl_vector_set(v86, 2, 3.0 );
+  double c_r85;
+  c_r85 = gsl_vector_max(v86);
+  puts("def test_gsl_vector_max()");
+  puts("  v86 = GSL_Vector.create(3)");
+  puts("  v86.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_max(v86)");
+  printf("  assert_in_delta r_r1, %.15g, EPSILON\n", c_r85);
+  puts("end");
+ /* c_test_gsl_vector_ispos */
+  gsl_vector * v91;
+  v91 = gsl_vector_alloc(3); 
+  gsl_vector_set(v91, 0, 1.0 ); gsl_vector_set(v91, 1, 2.0 ); gsl_vector_set(v91, 2, 3.0 );
+  int c_r90;
+  c_r90 = gsl_vector_ispos(v91);
+  puts("def test_gsl_vector_ispos()");
+  puts("  v91 = GSL_Vector.create(3)");
+  puts("  v91.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_ispos(v91)");
+  printf("  assert_in_delta r_r1, %lu, EPSILON\n", c_r90);
+  puts("end");
+ /* c_test_gsl_vector_min_index */
+  gsl_vector * v89;
+  v89 = gsl_vector_alloc(3); 
+  gsl_vector_set(v89, 0, 1.0 ); gsl_vector_set(v89, 1, 2.0 ); gsl_vector_set(v89, 2, 3.0 );
+  size_t c_r88;
+  c_r88 = gsl_vector_min_index(v89);
+  puts("def test_gsl_vector_min_index()");
+  puts("  v89 = GSL_Vector.create(3)");
+  puts("  v89.set( [1.0,2.0,3.0] )");
+  puts("  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_min_index(v89)");
+  printf("  assert_in_delta r_r1, %lu, EPSILON\n", c_r88);
+  puts("end");
+  puts("end");
+  return(0);
+}

data/test/vector_class_test.rb

@@ -0,0 +1,128 @@
+$: << File.join('..','lib')
+
+require 'rubygems'
+require 'ffi'
+
+require 'test/unit'
+require 'test/unit/autorunner'
+
+require 'gsl4r'
+require 'gsl4r/vector'
+
+include GSL4r::Vector
+include FFI
+
+class VectorTests < Test::Unit::TestCase
+
+  SIZE=50
+
+  def test_gsl_vector_alloc_and_free()
+    assert_nothing_raised do
+      vecptr = MemoryPointer.new :pointer
+      vecptr = GSL4r::Vector::Methods::gsl_vector_alloc( SIZE )
+      GSL4r::Vector::Methods::gsl_vector_free( vecptr )
+    end
+  end
+
+  # freed by leaving scope
+  def test_gsl_vector_calloc()
+    assert_nothing_raised do
+      vecptr = MemoryPointer.new :pointer
+      vecptr = GSL4r::Vector::Methods::gsl_vector_calloc( SIZE )
+      vec = GSL_Vector.new( vecptr )
+
+      v = vec.values
+      assert v.length == SIZE
+      v.each { |i| assert i == 0 }
+
+    end
+  end
+
+  def test_gsl_vector_new()
+    assert_nothing_raised do
+      vec = GSL_Vector.create( SIZE )
+
+      assert vec.length == SIZE
+    end
+  end
+
+  class ::GSL4r::Vector::GSL_Vector_Monitor < ::GSL4r::Vector::GSL_Vector
+    def self.release(ptr)
+      super
+      $release_count = $release_count + 1
+    end
+  end
+
+  class ::GSL4r::Vector::GSL_Vector_Cast_Monitor < ::GSL4r::Vector::GSL_Vector_Cast
+    def self.release(ptr)
+      super
+      $release_count = $release_count + 1
+    end
+  end
+
+  # tests using a cast object to retain a reference
+  # to a memory pointer, without freeing it when the
+  # object is garbage collected.
+  def test_gsl_vector_cast()
+    $release_count = 0
+
+    assert_nothing_raised do
+      vecptr = MemoryPointer.new :pointer
+      vecptr = GSL4r::Vector::Methods::gsl_vector_alloc( SIZE )
+
+      v1=GSL4r::Vector::GSL_Vector_Monitor.new( vecptr )
+      v2=GSL4r::Vector::GSL_Vector_Cast_Monitor.new( vecptr )
+
+      assert v1.length == v2.length
+    end
+
+    v1=nil
+    v2=nil
+    GC.start
+
+    assert $release_count == 1
+
+  end
+
+  def test_gsl_vector_set_values()
+    assert_nothing_raised do
+      vec = GSL_Vector.create( SIZE )
+
+      vec.set_with_array( (0..SIZE-1).to_a.collect { |i| i=i*i } )
+      v = vec.values
+      (0..SIZE-1).each { |i| assert v[i] == i*i }
+
+    end
+  end
+
+  def test_gsl_vector_set_all()
+    assert_nothing_raised do
+      vec = GSL_Vector.create( SIZE )
+
+      vec.set_all( 5 )
+      v = vec.values
+      s = 0.0
+      (0..SIZE-1).each { |i| s = s + v[i] }
+
+      assert (SIZE*5) == s
+    end
+  end
+
+  def test_gsl_vector_subvector()
+    assert_nothing_raised do
+      vec = GSL_Vector.create( SIZE )
+
+      a = Array.new( SIZE )
+      vec.set_with_array( (0..SIZE-1).to_a.collect { |i| i=i+1.0 } )
+
+      subvec = vec.subvector( 10, 5 )
+
+      v = subvec.values
+
+      test_a = (0..4).to_a.collect { |i| i=i+11.0 }
+
+      assert v == test_a
+    end
+  end
+
+end

data/test/vector_test.rb

@@ -0,0 +1,104 @@
+$: << File.join('..','lib')
+require 'test/unit'
+require 'test/unit/autorunner'
+require 'gsl4r/vector'
+include GSL4r::Vector
+class VectorTests < Test::Unit::TestCase
+  EPSILON = 5.0e-15
+def test_gsl_vector_sub()
+  v76 = GSL_Vector.create(3)
+  v76.set( [1.0,2.0,3.0] )
+  v77 = GSL_Vector.create(3)
+  v77.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_sub(v76,v77)
+  assert v76.values == [0,0,0]
+end
+def test_gsl_vector_scale()
+  v82 = GSL_Vector.create(3)
+  v82.set( [1.0,2.0,3.0] )
+
+  v83 = 2.0
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_scale(v82,v83)
+  assert v82.values == [2,4,6]
+end
+def test_gsl_vector_max_index()
+  v87 = GSL_Vector.create(3)
+  v87.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_max_index(v87)
+  assert_in_delta r_r1, 2, EPSILON
+end
+def test_gsl_vector_isneg()
+  v92 = GSL_Vector.create(3)
+  v92.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_isneg(v92)
+  assert_in_delta r_r1, 0, EPSILON
+end
+def test_gsl_vector_mul()
+  v78 = GSL_Vector.create(3)
+  v78.set( [1.0,2.0,3.0] )
+  v79 = GSL_Vector.create(3)
+  v79.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_mul(v78,v79)
+  assert v78.values == [1,4,9]
+end
+def test_gsl_vector_add_constant()
+  v84 = GSL_Vector.create(3)
+  v84.set( [1.0,2.0,3.0] )
+
+  v85 = 2.0
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_add_constant(v84,v85)
+  assert v84.values == [3,4,5]
+end
+def test_gsl_vector_isnull()
+  v90 = GSL_Vector.create(3)
+  v90.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_isnull(v90)
+  assert_in_delta r_r1, 0, EPSILON
+end
+def test_gsl_vector_min()
+  v88 = GSL_Vector.create(3)
+  v88.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_min(v88)
+  assert_in_delta r_r1, 1, EPSILON
+end
+def test_gsl_vector_isnonneg()
+  v93 = GSL_Vector.create(3)
+  v93.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_isnonneg(v93)
+  assert_in_delta r_r1, 1, EPSILON
+end
+def test_gsl_vector_add()
+  v74 = GSL_Vector.create(3)
+  v74.set( [1.0,2.0,3.0] )
+  v75 = GSL_Vector.create(3)
+  v75.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_add(v74,v75)
+  assert v74.values == [2,4,6]
+end
+def test_gsl_vector_div()
+  v80 = GSL_Vector.create(3)
+  v80.set( [1.0,2.0,3.0] )
+  v81 = GSL_Vector.create(3)
+  v81.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_div(v80,v81)
+  assert v80.values == [1,1,1]
+end
+def test_gsl_vector_max()
+  v86 = GSL_Vector.create(3)
+  v86.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_max(v86)
+  assert_in_delta r_r1, 3, EPSILON
+end
+def test_gsl_vector_ispos()
+  v91 = GSL_Vector.create(3)
+  v91.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_ispos(v91)
+  assert_in_delta r_r1, 1, EPSILON
+end
+def test_gsl_vector_min_index()
+  v89 = GSL_Vector.create(3)
+  v89.set( [1.0,2.0,3.0] )
+  r_r1 = ::GSL4r::Vector::Methods::gsl_vector_min_index(v89)
+  assert_in_delta r_r1, 0, EPSILON
+end
+end

metadata

@@ -1,7 +1,12 @@
 --- !ruby/object:Gem::Specification 
 name: gsl4r
 version: !ruby/object:Gem::Version 
-  version: 0.0.2
+  prerelease: false
+  segments: 
+  - 0
+  - 0
+  - 3
+  version: 0.0.3
 platform: ruby
 authors: 
 - Colby Gutierrez-Kraybill
@@ -9,7 +14,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2010-03-01 00:00:00 -08:00
+date: 2010-03-22 00:00:00 -07:00
 default_executable: 
 dependencies: []
 
@@ -32,6 +37,7 @@ files:
 - Rakefile
 - TODO
 - changelog
+- lib/b_t.rb
 - lib/gsl4r/block.rb
 - lib/gsl4r/complex.rb
 - lib/gsl4r/const.rb
@@ -46,10 +52,15 @@ files:
 - lib/gsl4r/vector.rb
 - lib/gsl4r.rb
 - lib/t.rb
+- test/block_class_test.rb
 - test/complex_test.rb
 - test/gsl_complex_tests_gen
 - test/gsl_complex_tests_gen.c
+- test/gsl_vector_tests_gen
+- test/gsl_vector_tests_gen.c
 - test/README
+- test/vector_class_test.rb
+- test/vector_test.rb
 has_rdoc: true
 homepage: http://gsl4r.rubyforge.org
 licenses: []
@@ -63,18 +74,20 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements: 
   - - ">="
     - !ruby/object:Gem::Version 
+      segments: 
+      - 0
       version: "0"
-  version: 
 required_rubygems_version: !ruby/object:Gem::Requirement 
   requirements: 
   - - ">="
     - !ruby/object:Gem::Version 
+      segments: 
+      - 0
       version: "0"
-  version: 
 requirements: 
 - GNU Scientific Library, version 1.13 or greater
 rubyforge_project: gsl4r
-rubygems_version: 1.3.5
+rubygems_version: 1.3.6
 signing_key: 
 specification_version: 3
 summary: GSL4r, ruby FFI wrappers around GNU Scientific Library