tioga 1.4

data/split/Dtable/dtable_intern.h

@@ -0,0 +1,144 @@
+/* dtable.h */
+/*
+   Copyright (C) 2004  Bill Paxton
+
+   This file is part of Dtable.
+
+   Dtable is free software; you can redistribute it and/or modify
+   it under the terms of the GNU General Library Public License as published
+   by the Free Software Foundation; either version 2 of the License, or
+   (at your option) any later version.
+
+   Dtable is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU Library General Public License for more details.
+
+   You should have received a copy of the GNU Library General Public License
+   along with Dtable; if not, write to the Free Software
+   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+*/
+
+#ifndef __dtable_H__
+#define __dtable_H__
+
+#include <namespace.h>   
+#include <ruby.h>
+#include "dvector.h"
+
+/*======================================================================*/
+
+
+PUBLIC void Init_Dtable();
+PRIVATE VALUE Read_Dtable(VALUE dest, char *filename, int skip_lines);
+PRIVATE double **Dtable_Ptr(VALUE dtable, long *num_cols, long *num_rows);
+
+PRIVATE bool Is_Dtable(VALUE obj);
+
+PRIVATE VALUE dtable_apply_math_op(VALUE source, double (*op)(double));
+PRIVATE VALUE dtable_apply_math_op2(VALUE ary1, VALUE ary2, double (*op)(double, double));
+PRIVATE VALUE dtable_apply_math_op_bang(VALUE ary, double (*op)(double));
+PRIVATE VALUE dtable_apply_math_op2_bang(VALUE ary1, VALUE ary2, double (*op)(double, double));
+
+PRIVATE VALUE dtable_dup(VALUE ary);
+PRIVATE VALUE dtable_transpose(VALUE ary);
+PRIVATE VALUE dtable_reverse_rows(VALUE ary);
+PRIVATE VALUE dtable_reverse_cols(VALUE ary);
+PRIVATE VALUE dtable_rotate_cw90(VALUE ary);
+PRIVATE VALUE dtable_rotate_ccw90(VALUE ary);
+PRIVATE VALUE dtable_num_cols(VALUE ary);
+PRIVATE VALUE dtable_num_rows(VALUE ary);
+PRIVATE VALUE dtable_min(VALUE ary);
+PRIVATE VALUE dtable_max(VALUE ary);
+PRIVATE VALUE dtable_minmax(VALUE ary);
+PRIVATE VALUE dtable_row(VALUE ary, VALUE row_num);
+PRIVATE VALUE dtable_column(VALUE ary, VALUE column_num);
+PRIVATE VALUE dtable_set_row(VALUE ary, VALUE row_num, VALUE dvec);
+PRIVATE VALUE dtable_set_column(VALUE ary, VALUE col_num, VALUE dvec);
+PRIVATE VALUE dtable_clear(VALUE ary, VALUE val);
+PRIVATE VALUE dtable_set(VALUE ary, VALUE val);
+PRIVATE VALUE dtable_neg(VALUE ary);
+PRIVATE VALUE dtable_abs(VALUE ary);
+PRIVATE VALUE dtable_sin(VALUE ary);
+PRIVATE VALUE dtable_cos(VALUE ary);
+PRIVATE VALUE dtable_tan(VALUE ary);
+PRIVATE VALUE dtable_asin(VALUE ary);
+PRIVATE VALUE dtable_acos(VALUE ary);
+PRIVATE VALUE dtable_atan(VALUE ary);
+PRIVATE VALUE dtable_sinh(VALUE ary);
+PRIVATE VALUE dtable_cosh(VALUE ary);
+PRIVATE VALUE dtable_tanh(VALUE ary);
+PRIVATE VALUE dtable_asinh(VALUE ary);
+PRIVATE VALUE dtable_acosh(VALUE ary);
+PRIVATE VALUE dtable_atanh(VALUE ary);
+PRIVATE VALUE dtable_ceil(VALUE ary);
+PRIVATE VALUE dtable_floor(VALUE ary);
+PRIVATE VALUE dtable_round(VALUE ary);
+PRIVATE VALUE dtable_exp(VALUE ary);
+PRIVATE VALUE dtable_exp10(VALUE ary);
+PRIVATE VALUE dtable_log(VALUE ary);
+PRIVATE VALUE dtable_log10(VALUE ary);
+PRIVATE VALUE dtable_inv(VALUE ary);
+PRIVATE VALUE dtable_sqrt(VALUE ary);
+PRIVATE VALUE dtable_neg_bang(VALUE ary);
+PRIVATE VALUE dtable_abs_bang(VALUE ary);
+PRIVATE VALUE dtable_sin_bang(VALUE ary);
+PRIVATE VALUE dtable_cos_bang(VALUE ary);
+PRIVATE VALUE dtable_tan_bang(VALUE ary);
+PRIVATE VALUE dtable_asin_bang(VALUE ary);
+PRIVATE VALUE dtable_acos_bang(VALUE ary);
+PRIVATE VALUE dtable_atan_bang(VALUE ary);
+PRIVATE VALUE dtable_sinh_bang(VALUE ary);
+PRIVATE VALUE dtable_cosh_bang(VALUE ary);
+PRIVATE VALUE dtable_tanh_bang(VALUE ary);
+PRIVATE VALUE dtable_asinh_bang(VALUE ary);
+PRIVATE VALUE dtable_acosh_bang(VALUE ary);
+PRIVATE VALUE dtable_atanh_bang(VALUE ary);
+PRIVATE VALUE dtable_ceil_bang(VALUE ary);
+PRIVATE VALUE dtable_floor_bang(VALUE ary);
+PRIVATE VALUE dtable_round_bang(VALUE ary);
+PRIVATE VALUE dtable_exp_bang(VALUE ary);
+PRIVATE VALUE dtable_exp10_bang(VALUE ary);
+PRIVATE VALUE dtable_log_bang(VALUE ary);
+PRIVATE VALUE dtable_log10_bang(VALUE ary);
+PRIVATE VALUE dtable_inv_bang(VALUE ary);
+PRIVATE VALUE dtable_sqrt_bang(VALUE ary);
+PRIVATE VALUE dtable_trim(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_log(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_log10(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_inv(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_asin(VALUE ary);
+PRIVATE VALUE dtable_safe_acos(VALUE ary);
+PRIVATE VALUE dtable_safe_sqrt(VALUE ary);
+PRIVATE VALUE dtable_atan2_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_modulo_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_remainder_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_trim_bang(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_pow_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_as_exponent_of_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_safe_log_bang(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_log10_bang(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_inv_bang(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_safe_sqrt_bang(VALUE ary);
+PRIVATE VALUE dtable_safe_asin_bang(VALUE ary);
+PRIVATE VALUE dtable_safe_acos_bang(VALUE ary);
+PRIVATE VALUE dtable_add(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_sub(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_mul(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_div(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_mod(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_remainder(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_pow(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_as_exponent_of(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_atan2(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_add_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_sub_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_mul_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_div_bang(VALUE ary, VALUE arg);
+PRIVATE VALUE dtable_read(int argc, VALUE *argv, VALUE self);
+PRIVATE VALUE dtable_at(VALUE ary, VALUE xloc, VALUE yloc);
+PRIVATE VALUE dtable_aset(VALUE ary, VALUE xloc, VALUE yloc, VALUE val);
+PRIVATE void dtable_store(VALUE ary, long i, long j, double v);
+
+#endif   /* __dtable_H__ */
+

data/split/Dtable/dvector.h

@@ -0,0 +1,61 @@
+/* Dvector.h */
+/*
+   Copyright (C) 2005  Bill Paxton
+
+   Dvector is free software; you can redistribute it and/or modify
+   it under the terms of the GNU General Library Public License as published
+   by the Free Software Foundation; either version 2 of the License, or
+   (at your option) any later version.
+
+   Dvector is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU Library General Public License for more details.
+
+   You should have received a copy of the GNU Library General Public License
+   along with Dvector; if not, write to the Free Software
+   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+*/
+
+#ifndef __Dvector_H__
+#define __Dvector_H__
+
+/* this file has been heavily modified by Vincent Fourmond to take care
+   of the 'RCR330' scheme of exporting symbols 
+*/
+
+#include <symbols.h>
+#include <stdbool.h>
+
+/*======================================================================*/
+
+/* functions for handling Dvectors: */
+
+DECLARE_SYMBOL(double *, Dvector_Data_for_Read, (VALUE dvector, long *len_ptr));
+ /* returns pointer to the dvector's data (which may be shared) */
+DECLARE_SYMBOL(double *, Dvector_Data_Copy, (VALUE dvector, long *len_ptr)); 
+/* like Dvector_Data_for_Read, but returns pointer to a copy of the data */
+DECLARE_SYMBOL(double *, Dvector_Data_for_Write, 
+		(VALUE dvector, long *len_ptr));
+DECLARE_SYMBOL(double *, Dvector_Data_Resize, (VALUE dvector, long new_len));
+DECLARE_SYMBOL(double *, Dvector_Data_Replace,
+	       (VALUE dvector, long len, double *data)); 
+/* copies the data into the dvector */
+DECLARE_SYMBOL(VALUE, Dvector_Create, (void));
+DECLARE_SYMBOL(void, Dvector_Store_Double, (VALUE ary, long idx, double val));
+/* pushes one element onto the vector */
+DECLARE_SYMBOL(void, Dvector_Push_Double, (VALUE ary, double val));
+
+/* functions for interpolation */
+DECLARE_SYMBOL(double, c_dvector_spline_interpolate,
+	       (double x, int n_pts_data, double *Xs, double *Ys, 
+		double *Bs, double *Cs, double *Ds));
+DECLARE_SYMBOL(double, c_dvector_linear_interpolate,
+	       (int num_pts, double *xs, double *ys, double x));
+DECLARE_SYMBOL(void, c_dvector_create_spline_interpolant,
+	       (int n_pts_data, double *Xs, double *Ys,
+		bool start_clamped, double start_slope, 
+		bool end_clamped, double end_slope,
+		double *As, double *Bs, double *Cs));
+#endif   /* __Dvector_H__ */
+

data/split/Dtable/extconf.rb

@@ -0,0 +1,4 @@
+# Dtable installation file
+require 'mkmf'
+
+create_makefile 'Dobjects/Dtable'

data/split/Dtable/include/dtable.h

@@ -0,0 +1,35 @@
+/* dtable.h */
+/*
+   Copyright (C) 2004  Bill Paxton
+
+   This file is part of Dtable.
+
+   Dtable is free software; you can redistribute it and/or modify
+   it under the terms of the GNU General Library Public License as published
+   by the Free Software Foundation; either version 2 of the License, or
+   (at your option) any later version.
+
+   Dtable is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU Library General Public License for more details.
+
+   You should have received a copy of the GNU Library General Public License
+   along with Dtable; if not, write to the Free Software
+   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+*/
+
+#ifndef __dtable_H__
+#define __dtable_H__
+
+#include <symbols.h>
+
+/*======================================================================*/
+
+
+DECLARE_SYMBOL(VALUE, Read_Dtable, 
+	       (VALUE dest, char *filename, int skip_lines));
+DECLARE_SYMBOL(double **, Dtable_Ptr, 
+	       (VALUE dtable, long *num_cols, long *num_rows));
+
+#endif

data/split/Dtable/lib/Dtable_extras.rb

@@ -0,0 +1,90 @@
+#  pre_ddata_extras.rb
+
+# add arithmetic operators for Float and Fixnum with Dtable
+
+class Float
+    alias :pre_ddata_add :+
+    def +(a)
+        if a.class == Dobjects::Dtable
+            a+self
+        else
+            pre_ddata_add(a)
+        end
+    end
+    alias :pre_ddata_sub :-
+    def -(a)
+        if a.class == Dobjects::Dtable
+            a.neg+self
+        else
+            pre_ddata_sub(a)
+        end
+    end
+    alias :pre_ddata_mult :*
+    def *(a)
+        if a.class == Dobjects::Dtable
+            a*self
+        else
+            pre_ddata_mult(a)
+        end
+    end
+    alias :pre_ddata_div :/
+    def /(a)
+        if a.class == Dobjects::Dtable
+            a.inv*self
+        else
+            pre_ddata_div(a)
+        end
+    end
+    alias :pre_ddata_pow :**
+    def **(a)
+        if a.class == Dobjects::Dtable
+            a.as_exponent_of(self)
+        else
+            pre_ddata_pow(a)
+        end
+    end
+end
+
+class Fixnum
+    alias :pre_ddata_add :+
+    def +(a)
+        if a.class == Dobjects::Dtable
+            a+self
+        else
+            pre_ddata_add(a)
+        end
+    end
+    alias :pre_ddata_sub :-
+    def -(a)
+        if a.class == Dobjects::Dtable
+            a.neg+self
+        else
+            pre_ddata_sub(a)
+        end
+    end
+    alias :pre_ddata_mult :*
+    def *(a)
+        if a.class == Dobjects::Dtable
+            a*self
+        else
+            pre_ddata_mult(a)
+        end
+    end
+    alias :pre_ddata_div :/
+    def /(a)
+        if a.class == Dobjects::Dtable
+            a.inv*self
+        else
+            pre_ddata_div(a)
+        end
+    end
+    alias :pre_ddata_pow :**
+    def **(a)
+        if a.class == Dobjects::Dtable
+            a.as_exponent_of(self)
+        else
+            pre_ddata_pow(a)
+        end
+    end
+end
+

data/split/Dtable/namespace.h

@@ -0,0 +1,47 @@
+/* namespace.h: an attempt at rationalizing shared objects
+   namespace use.
+  
+   Copyright (C) 2006  Vincent Fourmond
+
+   This program is free software; you can redistribute it and/or modify
+   it under the terms of the GNU General Library Public License as published
+   by the Free Software Foundation; either version 2 of the License, or
+   (at your option) any later version.
+   
+   This program is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU Library General Public License for more details.
+   
+   You should have received a copy of the GNU Library General Public License
+   along with this program; if not, write to the Free Software
+   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+*/
+
+#ifndef _NAMESPACE_H
+#define _NAMESPACE_H
+
+/* This header file provides two OS-specific macros for the definition of
+   extern symbols:
+
+   * PUBLIC, which has to be used to mark objects that will be used
+   outside the module
+
+   * PRIVATE, for symbols which are "extern" but intern to the module
+
+   Please don't add "extern" after the PRIVATE or PUBLIC declaration
+   as this would break compilation on Darwin.
+*/
+
+#ifdef __APPLE__
+# define PRIVATE __private_extern__
+# define PUBLIC 
+#elif __GNUC__ >= 4 /* we have the visibility attribute */
+# define PRIVATE __attribute__ ((visibility ("hidden"))) 
+# define PUBLIC __attribute__ ((visibility ("default"))) 
+#else /* not really good */
+# define PRIVATE 
+# define PUBLIC 
+#endif /* __APPLE__  and __GNU_C_ >= 4*/
+
+#endif

data/split/Dtable/safe_double.h

@@ -0,0 +1,104 @@
+/**********************************************************************
+
+   safe_double.h: small abstraction for storing double without endianess
+   problems
+   
+   Copyright (C) 2006 Vincent Fourmond
+
+   This program is free software; you can redistribute it and/or modify
+   it under the terms of the GNU General Library Public License as published
+   by the Free Software Foundation; either version 2 of the License, or
+   (at your option) any later version.
+   
+   This program is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU Library General Public License for more details.
+   
+   You should have received a copy of the GNU Library General Public License
+   along with this program; if not, write to the Free Software
+   Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301 USA
+   
+**********************************************************************/
+
+/* This file provides two functions: store_double and get_double, which
+   can be used in a 'safe' way to store doubles  and retrive them
+   in a hopefully platform-independent form. However, it just 
+   stores it without regards to that if ieee754.h isn't found...
+*/
+
+#ifndef _DOUBLE_H
+#define _DOUBLE_H
+
+#define STORE_LOWER_BYTE(a,p) do {\
+*(p++) = (a) & 0xFF; (a) >>= 8; }\
+while(0)
+#define STORE_UNSIGNED(a,p)   for(i = 0; i < 4; i++) STORE_LOWER_BYTE(a,p);
+#define GET_UNSIGNED(a,p) do { a = 0; for(i = 0; i < 4; i++) \
+(a) |= *(p++) << (i * 8); } while (0)
+
+#ifdef HAVE_IEEE754_H
+#include <ieee754.h>
+
+
+static inline void store_double(double a, unsigned char * p)
+{
+  unsigned int tmp;
+  int i;
+  union ieee754_double d;
+  d.d = a;
+  /* we store it with lower bytes firts */
+  tmp = d.ieee.mantissa1;
+  STORE_UNSIGNED(tmp, p);
+  tmp = d.ieee.negative << 31 |
+    d.ieee.exponent << 20 | 
+    d.ieee.mantissa0;
+  STORE_UNSIGNED(tmp, p);
+}
+
+
+static inline double get_double(const unsigned char * p)
+{
+  unsigned int tmp;
+  int i;
+  union ieee754_double d;
+  GET_UNSIGNED(tmp, p);
+  d.ieee.mantissa1 = tmp;
+  GET_UNSIGNED(tmp, p);
+  d.ieee.mantissa0 = tmp & 0xFFFFF;
+  d.ieee.exponent = (tmp >> 20) & 0x7FF;
+  d.ieee.negative = (tmp >> 31) & 0x1;
+  return d.d;
+}
+
+#else
+
+union basic_double{
+  double d;
+  struct {
+    unsigned int a:32;
+    unsigned int b:32;
+  } i;
+};
+
+static inline void store_double(double a, unsigned char * p)
+{
+  union basic_double d;
+  int i;
+  d.d = a;
+  STORE_UNSIGNED(d.i.a, p);
+  STORE_UNSIGNED(d.i.b, p);
+}
+
+static inline double get_double(const unsigned char * p)
+{
+  union basic_double d;
+  int i;
+  GET_UNSIGNED(d.i.a, p);
+  GET_UNSIGNED(d.i.b, p);
+  return d.d;
+}
+
+#endif /* HAVE_IEEE754_H */
+
+#endif /* _DOUBLE_H */