ruby-dcl 1.6.2

data/ChangeLog

@@ -0,0 +1,321 @@
+Thu Feb 24  2011  T Horinouchi
+	* ruby-dcl-1.6.2 released (cvs tag ruby-dcl-1_6_2)
+	  because of a packaging miss of 1.6.1
+Thu Feb 24  2011  T Horinouchi
+	* README: Updated (to leave the installation part to web page)
+Thu Feb 24  2011  T Koshiro
+	* add a dot file .keepme to an empty directory lib/ (for "cvs export")
+Thu Feb 24  2011  T Horinouchi
+	* ruby-dcl-1.6.1 released (cvs tag ruby-dcl-1_6_1)
+Thu Feb 24  2011  T Koshiro
+	* demo/grph1/scpack/scpkt7.rb, demo/grph1/sgpack/sgpk02.rb: debug.
+	* GenWrapper/dcl_narrayed_funcs.rb, GenWrapper/etc/p_header:
+	  correct prototype declarations of math1 functions:
+	  doublereal -> real
+	* GenWrapper/elim_ary_size.rb, GenWrapper/etc/Mk_proto,
+	  GenWrapper/etc/dcl_gen, GenWrapper/proto/Makefile:
+	  change a directory name: src/ -> ../
+	  (because CVS does not support symbolic links.)
+Wed Feb 23  2011  T Koshiro
+	* GenWrapper/proto/grph2_uepack.fp: added a prototype definition of
+	  uetonb; the change is made limited to #if DCLVER >= 534.
+	* demo/grph2/umpack/test10.rb, demo/grph2/udegpk/u2df09[bcef].rb:
+	  created (for DCL 5.3.4).
+	* dcl_obj2cary.c: To use StringValue, RSTRING_PTR, and RSTRING_LEN
+	  macros instead of rb_str2cstr for Ruby 1.8 or later.
+Fri Feb 18  2011  T Horinouchi
+	* ruby-dcl-1.6.1 preview
+Fri Feb 18  2011  T Horinouchi < T Koshiro [dennou-ruby: 003224] (2010-12-11)
+	* ruby-dcl-1.6.0-ruby1_9_2.patch
+	* Support ruby 1.9.2
+	* ruby-dcl-1.6.0-dcl_rb_GenWrapper.patch
+	* Moved a change in ruby-dcl-1.6.0, which was mistakenly
+	  coded directly in lib/dcl.rb, into GenWrapper:
+	  Added a new file dcl_rb_footing. dcl_rb_heading was modified too.
+Thu Jul 17  2008  T Horinouchi
+	* ruby-dcl-1.6.0 released
+Wed Jul 16  2008  T Horinouchi
+	* init.c.gtk,init.c.default : To use RARRAY_LEN and
+	  RARRAY_PTR macros
+Wed Jul 16  2008  T Horinouchi < M Otsuka [dennou-ruby: 002920]
+	* extconf.rb,dcl_obj2cary.c : applied a patch
+Tue Jul 15  2008  T Horinouchi
+	* dcl_obj2cary.c : introduced RARRAY_PTR and RARRAY_LEN
+	* dcl_rubydcloriginal.c : created -- defined uemrkz
+	* lib/dcl.rb : uemrkz support
+	* demo/rubydcloriginal/uemrkz1.rb : created
+	* GenWrapper/dcl_narrayed_funcs.rb : created -- generator
+	  of NArray versions of scalar functions such as stftrf.
+	  ( real,real -> real,real kind of subroutines in
+	  stpack, maplib, ctlib functions ) 
+	  NOTE: the remaining subprograms in stpack and maplib are 
+	  fully covered by *.fp, but prototypes of ctlib is 
+	  not made --- it's easy, you can directly use narray.
+	* dcl_narrayed_funcs.c : Automatically generated by
+	  GenWrapper/dcl_narrayed_funcs.rb -- DO NOT EDIT IT.
+	* GenWrapper/proto/grph1_stpack.fp : commented out prototype
+	  definitions of stftrf stitrf stfpr2 stipr2, which 
+	  are covered by dcl_narrayed_funcs.rb
+	* GenWrapper/proto/math1_maplib.fp : created to cover functions
+	  in MAPLIB, which are not treated by dcl_narrayed_funcs.rb
+	* init.c.default, init.c.gtk : to call the initializer in 
+	  dcl_rubydcloriginal.c and dcl_narrayed_funcs.c
+Mon Jul 14  2008  T Horinouchi < T Koshiro [dennou-ruby:002958]
+	* Debug in swp*
+Tue Dec  5  2006  T Horinouchi
+ 	* ruby-dcl-1.5.3 re-released (sorry!)
+	* init.c.gtk,init.c.default: debug (insufficient alloc length)
+Fri Dec  1  2006  T Horinouchi
+ 	* ruby-dcl-1.5.3 released
+	* extconf.rb, init.c.*: debug for cdcl 5.3.2
+Thu Nov 30  2006  T Horinouchi
+        * Re-generated C source files and dcl.rb
+	  % cd GenWrapper && make
+	* Changed to deal with DCL minor versions such as 5.3.2 etc
+	  in adition to 5.3 etc. Changed files are:
+	  extconf.rb, init.c.gtk, init.c.default,
+	  GenWrapper/proto/grph1_sgpack.fp, GenWrapper/proto/grph1_stpack.fp
+	  GenWrapper/proto/grph1_swpack.fp, GenWrapper/proto/grph2_grpack.fp,
+	  GenWrapper/proto/grph2_umpack.fp, GenWrapper/proto/math1_gt2dlib.fp
+        * Patch ruby-dcl-1.5.2_arg.patch from [dennou-ruby:002548]
+          (by S Nishizawa) was adapted in a very different way
+	  * init.c.gtk, init.c.default: added dcl_init_arg instead of 
+	    to math1_oslib.c. Also, changes are made limited to 
+	    #if DCLVER >= 532.
+	  * GenWrapper/dcl_rb_heading: added 
+	    init_arg($0,ARGV) if DCLVERSION >= "5.3.2"
+Thu Nov 30  2006  T Horinouchi < T Koshiro
+	* GenWrapper/dcl_rb_heading: /cygwin|mingw/ --> /cygwin|mingw|darwin/
+Wed Nov 17  2005  T Horinouchi
+	* ruby-dcl-1.5.2 released
+Mon Nov 14  2005  T Horinouchi < S Nishizawa 
+	(patch, Message-Id: <20051111204153.c17a86f0.seiya@kugi.kyoto-u.ac.jp>)
+	* grph1_zgpack.c: added methods zgqdrw and zgqpmp to enable us to
+	  extract Gtk objects made by DCL).
+Tue Aug  2 17:24:48 JST 2005
+        * demo/rakuraku/kihon/kihonc.rb: debug
+Fri May 13  2005  T Horinouchi
+	* ruby-dcl-1.5.1 released
+Mon Apr 11 2005  Takeshi Horinouchi 
+	* init.c.*: changed global variable mNumRu to file locale.
+Thu Mar 31  2005  T Horinouchi
+	* GenWrapper/etc/p_header: added extern
+Tue Dec 28  2004  T Horinouchi
+	* GenWrapper/dcl_rb_heading: added glib2 requirement for windows 
+	  (cygwin)
+
+Thu Nov 11  2004  T Horinouchi
+	* ruby-dcl-1.5.0 released
+	* extconf.rb: added messaging on DCL version
+	* demo/grph2/umpack/test09.rb,demo/grph2/ulpack/ulpk0[12]n.f,
+	  demo/math1/gt2dlib/gt2d01.rb: created (for DCL 5.3).
+	* demo/grph2/g2pack/g2pk0[12].rb: very minor update2]
+Tue Nov  9  2004  T Horinouchi
+	* Revisions are made to support both dcl-5.2 and dcl-5.3:
+          * demo/grph2/g2pack/g2pk0[12].rb: created (equivalent to
+	    g2pk0[12].f of the Fortran version)
+	  * GenWrapper/proto/:
+	    Makefile: minor improvement;  math1_gt2dlib.fp: created
+	    grph1_stpack.fp, grph1_swpack.fp, grph2_grpack.fp,
+	    grph2_umpack.fp: updated for dcl-5.3
+	  * GenWrapper/pparse.rb: to retain macros starting with /^#/.
+	  * GenWrapper/cproto2init.rb: deal with /^#if *DCLVER/ .. /^#endif/
+	  * GenWrapper/{variable.rb,prototype.rb,cproto2init.rb}: 
+	    updated for ruby 1.8 (and 1.6 as before).
+	  * extconf.rb: to add -DDCLVER=52 (or 53 ..) in $CFLAGS
+	  * init.c.{default|gtk}: added new constant NumRu::DCL::DCLVERSION
+	    which is either "5.2" or "5.3", currently.
+Thu Aug 04  2004  T Horinouchi < D Tsukahara
+	* lib/dcl.rb: debug for SHTLIB (Spherical Harmonics Transforpmation)
+	* GenWrapper/proto/math2_shtlib.fp: debug for SHTLIB
+Thu Jul 15  2004  T Horinouchi < T Koshiro
+	* ruby-dcl-1.4.1 released
+	* lib/dcl.rb: dependent library updated only for cygwin & mingw.
+Tue Jul 13  2004  T Horinouchi < T Koshiro
+	* extconf.rb: further minor update for cygwin & mingw.
+Tue Jun 29  2004  T Horinouchi < T Koshiro
+	* extconf.rb: debug for cygwin & mingw.
+Fri Mar  5  2004  T Horinouchi < S Nishizawa
+	* ruby-dcl-1.4.0 released
+	* extconf.rb,grph1_zgpack.c.org: modified to support gtk2
+	  (dcl-5.2-C.3-gtk2-5.tar.gz).
+	  gtk1 support stil remains, because dcl-5.2-C.3 supports gtk1.
+Mon Jan 19  2004  T Horinouchi
+	* ruby-dcl-1.3.2 released
+	* Debug in GenWrapper/proto/math2_shtlib.fp
+	  --> cd GenWrapper & make; Also did make dcl_rb and editting
+	      of lib/dcl.rb by hand. Especially needed was to avoid
+	      size checking of the "work" or "q" array initialized
+	      by shtint because not all info on its size is provided
+	      as array arguments (im, jm, or mm).
+Wed Dec 10  2003  T Horinouchi < T Koshiro
+	* ruby-dcl-1.3.1 released
+	* dcl_obj2cary.c: minor debug for Cygwin (EXTERN -> extern)
+Wed Oct  1  2003  T Horinouchi < S Nishizawa
+	* demo/rakuraku/map3d/map3d1.rb: debug 
+Tue Sep 23 2003  Takeshi Horinouchi < T Koshiro
+	* ruby-dcl-1.3.0 released
+	* dcl_obj2cary.c: support NArrayMiss
+	* lib/dcl.rb: modified not to exclude NArrayMiss
+	* GenWrapper/variable.c: modified not to exclude NArrayMiss
+	* GenWrapper/old/: eliminated
+	* extconf.rb: debug for windows (by T Koshiro)
+Wed Aug 27 2003  Takeshi Horinouchi
+	* ruby-dcl-1.2.2 released
+	* Verified that demo programs work fine with Ruby 1.8.0. (not all
+	  but so far all work fine)
+	* lib/dcl.rb: applied the gotoken patch (to eliminate unneccesary 
+	  spaces)
+	* init.c.default,init.c.gtk: Init_dcl: No longer asks whether NumRu
+	  has been defined -- ok for Ruby 1.6.5 or so or later
+	* extconf.rb: do not edit "install:" anymore (for Ruby 1.8; works
+	  with 1.6 too.)
+	
+Tue Jun 10  2003  T Horinouchi
+	* ruby-dcl-1.2.1 released
+
+Tue Jun 10  2003  T Horinouchi thanks to S Nishizawa
+	* gtk handling by S Nishizawa was merged in.
+	  added: ./grph1_zgpack.c.org ./init.c.gtk
+	  modified: extconf.rb
+	  (these are based on ruby-dcl-1.2.0-gtk-5.patch with some
+	  modification --- file handling in extconf.rb and extern->EXTERN
+	  in init.c.gtk)
+
+Tue Jun 10  2003  T Horinouchi
+	* extconf.rb: added have_library("narray") (for windows)
+	* init.c: extern -> EXTERN (for windows)
+	* lib/dcl.rb: added (gl|sg|us|..)p(set|stx|get)
+	
+Sun Jun  8  2003  T Horinouchi
+	* directory structure reorganization -- to re-include the wrapper
+	  generator etc in the tar ball. This time, the source generator
+	  is put in the subdirectory GenWrapper, and ordinary users do not
+	  have to generate the C sources.
+	* demo/ minor debug to add "DCL::" to where forgotten
+	* GenWrapper/proto/grph2_udpack.fp activated udlabl
+
+Tue May 13  2003  T Horinouchi
+	* ruby-dcl-1.2.0 released
+	* directory structure reorganization: separeted the C source
+	  generation from the distributed package.
+	   - renamed src/ as ruby-dcl-1.1.3/
+	   - move some files into (the old) src: 
+	     ChangeLog (this file)
+	   - Wrote the new README file.`
+	
+Tue May 21  2002  T Horinouchi
+	* ruby-dcl-1.1.2 released
+	* elim_ary_size.rb: make old methods (such as __udcntr) private
+	* demo:  minor refinement of some demo programs
+
+Mon May 20  2002  S Nishizawa
+	* demo/*/*/*.rb: include only NumRu instead of NumRu::DCL
+
+Mon May 20  2002  T Horinouchi
+	* demo/rakuraku/map3d/map3d7.rb: made from map3d6.rb (below)
+	  -- wire frame with painting.  great!
+	* demo/rakuraku/map3d/map3d6.rb: written based on map3d6.f
+	  in dcl-5.2 (FORTRAN). Wire frame!
+
+Mon May 20  2002  T Horinouchi
+	* ruby-dcl-1.1.1 released
+	* elim_ary_size.rb: debug: added calls of "module_fucntion" to the
+	  top of the source file to generate (src/lib/dcl.rb).
+	* demo/gokuraku/*/*.rb: modified to use array expressions
+	  (note: u2d/u2d[14].rb were rewritten only half way because of
+	  a bug in NArray to support a[0..2,0..3]=b[0..2,1] type
+	  substitution)
+	
+Sun May 19  2002  T Horinouchi
+	* demo/*/*/*.rb: downcased most constants (downcase_demos.rb)
+	
+Sat May 18  2002  T Horinouchi
+	* ruby-dcl-1.1.0 released
+	* src/dcl_cary2obj.c, variable.rb, etc/p_header:
+	  modified to return NArray instead of Array
+	  (it has long been supposed to do so, but it never worked)
+	
+Fri May 17  2002  T Horinouchi
+	* proto/misc1_(date|time)lib added protoypes of C functions,
+	  but they are commented out so far. (anyway, these two libraries
+	  are not very needed since RUby has date and time libs)
+
+Mon May 13-15  2002  T Horinouchi
+	* MAJOR CHANGE: Arguments to represent array sizes are now
+	  eliminated from all the methods. This is done as follows:
+	  the original ruby-dcl library is renamed from dcl to dcl_raw,
+	  and a new ruby source code, dcl.rb, is generated from the
+	  prototype definitions in proto/*.fp by using the filter
+	  elim_ary_size.rb (new).
+	* src/extconf.rb: rename dcl -> dcl_raw (dcl.so -> dcl_raw.so)
+	                  minor change to better handle the numru
+			  directory
+	* src/init.c: rename Init_dcl -> Init_dcl_raw
+	* elim_ary_size.rb: created (filter to eliminate array size arguments)
+	* Makefile
+	
+Thu Dec 28  2001  T Horinouchi
+	* ruby-dcl-1.0.0 released
+	* proto_old/,demo_new_notworking/: moved to old/
+	* test/u2df01.rb: eliminated, since it is obsolete
+	* README: install instruction updated
+	* src/extconf.rb: set $sitearchdir as a default include path for
+	NArray headers. Also, made possible to alter it with
+	--with-narray-include=
+Thu Dec 20  2001  T Horinouchi
+	* proto/*.fp: Commented out *P(SET|GET|PSTX) because the parameter
+	type they accept is not unique. Use instead *[RLC](SET|GET|STX)
+	(automatically done with ~horinout/davis/scripts/dclproto_comout.rb)
+
+Thu Dec 13  2001  T Horinouchi
+	* src/init.c: the Dec 06 modification was made more robust
+
+Thu Dec 06  2001  T Horinouchi
+	* ruby-dcl-1.beta.7 released
+	* src/init.c: rename mNum as mNumRu; initialize it only at first time
+
+Thu Nov 06  2001  Naoki Kawanabe
+        
+	* ruby-dcl-1.beta.6 released
+	* modified bug in make
+	
+Thu Nov 03  2001  Naoki Kawanabe & GOTO kentaro(gotoken@notwork.org)
+	* ruby-dcl-1.beta.5 released
+	* pparse.rb: does not depend on cpp(1).
+	* install directry changed (require 'dcl'-> require 'numru/dcl')
+	* added module 'NumRu'
+	* changed module name (DCL -> NumRu::DCL)
+Thu Oct 23  2001  Takeshi Horinouchi
+
+	* ruby-dcl-1.beta.4 released
+	* src/extconf.rb: now install: is made equiv to site-install: in
+	  src/Makefile (as well as in ./Makefile). Plus, the program was a 
+	  bit tidied up (no essential change).
+	* Makefile,proto/Makefile: add target distclean
+	* rename: demo_new --> demo_new_notworking
+	          proto --> proto_old
+		  proto_new --> proto
+	  edited Makefile etc accordingly (esp. for proto_new proto)
+	* rbdcl is renamed to RubyDCL
+	
+Thu Sep 28  2001  Naoki Kawanbe
+
+	* rbdcl-1.beta.3 released
+	* dcl_obj2cary.c,dcl_cary2obj.c: include libtinyf2c.h insead of f2c.h
+	
+Thu Sep 13  2001  Takeshi Horinouchi
+
+	* rbdcl-1.beta.2 released
+
+Thu Sep 12  2001  Takeshi Horinouchi
+
+	* updated: README
+
+	* debug: src/dcl_obj2cary.c
+	  (integer coarcing / messages)
+
+Thu Sep 12  2001  Takeshi Horinouchi
+
+	* rbdcl-1.beta.1 released
+

data/GenWrapper/Makefile

@@ -0,0 +1,14 @@
+all: source dcl_rb
+
+source:
+	cd proto ; \
+	$(MAKE) ; \
+	cd ../
+
+dcl_rb:
+	ruby elim_ary_size.rb proto/*.fp
+
+clean: clean.local
+
+clean.local:
+	$(RM) *~

data/GenWrapper/cproto2init.rb

@@ -0,0 +1,41 @@
+# Program: cproto2init.rb
+#          2000/11/20 Takeshi Horinouchi
+# Usage:
+#
+#     % ruby cproto2init.rb cprotofile
+#
+#  Here, cprotofile is the name of the C program obtained by,
+#  for example,
+#
+#     ruby proto2c.rb grph2_grpack.fp > grph2_grpack.c
+#
+#  The package name is determined from the file name such as
+#  'grph2_grpack'.
+
+package =  $*[0].gsub(/.*\//,'').sub(/\..*/,'')    # remove directory name
+                                                  # and suffix
+print("void\n","init_",package,"(mDCL)\n","VALUE mDCL;\n","{\n")
+
+macro_dclver = false 
+
+while (line=gets)
+  case line
+  when /^static VALUE/
+    funcdef = gets.chop
+    func = funcdef.sub(/\(.*$/,'')
+    nargs = funcdef.delete('^,').length
+    method = func.sub(/dcl_/,'')
+    print('    rb_define_module_function(mDCL, "',method,
+	   '", ',func,", ",nargs,");\n")
+  when  /^#if *DCLVER/
+    macro_dclver = true
+    print line
+  when /^#endif/
+    if macro_dclver
+      print line
+      macro_dclver = false
+    end
+  end
+end
+
+print "}\n"

data/GenWrapper/dcl_narrayed_funcs.rb

@@ -0,0 +1,181 @@
+# -*- C -*-
+# USAGE: ruby dcl_narray_funcs.rb
+# --> creates ../dcl_narray_funcs.c
+
+require "erb"
+
+erbsrc = <<EOS
+<% rr2rr_funcs = %w[ stftrf stitrf stfpr2 stipr2 mpfcyl mpicyl mpfmer mpimer mpfmwd mpimwd mpfmwl mpimwl mpfhmr mpihmr mpfek6 mpiek6 mpfktd mpiktd mpfcon mpicon mpfcoa mpicoa mpfcoc mpicoc mpfbon mpibon mpfotg mpiotg mpfpst mpipst mpfazm mpiazm mpfaza mpiaza ct2pc ct2cp ct2ec ct2bc ct2hc ct2ch ]
+%>
+/*
+ * Original methods of Ruby-DCL
+ *
+ * Alert: DO NOT EDIT THIS FILE!!!
+ *   This file was automatically generated by GenWrapper/dcl_narray_funcs.rb 
+ */
+
+#include <stdio.h>
+#include "ruby.h"
+#include "libtinyf2c.h"
+#include "narray.h"
+
+/* for compatibility with ruby 1.6 */
+#ifndef RARRAY_LEN
+#define RARRAY_LEN(a) (RARRAY(a)->len)
+#endif
+
+#define DFLT_SIZE 32
+
+extern char    *dcl_obj2ccharary(VALUE, int, int);
+extern integer *dcl_obj2cintegerary(VALUE);
+extern real    *dcl_obj2crealary(VALUE);
+extern complex *dcl_obj2ccomplexary(VALUE);
+extern logical *dcl_obj2clogicalary(VALUE);
+
+extern VALUE dcl_ccharary2obj(char *, int, int);
+extern VALUE dcl_cintegerary2obj(integer *, int, int, int *);
+extern VALUE dcl_crealary2obj(real *, int, int, int *);
+extern VALUE dcl_ccomplexary2obj(complex *, int, char *);
+extern VALUE dcl_clogicalary2obj(logical *, int, int, int *);
+
+extern void dcl_freeccharary(char *);
+extern void dcl_freecintegerary(integer *);
+extern void dcl_freecrealary(real *);
+extern void dcl_freeccomplexary(complex *);
+extern void dcl_freeclogicalary(logical *);
+
+/* for functions which return real */
+/* fnclib */
+extern real rd2r_(real *);
+extern real rr2d_(real *);
+extern real rexp_(real *, integer *, integer *);
+extern real rfpi_(void);
+extern real rmod_(real *, real *);
+/* gnmlib */
+extern real rgnlt_(real *);
+extern real rgnle_(real *);
+extern real rgngt_(real *);
+extern real rgnge_(real *);
+/* rfalib */
+extern real rmax_(real *, integer *, integer *);
+extern real rmin_(real *, integer *, integer *);
+extern real rsum_(real *, integer *, integer *);
+extern real rave_(real *, integer *, integer *);
+extern real rvar_(real *, integer *, integer *);
+extern real rstd_(real *, integer *, integer *);
+extern real rrms_(real *, integer *, integer *);
+extern real ramp_(real *, integer *, integer *);
+/* rfblib */
+extern real rprd_(real *, real *, integer *, integer *, integer *);
+extern real rcov_(real *, real *, integer *, integer *, integer *);
+extern real rcor_(real *, real *, integer *, integer *, integer *);
+
+<% rr2rr_funcs.each do |func| %>
+extern int <%= func%>_(real*,real*,real*,real*);
+<% end %>
+
+extern VALUE mDCL;
+
+/* driver of real*2->real*2 functions such as STFTRF
+ */
+static VALUE
+dcl_rr2rr_num(obj, ux, uy, func)
+    VALUE obj, ux, uy;
+    int (*func)(real*,real*,real*,real*);
+{
+    real i_ux;
+    real i_uy;
+    real o_vx;
+    real o_vy;
+    VALUE vx;
+    VALUE vy;
+
+    if (TYPE(ux) != T_FLOAT) {
+      ux = rb_funcall(ux, rb_intern("to_f"), 0);
+    }
+    if (TYPE(uy) != T_FLOAT) {
+      uy = rb_funcall(uy, rb_intern("to_f"), 0);
+    }
+
+    i_ux = (real)NUM2DBL(ux);
+    i_uy = (real)NUM2DBL(uy);
+
+
+    (*func)(&i_ux, &i_uy, &o_vx, &o_vy);
+
+    vx = rb_float_new((double)o_vx);
+    vy = rb_float_new((double)o_vy);
+
+
+    return rb_ary_new3(2, vx, vy);
+
+}
+
+/* driver of real*2->real*2 functions such as STFTRF
+ */
+static VALUE
+dcl_rr2rr_na(obj, p, q, func)
+    VALUE obj, p, q;
+    int (*func)(real*,real*,real*,real*);
+{
+    VALUE r, s;  // return value
+    real *i_p, *i_q;
+    real *i_r, *i_s;
+    size_t len,j,lenq;
+
+    i_p = dcl_obj2crealary(p);
+    i_q = dcl_obj2crealary(q);
+    len  = (TYPE(p)==T_ARRAY) ? RARRAY_LEN(p) : NA_TOTAL(p) ;
+    lenq = (TYPE(q)==T_ARRAY) ? RARRAY_LEN(q) : NA_TOTAL(q) ;
+    if (len != lenq) {
+         rb_raise(rb_eArgError, "lengths of the 1st & 2nd args do not agree");
+    }
+
+    r = na_make_object(NA_SFLOAT, 1, &len, cNArray);
+    s = na_make_object(NA_SFLOAT, 1, &len, cNArray);
+    i_r = NA_PTR_TYPE(r,real*);
+    i_s = NA_PTR_TYPE(s,real*);
+
+    for(j=0; j<len; j++){
+	(*func)(i_p+j, i_q+j, i_r+j, i_s+j);
+    }
+
+    return rb_ary_new3(2, r, s);
+}
+
+static VALUE
+dcl_rr2rr(obj, p, q, func)
+    VALUE obj, p, q;
+    int (*func)(real*,real*,real*,real*);
+{
+    switch (TYPE(p)){
+    case T_DATA:
+    case T_ARRAY:
+	return( dcl_rr2rr_na(obj, p, q, func) );
+    default:
+	return( dcl_rr2rr_num(obj, p, q, func) );
+    }
+}
+
+
+<% rr2rr_funcs.each do |func| %>
+static VALUE
+dcl_<%= func%>(obj, p, q)
+    VALUE obj, p, q;
+{   return( dcl_rr2rr(obj, p, q, &<%= func%>_) );  }
+<% end %>
+
+void
+init_narrayed_funcs(mDCL)
+VALUE mDCL;
+{
+<% rr2rr_funcs.each do |func| %>
+    rb_define_module_function(mDCL, "<%= func%>", dcl_<%= func%>, 2);
+<% end %>
+}
+
+EOS
+
+File.open('../dcl_narrayed_funcs.c','w'){|f| f.print(ERB.new(erbsrc).result)}
+
+