simple-proj-carray 1.0.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 9e86e595e2cb9e35de18e0c3ac4ad6a5a3125c672461cc7ba8c058825ede0000
+  data.tar.gz: 9b19b6cbb2ca770794a134593cc10b4a3afaa206a1d5e0f0c6f51ed7935913ad
+SHA512:
+  metadata.gz: 147a7445f08a219ac6dbe5ae90c2890b922efad0852b19a06aba66bae78f13c17378c024f279e9570884b689c6c0a0399d257cb22b25ac33d9956421c62d0f7c
+  data.tar.gz: 2473c957e5d1c5727b2596257107e782dfeea4470e68bf224be34f53a1d75dad0a84d01041c33b5caf663eb1b57b3400c493f13b003a9de13e0dec8cb2b950bf

data/Rakefile

@@ -0,0 +1,8 @@
+GEMSPEC = "simple-proj-carray.gemspec"
+
+task :install do
+  spec = eval File.read(GEMSPEC)
+  system %{
+    gem build #{GEMSPEC}; gem install #{spec.full_name}.gem
+  }
+end

data/ext/extconf.rb

@@ -0,0 +1,23 @@
+require "mkmf"
+require 'carray/mkmf'
+
+if defined?(Gem)
+  if Gem::VERSION >= "1.7.0"    
+    rb_proj_include = [$sitearchdir, *Gem::Specification.find_all_by_name("simple-proj").map { |s| File.join(s.gem_dir,"ext") }]
+  else
+    rb_proj_include = [$sitearchdir, *Gem.all_load_paths.reverse.grep(/simple-proj/)]
+  end
+else
+  rb_proj_include = [$sitearchdir]
+end
+
+dir_config("proj", possible_includes, possible_libs)
+
+$CFLAGS += rb_proj_include.map{|s| " -I #{s}"}.join(" ")
+
+if have_header("rb_proj.h") 
+  have_carray()
+  create_makefile("simple_proj_carray")
+end
+
+

data/ext/rb_simple_proj_carray.c

@@ -0,0 +1,517 @@
+#include "ruby.h"
+#include "rb_proj.h"
+#include "carray.h"
+
+static VALUE
+rb_proj_forward_ca_2 (VALUE self, volatile VALUE vlon, volatile VALUE vlat, 
+                                  volatile VALUE vx, volatile VALUE vy)
+{
+  Proj      *proj;
+  PJ_COORD   data_in, data_out;
+  
+  CArray    *clon, *clat, *cx, *cy;
+  double    *p1, *p2, *p3, *p4;
+  ca_size_t  s1,  s2,  s3,  s4;
+  int8_t    *m, *mp;
+  ca_size_t  count, err_count, i;
+  
+  int        angular_input;
+
+  Data_Get_Struct(self, Proj, proj);
+
+  if ( ! proj->is_src_latlong ) {
+    rb_raise(rb_eRuntimeError, "requires latlong src crs. use #transform_forward instead of #forward.");
+  }
+
+  angular_input = proj_angular_input(proj->ref, PJ_FWD);
+
+  clon = ca_wrap_readonly(vlon, CA_DOUBLE);
+  clat = ca_wrap_readonly(vlat, CA_DOUBLE);
+  cx   = ca_wrap_writable(vx,   CA_DOUBLE);
+  cy   = ca_wrap_writable(vy,   CA_DOUBLE);
+  
+  ca_attach_n(4, clon, clat, cx, cy);
+
+  count = ca_set_iterator(4, 
+                          clon, &p1, &s1,
+                          clat, &p2, &s2,
+                          cx,   &p3, &s3,
+                          cy,   &p4, &s4);
+
+  m = ca_allocate_mask_iterator(2, clon, clat);
+
+  mp = m;
+  err_count = 0;
+  for (i=0; i<count; i++) {
+    if ( *mp ) {
+      err_count += 1;
+      *p3 = 0.0/0.0;
+      *p4 = 0.0/0.0;
+    }
+    else {
+      if ( angular_input == 1 ) {
+        data_in.lp.lam = proj_torad(*p1);
+        data_in.lp.phi = proj_torad(*p2);
+      }
+      else {
+        data_in.xy.x = (*p1);
+        data_in.xy.y = (*p2);        
+      }
+      data_out = proj_trans(proj->ref, PJ_FWD, data_in);
+      if ( data_out.xy.x == HUGE_VAL ) {
+        err_count += 1;
+        *p3 = 0.0/0.0;
+        *p4 = 0.0/0.0;
+      }
+      else {
+        *p3 = data_out.xy.x;
+        *p4 = data_out.xy.y;          
+      }
+    }
+    p1+=s1; p2+=s2; p3+=s3; p4+=s4; mp++;
+  }
+
+  free(m);
+
+  ca_sync_n(2, cx, cy);
+  ca_detach_n(4, clon, clat, cx, cy);
+
+  return SIZE2NUM(err_count);
+}
+
+static VALUE
+rb_proj_forward_ca_3 (VALUE self, volatile VALUE vlon, volatile VALUE vlat, volatile VALUE vzi,
+                                  volatile VALUE vx, volatile VALUE vy, volatile VALUE vzo)
+{
+  Proj      *proj;
+  PJ_COORD   data_in, data_out;
+  
+  CArray    *clon, *clat, *czi, *cx, *cy, *czo;
+  double    *p1, *p2, *p3, *p4, *p5, *p6;
+  ca_size_t  s1,  s2,  s3,  s4,  s5,  s6;
+  int8_t    *m, *mp;
+  ca_size_t  count, err_count, i;
+  
+  int        angular_input;
+
+  Data_Get_Struct(self, Proj, proj);
+
+  if ( ! proj->is_src_latlong ) {
+    rb_raise(rb_eRuntimeError, "requires latlong src crs. use #transform_forward instead of #forward.");
+  }
+
+  angular_input = proj_angular_input(proj->ref, PJ_FWD);
+
+  clon = ca_wrap_readonly(vlon, CA_DOUBLE);
+  clat = ca_wrap_readonly(vlat, CA_DOUBLE);
+  czi  = ca_wrap_readonly(vzi,  CA_DOUBLE);
+  cx   = ca_wrap_writable(vx,   CA_DOUBLE);
+  cy   = ca_wrap_writable(vy,   CA_DOUBLE);
+  czo  = ca_wrap_writable(vzo,  CA_DOUBLE);
+  
+  ca_attach_n(6, clon, clat, czi, cx, cy, czo);
+
+  count = ca_set_iterator(6, 
+                          clon, &p1, &s1,
+                          clat, &p2, &s2,
+                          czi,  &p3, &s3,
+                          cx,   &p4, &s4,
+                          cy,   &p5, &s5,
+                          czo,  &p6, &s6);
+
+  m = ca_allocate_mask_iterator(3, clon, clat, czi);
+
+  mp = m;
+  err_count = 0;
+  for (i=0; i<count; i++) {
+    if ( *mp ) {
+      err_count += 1;
+      *p4 = 0.0/0.0;
+      *p5 = 0.0/0.0;
+      *p6 = 0.0/0.0;
+    }
+    else {
+      if ( proj->is_src_latlong == 1 ) {
+        data_in.lpz.lam = proj_torad(*p1);
+        data_in.lpz.phi = proj_torad(*p2);
+        data_in.lpz.z   = (*p3);
+      }
+      else {
+        data_in.xyz.x = (*p1);
+        data_in.xyz.y = (*p2);        
+        data_in.xyz.z = (*p3);
+      }
+      data_out = proj_trans(proj->ref, PJ_FWD, data_in);
+      if ( data_out.xyz.x == HUGE_VAL ) {
+        err_count += 1;
+        *p4 = 0.0/0.0;
+        *p5 = 0.0/0.0;
+        *p6 = 0.0/0.0;
+      }
+      else {
+        *p4 = data_out.xyz.x;
+        *p5 = data_out.xyz.y;
+        *p6 = data_out.xyz.z;
+      }
+    }
+    p1+=s1; p2+=s2; p3+=s3; p4+=s4; p5+=s5; p5+=s5; mp++;
+  }
+
+  free(m);
+
+  ca_sync_n(3, cx, cy, czo);
+  ca_detach_n(6, clon, clat, czi, cx, cy, czo);
+
+  return SIZE2NUM(err_count);
+}
+
+static VALUE
+rb_proj_inverse_ca_2 (VALUE self, volatile VALUE vx, volatile VALUE vy, 
+                                  volatile VALUE vlon, volatile VALUE vlat)
+{
+  Proj      *proj;
+  PJ_COORD   data_in, data_out;
+
+  CArray    *cx, *cy, *clon, *clat;
+  double    *p1, *p2, *p3, *p4;
+  ca_size_t  s1,  s2,  s3,  s4;
+  int8_t    *m, *mp;
+  ca_size_t  count, err_count, i;
+
+  int        angular_output;
+  
+
+  Data_Get_Struct(self, Proj, proj);
+
+  if ( ! proj->is_src_latlong ) {
+    rb_raise(rb_eRuntimeError, "requires latlong src crs. use #transform_inverse instead of #inverse.");
+  }
+
+  angular_output = proj_angular_output(proj->ref, PJ_INV);
+
+  cx   = ca_wrap_readonly(vx,   CA_DOUBLE);
+  cy   = ca_wrap_readonly(vy,   CA_DOUBLE);
+  clon = ca_wrap_writable(vlon, CA_DOUBLE);
+  clat = ca_wrap_writable(vlat, CA_DOUBLE);
+
+  ca_attach_n(4, cx, cy, clon, clat);
+
+  count = ca_set_iterator(4, 
+                          cx,   &p1, &s1,
+                          cy,   &p2, &s2,
+                          clon, &p3, &s3,
+                          clat, &p4, &s4);
+
+  m = ca_allocate_mask_iterator(2, clon, clat);
+  
+  mp = m;
+  err_count = 0;
+  for (i=0; i<count; i++) {
+    if ( *mp ) {
+      err_count += 1;
+      *p3 = 0.0/0.0;
+      *p4 = 0.0/0.0;
+    }
+    else {
+      data_in.xy.x = (*p1);
+      data_in.xy.y = (*p2);        
+      data_out = proj_trans(proj->ref, PJ_INV, data_in);
+      if ( data_out.lp.lam == HUGE_VAL ) {
+        err_count += 1;
+        *p3 = 0.0/0.0;
+        *p4 = 0.0/0.0;
+      }
+      else {
+        if ( angular_output == 1 ) {
+          *p3 = proj_todeg(data_out.lp.lam);
+          *p4 = proj_todeg(data_out.lp.phi);
+        }
+        else {
+          *p3 = data_out.xy.x;
+          *p4 = data_out.xy.y;          
+        }
+      }
+    }
+    p1+=s1; p2+=s2; p3+=s3; p4+=s4;
+    mp++;
+  }
+
+  free(m);
+
+  ca_sync_n(2, clon, clat);
+  ca_detach_n(4, cx, cy, clon, clat);
+
+  return SIZE2NUM(err_count);
+}
+
+
+static VALUE
+rb_proj_inverse_ca_3 (VALUE self, volatile VALUE vx, volatile VALUE vy, volatile VALUE vzi,
+                                  volatile VALUE vlon, volatile VALUE vlat, volatile VALUE vzo)
+{
+  Proj      *proj;
+  PJ_COORD   data_in, data_out;
+
+  CArray    *cx, *cy, *czi, *clon, *clat, *czo;
+  double    *p1, *p2, *p3, *p4, *p5, *p6;
+  ca_size_t  s1,  s2,  s3,  s4,  s5,  s6;
+  int8_t    *m, *mp;
+  ca_size_t  count, err_count, i;
+
+  int        angular_output;
+
+  Data_Get_Struct(self, Proj, proj);
+
+  if ( ! proj->is_src_latlong ) {
+    rb_raise(rb_eRuntimeError, "requires latlong src crs. use #transform_inverse instead of #inverse.");
+  }
+
+  angular_output = proj_angular_output(proj->ref, PJ_INV);
+
+  cx   = ca_wrap_readonly(vx,   CA_DOUBLE);
+  cy   = ca_wrap_readonly(vy,   CA_DOUBLE);
+  czi  = ca_wrap_readonly(vzi,  CA_DOUBLE);
+  clon = ca_wrap_writable(vlon, CA_DOUBLE);
+  clat = ca_wrap_writable(vlat, CA_DOUBLE);
+  czo  = ca_wrap_writable(vzo,  CA_DOUBLE);
+
+  ca_attach_n(6, cx, cy, czi, clon, clat, czo);
+
+  count = ca_set_iterator(6, 
+                          cx,   &p1, &s1,
+                          cy,   &p2, &s2,
+                          czi,  &p3, &s3,
+                          clon, &p4, &s4,
+                          clat, &p5, &s5,
+                          czo,  &p6, &s6);
+
+  m = ca_allocate_mask_iterator(3, clon, clat, czi);
+  
+  mp = m;
+  err_count = 0;
+  for (i=0; i<count; i++) {
+    if ( *mp ) {
+      err_count += 1;
+      *p4 = 0.0/0.0;
+      *p5 = 0.0/0.0;
+      *p6 = 0.0/0.0;
+    }
+    else {
+      data_in.xyz.x = (*p1);
+      data_in.xyz.y = (*p2);
+      data_in.xyz.z = (*p3);
+      data_out = proj_trans(proj->ref, PJ_INV, data_in);
+      if ( data_out.lpz.lam == HUGE_VAL ) {
+        err_count += 1;
+        *p4 = 0.0/0.0;
+        *p5 = 0.0/0.0;
+        *p6 = 0.0/0.0;
+      }
+      else {
+        if ( angular_output == 1 ) {
+          *p4 = proj_todeg(data_out.lpz.lam);
+          *p5 = proj_todeg(data_out.lpz.phi);
+          *p6 = data_out.lpz.z;
+        }
+        else {
+          *p4 = data_out.xyz.x;
+          *p5 = data_out.xyz.y;
+          *p6 = data_out.xyz.z;          
+        }
+      }
+    }
+    p1+=s1; p2+=s2; p3+=s3; p4+=s4;
+    mp++;
+  }
+
+  free(m);
+
+  ca_sync_n(3, clon, clat, czo);
+  ca_detach_n(6, cx, cy, czi, clon, clat, czo);
+
+  return SIZE2NUM(err_count);
+}
+
+static VALUE
+rb_proj_transform_ca_2_i (VALUE self, VALUE vx, VALUE vy,
+                                      VALUE vxo, VALUE vyo, PJ_DIRECTION direction)
+{
+  Proj *trans;
+  PJ_COORD c_in, c_out;
+  CArray     *cx, *cy, *cxo, *cyo;
+  double     *p1, *p2, *p4,  *p5;
+  ca_size_t   s1,  s2,  s4,   s5;
+  int8_t  *m, *mp;
+  ca_size_t     count, err_count, i;
+
+  Data_Get_Struct(self, Proj, trans);
+
+  cx   = ca_wrap_readonly(vx, CA_DOUBLE);
+  cy   = ca_wrap_readonly(vy, CA_DOUBLE);
+  cxo  = ca_wrap_writable(vxo, CA_DOUBLE);
+  cyo  = ca_wrap_writable(vyo, CA_DOUBLE);
+
+  count = ca_get_loop_count(2, cx, cy);
+  if ( count == -1 ) {
+    rb_raise(rb_eRuntimeError, "invalid data_num");
+  }
+  
+  ca_attach_n(4, cx, cy, cxo, cyo);
+
+  ca_set_iterator(4, 
+                  cx, &p1, &s1,
+                  cy, &p2, &s2,
+                  cxo, &p4, &s4,
+                  cyo, &p5, &s5);
+
+  m = ca_allocate_mask_iterator(2, cx, cy);
+
+  mp = m;
+  err_count = 0;
+  for (i=0; i<count; i++) {
+    if ( *mp ) {
+      err_count += 1;
+      *p4 = 0.0/0.0;
+      *p5 = 0.0/0.0;
+    }
+    else {
+      c_in.xyz.x = (*p1);
+      c_in.xyz.y = (*p2);
+      c_out = proj_trans(trans->ref, PJ_FWD, c_in);
+      (*p4) = c_out.xyz.x;
+      (*p5) = c_out.xyz.y;
+      if ( c_out.xyz.x == HUGE_VAL ) {
+	      err_count += 1;
+      }
+    }
+    p1+=s1; p2+=s2;  
+    p4+=s4; p5+=s5; 
+    mp++;
+  }
+
+  free(m);
+
+  ca_sync_n(2, cxo, cyo);
+  ca_detach_n(4, cx, cy, cxo, cyo);
+
+  return SIZE2NUM(err_count);
+}
+
+static VALUE
+rb_proj_transform_forward_ca_2 (VALUE self, VALUE vx, VALUE vy,  
+                                            VALUE vxo, VALUE vyo)
+{
+  return rb_proj_transform_ca_2_i(self, vx, vy, vxo, vyo, PJ_FWD);
+}
+
+static VALUE
+rb_proj_transform_inverse_ca_2 (VALUE self, VALUE vx, VALUE vy, 
+                                            VALUE vxo, VALUE vyo)
+{
+  return rb_proj_transform_ca_2_i(self, vx, vy, vxo, vyo, PJ_INV);
+}
+
+static VALUE
+rb_proj_transform_ca_3_i (VALUE self, VALUE vx, VALUE vy, VALUE vz, 
+                                      VALUE vxo, VALUE vyo, VALUE vzo, PJ_DIRECTION direction)
+{
+  Proj *trans;
+  PJ_COORD c_in, c_out;
+  CArray     *cx, *cy, *cz, *cxo, *cyo, *czo;
+  double     *p1, *p2, *p3,  *p4,  *p5,  *p6;
+  ca_size_t   s1,  s2,  s3,   s4,   s5,   s6;
+  int8_t  *m, *mp;
+  ca_size_t     count, err_count, i;
+
+  Data_Get_Struct(self, Proj, trans);
+
+  cx   = ca_wrap_readonly(vx, CA_DOUBLE);
+  cy   = ca_wrap_readonly(vy, CA_DOUBLE);
+  cz   = ca_wrap_readonly(vz, CA_DOUBLE);
+  cxo  = ca_wrap_writable(vxo, CA_DOUBLE);
+  cyo  = ca_wrap_writable(vyo, CA_DOUBLE);
+  czo  = ca_wrap_writable(vzo, CA_DOUBLE);
+
+  count = ca_get_loop_count(3, cx, cy, cz);
+  if ( count == -1 ) {
+    rb_raise(rb_eRuntimeError, "invalid data_num");
+  }
+  
+  ca_attach_n(6, cx, cy, cz, cxo, cyo, czo);
+
+  ca_set_iterator(6, 
+                  cx, &p1, &s1,
+                  cy, &p2, &s2,
+                  cz, &p3, &s3,
+                  cxo, &p4, &s4,
+                  cyo, &p5, &s5,
+                  czo, &p6, &s6);
+
+  m = ca_allocate_mask_iterator(3, cx, cy, cz);
+
+  mp = m;
+  err_count = 0;
+  for (i=0; i<count; i++) {
+    if ( *mp ) {
+      err_count += 1;
+      *p4 = 0.0/0.0;
+      *p5 = 0.0/0.0;
+      *p6 = 0.0/0.0;
+    }
+    else {
+      c_in.xyz.x = (*p1);
+      c_in.xyz.y = (*p2);
+      c_in.xyz.z = (*p3);        
+      c_out = proj_trans(trans->ref, PJ_FWD, c_in);
+      (*p4) = c_out.xyz.x;
+      (*p5) = c_out.xyz.y;
+      (*p6) = c_out.xyz.z;        
+/*      if ( status < 0 ) {
+	      err_count += 1;
+      }
+*/
+    }
+    p1+=s1; p2+=s2; p3+=s3; 
+    p4+=s4; p5+=s5; p6+=s6;
+    mp++;
+  }
+
+  free(m);
+
+  ca_sync_n(3, cxo, cyo, czo);
+  ca_detach_n(6, cx, cy, cz, cxo, cyo, czo);
+
+  return SIZE2NUM(err_count);
+}
+
+static VALUE
+rb_proj_transform_forward_ca_3 (VALUE self, VALUE vx, VALUE vy, VALUE vz, 
+                                            VALUE vxo, VALUE vyo, VALUE vzo)
+{
+  return rb_proj_transform_ca_3_i(self, vx, vy, vz, vxo, vyo, vzo, PJ_FWD);
+}
+
+static VALUE
+rb_proj_transform_inverse_ca_3 (VALUE self, VALUE vx, VALUE vy, VALUE vz, 
+                                            VALUE vxo, VALUE vyo, VALUE vzo)
+{
+  return rb_proj_transform_ca_3_i(self, vx, vy, vz, vxo, vyo, vzo, PJ_INV);
+}
+
+void
+Init_simple_proj_carray ()
+{
+  rb_cProj = rb_define_class("PROJ", rb_cObject);
+
+  rb_define_private_method(rb_cProj, "_forward_ca_2", rb_proj_forward_ca_2, 4);
+  rb_define_private_method(rb_cProj, "_forward_ca_3", rb_proj_forward_ca_3, 6);
+  rb_define_private_method(rb_cProj, "_inverse_ca_2", rb_proj_inverse_ca_2, 4);
+  rb_define_private_method(rb_cProj, "_inverse_ca_3", rb_proj_inverse_ca_3, 6);
+
+  rb_define_private_method(rb_cProj, "_transform_forward_ca_2", rb_proj_transform_forward_ca_2, 4);
+  rb_define_private_method(rb_cProj, "_transform_forward_ca_3", rb_proj_transform_forward_ca_3, 6);
+
+  rb_define_private_method(rb_cProj, "_transform_inverse_ca_2", rb_proj_transform_inverse_ca_2, 4);
+  rb_define_private_method(rb_cProj, "_transform_inverse_ca_3", rb_proj_transform_inverse_ca_3, 6);
+
+}

data/lib/simple-proj-carray.rb

@@ -0,0 +1,104 @@
+require 'simple-proj'
+require 'carray'
+require 'simple_proj_carray.so'
+
+class PROJ
+
+  def forward_carray (lon, lat, z = nil)
+    if z.nil?
+      argv = [lon, lat].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      lon, lat = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      x    = ref.template
+      y    = ref.template
+      _forward_ca_2(lon, lat, x, y)
+      return x, y
+    else
+      argv = [lon, lat, z].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      lat, lon, z = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      x    = ref.template
+      y    = ref.template
+      zo   = ref.template
+      _forward_ca_3(lon, lat, z, x, y, zo)
+      return x, y, zo
+    end
+  end
+  
+  alias forward_lonlat_carray forward_carray
+  
+  def forward_latlon_carray (lat, lon, z = nil)
+    return forward_carray(lon, lat, z)
+  end
+  
+  def inverse_carray (x, y, z = nil)
+    if z.nil?
+      argv = [x, y].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      x, y = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      lon    = ref.template
+      lat    = ref.template
+      _inverse_ca_2(x, y, lon, lat)
+      return lon, lat
+    else
+      argv = [x, y, z].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      x, y, z = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      lon  = ref.template
+      lat  = ref.template
+      zo   = ref.template
+      _inverse_ca_3(x, y, z, lon, lat, zo)
+      return lon, lat, zo
+    end
+  end
+
+  alias inverse_lonlat_carray inverse_carray
+  
+  def inverse_lonlat_carray (lat, lon, z = nil)
+    return inverse_carray(lon, lat, z)
+  end
+
+  def transform_forward_carray (x, y, z = nil)
+    if z.nil?
+      argv = [x, y].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      x, y = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      xo   = ref.template
+      yo   = ref.template
+      p x,y
+      _transform_forward_ca_2(x, y, xo, yo)
+      return xo, yo
+    else
+      argv = [x, y, z].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      x, y, z = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      xo   = ref.template
+      yo   = ref.template
+      zo   = ref.template
+      _transform_forward_ca_3(x, y, z, xo, yo, zo)
+      return xo, yo, zo
+    end
+  end
+
+  def transform_inverse_carray (x, y, z = nil)
+    if z.nil?
+      argv = [x, y].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      x, y = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      xo  = ref.template
+      yo  = ref.template
+      _transform_inverse_ca_2(x, y, xo, yo)
+      return xo, yo
+    else
+      argv = [x, y, z].map{|a| a.is_a?(CArray) ? a : CA_DOUBLE(a) }
+      x, y, z = *argv
+      ref  = argv[argv.map(&:size).to_ca.max_addr]
+      xo   = ref.template
+      yo   = ref.template
+      zo   = ref.template
+      _transform_inverse_ca_3(x, y, z, xo, yo, zo)
+      return xo, yo, zo
+    end
+  end
+
+end

data/simple-proj-carray.gemspec

@@ -0,0 +1,27 @@
+
+Gem::Specification::new do |s|
+  version = "1.0.0"
+
+  files = Dir.glob("**/*") - [ 
+                               Dir.glob("simple-proj-carray-*.gem"), 
+                               Dir.glob("doc/**/*"),
+                               Dir.glob("test/**/*"),
+                             ].flatten
+
+  s.platform    = Gem::Platform::RUBY
+  s.name        = "simple-proj-carray"
+  s.summary     = "An extension library for PROJ for carray"
+  s.description = <<-HERE
+    An extension library for PROJ for carray
+  HERE
+  s.version     = version
+  s.licenses    = ['MIT']
+  s.author      = "Hiroki Motoyoshi"
+  s.email       = ""
+  s.homepage    = 'https://github.com/himotoyoshi/simple-proj-carray'
+  s.files       = files
+  s.extensions  = [ "ext/extconf.rb" ]
+  s.required_ruby_version = ">= 1.8.1"
+  s.add_runtime_dependency 'simple-proj', '~> 1.0'
+  s.add_runtime_dependency 'carray', '~> 1.3'
+end

metadata

@@ -0,0 +1,77 @@
+--- !ruby/object:Gem::Specification
+name: simple-proj-carray
+version: !ruby/object:Gem::Version
+  version: 1.0.0
+platform: ruby
+authors:
+- Hiroki Motoyoshi
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2020-09-16 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: simple-proj
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.0'
+- !ruby/object:Gem::Dependency
+  name: carray
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.3'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.3'
+description: "    An extension library for PROJ for carray\n"
+email: ''
+executables: []
+extensions:
+- ext/extconf.rb
+extra_rdoc_files: []
+files:
+- README.md
+- Rakefile
+- ext/extconf.rb
+- ext/rb_simple_proj_carray.c
+- lib/simple-proj-carray.rb
+- simple-proj-carray.gemspec
+homepage: https://github.com/himotoyoshi/simple-proj-carray
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 1.8.1
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.1.2
+signing_key: 
+specification_version: 4
+summary: An extension library for PROJ for carray
+test_files: []