numo-narray-alt 0.9.4 → 0.9.6

data/ext/numo/narray/index.c

@@ -71,7 +71,10 @@ static void na_index_set_step(na_index_arg_t* q, int i, size_t n, size_t beg, ss
 }
 
 static void na_index_set_scalar(na_index_arg_t* q, int i, ssize_t size, ssize_t x) {
-  if (x < -size || x >= size) rb_raise(rb_eRangeError, "array index (%" SZF "d) is out of array size (%" SZF "d)", x, size);
+  if (x < -size || x >= size)
+    rb_raise(
+      rb_eRangeError, "array index (%" SZF "d) is out of array size (%" SZF "d)", x, size
+    );
   if (x < 0) x += size;
   q->n = 1;
   q->beg = x;
@@ -172,7 +175,8 @@ static void na_parse_narray_index(VALUE a, int orig_dim, ssize_t size, na_index_
   q->orig_dim = orig_dim;
 }
 
-static void na_parse_range(VALUE range, ssize_t step, int orig_dim, ssize_t size, na_index_arg_t* q) {
+static void
+na_parse_range(VALUE range, ssize_t step, int orig_dim, ssize_t size, na_index_arg_t* q) {
   int n;
   ssize_t beg, end, beg_orig, end_orig;
   const char *dot = "..", *edot = "...";
@@ -192,7 +196,9 @@ static void na_parse_range(VALUE range, ssize_t step, int orig_dim, ssize_t size
       dot = edot;
     }
     if (beg < 0 || beg >= size) {
-      rb_raise(rb_eRangeError, "%" SZF "d%s is out of range for size=%" SZF "d", beg_orig, dot, size);
+      rb_raise(
+        rb_eRangeError, "%" SZF "d%s is out of range for size=%" SZF "d", beg_orig, dot, size
+      );
     }
   } else {
     end = end_orig = NUM2SSIZET(x.end);
@@ -204,7 +210,10 @@ static void na_parse_range(VALUE range, ssize_t step, int orig_dim, ssize_t size
       dot = edot;
     }
     if (beg < 0 || beg >= size || end < 0 || end >= size) {
-      rb_raise(rb_eRangeError, "%" SZF "d%s%" SZF "d is out of range for size=%" SZF "d", beg_orig, dot, end_orig, size);
+      rb_raise(
+        rb_eRangeError, "%" SZF "d%s%" SZF "d is out of range for size=%" SZF "d", beg_orig,
+        dot, end_orig, size
+      );
     }
   }
 #else
@@ -224,7 +233,10 @@ static void na_parse_range(VALUE range, ssize_t step, int orig_dim, ssize_t size
     dot = edot;
   }
   if (beg < 0 || beg >= size || end < 0 || end >= size) {
-    rb_raise(rb_eRangeError, "%" SZF "d%s%" SZF "d is out of range for size=%" SZF "d", beg_orig, dot, end_orig, size);
+    rb_raise(
+      rb_eRangeError, "%" SZF "d%s%" SZF "d is out of range for size=%" SZF "d", beg_orig, dot,
+      end_orig, size
+    );
   }
 #endif
   n = (int)((end - beg) / step + 1);
@@ -342,7 +354,8 @@ static void na_index_parse_each(volatile VALUE a, ssize_t size, int i, na_index_
   }
 }
 
-static void na_at_parse_each(volatile VALUE a, ssize_t size, int i, VALUE* idx, ssize_t stride) {
+static void
+na_at_parse_each(volatile VALUE a, ssize_t size, int i, VALUE* idx, ssize_t stride) {
   na_index_arg_t q;
   size_t n, k;
   ssize_t* index;
@@ -444,7 +457,6 @@ static size_t na_index_parse_args(VALUE args, narray_t* na, na_index_arg_t* q, i
     // rest (ellipsis) dimension
     if (v == Qfalse) {
       for (l = ndim - (nidx - 1); l > 0; l--) {
-        // printf("i=%d j=%d k=%d l=%d ndim=%d nidx=%d\n",i,j,k,l,ndim,nidx);
         na_index_parse_each(Qtrue, na->shape[k], k, &q[j]);
         if (q[j].n > 1) {
           total *= q[j].n;
@@ -481,8 +493,10 @@ static void na_get_strides_nadata(const narray_data_t* na, ssize_t* strides, ssi
   }
 }
 
-static void na_index_aref_nadata(narray_data_t* na1, narray_view_t* na2, na_index_arg_t* q, ssize_t elmsz, int ndim,
-                                 int keep_dim) {
+static void na_index_aref_nadata(
+  narray_data_t* na1, narray_view_t* na2, na_index_arg_t* q, ssize_t elmsz, int ndim,
+  int keep_dim
+) {
   int i, j;
   ssize_t size, k, total = 1;
   ssize_t stride1;
@@ -531,8 +545,10 @@ static void na_index_aref_nadata(narray_data_t* na1, narray_view_t* na2, na_inde
   na2->base.size = total;
 }
 
-static void na_index_aref_naview(narray_view_t* na1, narray_view_t* na2, na_index_arg_t* q, ssize_t elmsz, int ndim,
-                                 int keep_dim) {
+static void na_index_aref_naview(
+  narray_view_t* na1, narray_view_t* na2, na_index_arg_t* q, ssize_t elmsz, int ndim,
+  int keep_dim
+) {
   int i, j;
   ssize_t total = 1;
 
@@ -779,10 +795,12 @@ static int check_index_count(int argc, int na_ndim, int count_new, int count_res
   case 0:
     if (argc == 1 && count_new == 0) return 1;
     if (argc == result_nd) return result_nd;
-    rb_raise(rb_eIndexError,
-             "# of index(=%i) should be "
-             "equal to ndim(=%i) or 1",
-             argc, na_ndim);
+    rb_raise(
+      rb_eIndexError,
+      "# of index(=%i) should be "
+      "equal to ndim(=%i) or 1",
+      argc, na_ndim
+    );
     break;
   case 1:
     if (argc - 1 <= result_nd) return result_nd;
@@ -794,7 +812,9 @@ static int check_index_count(int argc, int na_ndim, int count_new, int count_res
   return -1;
 }
 
-int na_get_result_dimension(VALUE self, int argc, VALUE* argv, ssize_t stride, size_t* pos_idx) {
+int na_get_result_dimension(
+  VALUE self, int argc, VALUE* argv, ssize_t stride, size_t* pos_idx
+) {
   int i, j;
   int count_new = 0;
   int count_rest = 0;
@@ -890,10 +910,12 @@ int na_get_result_dimension(VALUE self, int argc, VALUE* argv, ssize_t stride, s
       return 0;
     }
   }
-  rb_raise(rb_eIndexError,
-           "# of index(=%i) should be "
-           "equal to ndim(=%i) or 1",
-           argc, na->ndim);
+  rb_raise(
+    rb_eIndexError,
+    "# of index(=%i) should be "
+    "equal to ndim(=%i) or 1",
+    argc, na->ndim
+  );
   return -1;
 }
 
@@ -936,11 +958,12 @@ static VALUE na_slice(int argc, VALUE* argv, VALUE self) {
 
 /*
   Multi-dimensional element reference.
-  Returns an element at `dim0`, `dim1`, ... are Numeric indices for each dimension, or returns a NArray View as a sliced array
-  if `dim0`, `dim1`, ... includes other than Numeric index, e.g., Range or Array or true.
+  Returns an element at `dim0`, `dim1`, ... are Numeric indices for each dimension, or returns a
+  NArray View as a sliced array if `dim0`, `dim1`, ... includes other than Numeric index, e.g.,
+  Range or Array or true.
   @overload [](dim0,...,dimL)
-  @param [Numeric,Range,Array,Numo::Int32,Numo::Int64,Numo::Bit,TrueClass,FalseClass,Symbol] dim0,...,dimL  multi-dimensional
-  indices.
+  @param [Numeric,Range,Array,Numo::Int32,Numo::Int64,Numo::Bit,TrueClass,FalseClass,Symbol]
+  dim0,...,dimL  multi-dimensional indices.
   @return [Numeric,Numo::NArray] an element or NArray view.
   @see #[]=
   @see #at
@@ -980,8 +1003,8 @@ static VALUE na_slice(int argc, VALUE* argv, VALUE self) {
   Replace element(s) at `dim0`, `dim1`, ... .
   Broadcasting mechanism is applied.
   @overload []=(dim0,...,dimL,val)
-  @param [Numeric,Range,Array,Numo::Int32,Numo::Int64,Numo::Bit,TrueClass,FalseClass,Symbol] dim0,...,dimL  multi-dimensional
-  indices.
+  @param [Numeric,Range,Array,Numo::Int32,Numo::Int64,Numo::Bit,TrueClass,FalseClass,Symbol]
+  dim0,...,dimL  multi-dimensional indices.
   @param [Numeric,Numo::NArray,Array] val  Value(s) to be set to self.
   @return [Numeric,Numo::NArray,Array] returns `val` (last argument).
   @see #[]
@@ -1022,8 +1045,8 @@ static VALUE na_slice(int argc, VALUE* argv, VALUE self) {
   Similar to numpy's tuple indexing, i.e., `a[[1,2,..],[3,4,..]]`
   Same as Numo::NArray#[] for one-dimensional NArray.
   @overload at(dim0,...,dimL)
-  @param [Range,Array,Numo::Int32,Numo::Int64] dim0,...,dimL  multi-dimensional index arrays.
-  @return [Numo::NArray] one-dimensional NArray view.
+    @param [Range,Array,Numo::Int32,Numo::Int64] dim0,...,dimL  multi-dimensional index arrays.
+    @return [Numo::NArray] one-dimensional NArray view.
   @see #[]
 
   @example

data/ext/numo/narray/kwargs.c

@@ -15,10 +15,10 @@ struct RBasicRaw {
   VALUE klass;
 };
 
-#define RBASIC_SET_CLASS(obj, cls)                                                                                             \
-  do {                                                                                                                         \
-    VALUE _obj_ = (obj);                                                                                                       \
-    RB_OBJ_WRITE(_obj_, &((struct RBasicRaw*)(_obj_))->klass, cls);                                                            \
+#define RBASIC_SET_CLASS(obj, cls)                                                             \
+  do {                                                                                         \
+    VALUE _obj_ = (obj);                                                                       \
+    RB_OBJ_WRITE(_obj_, &((struct RBasicRaw*)(_obj_))->klass, cls);                            \
   } while (0)
 
 /* from class.c */
@@ -69,7 +69,7 @@ static int separate_symbol(st_data_t key, st_data_t value, st_data_t arg) {
 
 VALUE
 rb_extract_keywords(VALUE* orighash) {
-  VALUE parthash[2] = {0, 0};
+  VALUE parthash[2] = { 0, 0 };
   VALUE hash = *orighash;
 
   if (RHASH_EMPTY_P(hash)) {
@@ -84,15 +84,17 @@ rb_extract_keywords(VALUE* orighash) {
   return parthash[0];
 }
 
-int rb_get_kwargs(VALUE keyword_hash, const ID* table, int required, int optional, VALUE* values) {
+int rb_get_kwargs(
+  VALUE keyword_hash, const ID* table, int required, int optional, VALUE* values
+) {
   int i = 0, j;
   int rest = 0;
   VALUE missing = Qnil;
   st_data_t key;
 
-#define extract_kwarg(keyword, val)                                                                                            \
-  (key = (st_data_t)(keyword),                                                                                                 \
-   values ? st_delete(rb_hash_tbl_raw(keyword_hash), &key, (val)) : st_lookup(rb_hash_tbl_raw(keyword_hash), key, (val)))
+#define extract_kwarg(keyword, val)                                                            \
+  (key = (st_data_t)(keyword), values ? st_delete(rb_hash_tbl_raw(keyword_hash), &key, (val))  \
+                                      : st_lookup(rb_hash_tbl_raw(keyword_hash), key, (val)))
 
   if (NIL_P(keyword_hash)) keyword_hash = 0;
 

data/ext/numo/narray/math.c

@@ -40,8 +40,10 @@ static VALUE nary_math_cast2(VALUE type1, VALUE type2) {
   if (RTEST(rb_class_inherited_p(type2, cNArray))) {
     return nary_type_s_upcast(type2, type1);
   }
-  if (RTEST(rb_class_inherited_p(type1, rb_cNumeric)) && RTEST(rb_class_inherited_p(type2, rb_cNumeric))) {
-    if (RTEST(rb_class_inherited_p(type1, rb_cComplex)) || RTEST(rb_class_inherited_p(type2, rb_cComplex))) {
+  if (RTEST(rb_class_inherited_p(type1, rb_cNumeric)) &&
+      RTEST(rb_class_inherited_p(type2, rb_cNumeric))) {
+    if (RTEST(rb_class_inherited_p(type1, rb_cComplex)) ||
+        RTEST(rb_class_inherited_p(type2, rb_cComplex))) {
       return rb_cComplex;
     }
     return rb_cFloat;
@@ -70,9 +72,9 @@ static VALUE nary_mathcast(int argc, VALUE* argv) {
   Dispatches method to Math module of upcasted type,
   eg, Numo::DFloat::Math.
   @overload method_missing(name,x,...)
-  @param [Symbol] name  method name.
-  @param [NArray,Numeric] x  input array.
-  @return [NArray] result.
+    @param [Symbol] name  method name.
+    @param [NArray,Numeric] x  input array.
+    @return [NArray] result.
 */
 static VALUE nary_math_method_missing(int argc, VALUE* argv, VALUE mod) {
   VALUE type, ans, typemod, hash;
@@ -99,11 +101,15 @@ static VALUE nary_math_method_missing(int argc, VALUE* argv, VALUE mod) {
 void Init_nary_math(void) {
   VALUE hCast;
 
+  /**
+   * Document-module: Numo::NMath
+   *
+   * This module provides mathematical functions for NArray.
+   */
   numo_mNMath = rb_define_module_under(mNumo, "NMath");
   rb_define_singleton_method(numo_mNMath, "method_missing", nary_math_method_missing, -1);
 
   hCast = rb_hash_new();
-  rb_define_const(numo_mNMath, "DISPATCH", hCast);
   rb_hash_aset(hCast, numo_cInt64, numo_mDFloatMath);
   rb_hash_aset(hCast, numo_cInt32, numo_mDFloatMath);
   rb_hash_aset(hCast, numo_cInt16, numo_mDFloatMath);
@@ -125,6 +131,8 @@ void Init_nary_math(void) {
 #endif
   rb_hash_aset(hCast, rb_cFloat, rb_mMath);
   rb_hash_aset(hCast, rb_cComplex, numo_mDComplexMath);
+  /* Dispatch table representing the corresponding Math module. */
+  rb_define_const(numo_mNMath, "DISPATCH", hCast);
 
   id_send = rb_intern("send");
   id_UPCAST = rb_intern("UPCAST");

data/ext/numo/narray/narray.c

@@ -321,9 +321,9 @@ static void na_setup(VALUE self, int ndim, size_t* shape) {
 /*
   @overload initialize(shape)
   @overload initialize(size0, size1, ...)
-  @param [Array] shape (array of sizes along each dimension)
-  @param [Integer] sizeN (size along Nth-dimension)
-  @return [Numo::NArray] unallocated narray.
+    @param [Array] shape (array of sizes along each dimension)
+    @param [Integer] sizeN (size along Nth-dimension)
+    @return [Numo::NArray] unallocated narray.
 
   Constructs an instance of NArray class using the given
   and <i>shape</i> or <i>sizes</i>.
@@ -398,8 +398,8 @@ nary_view_new(VALUE klass, int ndim, size_t* shape) {
   Replaces the contents of self with the contents of other narray.
   Used in dup and clone method.
   @overload initialize_copy(other)
-  @param [Numo::NArray] other
-  @return [Numo::NArray] self
+    @param [Numo::NArray] other
+    @return [Numo::NArray] self
  */
 static VALUE na_initialize_copy(VALUE self, VALUE orig) {
   narray_t* na;
@@ -459,10 +459,10 @@ static VALUE na_s_ones(int argc, VALUE* argv, VALUE klass) {
   but for typed NArray subclasses, e.g., DFloat, Int64.
 
   @overload linspace(x1, x2, [n])
-  @param [Numeric] x1   The start value
-  @param [Numeric] x2   The end value
-  @param [Integer] n    The number of elements. (default is 100).
-  @return [Numo::NArray]  result array.
+    @param [Numeric] x1   The start value
+    @param [Numeric] x2   The end value
+    @param [Integer] n    The number of elements. (default is 100).
+    @return [Numo::NArray]  result array.
 
   @example
     a = Numo::DFloat.linspace(-5,5,7)
@@ -495,11 +495,11 @@ static VALUE na_s_linspace(int argc, VALUE* argv, VALUE klass) {
   i.e., DFloat, SFloat, DComplex, and SComplex.
 
   @overload logspace(a, b, [n, base])
-  @param [Numeric] a  The start value
-  @param [Numeric] b  The end value
-  @param [Integer] n  The number of elements. (default is 50)
-  @param [Numeric] base  The base of log space. (default is 10)
-  @return [Numo::NArray]  result array.
+    @param [Numeric] a  The start value
+    @param [Numeric] b  The end value
+    @param [Integer] n  The number of elements. (default is 50)
+    @param [Numeric] base  The base of log space. (default is 10)
+    @return [Numo::NArray]  result array.
 
   @example
     Numo::DFloat.logspace(4,0,5,2)
@@ -535,8 +535,8 @@ static VALUE na_s_logspace(int argc, VALUE* argv, VALUE klass) {
 /*
   Returns a NArray with shape=(n,n) whose diagonal elements are 1, otherwise 0.
   @overload  eye(n)
-  @param [Integer] n  Size of NArray. Creates 2-D NArray with shape=(n,n)
-  @return [Numo::NArray]  created NArray.
+    @param [Integer] n  Size of NArray. Creates 2-D NArray with shape=(n,n)
+    @return [Numo::NArray]  created NArray.
   @example
     a = Numo::DFloat.eye(3)
     # => Numo::DFloat#shape=[3,3]
@@ -857,7 +857,8 @@ void na_copy_flags(VALUE src, VALUE dst) {
   na2->flag[0] = na1->flag[0];
   // na2->flag[1] = NA_FL1_INIT;
 
-  RBASIC(dst)->flags |= (RBASIC(src)->flags) & (FL_USER1 | FL_USER2 | FL_USER3 | FL_USER4 | FL_USER5 | FL_USER6 | FL_USER7);
+  RBASIC(dst)->flags |= (RBASIC(src)->flags) & (FL_USER1 | FL_USER2 | FL_USER3 | FL_USER4 |
+                                                FL_USER5 | FL_USER6 | FL_USER7);
 }
 
 // fix name, ex, allow_stride_for_flatten_view
@@ -1021,8 +1022,9 @@ na_make_view(VALUE self) {
  *
  *  Expand the shape of an array. Insert a new axis with size=1
  *  at a given dimension.
- *  @param [Integer] dim  dimension at which new axis is inserted.
- *  @return [Numo::NArray]  result narray view.
+ *  @overload expand_dims(dim)
+ *    @param [Integer] dim  dimension at which new axis is inserted.
+ *    @return [Numo::NArray]  result narray view.
  */
 static VALUE na_expand_dims(VALUE self, VALUE vdim) {
   int i, j, nd, dim;
@@ -1188,8 +1190,8 @@ numo_na_upcast(VALUE type1, VALUE type2) {
   Note that NArray has distinct UPCAST mechanism.
   Coerce is used for operation between non-NArray and NArray.
   @overload coerce(other)
-  @param [Object] other  numeric object.
-  @return [Array]  NArray-casted [other,self]
+    @param [Object] other  numeric object.
+    @return [Array]  NArray-casted [other,self]
 */
 static VALUE nary_coerce(VALUE x, VALUE y) {
   VALUE type;
@@ -1226,9 +1228,9 @@ static VALUE nary_s_byte_size(VALUE type) {
 /*
   Returns a new 1-D array initialized from binary raw data in a string.
   @overload from_binary(string,[shape])
-  @param [String] string  Binary raw data.
-  @param [Array] shape  array of integers representing array shape.
-  @return [Numo::NArray] NArray containing binary data.
+    @param [String] string  Binary raw data.
+    @param [Array] shape  array of integers representing array shape.
+    @return [Numo::NArray] NArray containing binary data.
  */
 static VALUE nary_s_from_binary(int argc, VALUE* argv, VALUE type) {
   size_t len, str_len, byte_size;
@@ -1301,9 +1303,9 @@ static VALUE nary_s_from_binary(int argc, VALUE* argv, VALUE type) {
 /*
   Returns a new 1-D array initialized from binary raw data in a string.
   @overload store_binary(string,[offset])
-  @param [String] string  Binary raw data.
-  @param [Integer] (optional) offset  Byte offset in string.
-  @return [Integer] stored length.
+    @param [String] string  Binary raw data.
+    @param [Integer] (optional) offset  Byte offset in string.
+    @return [Integer] stored length.
  */
 static VALUE nary_store_binary(int argc, VALUE* argv, VALUE self) {
   size_t size, str_len, byte_size, offset;
@@ -1350,7 +1352,7 @@ static VALUE nary_store_binary(int argc, VALUE* argv, VALUE self) {
 /*
   Returns string containing the raw data bytes in NArray.
   @overload to_binary()
-  @return [String] String object containing binary raw data.
+    @return [String] String object containing binary raw data.
  */
 static VALUE nary_to_binary(VALUE self) {
   size_t len, offset = 0;
@@ -1376,7 +1378,7 @@ static VALUE nary_to_binary(VALUE self) {
 /*
   Dump marshal data.
   @overload marshal_dump()
-  @return [Array] Array containing marshal data.
+    @return [Array] Array containing marshal data.
  */
 static VALUE nary_marshal_dump(VALUE self) {
   VALUE a;
@@ -1410,8 +1412,8 @@ static VALUE na_inplace(VALUE self);
 /*
   Load marshal data.
   @overload marshal_load(data)
-  @param [Array] Array containing marshal data.
-  @return [nil]
+    @param [Array] Array containing marshal data.
+    @return [nil]
  */
 static VALUE nary_marshal_load(VALUE self, VALUE a) {
   VALUE v;
@@ -1423,10 +1425,12 @@ static VALUE nary_marshal_load(VALUE self, VALUE a) {
     rb_raise(rb_eArgError, "marshal array size should be 4");
   }
   if (RARRAY_AREF(a, 0) != INT2FIX(1)) {
-    rb_raise(rb_eArgError,
-             "NArray marshal version %d is not supported "
-             "(only version 1)",
-             NUM2INT(RARRAY_AREF(a, 0)));
+    rb_raise(
+      rb_eArgError,
+      "NArray marshal version %d is not supported "
+      "(only version 1)",
+      NUM2INT(RARRAY_AREF(a, 0))
+    );
   }
   na_initialize(self, RARRAY_AREF(a, 1));
   NA_FL0_SET(self, FIX2INT(RARRAY_AREF(a, 2)));
@@ -1458,8 +1462,8 @@ static VALUE nary_marshal_load(VALUE self, VALUE a) {
 /*
   Cast self to another NArray datatype.
   @overload cast_to(datatype)
-  @param [Class] datatype NArray datatype.
-  @return [Numo::NArray]
+    @param [Class] datatype NArray datatype.
+    @return [Numo::NArray]
  */
 static VALUE nary_cast_to(VALUE obj, VALUE type) {
   return rb_funcall(type, id_cast, 1, obj);
@@ -1538,7 +1542,6 @@ static VALUE na_get_reduce_flag_from_axes(VALUE na_obj, VALUE axes) {
   narg = RARRAY_LEN(axes);
   for (i = 0; i < narg; i++) {
     v = RARRAY_AREF(axes, i);
-    // printf("argv[%d]=",i);rb_p(v);
     if (TYPE(v) == T_FIXNUM) {
       beg = FIX2INT(v);
       if (beg < 0) beg += ndim;
@@ -1547,7 +1550,6 @@ static VALUE na_get_reduce_flag_from_axes(VALUE na_obj, VALUE axes) {
       }
       len = 1;
       step = 0;
-      // printf("beg=%d step=%d len=%d\n",beg,step,len);
     } else if (rb_obj_is_kind_of(v, rb_cRange)
 #ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
                || rb_obj_is_kind_of(v, rb_cArithSeq)
@@ -1611,12 +1613,14 @@ nary_reduce_options(VALUE axes, VALUE* opts, int naryc, VALUE* naryv, ndfunc_t*
 }
 
 VALUE
-nary_reduce_dimension(int argc, VALUE* argv, int naryc, VALUE* naryv, ndfunc_t* ndf, na_iter_func_t iter_nan) {
+nary_reduce_dimension(
+  int argc, VALUE* argv, int naryc, VALUE* naryv, ndfunc_t* ndf, na_iter_func_t iter_nan
+) {
   long narg;
   VALUE axes;
   VALUE kw_hash = Qnil;
-  ID kw_table[3] = {id_axis, id_keepdims, id_nan};
-  VALUE opts[3] = {Qundef, Qundef, Qundef};
+  ID kw_table[3] = { id_axis, id_keepdims, id_nan };
+  VALUE opts[3] = { Qundef, Qundef, Qundef };
 
   narg = rb_scan_args(argc, argv, "*:", &axes, &kw_hash);
   rb_get_kwargs(kw_hash, kw_table, 0, 3, opts);
@@ -1727,7 +1731,7 @@ static VALUE na_profile_set(VALUE mod, VALUE val) {
 /*
   Returns the number of rows used for NArray#inspect
   @overload inspect_rows
-  @return [Integer or nil]  the number of rows.
+    @return [Integer or nil]  the number of rows.
 */
 static VALUE na_inspect_rows(VALUE mod) {
   if (numo_na_inspect_rows > 0) {
@@ -1740,8 +1744,8 @@ static VALUE na_inspect_rows(VALUE mod) {
 /*
   Set the number of rows used for NArray#inspect
   @overload inspect_rows=(rows)
-  @param [Integer or nil] rows  the number of rows
-  @return [nil]
+    @param [Integer or nil] rows  the number of rows
+    @return [nil]
 */
 static VALUE na_inspect_rows_set(VALUE mod, VALUE num) {
   if (RTEST(num)) {
@@ -1755,7 +1759,7 @@ static VALUE na_inspect_rows_set(VALUE mod, VALUE num) {
 /*
   Returns the number of cols used for NArray#inspect
   @overload inspect_cols
-  @return [Integer or nil]  the number of cols.
+    @return [Integer or nil]  the number of cols.
 */
 static VALUE na_inspect_cols(VALUE mod) {
   if (numo_na_inspect_cols > 0) {
@@ -1768,8 +1772,8 @@ static VALUE na_inspect_cols(VALUE mod) {
 /*
   Set the number of cols used for NArray#inspect
   @overload inspect_cols=(cols)
-  @param [Integer or nil] cols  the number of cols
-  @return [nil]
+    @param [Integer or nil] cols  the number of cols
+    @return [nil]
 */
 static VALUE na_inspect_cols_set(VALUE mod, VALUE num) {
   if (RTEST(num)) {
@@ -1784,8 +1788,8 @@ static VALUE na_inspect_cols_set(VALUE mod, VALUE num) {
   Equality of self and other in view of numerical array.
   i.e., both arrays have same shape and corresponding elements are equal.
   @overload == other
-  @param [Object] other
-  @return [Boolean] true if self and other is equal.
+    @param [Object] other
+    @return [Boolean] true if self and other is equal.
 */
 static VALUE na_equal(VALUE self, volatile VALUE other) {
   volatile VALUE vbool;
@@ -1820,16 +1824,21 @@ void Init_narray(void) {
   rb_ext_ractor_safe(true);
 #endif
 
+  /**
+   * Document-module: Numo
+   *
+   * Ruby/Numo (NUmerical MOdules)
+   */
   mNumo = rb_define_module("Numo");
 
-  /*
-    Document-class: Numo::NArray
-
-    Numo::NArray is the abstract super class for
-    Numerical N-dimensional Array in the Ruby/Numo module.
-    Use Typed Subclasses of NArray (Numo::DFloat, Int32, etc)
-    to create data array instances.
-  */
+  /**
+   * Document-class: Numo::NArray
+   *
+   * Numo::NArray is the abstract super class for
+   * Numerical N-dimensional Array in the Ruby/Numo module.
+   * Use Typed Subclasses of NArray (Numo::DFloat, Int32, etc)
+   * to create data array instances.
+   */
   cNArray = rb_define_class_under(mNumo, "NArray", rb_cObject);
 
 #ifndef HAVE_RB_CCOMPLEX
@@ -1840,12 +1849,38 @@ void Init_narray(void) {
   rb_cArithSeq = rb_path2class("Enumerator::ArithmeticSequence");
 #endif
 
+  /* The version of Numo::NArray Alternative. */
   rb_define_const(cNArray, "VERSION", rb_str_new2(NARRAY_VERSION));
 
+  /**
+   * Document-class: Numo::NArray::CastError
+   *
+   * Exception raised when type casting is not possible.
+   */
   nary_eCastError = rb_define_class_under(cNArray, "CastError", rb_eStandardError);
+  /**
+   * Document-class: Numo::NArray::ShapeError
+   *
+   * Exception raised when shape is invalid.
+   */
   nary_eShapeError = rb_define_class_under(cNArray, "ShapeError", rb_eStandardError);
+  /**
+   * Document-class: Numo::NArray::OperationError
+   *
+   * Exception raised when operation is not appropriate.
+   */
   nary_eOperationError = rb_define_class_under(cNArray, "OperationError", rb_eStandardError);
+  /**
+   * Document-class: Numo::NArray::DimensionError
+   *
+   * Exception raised when dimension is invalid.
+   */
   nary_eDimensionError = rb_define_class_under(cNArray, "DimensionError", rb_eStandardError);
+  /**
+   * Document-class: Numo::NArray::ValueError
+   *
+   * Exception raised when occurred value error.
+   */
   nary_eValueError = rb_define_class_under(cNArray, "ValueError", rb_eStandardError);
 
   rb_define_singleton_method(cNArray, "debug=", na_debug_set, 1);