carray 1.3.7 → 1.4.0

data/ext/fortio/ruby_fortio.c

@@ -1,182 +0,0 @@
-/* ---------------------------------------------------------------------------
-
-  carray/carray_fortio.c
-
-  This file is part of Ruby/CArray extension library.
-  You can redistribute it and/or modify it under the terms of
-  the Ruby Licence.
-
-  Copyright (C) 2005 Hiroki Motoyoshi
-
----------------------------------------------------------------------------- */
-
-#include "ruby.h"
-#include <math.h>
-
-static int
-read_F (char *buf, int buflen, int scale, int length, int prec, double *a)
-{
-  char fmt[128];
-  int status;
-  if ( ! strchr(buf, '.') ) {
-    int tmp = 0;
-    snprintf(fmt, 128, "%%%i", length);
-    status = sscanf(buf, fmt, &tmp);
-    *a = tmp * pow(10, -prec);
-  }
-  else {
-    snprintf(fmt, 128, "%%%i.%if", length, prec);     
-    status = sscanf(buf, fmt, a);
-  }
-  return status;
-}
-
-static VALUE
-rb_ff_read_F (VALUE mod, VALUE vbuffer, VALUE vscale, VALUE vlength, VALUE vprec)
-{
-  double val;
-  int status;
-  status = read_F(StringValuePtr(vbuffer), (int) RSTRING_LEN(vbuffer),
-                  NUM2INT(vscale), NUM2INT(vlength), NUM2INT(vprec), &val);
-  if ( ! status ) {
-    rb_raise(rb_eRuntimeError, "invalid string for F descriptor");
-  }
-  return rb_float_new(val);
-}
-
-static int
-write_F (char *buf, int buflen, int plus, int scale, int length, int prec, double a)
-{
-  char fmt[128];
-  if ( plus ) {
-    snprintf(fmt, 128, "%%+%i.%if", length, prec);
-  }
-  else {
-    snprintf(fmt, 128, "%%%i.%if", length, prec);
-  }
-  if ( scale == 0 ) {
-    return snprintf(buf, buflen, fmt, a);
-  }
-  else {
-    return snprintf(buf, buflen, fmt, a * pow(10, scale));
-  }
-}
-
-static VALUE
-rb_ff_write_F (VALUE mod, VALUE vsign, VALUE vscale, VALUE vlength, VALUE vprec, VALUE va)
-{
-  int length = NUM2INT(vlength);
-  char buf[256];
-  if ( length > 255 ) {
-    rb_raise(rb_eRuntimeError, "too long decimal format for FortranFormat");
-  }
-  write_F(buf, 255, 
-             RTEST(vsign), NUM2INT(vscale), length, NUM2INT(vprec),  NUM2DBL(va));
-  return rb_str_new2(buf);
-}
-
-
-static int
-write_E (char *buf, int buflen, int plus, int scale, int length, int prec, int iexp, double a)
-{
-  char fmt[128];
-  int prec0;
-  int e, sign = 1;
-  prec0 = ( 1 - scale < 0 ) ? prec - scale + 1 : prec;
-  sign = 1;
-  if ( a < 0 ) {
-    e = floor(log10(-a))+1;
-    sign = -1;
-  }
-  if ( a > 0 ) {
-    e = floor(log10(a))+1;      
-  }
-  else {
-    e = 0;
-  }
-  if ( prec0 == 0 ) {
-    if ( plus ) {
-      snprintf(fmt, 128, "%%+%i.0f.E%%+0%ii", length-iexp-3, iexp+1);
-    }
-    else {
-      snprintf(fmt, 128, "%%%i.0f.E%%+0%ii", length-iexp-3, iexp+1);
-    }
-  }
-  else {
-    if ( plus ) {
-      snprintf(fmt, 128, "%%+%i.%ifE%%+0%ii", length-iexp-2, prec0, iexp+1);
-    }
-    else {
-      snprintf(fmt, 128, "%%%i.%ifE%%+0%ii", length-iexp-2, prec0, iexp+1);
-    }
-  }
-  return snprintf(buf, buflen, fmt, sign * a * pow(10,scale-e), e-scale);
-}
-
-static VALUE
-rb_ff_write_E (VALUE mod, VALUE vsign, VALUE vscale, VALUE vlength, VALUE vprec, VALUE vexp, VALUE va)
-{
-  int length = NUM2INT(vlength);
-  int iexp   = NIL_P(vexp) ? 2 : NUM2INT(vexp);
-  char buf[256];
-  if ( length > 255 ) {
-    rb_raise(rb_eRuntimeError, "too long decimal format for FortranFormat");
-  }
-  write_E(buf, 255, 
-             RTEST(vsign), NUM2INT(vscale), length, NUM2INT(vprec), iexp, NUM2DBL(va));
-  return rb_str_new2(buf);
-}
-
-static int
-write_GF (char *buf, int buflen, int plus, int length, int prec, int iexp, double a)
-{
-  int e;
-  if ( a < 0 ) {
-    e = floor(log10(-a))+1;
-  }
-  if ( a > 0 ) {
-    e = floor(log10(a))+1;      
-  }
-  else {
-    e = 0;
-  }
-  write_E(buf, 255, plus, e, length, prec-1, iexp, a);
-  memset(buf+(length-iexp-2), ' ', iexp+2);
-  return 1;
-}
-
-static VALUE
-rb_ff_write_G (VALUE mod, VALUE vsign, VALUE vscale, VALUE vlength, VALUE vprec, VALUE vexp, VALUE va)
-{
-  int length = NUM2INT(vlength);
-  int iprec  = NUM2INT(vprec);
-  int iexp   = NIL_P(vexp) ? 2 : NUM2INT(vexp);
-  double val = NUM2DBL(va);
-  char buf[256];
-  if ( length > 255 ) {
-    rb_raise(rb_eRuntimeError, "too long decimal format for FortranFormat");
-  }
-  if ( ( fabs(val) <= 0.1 ) ||
-       ( fabs(val) >= pow(10, iprec) ) ) {
-    write_E(buf, 255, 
-            RTEST(vsign), NUM2INT(vscale), length, iprec, iexp, val);
-  }
-  else {
-    write_GF(buf, 255, RTEST(vsign), length, iprec, iexp, val);
-  }
-  return rb_str_new2(buf);
-}
-
-void
-Init_fortio_ext ()
-{
-  VALUE rb_cFF = rb_define_class("FortranFormat", rb_cObject);
-  rb_define_singleton_method(rb_cFF, "read_F",  rb_ff_read_F,  4);
-  rb_define_singleton_method(rb_cFF, "write_F", rb_ff_write_F, 5);
-  rb_define_singleton_method(rb_cFF, "write_E", rb_ff_write_E, 6);
-  rb_define_singleton_method(rb_cFF, "write_G", rb_ff_write_G, 6);
-}
-
-
-
-

data/ext/fortio/test/test_H.rb

@@ -1,5 +0,0 @@
-require "fortio"
-
-p FortranFormat.new("((5(5Hhello)$))")
-
-puts fortran_format("(2(5(5Hhello))$)")

data/ext/fortio/test/test_T.rb

@@ -1,7 +0,0 @@
-require "fortio"
-
-fortran_format_write(STDOUT, "(TR12,A1)", "a")
-fortran_format_write(STDOUT, "(TR12,A1,TL5,A1,TL3,A1)", "a", "b", "c")
-fortran_format_write(STDOUT, "(TR12,A1,T5,A1,T3,A1)", "a", "b", "c")
-
-p FortranFormat.new("(A12)")

data/ext/fortio/test/test_fortran_format.rb

@@ -1,86 +0,0 @@
-require 'fortio'
-
-include Math
-
-#F77 = "ifort"
-F77 = "gfortran"
-#F77 = "g95"
-#F77 = "g77"
-
-def f77_compile_and_run (source)
-  output = nil
-  begin
-    open("test.f", "w") { |io|
-      io.write(source)
-    }
-    system("#{F77} test.f")
-    output = `./a.out`.chomp
-  ensure
-    File.unlink("test.f")
-    File.unlink("a.out")
-  end
-end
-
-[
- ["F12.3",    10.0,        nil],
- ["2F12.3",   [10.0,11.0], nil],
- ["E12.2",    10.0,        nil],
- ["2E12.2",   [10.0,11.0], nil],
- ["2PE12.2",  10.0,        nil],
- ["2P2E12.2", [10.0,11.0], nil],
- ["2G12.2",   [10.0,11.0], nil],
- ["2PG12.2",  10.0,        nil],
- ["2PG12.2E3",10.0,        nil],
- ["2PG12.2",  1234567890.0,        nil],
- ["2PG12.2E3",1234567890.0,        nil],
- # ---
- ["F15.8",    PI, "#{PI}D0"],
- ["E15.8",    PI, "#{PI}D0"],
- ["2PE15.8",  PI, "#{PI}D0"],
- # ---
- ["F10.8",    PI, "#{PI}D0"],
- ["E10.8",    PI, "#{PI}D0"],
- ["2PE10.8",  PI, "#{PI}D0"],
- # ---
- ["I4.3",     123456,     nil],
- ["I6.5",     123,        nil],
- ["I6.5",     123,        nil],
- ["I7.6",     123456,     nil],
- ["I6",       123456,     nil],
- # ---
- ["SPI6",      123456,    nil],
- ["SPI7",      123456,    nil],
- ["SPF12.4",   123456.0,  nil],
- ["SPE12.4E3", 123456.0,  nil],
- ["SPG12.4E3", 123456.0,  nil],
- ["SPG12.4E3", 123.0,     nil],
- ["TR2SPSG12.4E3", 123.0, nil],
- # ---
- ["L1",        true,     ".TRUE."],
- ["L1",        false,    ".FALSE."],
- ["L2",        true,     ".TRUE."],
- ["L2",        false,    ".FALSE."],
- ["L4",        true,     ".TRUE."],
- ["L4",        false,    ".FALSE."],
- # ---
- ["A1",        "a",     "'a'"],
- ["A2",        "ab",    "'ab'"],
- ["5Hhello",   nil,  nil],
- # ---
- ["/",        "\n",     ""],
- ["2/",        "\n\n",     ""],
-].each do |fmt, rval, fval|
-  fort_out = f77_compile_and_run %{
-      write(6, '(#{fmt})') #{[fval||rval||[]].flatten.join(',')}
-      end
-  }
-  ruby_out = fortran_format(fmt+"$", *rval)
-  puts "-------------------------------------------------------------"
-  printf("format         : %s\n", fmt)
-  printf("fortran input  : %s\n", [fval||rval].flatten.join(','))
-  printf("ruby input     : %s\n", [rval].flatten.join(","))
-  printf("fortran output : %s\n", fort_out.inspect)
-  printf("ruby output    : %s\n", ruby_out.inspect)
-  printf("comparison     : %s\n", fort_out == ruby_out)
-end
-

data/ext/fortio/test/test_namelist.rb

@@ -1,25 +0,0 @@
-require "fortio"
-require "pp"
-
-text =<<END
-aaaa bbbb cccc
-dddd eeee ffff
-&name
- c=30.d0,
- a=10.d0,
- b=20.d0,
- d(1:3) = 1,2,3
- &end
-&foo
- c=10.d+32,
- a=50.d0,
- b=90.d0,
- d(1:5) = 3*4,,3
- e="hello""world!"
- &end
-END
-
-pp nml = fortran_namelist_read(text, "name", {:d=>CArray.int(3)})
-pp nml = fortran_namelist_read(text, "foo", {:d=>CArray.int(5)})
-
-puts fortran_namelist("bar", nml)

data/ext/fortio/test/test_namelist_write.rb

@@ -1,10 +0,0 @@
-require "fortio"
-
-a = {
-  :array => [1,2,3],
-  :x => 1,
-  :y => 2.2
-}
-
-fortran_namelist_write(STDOUT, "a", a)
-fortran_namelist_write(STDOUT, "a", a)

data/ext/fortio/test/test_sequential.rb

@@ -1,13 +0,0 @@
-require "fortio"
-
-fs = FortranSequential.open("test.dat", "w", :endian=>"big")
-rec = fs.record
-rec.write("d4l", [1.2, 2008, 6, 21, 9])
-fs.close
-
-fs = FortranSequential.open("test.dat", "r", :endian=>"big")
-rec = fs.record
-p rec.read("d4l")
-fs.close
-
-File.unlink("test.dat")

data/ext/fortio/test/test_sequential2.rb

@@ -1,13 +0,0 @@
-require "fortio"
-
-a = CArray.float(3).seq!
-
-fortran_sequential_open("test.dat", "w") { |io|
-  io.write(a)
-}
-
-fortran_sequential_open("test.dat") { |io|
-  p io.read.unpack("e3")
-}
-
-File.unlink("test.dat")

data/ext/fortio/work/test.rb

@@ -1,10 +0,0 @@
-require "carray"
-
-text =<<END
-  1  2  3  4
-  5  6  7  8
-END
-
-ca = CArray.from_fortran_format("2I3,2F3.0/2F3.0,2I3", text)
-
-ca.to_fortran_format("2I3,1P,2E10.2E3", STDOUT)

data/ext/fortio/work/test_e.rb

@@ -1,19 +0,0 @@
-require "fortio"
-
-p f = FortranFormat("4(2P,3E6.2E4,2I2.1)$")
-p f.count_args
-p fortran_format("2P,E10.2E4,$", 0.000012345)
-p f = FortranFormat("2(E10.2E4,',')")
-p fortran_format("2(E10.2E4,',')", 0.000012345, 2.0)
-p f = FortranFormat('2(-2PG15.5 4X2(I2)/)')
-puts fortran_format(f, 3.1415926,2,4, 2.71828,3,5)
-
-p f = FortranFormat('3P,4(G15.5/)')
-p f.count_args
-puts fortran_format(f, 0.1, 0.11, 10**5-1, 10**5)
-
-p f = FortranFormat('A10,"""\'|\'"""')
-p f.count_args
-puts fortran_format(f, "hello")
-
-puts fortran_format_read("12345", "F5.2")

data/ext/fortio/work/test_ep.rb

@@ -1,10 +0,0 @@
-require "fortio"
-
-fmt200 = FortranFormat.new("(3pe13.6,-2p2e13.6)")
-
-p fortran_format_write("", fmt200, *[0.000123456]*3)
-
-p text = fortran_format_write("", "3pf10.5", 3)
-p fortran_format_read(text, "3pf10.5")
-
-p fortran_format_write("", "g13.5", 893)

data/ext/fortio/work/test_parse.rb

@@ -1,12 +0,0 @@
-require "carray"
-require "fortio"
-
-head = "    1440   721    0.00   90.00  0.2500 31 0.18712E+00 0.44215E-01 0.35660E-01 0.30406E-01 0.34569E-01 0.17713E-01 0.25666E-01 0.26762E-01 0.10912E-01 0.10623E-01 0.62775E-02 0.66029E-02 0.11870E+00 0.47465E-02 0.86580E-01 0.33813E-02 0.10799E+00 0.22265E-02 0.76236E-01 0.65007E-01 0.51308E-01 0.58308E-01 0.43458E-01 0.21357E-01 0.89651E-02 0.10000E-01 0.10000E-01 0.10000E-01 0.10000E-01 0.10000E-02 0.10000E+01"
-
-ffp = FortranFormat.new('2X,2i6,2f8.2,f8.4,i3,31e12.5')
-p list = ffp.read(head)
-ffp.write(test="", *list)
-
-p head 
-p test
-p head == test

data/ext/mathfunc/lib/autoload/autoload_math_mathfunc.rb

@@ -1 +0,0 @@
-require "carray/math/mathfunc"

data/ext/mathfunc/lib/math/mathfunc.rb

@@ -1,15 +0,0 @@
-# ----------------------------------------------------------------------------
-#
-#  lib/carray/graphics/imagemap.rb
-#
-#  This file is part of Ruby/CArray extension library.
-#  You can redistribute it and/or modify it under the terms of
-#  the Ruby Licence.
-#
-#  Copyright (C) 2005 Hiroki Motoyoshi
-#
-# ----------------------------------------------------------------------------
-
-require "carray/carray_mathfunc"
-
-

data/ext/narray/README

@@ -1,22 +0,0 @@
-
-This directory contains the plugin library used as the interface to NArray
-and NArrayMiss library from Ruby/CArray. 
-
----
-NArray Library (http://narray.rubyforge.org/) by Masahiro Tanaka.
-
-  NArray is an Numerical N-dimensional Array class. Supported element types 
-  are 1/2/4-byte Integer, single/double-precision Real/Complex, and Ruby 
-  Object. This extension library incorporates fast calculation and easy
-  manipulation of large numerical arrays into the Ruby language. NArray 
-  has features similar to NumPy, but NArray has vector and matrix subclasses.
-
-  (License: Ruby License)
-
----
-NArrayMiss Library (http://ruby.gfd-dennou.org/index.htm) by Seiya Nishizawa.
-
-  NArrayMiss is a additional class processing of missing value with to NArray
-  for Ruby.
-
-  (License: GNU Lesser General Public License (LGPL) version 2.)