carray 1.3.5 → 1.5.2

data/{carray_math.rb → ext/carray_math.rb}

@@ -10,7 +10,7 @@
 #
 # ----------------------------------------------------------------------------
 
-require './mkmath.rb'
+require_relative './mkmath'
 
 HEADERS << <<HERE_END
 /* ---------------------------------------------------------------------------
@@ -18,12 +18,10 @@ HEADERS << <
   carray_math.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.
 
-  This file is automatically generated from carray_math.rb.
 
-  Copyright (C) 2005 Hiroki Motoyoshi
+  This file is automatically generated from carray_math.rb.
+  Copyright (C) 2005-2020 Hiroki Motoyoshi
 
 ---------------------------------------------------------------------------- */
 
@@ -228,28 +226,28 @@ binop("pmin", "pmin",
       CMPLX_TYPES => nil,
       OBJ_TYPES =>'(#3) = rb_funcall(rb_assoc_new((#1),(#2)), rb_intern("min"), 0);')
 
-binop("add", "add",
+binop("+", "add",
       ALL_TYPES =>"(#3) = (#1) + (#2);",
       CMPLX_TYPES => HAVE_COMPLEX ? "(#3) = (#1) + (#2);" : nil,
       OBJ_TYPES =>'(#3) = rb_funcall((#1), id_plus, 1, (#2));')
 
-binop("sub", "sub",
+binop("-", "sub",
       ALL_TYPES => "(#3) = (#1) - (#2);",
       CMPLX_TYPES => HAVE_COMPLEX ? "(#3) = (#1) - (#2);" : nil,
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_minus, 1, (#2));')
 
-binop("mul", "mul",
+binop("*", "mul",
       ALL_TYPES => "(#3) = (#1) * (#2);",
       CMPLX_TYPES => HAVE_COMPLEX ? "(#3) = (#1) * (#2);" : nil,
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_star, 1, (#2));')
 
-binop("div", "div",
+binop("/", "div",
       INT_TYPES => "if ((#2)==0) {ca_zerodiv();}; (#3) = (#1) / (#2);",
       FLOAT_TYPES => "(#3) = (#1) / (#2);",
       CMPLX_TYPES => HAVE_COMPLEX ? "(#3) = (#1) / (#2);" : nil,
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_slash, 1, (#2));')
 
-binop("quo_i", nil,
+binop("quo_i", "quo_i",
       OBJ_TYPES => '(#3) = rb_funcall((#1), rb_intern("quo"), 1, (#2));')
 
 binop("rcp_mul", "rcp_mul",
@@ -258,21 +256,21 @@ binop("rcp_mul", "rcp_mul",
       CMPLX_TYPES => HAVE_COMPLEX ? "(#3) = (#2) / (#1);" : nil,
       OBJ_TYPES => '(#3) = rb_funcall((#2), id_slash, 1, (#1));')
 
-binop("mod", "mod",
+binop("%", "mod",
       INT_TYPES => "if ((#2)==0) {ca_zerodiv();}; (#3) = (#1) % (#2);",
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_percent, 1, (#2));')
 
-binop("bit_and", "bit_and_i",
+binop("&", "bit_and_i",
       BOOL_TYPES => "(#3) = (#1) & (#2);",
       INT_TYPES => "(#3) = (#1) & (#2);",
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_and, 1, (#2));')
 
-binop("bit_or", "bit_or_i",
+binop("|", "bit_or_i",
       BOOL_TYPES => "(#3) = (#1) | (#2);",
       INT_TYPES => "(#3) = (#1) | (#2);",
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_or, 1, (#2));')
 
-binop("bit_xor", "bit_xor_i",
+binop("^", "bit_xor_i",
       BOOL_TYPES => "(#3) = ((#1) != (#2)) ? 1 : 0;",
       INT_TYPES => "(#3) = (#1) ^ (#2);",
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_xor, 1, (#2));')
@@ -329,11 +327,11 @@ METHODS << %{
   rb_define_method(rb_cCArray, "bit_xor", rb_ca_bit_xor, 1);
 }
 
-binop("bit_lshift", "bit_lshift",
+binop("<<", "bit_lshift",
       INT_TYPES => "(#3) = (#1) << (#2);",
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_lshift, 1, (#2));')
 
-binop("bit_rshift", "bit_rshift",
+binop(">>", "bit_rshift",
       INT_TYPES => "(#3) = (#1) >> (#2);",
       OBJ_TYPES => '(#3) = rb_funcall((#1), id_rshift, 1, (#2));')
 
@@ -461,16 +459,16 @@ METHODS << %{
   rb_define_method(rb_cCArray, "xor", rb_ca_xor, 1);
 }
 
-alias_op("+", "add")
-alias_op("-", "sub")
-alias_op("*", "mul")
-alias_op("/", "div")
-alias_op("%", "mod")
-alias_op("&", "bit_and")
-alias_op("|", "bit_or")
-alias_op("^", "bit_xor")
-alias_op("<<", "bit_lshift")
-alias_op(">>", "bit_rshift")
+alias_op("add", "+")
+alias_op("sub", "-")
+alias_op("mul", "*")
+alias_op("div", "/")
+alias_op("mod", "%")
+alias_op("bit_and", "&")
+alias_op("bit_or", "|")
+alias_op("bit_xor", "^")
+alias_op("bit_lshift", "<<")
+alias_op("bit_rshift", ">>")
 
 #alias_op("==", "eq")
 #alias_op("!=", "ne")
@@ -697,7 +695,7 @@ rb_ca_ipower (VALUE self, VALUE other)
   /* unresolved unbound repeat array generates unbound repeat array again */
   if ( ca->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
     CAUnboundRepeat *cx = (CAUnboundRepeat *) ca;
-    obj = rb_ca_ubrep_new(obj, cx->rep_rank, cx->rep_dim);
+    obj = rb_ca_ubrep_new(rb_ca_ubrep_shave(self, obj), cx->rep_ndim, cx->rep_dim);
   }
 
   return obj;
@@ -775,8 +773,8 @@ static VALUE rb_ca_pow (VALUE self, VALUE other)
 
     /* unresolved unbound repeat array generates unbound repeat array again */
     if ( ca->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
-      CAUnboundRepeat *cx = (CAUnboundRepeat *) ca;
-      obj = rb_ca_ubrep_new(obj, cx->rep_rank, cx->rep_dim);
+      CAUnboundRepeat *cx = (CAUnboundRepeat *) ca;      
+      obj = rb_ca_ubrep_new(rb_ca_ubrep_shave(self, obj), cx->rep_ndim, cx->rep_dim);
     }
 
     return obj;

data/ext/{mathfunc/carray_mathfunc.c → carray_mathfunc.c}

@@ -3,10 +3,8 @@
   carray/mathfunc/carray_mathfunc.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
+  Copyright (C) 2005-2020 Hiroki Motoyoshi
 
 ---------------------------------------------------------------------------- */
 

data/{carray_numeric.c → ext/carray_numeric.c}

@@ -3,10 +3,8 @@
   carray_numeric.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
+  Copyright (C) 2005-2020 Hiroki Motoyoshi
 
 ---------------------------------------------------------------------------- */
 
@@ -168,7 +166,7 @@ rb_ca_arg (VALUE self)
 
   Data_Get_Struct(self, CArray, ca);
 
-  co = carray_new(CA_FLOAT64, ca->rank, ca->dim, 0, NULL);
+  co = carray_new(CA_FLOAT64, ca->ndim, ca->dim, 0, NULL);
   obj = ca_wrap_struct(co);
 
   if ( ca_has_mask(ca) ) {
@@ -220,6 +218,15 @@ Init_carray_numeric ()
   rb_define_alias(rb_cFalseClass, "__and__", "&");
   rb_define_alias(rb_cFalseClass, "__xor__", "^");
 
+#if RUBY_VERSION_CODE >= 240
+
+  rb_define_alias(rb_cInteger, "__or__", "|");
+  rb_define_alias(rb_cInteger, "__and__", "&");
+  rb_define_alias(rb_cInteger, "__xor__", "^");
+  rb_define_alias(rb_cInteger, "__lshift__", "<<");
+  rb_define_alias(rb_cInteger, "__rshift__", ">>");
+
+#else
   rb_define_alias(rb_cFixnum, "__or__", "|");
   rb_define_alias(rb_cFixnum, "__and__", "&");
   rb_define_alias(rb_cFixnum, "__xor__", "^");
@@ -231,6 +238,7 @@ Init_carray_numeric ()
   rb_define_alias(rb_cBignum, "__xor__", "^");
   rb_define_alias(rb_cBignum, "__lshift__", "<<");
   rb_define_alias(rb_cBignum, "__rshift__", ">>");
+#endif
 
   rb_define_method(rb_cTrueClass, "|", rb_hack_or, 1);
   rb_define_method(rb_cTrueClass, "&", rb_hack_and, 1);
@@ -242,6 +250,13 @@ Init_carray_numeric ()
   rb_define_method(rb_cFalseClass, "^", rb_hack_xor, 1);
   rb_define_method(rb_cFalseClass, "*", rb_hack_star, 1);
 
+#if RUBY_VERSION_CODE >= 240
+  rb_define_method(rb_cInteger, "|", rb_hack_or, 1);
+  rb_define_method(rb_cInteger, "&", rb_hack_and, 1);
+  rb_define_method(rb_cInteger, "^", rb_hack_xor, 1);
+  rb_define_method(rb_cInteger, "<<", rb_hack_lshift, 1);
+  rb_define_method(rb_cInteger, ">>", rb_hack_rshift, 1);
+#else
   rb_define_method(rb_cFixnum, "|", rb_hack_or, 1);
   rb_define_method(rb_cFixnum, "&", rb_hack_and, 1);
   rb_define_method(rb_cFixnum, "^", rb_hack_xor, 1);
@@ -253,6 +268,7 @@ Init_carray_numeric ()
   rb_define_method(rb_cBignum, "^", rb_hack_xor, 1);
   rb_define_method(rb_cBignum, "<<", rb_hack_lshift, 1);
   rb_define_method(rb_cBignum, ">>", rb_hack_rshift, 1);
+#endif
 
 #ifdef HAVE_COMPLEX_H
   rb_define_method(rb_cCArray, "arg", rb_ca_arg, 0);

data/{carray_operator.c → ext/carray_operator.c}

@@ -3,10 +3,8 @@
   carray_operator.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
+  Copyright (C) 2005-2020 Hiroki Motoyoshi
 
 ---------------------------------------------------------------------------- */
 
@@ -28,7 +26,6 @@ ca_zerodiv ()
   rb_raise(rb_eZeroDivError, "divided by 0");
 }
 
-
 VALUE
 rb_ca_call_monop (VALUE self, ca_monop_func_t func[])
 {
@@ -50,7 +47,7 @@ rb_ca_call_monop (VALUE self, ca_monop_func_t func[])
   ca_attach(ca1);
   ca_copy_mask_overlay(ca2, ca2->elements, 1, ca1);
   func[ca1->data_type](ca1->elements,
-                       ( ca2->mask ) ? ca2->mask->ptr : NULL,
+                       ( ca2->mask ) ? (boolean8_t *)ca2->mask->ptr : NULL,
                        ca1->ptr, 1,
                        ca2->ptr, 1);
   ca_detach(ca1);
@@ -58,7 +55,7 @@ rb_ca_call_monop (VALUE self, ca_monop_func_t func[])
   /* unresolved unbound repeat array generates unbound repeat array again */
   if ( ca1->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
     CAUnboundRepeat *cx = (CAUnboundRepeat *) ca1;
-    out = rb_ca_ubrep_new(out, cx->rep_rank, cx->rep_dim);
+    out = rb_ca_ubrep_new(rb_ca_ubrep_shave(self, out), cx->rep_ndim, cx->rep_dim);
   }
 
   return out;
@@ -75,7 +72,7 @@ rb_ca_call_monop_bang (VALUE self, ca_monop_func_t func[])
 
   ca_attach(ca1);
   func[ca1->data_type](ca1->elements,
-                       ( ca1->mask ) ? ca1->mask->ptr : NULL,
+                       ( ca1->mask ) ? (boolean8_t *)ca1->mask->ptr : NULL,
                        ca1->ptr, 1,
                        ca1->ptr, 1);
   ca_sync(ca1);
@@ -84,6 +81,29 @@ rb_ca_call_monop_bang (VALUE self, ca_monop_func_t func[])
   return self;
 }
 
+int
+rb_ca_test_castable (VALUE other)
+{
+  volatile VALUE retval;
+  if ( rb_respond_to(other, rb_intern("castable_to_carray?")) ) {
+    retval = rb_funcall(other, rb_intern("castable_to_carray?"), 0);
+    return RTEST(retval); 
+  }
+  else {
+    return 1;
+  }
+}
+
+VALUE 
+rb_ca_binop_pass_to_other (VALUE self, VALUE other, ID method)
+{
+  volatile VALUE pair;
+  pair = rb_funcall(other, rb_intern("coerce"), 1, self);
+  self  = rb_ary_entry(pair, 0);
+  other = rb_ary_entry(pair, 1);
+  return rb_funcall(self, method, 1, other);
+}
+
 VALUE
 rb_ca_call_binop (volatile VALUE self, volatile VALUE other,
                                          ca_binop_func_t func[])
@@ -111,7 +131,8 @@ rb_ca_call_binop (volatile VALUE self, volatile VALUE other,
       out = ca_wrap_struct(ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, 0,
                            ca2->ptr, 0,
                            ca3->ptr, 0);
@@ -126,7 +147,8 @@ rb_ca_call_binop (volatile VALUE self, volatile VALUE other,
       out = ca_wrap_struct(ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca2->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca2->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, 0,
                            ca2->ptr, 1,
                            ca3->ptr, 1);
@@ -143,7 +165,8 @@ rb_ca_call_binop (volatile VALUE self, volatile VALUE other,
       out = ca_wrap_struct(ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, 1,
                            ca2->ptr, 0,
                            ca3->ptr, 1);
@@ -162,7 +185,8 @@ rb_ca_call_binop (volatile VALUE self, volatile VALUE other,
       out = ca_wrap_struct(ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, 1,
                            ca2->ptr, 1,
                            ca3->ptr, 1);
@@ -174,13 +198,13 @@ rb_ca_call_binop (volatile VALUE self, volatile VALUE other,
   /* unresolved unbound repeat array generates unbound repeat array again */
   if ( ca1->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
     CAUnboundRepeat *cx = (CAUnboundRepeat *) ca1;
-    out = rb_ca_ubrep_new(out, cx->rep_rank, cx->rep_dim);
+    out = rb_ca_ubrep_new(rb_ca_ubrep_shave(self, out), cx->rep_ndim, cx->rep_dim);
   }
 
   /* unresolved unbound repeat array generates unbound repeat array again */
   if ( ca2->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
     CAUnboundRepeat *cx = (CAUnboundRepeat *) ca2;
-    out = rb_ca_ubrep_new(out, cx->rep_rank, cx->rep_dim);
+    out = rb_ca_ubrep_new(rb_ca_ubrep_shave(other, out), cx->rep_ndim, cx->rep_dim);
   }
 
   return out;
@@ -205,7 +229,8 @@ rb_ca_call_binop_bang (VALUE self, VALUE other, ca_binop_func_t func[])
   if ( rb_obj_is_cscalar(self) ) {
     if ( rb_obj_is_cscalar(other) ) { /* scalar vs scalar */
       ca_copy_mask_overlay(ca1, ca1->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca1->mask ) ? ca1->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca1->mask ) ? (boolean8_t *) ca1->mask->ptr : NULL,
                            ca1->ptr, 0,
                            ca2->ptr, 0,
                            ca1->ptr, 0);
@@ -217,7 +242,8 @@ rb_ca_call_binop_bang (VALUE self, VALUE other, ca_binop_func_t func[])
       }
 
       ca_copy_mask_overlay(ca1, ca1->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca1->mask ) ? ca1->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca1->mask ) ? (boolean8_t *) ca1->mask->ptr : NULL,
                            ca1->ptr, 0,
                            ca2->ptr, 0,
                            ca1->ptr, 0);
@@ -226,7 +252,8 @@ rb_ca_call_binop_bang (VALUE self, VALUE other, ca_binop_func_t func[])
   else {
     if ( rb_obj_is_cscalar(other) ) { /* array vs scalar */
       ca_copy_mask_overlay(ca1, ca1->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca1->mask ) ? ca1->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca1->mask ) ? (boolean8_t *) ca1->mask->ptr : NULL,
                            ca1->ptr, 1,
                            ca2->ptr, 0,
                            ca1->ptr, 1);
@@ -238,7 +265,8 @@ rb_ca_call_binop_bang (VALUE self, VALUE other, ca_binop_func_t func[])
       }
 
       ca_copy_mask_overlay(ca1, ca1->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca1->mask ) ? ca1->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca1->mask ) ? (boolean8_t *) ca1->mask->ptr : NULL,
                            ca1->ptr, 1,
                            ca2->ptr, 1,
                            ca1->ptr, 1);
@@ -264,22 +292,23 @@ rb_ca_call_moncmp (VALUE self, ca_moncmp_func_t func[])
     out = rb_cscalar_new(CA_BOOLEAN, 0, NULL);
   }
   else {
-    out = rb_carray_new(CA_BOOLEAN, ca1->rank, ca1->dim, 0, NULL);
+    out = rb_carray_new(CA_BOOLEAN, ca1->ndim, ca1->dim, 0, NULL);
   }
 
   Data_Get_Struct(out, CArray, ca2);
 
   ca_attach(ca1);
   ca_copy_mask_overlay(ca2, ca2->elements, 1, ca1);
-  func[ca1->data_type](ca1->elements, ( ca2->mask ) ? ca2->mask->ptr : NULL,
+  func[ca1->data_type](ca1->elements, 
+                       ( ca2->mask ) ? (boolean8_t *) ca2->mask->ptr : NULL,
                        ca1->ptr, 1,
-                       ca2->ptr, 1);
+                       (boolean8_t *) ca2->ptr, 1);
   ca_detach(ca1);
 
   /* unresolved unbound repeat array generates unbound repeat array again */
   if ( ca1->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
     CAUnboundRepeat *cx = (CAUnboundRepeat *) ca1;
-    out = rb_ca_ubrep_new(out, cx->rep_rank, cx->rep_dim);
+    out = rb_ca_ubrep_new(rb_ca_ubrep_shave(self, out), cx->rep_ndim, cx->rep_dim);
   }
 
   return out;
@@ -298,10 +327,10 @@ rb_ca_call_bincmp (volatile VALUE self, volatile VALUE other,
 
   /* check for comparison with CA_UNDEF */
   if ( other == CA_UNDEF ) {
-    if ( func == ca_bincmp_eq ) {      /* a.eq(UNDEF) -> a.is_masked */
+    if ( (ca_bincmp_func_t) func == (ca_bincmp_func_t) ca_bincmp_eq ) {  /* a.eq(UNDEF) -> a.is_masked */
       return rb_ca_is_masked(self);
     }
-    else if ( func == ca_bincmp_ne ) { /* a.ne(UNDEF) -> a.is_not_masked */
+    else if ( (ca_bincmp_func_t) func == (ca_bincmp_func_t) ca_bincmp_ne ) { /* a.ne(UNDEF) -> a.is_not_masked */
       return rb_ca_is_not_masked(self);
     }
     else {
@@ -324,17 +353,19 @@ rb_ca_call_bincmp (volatile VALUE self, volatile VALUE other,
       Data_Get_Struct(out, CArray, ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, ca1->bytes, 0,
                            ca2->ptr, ca2->bytes, 0,
                            ca3->ptr, ca3->bytes, 0);
     }
     else {                                          /* scalar vs array */
-      out = rb_carray_new(CA_BOOLEAN, ca2->rank, ca2->dim, 0, NULL);
+      out = rb_carray_new(CA_BOOLEAN, ca2->ndim, ca2->dim, 0, NULL);
       Data_Get_Struct(out, CArray, ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca2->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca2->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, ca1->bytes, 0,
                            ca2->ptr, ca2->bytes, 1,
                            ca3->ptr, ca3->bytes, 1);
@@ -342,11 +373,12 @@ rb_ca_call_bincmp (volatile VALUE self, volatile VALUE other,
   }
   else {
     if ( rb_obj_is_cscalar(other) ) {  /* array vs scalar */
-      out = rb_carray_new(CA_BOOLEAN, ca1->rank, ca1->dim, 0, NULL);
+      out = rb_carray_new(CA_BOOLEAN, ca1->ndim, ca1->dim, 0, NULL);
       Data_Get_Struct(out, CArray, ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, ca1->bytes, 1,
                            ca2->ptr, ca2->bytes, 0,
                            ca3->ptr, ca3->bytes, 1);
@@ -356,11 +388,12 @@ rb_ca_call_bincmp (volatile VALUE self, volatile VALUE other,
         rb_raise(rb_eRuntimeError, "elements mismatch in bincmp (%lld <-> %lld)",
                                  (ca_size_t) ca1->elements, (ca_size_t) ca2->elements);
       }
-      out = rb_carray_new(CA_BOOLEAN, ca1->rank, ca1->dim, 0, NULL);
+      out = rb_carray_new(CA_BOOLEAN, ca1->ndim, ca1->dim, 0, NULL);
       Data_Get_Struct(out, CArray, ca3);
 
       ca_copy_mask_overlay(ca3, ca3->elements, 2, ca1, ca2);
-      func[ca1->data_type](ca1->elements, ( ca3->mask ) ? ca3->mask->ptr : NULL,
+      func[ca1->data_type](ca1->elements, 
+                           ( ca3->mask ) ? (boolean8_t *) ca3->mask->ptr : NULL,
                            ca1->ptr, ca1->bytes, 1,
                            ca2->ptr, ca2->bytes, 1,
                            ca3->ptr, ca3->bytes, 1);
@@ -372,7 +405,7 @@ rb_ca_call_bincmp (volatile VALUE self, volatile VALUE other,
   /* unresolved unbound repeat array generates unbound repeat array again */
   if ( ca1->obj_type == CA_OBJ_UNBOUND_REPEAT ) {
     CAUnboundRepeat *cx = (CAUnboundRepeat *) ca1;
-    out = rb_ca_ubrep_new(out, cx->rep_rank, cx->rep_dim);
+    out = rb_ca_ubrep_new(rb_ca_ubrep_shave(self, out), cx->rep_ndim, cx->rep_dim);
   }
 
   return out;
@@ -435,11 +468,9 @@ ca_math_call (VALUE mod, VALUE arg, ID id)
   }
 }
 
-/* rdoc:
-  class CArray
-    def coerce (other)
-    end
-  end
+/* @overload coerece (other)
+
+[TBD]
 */
 
 static VALUE
@@ -501,6 +532,11 @@ rb_ca_coerce (VALUE self, VALUE other)
     } \
   }
 
+/* @overload mul_add (weight, min_count=nil, fill_value=nil)
+
+[TBD]
+*/
+
 static VALUE
 rb_ca_mul_add (int argc, VALUE *argv, volatile VALUE self)
 {