ruby-gsl-ng 0.2.2 → 0.2.3

data/History.txt

@@ -1,3 +1,24 @@
+=== 0.2.3
+* FIXES:
+  * HUGE memory leak, no GSL object was ever freed.
+  * Remove -march flag for MAC
+* LOTS of improvements:
+  * Added map_array: a map() version that actually returns an Array instead of a Vector/Matrix
+  * Now you use the & idiom for map methods and such (like map_array(&:to_s))
+  * Added Matrix#each_row / Matrix#each_column (which pass row and column matrices to the block)
+  * Added Matrix#row_vecview / Matrix#column_vecview (construct a row and column Matrix::View)
+  * Added Matrix#each_vec_row / Matrix#each_vec_column (same as the previous ones, but the block will receive a Vector::View instead)
+  * Added x, y, z and w accessors for Vector (this will probably go into Vector2, Vector3 and Vector4 in the future)
+  * "negate" operator
+  * Documentation is now YARD based. This will significantly improve the documentation readability.
+  * Added Vector#map_index! and Vector#wrap!
+  * Initial support for Random Number Generators (only Uniform and Gaussian distributions for the moment)
+  * Added Vector#<, #>, #<= and #>=
+  * Added Matrix#map_with_index!
+  * No more #fast_each and stuff. Now #each IS fast. This is accomplished by writing actual Ruby-C code instead of using FFI for
+    these functions.
+  * Initial support for GSL's special functions
+
 === 0.2.1
 * Fixed for FFI 0.6.0
 

data/Manifest

@@ -8,15 +8,24 @@ lib/gslng.rb
 lib/gslng/backend.rb
 lib/gslng/backend_components/error_handling.rb
 lib/gslng/backend_components/matrix.rb
+lib/gslng/backend_components/rng.rb
+lib/gslng/backend_components/special.rb
 lib/gslng/backend_components/vector.rb
 lib/gslng/finalizer.rb
 lib/gslng/matrix.rb
 lib/gslng/matrix_view.rb
+lib/gslng/rng/gaussian.rb
+lib/gslng/rng/rng.rb
+lib/gslng/rng/uniform.rb
+lib/gslng/special.rb
 lib/gslng/vector.rb
 lib/gslng/vector_view.rb
 ruby-gsl-ng.gemspec
 test/benchmark.rb
+test/benchmark_results
 test/matrix_test.rb
+test/rng_test.rb
 test/test_gsl.rb
+test/test_special.rb
 test/vector_test.rb
 Manifest

data/Rakefile

@@ -7,14 +7,18 @@ require 'yard'
 Echoe.new('ruby-gsl-ng') do |p|
   p.author = 'v01d'
   p.summary = "Ruby Object Oriented Graph LIbrary"
-  p.url = "http://github.com/v01d/roogli"
-  p.version = "0.2.2"
+  p.url = "http://github.com/v01d/ruby-gsl-ng"
+  p.version = "0.2.3"
   p.dependencies = ['yard', 'ffi']
 #  p.eval = proc { s.has_rdoc = 'yard' }
 end
 
-YARD::Rake::YardocTask.new do |t|
-  t.files   = ['lib/**/*.rb']   # optional
-  t.options = ['--verbose','--no-private'] # optional
+Rake::TaskManager.class_eval do
+  def remove_task(task)
+    @tasks.delete(task.to_s)
+  end
 end
 
+Rake.application.remove_task(:docs)
+YARD::Rake::YardocTask.new(:docs) {|t| t.options = ['--verbose','--no-private','--hide-void']}
+

data/ext/extconf.rb

@@ -1,6 +1,6 @@
 #!/usr/bin/ruby
 require 'mkmf'
 gsl_vars = pkg_config('gsl') or raise 'GSL not found!'
-extra_flags='-O3 -march=native'
-with_cppflags("#{$CPPFLAGS} #{extra_flags}") { true }
+extra_flags='-O4 -DHAVE_INLINE'
+with_cppflags("#{extra_flags}") { true }
 create_makefile('gslng_extensions')

data/ext/gslng_extensions.cpp

@@ -4,43 +4,66 @@
  * This could actually a standard C library (not a Ruby extension), but it's easier this way.
  */
 
+#include <ruby.h>
 #include <gsl/gsl_vector.h>
 #include <gsl/gsl_matrix.h>
 
-extern "C" void Init_gslng_extensions(void) { }
+/************************* Vector functions *****************************/
 
-/**** Vector *****/
-
-// For Vector::map!
-typedef double (*gsl_vector_callback_t)(double);
-
-extern "C" void gsl_vector_map(gsl_vector* v, gsl_vector_callback_t callback) {
+static VALUE gsl_vector_map(VALUE self, VALUE ptr) {
+	gsl_vector* v = (gsl_vector*)FIX2ULONG(ptr);
 	for (size_t i = 0; i < v->size; i++)
-		*gsl_vector_ptr(v, i) = (*callback)(*gsl_vector_const_ptr(v, i));
+		gsl_vector_set(v, i, NUM2DBL(rb_yield(rb_float_new(gsl_vector_get(v, i)))));
+	
+	return self;
 }
 
-// For Vector::map_index!
-typedef double (*gsl_vector_index_callback_t)(size_t);
-
-extern "C" void gsl_vector_map_index(gsl_vector* v, gsl_vector_index_callback_t callback) {
-	for (size_t i = 0; i < v->size; i++)
-		*gsl_vector_ptr(v, i) = (*callback)(i);
+static VALUE gsl_vector_map_index(VALUE self, VALUE ptr) {
+	gsl_vector* v = (gsl_vector*)FIX2ULONG(ptr);
+	for (size_t i = 0; i < v->size; i++) {
+		VALUE vi = ULONG2NUM(i);
+		gsl_vector_set(v, i, NUM2DBL(rb_yield(vi)));
+	}
+	
+	return self;
 }
 
-// A fast "each" for cases where there's no expected return value from the block
-typedef void (*gsl_vector_each_callback_t)(double);
+static VALUE gsl_vector_each_with_index(VALUE self, VALUE ptr) {
+	gsl_vector* v = (gsl_vector*)FIX2ULONG(ptr);
+	for (size_t i = 0; i < v->size; i++) {
+		VALUE vi = ULONG2NUM(i);
+		rb_yield_values(2, rb_float_new(gsl_vector_get(v, i)), vi);
+	}
+	
+	return self;
+}
 
-extern "C" void gsl_vector_each(gsl_vector* v, gsl_vector_each_callback_t callback) {
+static VALUE gsl_vector_each(VALUE self, VALUE ptr) {
+	gsl_vector* v = (gsl_vector*)FIX2ULONG(ptr);
 	for (size_t i = 0; i < v->size; i++)
-		(*callback)(*gsl_vector_const_ptr(v, i));
+		rb_yield(rb_float_new(gsl_vector_get(v, i)));
+	
+	return self;
 }
 
-// A fast "each_with_index" for cases where there's no expected return value from the block
-typedef void (*gsl_vector_each_with_index_callback_t)(double, size_t);
+static VALUE gsl_vector_to_a(VALUE self, VALUE ptr) {
+	gsl_vector* v = (gsl_vector*)FIX2ULONG(ptr);
+	
+	VALUE array = rb_ary_new2(v->size);
+	for (size_t i = 0; i < v->size; i++)
+		rb_ary_store(array, i, rb_float_new(gsl_vector_get(v, i)));
+	
+	return array;
+}
 
-extern "C" void gsl_vector_each_with_index(gsl_vector* v, gsl_vector_each_with_index_callback_t callback) {
+static VALUE gsl_vector_from_array(VALUE self, VALUE ptr, VALUE array) {
+	gsl_vector* v = (gsl_vector*)FIX2ULONG(ptr);
+	if (v->size != RARRAY_LEN(array)) rb_raise(rb_eRuntimeError, "Sizes differ!");
+	
 	for (size_t i = 0; i < v->size; i++)
-		(*callback)(*gsl_vector_const_ptr(v, i), i);
+		gsl_vector_set(v, i, NUM2DBL(rb_ary_entry(array, i)));
+	
+	return self;
 }
 
 // Hide the view in a new vector (gsl_vector_subvector)
@@ -59,41 +82,98 @@ extern "C" gsl_vector* gsl_vector_subvector2(gsl_vector* v, size_t offset, size_
   return vector_view;
 }
 
-/***** Matrix *****/
-// For Matrix::map!
-typedef double (*gsl_matrix_callback_t)(double);
+/************************* Matrix functions *****************************/
 
-extern "C" void gsl_matrix_map(gsl_matrix* m, gsl_matrix_callback_t callback) {
+static VALUE gsl_matrix_map(VALUE self, VALUE ptr) {
+  gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
   size_t size1 = m->size1;
   size_t size2 = m->size2;
 
 	for (size_t i = 0; i < size1; i++)
     for (size_t j = 0; j < size2; j++)
-      *gsl_matrix_ptr(m, i, j) = (*callback)(*gsl_matrix_const_ptr(m, i, j));
-}
+      gsl_matrix_set(m, i, j, NUM2DBL(rb_yield(rb_float_new(gsl_matrix_get(m, i, j)))));
 
-// For Matrix::map_index!
-typedef double (*gsl_matrix_index_callback_t)(size_t, size_t);
+  return self;
+}
 
-extern "C" void gsl_matrix_map_index(gsl_matrix* m, gsl_matrix_index_callback_t callback) {
+static VALUE gsl_matrix_map_index(VALUE self, VALUE ptr) {
+  gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
   size_t size1 = m->size1;
   size_t size2 = m->size2;
 
 	for (size_t i = 0; i < size1; i++)
     for (size_t j = 0; j < size2; j++)
-  		*gsl_matrix_ptr(m, i, j) = (*callback)(i, j);
+      gsl_matrix_set(m, i, j, NUM2DBL(rb_yield_values(2, ULONG2NUM(i), ULONG2NUM(j))));
+
+  return self;
 }
 
-// A fast "each" for cases where there's no expected return value from the block
-typedef void (*gsl_matrix_each_callback_t)(double);
+static VALUE gsl_matrix_map_with_index(VALUE self, VALUE ptr) {
+  gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
+  size_t size1 = m->size1;
+  size_t size2 = m->size2;
 
-extern "C" void gsl_matrix_each(gsl_matrix* m, gsl_matrix_each_callback_t callback) {
+	for (size_t i = 0; i < size1; i++)
+    for (size_t j = 0; j < size2; j++)
+      gsl_matrix_set(m, i, j, NUM2DBL(rb_yield_values(3, rb_float_new(gsl_matrix_get(m, i, j)), ULONG2NUM(i), ULONG2NUM(j))));
+
+  return self;
+}
+
+static VALUE gsl_matrix_each(VALUE self, VALUE ptr) {
+  gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
   size_t size1 = m->size1;
   size_t size2 = m->size2;
 
 	for (size_t i = 0; i < size1; i++)
     for (size_t j = 0; j < size2; j++)
-  		(*callback)(*gsl_matrix_const_ptr(m, i, j));
+      rb_yield(rb_float_new(gsl_matrix_get(m, i, j)));
+
+  return self;
+}
+
+static VALUE gsl_matrix_each_with_index(VALUE self, VALUE ptr) {
+  gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
+  size_t size1 = m->size1;
+  size_t size2 = m->size2;
+
+	for (size_t i = 0; i < size1; i++) {
+    VALUE vi = ULONG2NUM(i);
+    for (size_t j = 0; j < size2; j++)
+      rb_yield_values(3, rb_float_new(gsl_matrix_get(m, i, j)), vi, ULONG2NUM(j));
+  }
+  
+  return self;
+}
+
+static VALUE gsl_matrix_to_a(VALUE self, VALUE ptr) {
+	gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
+
+	VALUE array = rb_ary_new2(m->size1);
+	for (size_t i = 0; i < m->size1; i++) {
+    VALUE row = rb_ary_new2(m->size2);
+    for (size_t j = 0; j < m->size2; j++) {
+      rb_ary_store(row, j, rb_float_new(gsl_matrix_get(m, i, j)));
+    }
+    rb_ary_store(array, i, row);
+  }
+
+	return array;
+}
+
+static VALUE gsl_matrix_from_array(VALUE self, VALUE ptr, VALUE array) {
+  gsl_matrix* m = (gsl_matrix*)FIX2ULONG(ptr);
+  if (m->size1 != RARRAY_LEN(array)) rb_raise(rb_eRuntimeError, "Sizes differ!");
+
+	for (size_t i = 0; i < m->size1; i++) {
+    VALUE row = rb_ary_entry(array, i);
+    if (m->size2 != RARRAY_LEN(row)) rb_raise(rb_eRuntimeError, "Sizes differ!");
+
+    for (size_t j = 0; j < m->size2; j++)
+      gsl_matrix_set(m, i, j, NUM2DBL(rb_ary_entry(row, j)));
+  }
+
+	return self;
 }
 
 // Hide the view in a new matrix (gsl_matrix_submatrix)
@@ -116,4 +196,29 @@ extern "C" gsl_vector* gsl_matrix_column_view(gsl_matrix* m_ptr, size_t column,
   gsl_vector* vector_view = gsl_vector_alloc(view.vector.size);
   *vector_view = view.vector;
   return vector_view;
+}
+
+
+/************************* Module initialization *****************************/
+
+extern "C" void Init_gslng_extensions(void) {
+	VALUE GSLng_module = rb_define_module("GSLng");
+	VALUE Backend_module = rb_funcall(GSLng_module, rb_intern("backend"), 0);
+
+	// vector
+	rb_define_module_function(Backend_module, "gsl_vector_map!", (VALUE(*)(ANYARGS))gsl_vector_map, 1);
+	rb_define_module_function(Backend_module, "gsl_vector_map_index!", (VALUE(*)(ANYARGS))gsl_vector_map_index, 1);
+	rb_define_module_function(Backend_module, "gsl_vector_each_with_index", (VALUE(*)(ANYARGS))gsl_vector_each_with_index, 1);
+	rb_define_module_function(Backend_module, "gsl_vector_each", (VALUE(*)(ANYARGS))gsl_vector_each, 1);
+	rb_define_module_function(Backend_module, "gsl_vector_to_a", (VALUE(*)(ANYARGS))gsl_vector_to_a, 1);
+	rb_define_module_function(Backend_module, "gsl_vector_from_array", (VALUE(*)(ANYARGS))gsl_vector_from_array, 2);
+
+  // matrix
+	rb_define_module_function(Backend_module, "gsl_matrix_map!", (VALUE(*)(ANYARGS))gsl_matrix_map, 1);
+	rb_define_module_function(Backend_module, "gsl_matrix_map_index!", (VALUE(*)(ANYARGS))gsl_matrix_map_index, 1);
+  rb_define_module_function(Backend_module, "gsl_matrix_map_with_index!", (VALUE(*)(ANYARGS))gsl_matrix_map_with_index, 1);
+	rb_define_module_function(Backend_module, "gsl_matrix_each_with_index", (VALUE(*)(ANYARGS))gsl_matrix_each_with_index, 1);
+	rb_define_module_function(Backend_module, "gsl_matrix_each", (VALUE(*)(ANYARGS))gsl_matrix_each, 1);
+	rb_define_module_function(Backend_module, "gsl_matrix_to_a", (VALUE(*)(ANYARGS))gsl_matrix_to_a, 1);
+	rb_define_module_function(Backend_module, "gsl_matrix_from_array", (VALUE(*)(ANYARGS))gsl_matrix_from_array, 2);
 }

data/lib/gslng.rb

@@ -1,8 +1,9 @@
-require 'gslng/backend'
-require 'gslng/finalizer'
-
-require 'gslng/vector'
-require 'gslng/vector_view'
-require 'gslng/matrix'
-require 'gslng/matrix_view'
-
+require 'gslng/backend'
+require 'gslng/finalizer'
+
+require 'gslng/vector'
+require 'gslng/vector_view'
+require 'gslng/matrix'
+require 'gslng/matrix_view'
+require 'gslng/rng/rng'
+require 'gslng/special'

data/lib/gslng/backend.rb

@@ -1,15 +1,15 @@
 require 'ffi'
 
 module GSLng
-	# Anonymous module: avoids exposing this internal module when doing "include GSLng" at the top-level.
-	# If ruby had "private" modules I wouldn't have to do this.
+  # Anonymous module: avoids exposing this internal module when doing "include GSLng" at the top-level.
+  # If ruby had "private" modules I wouldn't have to do this.
   @backend = Module.new do
     extend FFI::Library
     ffi_lib(FFI::CURRENT_PROCESS)
   end
 
   # Returns the internal backend module
-	def GSLng.backend # @private
+  def GSLng.backend # @private
     @backend
   end
 end
@@ -18,3 +18,5 @@ require 'gslng_extensions'
 require 'gslng/backend_components/vector'
 require 'gslng/backend_components/matrix'
 require 'gslng/backend_components/error_handling'
+require 'gslng/backend_components/rng'
+require 'gslng/backend_components/special'

data/lib/gslng/backend_components/error_handling.rb

@@ -1,12 +1,12 @@
-module GSLng
-	backend.instance_eval do
-    # All of this raises an exception when GSL registers an error
-		callback :error_handler_callback, [ :string, :string, :int, :int ], :void
-    attach_function :gsl_set_error_handler, [ :error_handler_callback ], :error_handler_callback
-
-    ErrorHandlerCallback = Proc.new {|reason, file, line, errno|
-      raise RuntimeError, "#{reason} (errno: #{errno})", caller[2..-1]
-    }
-    gsl_set_error_handler(ErrorHandlerCallback)
-	end
+module GSLng
+  backend.instance_eval do
+    # All of this raises an exception when GSL registers an error
+    callback :error_handler_callback, [ :string, :string, :int, :int ], :void
+    attach_function :gsl_set_error_handler, [ :error_handler_callback ], :error_handler_callback
+
+    ErrorHandlerCallback = Proc.new {|reason, file, line, errno|
+      raise RuntimeError, "#{reason} (errno: #{errno})", caller[2..-1]
+    }
+    gsl_set_error_handler(ErrorHandlerCallback)
+  end
 end

data/lib/gslng/backend_components/matrix.rb

@@ -1,75 +1,66 @@
-module GSLng
-	backend.instance_eval do
-		# memory handling
-		attach_function :gsl_matrix_alloc, [ :size_t, :size_t ], :pointer
-		attach_function :gsl_matrix_calloc, [ :size_t, :size_t ], :pointer
-		attach_function :gsl_matrix_free, [ :pointer ], :void
-
-    # initializing
-    attach_function :gsl_matrix_set_all, [ :pointer, :double ], :void
-    attach_function :gsl_matrix_set_zero, [ :pointer ], :void
-    attach_function :gsl_matrix_set_identity, [ :pointer ], :void
-
-    # copying
-    attach_function :gsl_matrix_memcpy, [ :pointer, :pointer ], :int
-
-    # operations
-    attach_function :gsl_matrix_add, [ :pointer, :pointer ], :int
-    attach_function :gsl_matrix_sub, [ :pointer, :pointer ], :int
-    attach_function :gsl_matrix_mul_elements, [ :pointer, :pointer ], :int
-    attach_function :gsl_matrix_div_elements, [ :pointer, :pointer ], :int
-    attach_function :gsl_matrix_scale, [ :pointer, :double ], :int
-    attach_function :gsl_matrix_add_constant, [ :pointer, :double ], :int
-
-		# copy rows/columns into
-		attach_function :gsl_matrix_get_row, [ :pointer, :pointer, :size_t ], :int
-		attach_function :gsl_matrix_get_col, [ :pointer, :pointer, :size_t ], :int
-		attach_function :gsl_matrix_set_row, [ :pointer, :size_t, :pointer ], :int
-		attach_function :gsl_matrix_set_col, [ :pointer, :size_t, :pointer ], :int
-
-		# element access
-    attach_function :gsl_matrix_get, [ :pointer, :size_t, :size_t ], :double
-    attach_function :gsl_matrix_set, [ :pointer, :size_t, :size_t, :double ], :void
-    
-    # properties
-    attach_function :gsl_matrix_isnull, [ :pointer ], :int
-    attach_function :gsl_matrix_ispos, [ :pointer ], :int
-    attach_function :gsl_matrix_isneg, [ :pointer ], :int
-    attach_function :gsl_matrix_isnonneg, [ :pointer ], :int
-
-    # max and min
-    attach_function :gsl_matrix_max, [ :pointer ], :double
-    attach_function :gsl_matrix_min, [ :pointer ], :double
-    attach_function :gsl_matrix_minmax, [ :pointer, :buffer_out, :buffer_out ], :void
-    attach_function :gsl_matrix_max_index, [ :pointer, :buffer_out, :buffer_out ], :void
-    attach_function :gsl_matrix_min_index, [ :pointer, :buffer_out, :buffer_out ], :void
-    attach_function :gsl_matrix_minmax_index, [ :pointer, :buffer_out, :buffer_out, :buffer_out, :buffer_out ], :void
-
-    # exchange elements
-    attach_function :gsl_matrix_swap_rows, [ :pointer, :size_t, :size_t ], :int
-    attach_function :gsl_matrix_swap_columns, [ :pointer, :size_t, :size_t ], :int
-    attach_function :gsl_matrix_swap_rowcol, [ :pointer, :size_t, :size_t ], :int
-    attach_function :gsl_matrix_transpose_memcpy, [ :pointer, :pointer ], :int
-    attach_function :gsl_matrix_transpose, [ :pointer ], :int
-		
-    # From local extension
-		callback :gsl_matrix_callback, [ :double ], :double
-		attach_function :gsl_matrix_map, [ :pointer, :gsl_matrix_callback ], :void
-
-		callback :gsl_matrix_index_callback, [ :size_t, :size_t ], :double
-		attach_function :gsl_matrix_map_index, [ :pointer, :gsl_matrix_index_callback ], :void
-
-		callback :gsl_matrix_each_callback, [ :double ], :void
-		attach_function :gsl_matrix_each, [ :pointer, :gsl_matrix_each_callback ], :void
-
-    # views
-    attach_function :gsl_matrix_submatrix2, [ :pointer, :size_t, :size_t, :size_t, :size_t ], :pointer
-    attach_function :gsl_matrix_row_view, [ :pointer, :size_t, :size_t, :size_t ], :pointer
-    attach_function :gsl_matrix_column_view, [ :pointer, :size_t, :size_t, :size_t ], :pointer
-
-    # BLAS interface
-    enum :cblas_transpose_t, [ :no_transpose, 111, :transpose, :conjugate_transpose ]
-    attach_function :gsl_blas_dgemv, [ :cblas_transpose_t, :double, :pointer, :pointer, :double, :pointer ], :int
-    attach_function :gsl_blas_dgemm, [ :cblas_transpose_t, :cblas_transpose_t, :double, :pointer, :pointer, :double, :pointer ], :int
-	end
-end
+module GSLng
+  backend.instance_eval do
+    # memory handling
+    attach_function :gsl_matrix_alloc, [ :size_t, :size_t ], :pointer
+    attach_function :gsl_matrix_calloc, [ :size_t, :size_t ], :pointer
+    attach_function :gsl_matrix_free, [ :pointer ], :void
+
+    # initializing
+    attach_function :gsl_matrix_set_all, [ :pointer, :double ], :void
+    attach_function :gsl_matrix_set_zero, [ :pointer ], :void
+    attach_function :gsl_matrix_set_identity, [ :pointer ], :void
+
+    # copying
+    attach_function :gsl_matrix_memcpy, [ :pointer, :pointer ], :int
+
+    # operations
+    attach_function :gsl_matrix_add, [ :pointer, :pointer ], :int
+    attach_function :gsl_matrix_sub, [ :pointer, :pointer ], :int
+    attach_function :gsl_matrix_mul_elements, [ :pointer, :pointer ], :int
+    attach_function :gsl_matrix_div_elements, [ :pointer, :pointer ], :int
+    attach_function :gsl_matrix_scale, [ :pointer, :double ], :int
+    attach_function :gsl_matrix_add_constant, [ :pointer, :double ], :int
+
+    # copy rows/columns into
+    attach_function :gsl_matrix_get_row, [ :pointer, :pointer, :size_t ], :int
+    attach_function :gsl_matrix_get_col, [ :pointer, :pointer, :size_t ], :int
+    attach_function :gsl_matrix_set_row, [ :pointer, :size_t, :pointer ], :int
+    attach_function :gsl_matrix_set_col, [ :pointer, :size_t, :pointer ], :int
+
+    # element access
+    attach_function :gsl_matrix_get, [ :pointer, :size_t, :size_t ], :double
+    attach_function :gsl_matrix_set, [ :pointer, :size_t, :size_t, :double ], :void
+    
+    # properties
+    attach_function :gsl_matrix_isnull, [ :pointer ], :int
+    attach_function :gsl_matrix_ispos, [ :pointer ], :int
+    attach_function :gsl_matrix_isneg, [ :pointer ], :int
+    attach_function :gsl_matrix_isnonneg, [ :pointer ], :int
+
+    # max and min
+    attach_function :gsl_matrix_max, [ :pointer ], :double
+    attach_function :gsl_matrix_min, [ :pointer ], :double
+    attach_function :gsl_matrix_minmax, [ :pointer, :buffer_out, :buffer_out ], :void
+    attach_function :gsl_matrix_max_index, [ :pointer, :buffer_out, :buffer_out ], :void
+    attach_function :gsl_matrix_min_index, [ :pointer, :buffer_out, :buffer_out ], :void
+    attach_function :gsl_matrix_minmax_index, [ :pointer, :buffer_out, :buffer_out, :buffer_out, :buffer_out ], :void
+
+    # exchange elements
+    attach_function :gsl_matrix_swap_rows, [ :pointer, :size_t, :size_t ], :int
+    attach_function :gsl_matrix_swap_columns, [ :pointer, :size_t, :size_t ], :int
+    attach_function :gsl_matrix_swap_rowcol, [ :pointer, :size_t, :size_t ], :int
+    attach_function :gsl_matrix_transpose_memcpy, [ :pointer, :pointer ], :int
+    attach_function :gsl_matrix_transpose, [ :pointer ], :int
+    
+    # From local extension    
+    # views
+    attach_function :gsl_matrix_submatrix2, [ :pointer, :size_t, :size_t, :size_t, :size_t ], :pointer
+    attach_function :gsl_matrix_row_view, [ :pointer, :size_t, :size_t, :size_t ], :pointer
+    attach_function :gsl_matrix_column_view, [ :pointer, :size_t, :size_t, :size_t ], :pointer
+
+    # BLAS interface
+    enum :cblas_transpose_t, [ :no_transpose, 111, :transpose, :conjugate_transpose ]
+    attach_function :gsl_blas_dgemv, [ :cblas_transpose_t, :double, :pointer, :pointer, :double, :pointer ], :int
+    attach_function :gsl_blas_dgemm, [ :cblas_transpose_t, :cblas_transpose_t, :double, :pointer, :pointer, :double, :pointer ], :int
+  end
+end