rb-gsl 1.15.3.1

data/include/rb_gsl_with_narray.h

@@ -0,0 +1,29 @@
+#ifndef RB_GSL_WITH_NARRAY_H
+#define RB_GSL_WITH_NARRAY_H
+
+#ifdef HAVE_NARRAY_H
+
+#include "narray.h"
+#include "gsl/gsl_vector.h"
+
+gsl_vector* make_cvector_from_narray(VALUE);
+void cvector_set_from_narray(gsl_vector*, VALUE);
+void carray_set_from_narray(double*, VALUE);
+
+VALUE rb_gsl_na_to_gsl_vector_view_method(VALUE na);
+VALUE rb_gsl_na_to_gsl_matrix(VALUE obj, VALUE nna);
+gsl_vector_view* na_to_gv_view(VALUE na);
+gsl_matrix_view* na_to_gm_view(VALUE nna);
+
+gsl_vector_int_view* na_to_gv_int_view(VALUE na);
+gsl_matrix_int_view* na_to_gm_int_view(VALUE nna);
+gsl_vector* na_to_gv(VALUE na);
+gsl_vector_int* na_to_gv_int(VALUE na);
+gsl_matrix* na_to_gm(VALUE nna);
+gsl_matrix_int* na_to_gm_int(VALUE nna);
+extern VALUE cNVector, cNMatrix;
+
+gsl_vector_complex* na_to_gv_complex(VALUE na);
+gsl_vector_complex_view* na_to_gv_complex_view(VALUE na);
+#endif // HAVE_NARRAY_H
+#endif // RB_GSL_WITH_NARRAY_H

data/include/templates_off.h

@@ -0,0 +1,87 @@
+#ifdef FUNCTION
+#undef FUNCTION
+#endif
+
+#ifdef CONCAT4
+#undef CONCAT4
+#endif
+
+#ifdef CONCAT4x
+#undef CONCAT4x
+#endif
+
+#ifdef CONCAT3
+#undef CONCAT3
+#endif
+
+#ifdef CONCAT3x
+#undef CONCAT3x
+#endif
+
+#ifdef CONCAT2
+#undef CONCAT2
+#endif
+
+#ifdef CONCAT2x
+#undef CONCAT2x
+#endif
+
+#ifdef GSL_TYPE
+#undef GSL_TYPE
+#endif
+
+#ifdef REAL_GSL_TYPE
+#undef REAL_GSL_TYPE
+#endif
+
+#ifdef QUALIFIED_GSL_TYPE
+#undef QUALIFIED_GSL_TYPE
+#endif
+
+#ifdef VIEW
+#undef VIEW
+#endif
+
+#ifdef REAL_VIEW
+#undef REAL_VIEW
+#endif
+
+#ifdef QUALIFIED_VIEW
+#undef QUALIFIED_VIEW
+#endif
+
+#ifdef QUALIFIED_REAL_GSL_TYPE
+#undef QUALIFIED_REAL_GSL_TYPE
+#endif
+
+#ifdef QUALIFIED_REAL_VIEW
+#undef QUALIFIED_REAL_VIEW
+#endif
+
+#ifdef USES_LONGDOUBLE
+#undef USES_LONGDOUBLE
+#endif
+
+#ifdef SHORT_REAL
+#undef SHORT_REAL
+#endif
+
+#ifndef USE_QUALIFIER
+#ifdef QUALIFIER
+#undef QUALIFIER
+#endif
+#endif
+
+#undef BASE
+#undef BASE_EPSILON
+#undef SHORT
+#undef ATOMIC
+#undef MULTIPLICITY
+#undef IN_FORMAT
+#undef OUT_FORMAT
+#undef ATOMIC_IO
+#undef ZERO
+#undef ONE
+#undef NAME
+#undef STRING
+#undef EXPAND

data/include/templates_on.h

@@ -0,0 +1,241 @@
+/* If BASE is undefined we use function names like gsl_name()
+   and assume that we are using doubles.
+
+   If BASE is defined we used function names like gsl_BASE_name()
+   and use BASE as the base datatype      */
+
+#if   defined(BASE_GSL_COMPLEX_LONG)
+#define BASE gsl_complex_long_double
+#define SHORT complex_long_double
+#define SHORT_REAL long_double
+#define ATOMIC long double
+#define USES_LONGDOUBLE 1
+#define MULTIPLICITY 2
+#define IN_FORMAT "%Lg"
+#define OUT_FORMAT "%Lg"
+#define ATOMIC_IO ATOMIC
+#define ZERO {{0.0L,0.0L}}
+#define ONE {{1.0L,0.0L}}
+#define BASE_EPSILON GSL_DBL_EPSILON
+
+#elif defined(BASE_GSL_COMPLEX)
+#define BASE gsl_complex
+#define SHORT complex
+#define SHORT_REAL
+#define ATOMIC double
+#define MULTIPLICITY 2
+#define IN_FORMAT "%lg"
+#define OUT_FORMAT "%g"
+#define ATOMIC_IO ATOMIC
+#define ZERO {{0.0,0.0}}
+#define ONE {{1.0,0.0}}
+#define BASE_EPSILON GSL_DBL_EPSILON
+
+#elif defined(BASE_GSL_COMPLEX_FLOAT)
+#define BASE gsl_complex_float
+#define SHORT complex_float
+#define SHORT_REAL float
+#define ATOMIC float
+#define MULTIPLICITY 2
+#define IN_FORMAT "%g"
+#define OUT_FORMAT "%g"
+#define ATOMIC_IO ATOMIC
+#define ZERO {{0.0F,0.0F}}
+#define ONE {{1.0F,0.0F}}
+#define BASE_EPSILON GSL_FLT_EPSILON
+
+#elif defined(BASE_LONG_DOUBLE)
+#define BASE long double
+#define SHORT long_double
+#define ATOMIC long double
+#define USES_LONGDOUBLE 1
+#define MULTIPLICITY 1
+#define IN_FORMAT "%Lg"
+#define OUT_FORMAT "%Lg"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0.0L
+#define ONE 1.0L
+#define BASE_EPSILON GSL_DBL_EPSILON
+
+#elif defined(BASE_DOUBLE)
+#define BASE double
+#define SHORT
+#define ATOMIC double
+#define MULTIPLICITY 1
+#define IN_FORMAT "%lg"
+#define OUT_FORMAT "%g"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0.0
+#define ONE 1.0
+#define BASE_EPSILON GSL_DBL_EPSILON
+
+#elif defined(BASE_FLOAT)
+#define BASE float
+#define SHORT float
+#define ATOMIC float
+#define MULTIPLICITY 1
+#define IN_FORMAT "%g"
+#define OUT_FORMAT "%g"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0.0F
+#define ONE 1.0F
+#define BASE_EPSILON GSL_FLT_EPSILON
+
+#elif defined(BASE_ULONG)
+#define BASE unsigned long
+#define SHORT ulong
+#define ATOMIC unsigned long
+#define MULTIPLICITY 1
+#define IN_FORMAT "%lu"
+#define OUT_FORMAT "%lu"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0UL
+#define ONE 1UL
+
+#elif defined(BASE_LONG)
+#define BASE long
+#define SHORT long
+#define ATOMIC long
+#define MULTIPLICITY 1
+#define IN_FORMAT "%ld"
+#define OUT_FORMAT "%ld"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0L
+#define ONE 1L
+
+#elif defined(BASE_UINT)
+#define BASE unsigned int
+#define SHORT uint
+#define ATOMIC unsigned int
+#define MULTIPLICITY 1
+#define IN_FORMAT "%u"
+#define OUT_FORMAT "%u"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0U
+#define ONE 1U
+
+#elif defined(BASE_INT)
+#define BASE int
+#define SHORT int
+#define ATOMIC int
+#define MULTIPLICITY 1
+#define IN_FORMAT "%d"
+#define OUT_FORMAT "%d"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0
+#define ONE 1
+
+#elif defined(BASE_USHORT)
+#define BASE unsigned short
+#define SHORT ushort
+#define ATOMIC unsigned short
+#define MULTIPLICITY 1
+#define IN_FORMAT "%hu"
+#define OUT_FORMAT "%hu"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0U
+#define ONE 1U
+
+#elif defined(BASE_SHORT)
+#define BASE short
+#define SHORT short
+#define ATOMIC short
+#define MULTIPLICITY 1
+#define IN_FORMAT "%hd"
+#define OUT_FORMAT "%hd"
+#define ATOMIC_IO ATOMIC
+#define ZERO 0
+#define ONE 1
+
+#elif defined(BASE_UCHAR)
+#define BASE unsigned char
+#define SHORT uchar
+#define ATOMIC unsigned char
+#define MULTIPLICITY 1
+#define IN_FORMAT "%u"
+#define OUT_FORMAT "%u"
+#define ATOMIC_IO unsigned int
+#define ZERO 0U
+#define ONE 1U
+
+#elif defined(BASE_CHAR)
+#define BASE char
+#define SHORT char
+#define ATOMIC char
+#define MULTIPLICITY 1
+#define IN_FORMAT "%d"
+#define OUT_FORMAT "%d"
+#define ATOMIC_IO int
+#define ZERO 0
+#define ONE 1
+
+#else
+#error unknown BASE_ directive in source.h
+#endif
+
+#define CONCAT2x(a,b) a ## _ ## b 
+#define CONCAT2(a,b) CONCAT2x(a,b)
+#define CONCAT3x(a,b,c) a ## _ ## b ## _ ## c
+#define CONCAT3(a,b,c) CONCAT3x(a,b,c)
+#define CONCAT4x(a,b,c,d) a ## _ ## b ## _ ## c ## _ ## d
+#define CONCAT4(a,b,c,d) CONCAT4x(a,b,c,d)
+
+#ifndef USE_QUALIFIER
+#define QUALIFIER
+#endif
+
+#ifdef USE_QUALIFIER
+#if defined(BASE_DOUBLE)
+#define FUNCTION(dir,name) CONCAT3(dir,QUALIFIER,name)
+#define GSL_TYPE(dir) dir
+#define VIEW(dir,name) CONCAT2(dir,name)
+#define QUALIFIED_GSL_TYPE(dir) QUALIFIER dir
+#define QUALIFIED_VIEW(dir,name) CONCAT3(dir,QUALIFIER,name)
+#else
+#define FUNCTION(a,c) CONCAT4(a,SHORT,QUALIFIER,c)
+#define GSL_TYPE(dir) CONCAT2(dir,SHORT)
+#define VIEW(dir,name) CONCAT3(dir,SHORT,name)
+#define QUALIFIED_GSL_TYPE(dir) QUALIFIER CONCAT2(dir,SHORT)
+#define QUALIFIED_VIEW(dir,name) CONCAT4(dir,SHORT,QUALIFIER,name)
+#endif
+#if defined(BASE_GSL_COMPLEX)
+#define REAL_GSL_TYPE(dir) dir
+#define REAL_VIEW(dir,name) CONCAT2(dir,name)
+#define QUALIFIED_REAL_GSL_TYPE(dir) QUALIFIER dir
+#define QUALIFIED_REAL_VIEW(dir,name) CONCAT3(dir,QUALIFIER,name)
+#else
+#define REAL_GSL_TYPE(dir) CONCAT2(dir,SHORT_REAL)
+#define REAL_VIEW(dir,name) CONCAT3(dir,SHORT_REAL,name)
+#define QUALIFIED_REAL_GSL_TYPE(dir) QUALIFIER CONCAT2(dir,SHORT_REAL)
+#define QUALIFIED_REAL_VIEW(dir,name) CONCAT4(dir,SHORT_REAL,QUALIFIER,name)
+#endif
+#else
+#if defined(BASE_DOUBLE)
+#define FUNCTION(dir,name) CONCAT2(dir,name)
+#define GSL_TYPE(dir) dir
+#define VIEW(dir,name) CONCAT2(dir,name)
+#define QUALIFIED_GSL_TYPE(dir) GSL_TYPE(dir)
+#define QUALIFIED_VIEW(dir,name) CONCAT2(dir,name)
+#else
+#define FUNCTION(a,c) CONCAT3(a,SHORT,c)
+#define GSL_TYPE(dir) CONCAT2(dir,SHORT)
+#define VIEW(dir,name) CONCAT3(dir,SHORT,name)
+#define QUALIFIED_GSL_TYPE(dir) GSL_TYPE(dir)
+#define QUALIFIED_VIEW(dir,name) CONCAT3(dir,SHORT,name)
+#endif
+#if defined(BASE_GSL_COMPLEX)
+#define REAL_GSL_TYPE(dir) dir
+#define REAL_VIEW(dir,name) CONCAT2(dir,name)
+#define QUALIFIED_REAL_GSL_TYPE(dir) dir
+#define QUALIFIED_REAL_VIEW(dir,name) CONCAT2(dir,name)
+#else
+#define REAL_GSL_TYPE(dir) CONCAT2(dir,SHORT_REAL)
+#define REAL_VIEW(dir,name) CONCAT3(dir,SHORT_REAL,name)
+#define QUALIFIED_REAL_GSL_TYPE(dir) CONCAT2(dir,SHORT_REAL)
+#define QUALIFIED_REAL_VIEW(dir,name) CONCAT3(dir,SHORT_REAL,name)
+#endif
+#endif
+
+#define STRING(x) #x
+#define EXPAND(x) STRING(x)
+#define NAME(x) EXPAND(GSL_TYPE(x))

data/lib/gsl.rb

@@ -0,0 +1,3 @@
+require('narray')
+require('rb_gsl')
+require('gsl/oper.rb')

data/lib/gsl/gnuplot.rb

@@ -0,0 +1,41 @@
+require("gnuplot")
+
+class Array
+  def to_gplot
+    if ( self[0].kind_of? Array ) then
+      tmp = self[0].zip( *self[1..-1] )
+      tmp.collect { |a| a.join(" ") }.join("\n") + "\ne"
+    elsif ( self[0].kind_of? Numeric ) then
+      s = ""
+      self.length.times { |i| s << "#{self[i]}\n" }
+      s
+    elsif ( self[0].kind_of? GSL::Vector ) then
+      tmp = self[0].zip( *self[1..-1] )
+      tmp.collect { |a| a.join(" ") }.join("\n") + "\ne"
+    else
+      self[0].zip( *self[1..-1] ).to_gplot
+    end
+  end
+  def to_gsplot
+    f = ""
+    
+    if ( self[0].kind_of? Array ) then
+      x = self[0]
+      y = self[1]
+      d = self[2]
+
+      x.each_with_index do |xv, i|
+        y.each_with_index do |yv, j|
+          f << [ xv, yv, d[i][j] ].join(" ") << "\n"
+        end
+        # f << "\n"
+      end
+    elsif ( self[0].kind_of? Numeric ) then
+      self.length.times do |i| f << "#{self[i]}\n" end
+    else
+      self[0].zip( *self[1..-1] ).to_gsplot
+    end
+    
+    f
+  end
+end

data/lib/gsl/oper.rb

@@ -0,0 +1,68 @@
+class Fixnum
+  alias :_orig_mul :*
+  alias :_orig_div :/
+  def *(other)
+    if other.kind_of?(GSL::Matrix) or other.kind_of?(GSL::Vector) or other.kind_of?(GSL::Matrix::Int) or other.kind_of?(GSL::Vector::Int) or other.kind_of?(GSL::Vector::Complex)or other.kind_of?(GSL::Matrix::Complex)
+      other.scale(self)
+    else
+      if GSL.have_tensor?
+        if other.kind_of?(GSL::Tensor) or other.kind_of?(GSL::Tensor::Int)
+          other.scale(self)
+        else
+          self._orig_mul(other)
+        end
+      else
+        self._orig_mul(other)
+      end
+    end
+  end
+
+  def /(other)
+    if other.kind_of?(GSL::Poly) or other.kind_of?(GSL::Poly::Int)
+      a = GSL::Poly[1]; a[0] = self
+      GSL::Rational.new(a, other)
+    elsif other.kind_of?(GSL::Vector::Col) 
+      other.scale(1.0/GSL::pow_2(other.dnrm2))
+    elsif other.kind_of?(GSL::Vector::Int::Col)
+      v = other.to_f
+      v.scale(1.0/GSL::pow_2(v.dnrm2))
+    else
+      self._orig_div(other)
+    end
+  end
+end
+
+class Float
+  alias :_orig_mul :*
+  alias :_orig_div :/
+
+  def *(other)
+    if other.kind_of?(GSL::Matrix) or other.kind_of?(GSL::Vector) or other.kind_of?(GSL::Matrix::Int) or other.kind_of?(GSL::Vector::Int) or other.kind_of?(GSL::Vector::Complex)or other.kind_of?(GSL::Matrix::Complex)
+      other.scale(self)
+    else
+      if GSL.have_tensor?
+        if other.kind_of?(GSL::Tensor) or other.kind_of?(GSL::Tensor::Int)
+          other.scale(self)
+        else
+          self._orig_mul(other)
+        end
+      else
+        self._orig_mul(other)
+      end
+    end
+  end
+
+  def /(other)
+    if other.kind_of?(GSL::Poly) or other.kind_of?(GSL::Poly::Int)
+      a = GSL::Poly[1]; a[0] = self
+      GSL::Rational.new(a, other)
+    elsif other.kind_of?(GSL::Vector::Col) 
+      other.scale(1.0/GSL::pow_2(other.dnrm2))
+    elsif other.kind_of?(GSL::Vector::Int::Col)
+      v = other.to_f
+      v.scale(1.0/GSL::pow_2(v.dnrm2))
+    else
+      self._orig_div(other)
+    end
+  end
+end