carray 3.0.0 → 3.0.1

data/ext/carray_kernels_monop.c

@@ -34,6 +34,7 @@ static inline VALUE BOOL2VAL (boolean8_t x) { return x ? Qtrue : Qfalse; }
 #endif
 
 #include "ca_op_powi.h"
+#include "ca_op_cmplx64.h"
 
 /* ===== monop zero ============================================ */
 
@@ -3166,7 +3167,7 @@ ca_monop_abs_i_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cabs((cmplx128_t)*p1); }
+        { (*p2) = cabsf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -3181,14 +3182,14 @@ ca_monop_abs_i_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cabs((cmplx128_t)*p1); }
+      { (*p2) = cabsf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cabs((cmplx128_t)*p1); }
+      { (*p2) = cabsf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -3211,7 +3212,7 @@ ca_monop_abs_i_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cabs((cmplx128_t)*p1); }
+        { (*p2) = cabs(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -3226,14 +3227,14 @@ ca_monop_abs_i_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cabs((cmplx128_t)*p1); }
+      { (*p2) = cabs(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cabs((cmplx128_t)*p1); }
+      { (*p2) = cabs(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -3772,18 +3773,18 @@ ca_monop_abs_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
   cmplx64_t  * __restrict__ q1 = (cmplx64_t  *) ptr1;
-  double * __restrict__ q2 = (double *) ptr2;
+  float * __restrict__ q2 = (float *) ptr2;
   ca_size_t k;
   if ( m ) {
     cmplx64_t  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     boolean8_t *pm;
     for (k=0; k<n; k++) {
       pm = m + k;
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cabs(*p1); }
+        { (*p2) = cabsf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -3793,19 +3794,19 @@ ca_monop_abs_cmplx64_t (ca_size_t n, boolean8_t *m,
        lane-parallel loads/stores.  Element-wise monop has no
        reassoc, so ε-close numeric contract is unaffected. */
     cmplx64_t  *p1;
-    double *p2;
+    float *p2;
     (void)p1; (void)p2;
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cabs(*p1); }
+      { (*p2) = cabsf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cabs(*p1); }
+      { (*p2) = cabsf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -3891,7 +3892,7 @@ int8_t ca_monop_abs_out_data_type[CA_NTYPE] = {
   CA_FLOAT32,
   CA_FLOAT64,
   -1,
-  CA_FLOAT64,
+  CA_FLOAT32,
   CA_FLOAT64,
   -1,
   -1,
@@ -4380,18 +4381,18 @@ ca_monop_abs2_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
   cmplx64_t  * __restrict__ q1 = (cmplx64_t  *) ptr1;
-  double * __restrict__ q2 = (double *) ptr2;
+  float * __restrict__ q2 = (float *) ptr2;
   ca_size_t k;
   if ( m ) {
     cmplx64_t  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     boolean8_t *pm;
     for (k=0; k<n; k++) {
       pm = m + k;
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { { double _r = creal(*p1); double _i = cimag(*p1); (*p2) = _r * _r + _i * _i; } }
+        { { float _r = crealf(*p1); float _i = cimagf(*p1); (*p2) = _r * _r + _i * _i; } }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -4401,19 +4402,19 @@ ca_monop_abs2_cmplx64_t (ca_size_t n, boolean8_t *m,
        lane-parallel loads/stores.  Element-wise monop has no
        reassoc, so ε-close numeric contract is unaffected. */
     cmplx64_t  *p1;
-    double *p2;
+    float *p2;
     (void)p1; (void)p2;
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { { double _r = creal(*p1); double _i = cimag(*p1); (*p2) = _r * _r + _i * _i; } }
+      { { float _r = crealf(*p1); float _i = cimagf(*p1); (*p2) = _r * _r + _i * _i; } }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     for (k=0; k<n; k++) {
-      { { double _r = creal(*p1); double _i = cimag(*p1); (*p2) = _r * _r + _i * _i; } }
+      { { float _r = crealf(*p1); float _i = cimagf(*p1); (*p2) = _r * _r + _i * _i; } }
       p1 += i1; p2 += i2;
     }
   }
@@ -4499,7 +4500,7 @@ int8_t ca_monop_abs2_out_data_type[CA_NTYPE] = {
   CA_FLOAT32,
   CA_FLOAT64,
   -1,
-  CA_FLOAT64,
+  CA_FLOAT32,
   CA_FLOAT64,
   -1,
   -1,
@@ -4999,7 +5000,7 @@ ca_monop_conj_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = conj(*p1); }
+        { (*p2) = conjf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -5014,14 +5015,14 @@ ca_monop_conj_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = conj(*p1); }
+      { (*p2) = conjf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = conj(*p1); }
+      { (*p2) = conjf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -5661,7 +5662,7 @@ ca_monop_sign_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { { double _m = cabs(*p1); (*p2) = (_m == 0.0) ? 0 : ((*p1) / _m); } }
+        { { float _m = cabsf(*p1); (*p2) = (_m == 0) ? 0 : ((*p1) / _m); } }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -5676,14 +5677,14 @@ ca_monop_sign_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { { double _m = cabs(*p1); (*p2) = (_m == 0.0) ? 0 : ((*p1) / _m); } }
+      { { float _m = cabsf(*p1); (*p2) = (_m == 0) ? 0 : ((*p1) / _m); } }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { { double _m = cabs(*p1); (*p2) = (_m == 0.0) ? 0 : ((*p1) / _m); } }
+      { { float _m = cabsf(*p1); (*p2) = (_m == 0) ? 0 : ((*p1) / _m); } }
       p1 += i1; p2 += i2;
     }
   }
@@ -5706,7 +5707,7 @@ ca_monop_sign_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { { double _m = cabs(*p1); (*p2) = (_m == 0.0) ? 0 : ((*p1) / _m); } }
+        { { double _m = cabs(*p1); (*p2) = (_m == 0) ? 0 : ((*p1) / _m); } }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -5721,14 +5722,14 @@ ca_monop_sign_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { { double _m = cabs(*p1); (*p2) = (_m == 0.0) ? 0 : ((*p1) / _m); } }
+      { { double _m = cabs(*p1); (*p2) = (_m == 0) ? 0 : ((*p1) / _m); } }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { { double _m = cabs(*p1); (*p2) = (_m == 0.0) ? 0 : ((*p1) / _m); } }
+      { { double _m = cabs(*p1); (*p2) = (_m == 0) ? 0 : ((*p1) / _m); } }
       p1 += i1; p2 += i2;
     }
   }
@@ -6132,11 +6133,11 @@ ca_monop_arg_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
   float  * __restrict__ q1 = (float  *) ptr1;
-  double * __restrict__ q2 = (double *) ptr2;
+  float * __restrict__ q2 = (float *) ptr2;
   ca_size_t k;
   if ( m ) {
     float  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     boolean8_t *pm;
     for (k=0; k<n; k++) {
       pm = m + k;
@@ -6153,7 +6154,7 @@ ca_monop_arg_float32_t (ca_size_t n, boolean8_t *m,
        lane-parallel loads/stores.  Element-wise monop has no
        reassoc, so ε-close numeric contract is unaffected. */
     float  *p1;
-    double *p2;
+    float *p2;
     (void)p1; (void)p2;
     #pragma omp simd
     for (k=0; k<n; k++) {
@@ -6163,7 +6164,7 @@ ca_monop_arg_float32_t (ca_size_t n, boolean8_t *m,
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     for (k=0; k<n; k++) {
       { (*p2) = carg((cmplx128_t)(*p1)); }
       p1 += i1; p2 += i2;
@@ -6222,18 +6223,18 @@ ca_monop_arg_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
   cmplx64_t  * __restrict__ q1 = (cmplx64_t  *) ptr1;
-  double * __restrict__ q2 = (double *) ptr2;
+  float * __restrict__ q2 = (float *) ptr2;
   ca_size_t k;
   if ( m ) {
     cmplx64_t  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     boolean8_t *pm;
     for (k=0; k<n; k++) {
       pm = m + k;
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = carg(*p1); }
+        { (*p2) = cargf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -6243,19 +6244,19 @@ ca_monop_arg_cmplx64_t (ca_size_t n, boolean8_t *m,
        lane-parallel loads/stores.  Element-wise monop has no
        reassoc, so ε-close numeric contract is unaffected. */
     cmplx64_t  *p1;
-    double *p2;
+    float *p2;
     (void)p1; (void)p2;
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = carg(*p1); }
+      { (*p2) = cargf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
-    double *p2 = q2;
+    float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = carg(*p1); }
+      { (*p2) = cargf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -6338,10 +6339,10 @@ int8_t ca_monop_arg_out_data_type[CA_NTYPE] = {
   CA_FLOAT64,
   CA_FLOAT64,
   CA_FLOAT64,
-  CA_FLOAT64,
+  CA_FLOAT32,
   CA_FLOAT64,
   -1,
-  CA_FLOAT64,
+  CA_FLOAT32,
   CA_FLOAT64,
   -1,
   -1,
@@ -6761,7 +6762,7 @@ ca_monop_imag_i_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cimag(*p1); }
+        { (*p2) = cimagf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -6776,14 +6777,14 @@ ca_monop_imag_i_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cimag(*p1); }
+      { (*p2) = cimagf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cimag(*p1); }
+      { (*p2) = cimagf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -6928,7 +6929,7 @@ ca_monop_arg_i_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = (float)carg((cmplx128_t)(*p1)); }
+        { (*p2) = cargf((cmplx64_t)(*p1)); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -6943,14 +6944,14 @@ ca_monop_arg_i_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = (float)carg((cmplx128_t)(*p1)); }
+      { (*p2) = cargf((cmplx64_t)(*p1)); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = (float)carg((cmplx128_t)(*p1)); }
+      { (*p2) = cargf((cmplx64_t)(*p1)); }
       p1 += i1; p2 += i2;
     }
   }
@@ -7018,7 +7019,7 @@ ca_monop_arg_i_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = (float)carg((cmplx128_t)(*p1)); }
+        { (*p2) = cargf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -7033,14 +7034,14 @@ ca_monop_arg_i_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = (float)carg((cmplx128_t)(*p1)); }
+      { (*p2) = cargf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = (float)carg((cmplx128_t)(*p1)); }
+      { (*p2) = cargf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -9695,7 +9696,7 @@ ca_monop_rcp_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = 1/(*p1); }
+        { (*p2) = op_crcp_cmplx64(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -9710,14 +9711,14 @@ ca_monop_rcp_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = 1/(*p1); }
+      { (*p2) = op_crcp_cmplx64(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = 1/(*p1); }
+      { (*p2) = op_crcp_cmplx64(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -9862,7 +9863,7 @@ ca_monop_sqrt_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = sqrt(*p1); }
+        { (*p2) = sqrtf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -9877,14 +9878,14 @@ ca_monop_sqrt_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = sqrt(*p1); }
+      { (*p2) = sqrtf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = sqrt(*p1); }
+      { (*p2) = sqrtf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -9952,7 +9953,7 @@ ca_monop_sqrt_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = csqrt(*p1); }
+        { (*p2) = csqrtf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -9967,14 +9968,14 @@ ca_monop_sqrt_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = csqrt(*p1); }
+      { (*p2) = csqrtf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = csqrt(*p1); }
+      { (*p2) = csqrtf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -10148,7 +10149,7 @@ ca_monop_exp_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = exp(*p1); }
+        { (*p2) = expf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10163,14 +10164,14 @@ ca_monop_exp_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = exp(*p1); }
+      { (*p2) = expf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = exp(*p1); }
+      { (*p2) = expf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -10238,7 +10239,7 @@ ca_monop_exp_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cexp(*p1); }
+        { (*p2) = cexpf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10253,14 +10254,14 @@ ca_monop_exp_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cexp(*p1); }
+      { (*p2) = cexpf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cexp(*p1); }
+      { (*p2) = cexpf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -10415,10 +10416,10 @@ VALUE rb_ca_exp_bang (VALUE self)
 VALUE rb_cmath_exp (VALUE mod, VALUE arg)
 { return ca_math_call(mod, arg, rb_intern("exp")); }
 
-/* ===== monop log ============================================ */
+/* ===== monop sin ============================================ */
 
 static void
-ca_monop_log_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_sin_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10434,7 +10435,7 @@ ca_monop_log_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = log(*p1); }
+        { (*p2) = sinf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10449,21 +10450,21 @@ ca_monop_log_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = log(*p1); }
+      { (*p2) = sinf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = log(*p1); }
+      { (*p2) = sinf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_log_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_sin_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10479,7 +10480,7 @@ ca_monop_log_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = log(*p1); }
+        { (*p2) = sin(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10494,21 +10495,21 @@ ca_monop_log_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = log(*p1); }
+      { (*p2) = sin(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = log(*p1); }
+      { (*p2) = sin(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_log_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_sin_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10524,7 +10525,7 @@ ca_monop_log_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = clog(*p1); }
+        { (*p2) = csinf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10539,21 +10540,21 @@ ca_monop_log_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = clog(*p1); }
+      { (*p2) = csinf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = clog(*p1); }
+      { (*p2) = csinf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_log_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_sin_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10569,7 +10570,7 @@ ca_monop_log_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = clog(*p1); }
+        { (*p2) = csin(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10584,21 +10585,21 @@ ca_monop_log_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = clog(*p1); }
+      { (*p2) = csin(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = clog(*p1); }
+      { (*p2) = csin(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_log_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_sin_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10617,9 +10618,9 @@ ca_monop_log_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(log(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(sin(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("log"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("sin"), 0);
   }
 }
  }
@@ -10640,9 +10641,9 @@ ca_monop_log_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(log(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(sin(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("log"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("sin"), 0);
   }
 }
  }
@@ -10655,9 +10656,9 @@ ca_monop_log_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(log(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(sin(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("log"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("sin"), 0);
   }
 }
  }
@@ -10667,7 +10668,7 @@ ca_monop_log_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_log[CA_NTYPE] = {
+ca_monop_sin[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -10678,33 +10679,33 @@ ca_monop_log[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_log_float32_t,
-  ca_monop_log_float64_t,
+  ca_monop_sin_float32_t,
+  ca_monop_sin_float64_t,
   ca_monop_not_implement,
-  ca_monop_log_cmplx64_t,
-  ca_monop_log_cmplx128_t,
+  ca_monop_sin_cmplx64_t,
+  ca_monop_sin_cmplx128_t,
   ca_monop_not_implement,
-  ca_monop_log_VALUE,
+  ca_monop_sin_VALUE,
 };
 
-VALUE rb_ca_log (VALUE self)
+VALUE rb_ca_sin (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_log);
+  return rb_ca_call_monop(self, ca_monop_sin);
 }
 
-VALUE rb_ca_log_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_log); }
+VALUE rb_ca_sin_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_sin); }
 
-VALUE rb_cmath_log (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("log")); }
+VALUE rb_cmath_sin (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("sin")); }
 
-/* ===== monop sin ============================================ */
+/* ===== monop cos ============================================ */
 
 static void
-ca_monop_sin_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_cos_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10720,7 +10721,7 @@ ca_monop_sin_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = sin(*p1); }
+        { (*p2) = cosf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10735,21 +10736,21 @@ ca_monop_sin_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = sin(*p1); }
+      { (*p2) = cosf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = sin(*p1); }
+      { (*p2) = cosf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_sin_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_cos_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10765,7 +10766,7 @@ ca_monop_sin_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = sin(*p1); }
+        { (*p2) = cos(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10780,21 +10781,21 @@ ca_monop_sin_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = sin(*p1); }
+      { (*p2) = cos(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = sin(*p1); }
+      { (*p2) = cos(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_sin_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_cos_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10810,7 +10811,7 @@ ca_monop_sin_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = csin(*p1); }
+        { (*p2) = ccosf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10825,21 +10826,21 @@ ca_monop_sin_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = csin(*p1); }
+      { (*p2) = ccosf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = csin(*p1); }
+      { (*p2) = ccosf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_sin_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_cos_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10855,7 +10856,7 @@ ca_monop_sin_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = csin(*p1); }
+        { (*p2) = ccos(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -10870,21 +10871,21 @@ ca_monop_sin_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = csin(*p1); }
+      { (*p2) = ccos(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = csin(*p1); }
+      { (*p2) = ccos(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_sin_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_cos_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -10903,9 +10904,9 @@ ca_monop_sin_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(sin(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(cos(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("sin"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("cos"), 0);
   }
 }
  }
@@ -10926,9 +10927,9 @@ ca_monop_sin_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(sin(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(cos(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("sin"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("cos"), 0);
   }
 }
  }
@@ -10941,9 +10942,9 @@ ca_monop_sin_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(sin(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(cos(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("sin"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("cos"), 0);
   }
 }
  }
@@ -10953,7 +10954,7 @@ ca_monop_sin_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_sin[CA_NTYPE] = {
+ca_monop_cos[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -10964,33 +10965,33 @@ ca_monop_sin[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_sin_float32_t,
-  ca_monop_sin_float64_t,
+  ca_monop_cos_float32_t,
+  ca_monop_cos_float64_t,
   ca_monop_not_implement,
-  ca_monop_sin_cmplx64_t,
-  ca_monop_sin_cmplx128_t,
+  ca_monop_cos_cmplx64_t,
+  ca_monop_cos_cmplx128_t,
   ca_monop_not_implement,
-  ca_monop_sin_VALUE,
+  ca_monop_cos_VALUE,
 };
 
-VALUE rb_ca_sin (VALUE self)
+VALUE rb_ca_cos (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_sin);
+  return rb_ca_call_monop(self, ca_monop_cos);
 }
 
-VALUE rb_ca_sin_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_sin); }
+VALUE rb_ca_cos_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_cos); }
 
-VALUE rb_cmath_sin (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("sin")); }
+VALUE rb_cmath_cos (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("cos")); }
 
-/* ===== monop cos ============================================ */
+/* ===== monop tan ============================================ */
 
 static void
-ca_monop_cos_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_tan_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11006,7 +11007,7 @@ ca_monop_cos_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cos(*p1); }
+        { (*p2) = tanf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11021,21 +11022,21 @@ ca_monop_cos_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cos(*p1); }
+      { (*p2) = tanf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cos(*p1); }
+      { (*p2) = tanf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_cos_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_tan_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11051,7 +11052,7 @@ ca_monop_cos_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cos(*p1); }
+        { (*p2) = tan(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11066,21 +11067,21 @@ ca_monop_cos_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cos(*p1); }
+      { (*p2) = tan(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cos(*p1); }
+      { (*p2) = tan(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_cos_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_tan_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11096,7 +11097,7 @@ ca_monop_cos_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = ccos(*p1); }
+        { (*p2) = ctanf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11111,21 +11112,21 @@ ca_monop_cos_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = ccos(*p1); }
+      { (*p2) = ctanf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = ccos(*p1); }
+      { (*p2) = ctanf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_cos_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_tan_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11141,7 +11142,7 @@ ca_monop_cos_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = ccos(*p1); }
+        { (*p2) = ctan(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11156,21 +11157,21 @@ ca_monop_cos_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = ccos(*p1); }
+      { (*p2) = ctan(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = ccos(*p1); }
+      { (*p2) = ctan(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_cos_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_tan_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11189,9 +11190,9 @@ ca_monop_cos_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(cos(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(tan(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("cos"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("tan"), 0);
   }
 }
  }
@@ -11212,9 +11213,9 @@ ca_monop_cos_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(cos(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(tan(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("cos"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("tan"), 0);
   }
 }
  }
@@ -11227,9 +11228,9 @@ ca_monop_cos_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(cos(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(tan(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("cos"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("tan"), 0);
   }
 }
  }
@@ -11239,8 +11240,7 @@ ca_monop_cos_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_cos[CA_NTYPE] = {
-  ca_monop_not_implement,
+ca_monop_tan[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -11250,33 +11250,34 @@ ca_monop_cos[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_cos_float32_t,
-  ca_monop_cos_float64_t,
   ca_monop_not_implement,
-  ca_monop_cos_cmplx64_t,
-  ca_monop_cos_cmplx128_t,
+  ca_monop_tan_float32_t,
+  ca_monop_tan_float64_t,
   ca_monop_not_implement,
-  ca_monop_cos_VALUE,
+  ca_monop_tan_cmplx64_t,
+  ca_monop_tan_cmplx128_t,
+  ca_monop_not_implement,
+  ca_monop_tan_VALUE,
 };
 
-VALUE rb_ca_cos (VALUE self)
+VALUE rb_ca_tan (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_cos);
+  return rb_ca_call_monop(self, ca_monop_tan);
 }
 
-VALUE rb_ca_cos_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_cos); }
+VALUE rb_ca_tan_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_tan); }
 
-VALUE rb_cmath_cos (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("cos")); }
+VALUE rb_cmath_tan (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("tan")); }
 
-/* ===== monop tan ============================================ */
+/* ===== monop asin ============================================ */
 
 static void
-ca_monop_tan_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_asin_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11292,7 +11293,7 @@ ca_monop_tan_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = tan(*p1); }
+        { (*p2) = asinf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11307,21 +11308,21 @@ ca_monop_tan_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = tan(*p1); }
+      { (*p2) = asinf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = tan(*p1); }
+      { (*p2) = asinf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_tan_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_asin_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11337,7 +11338,7 @@ ca_monop_tan_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = tan(*p1); }
+        { (*p2) = asin(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11352,21 +11353,21 @@ ca_monop_tan_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = tan(*p1); }
+      { (*p2) = asin(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = tan(*p1); }
+      { (*p2) = asin(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_tan_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_asin_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11382,7 +11383,7 @@ ca_monop_tan_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = ctan(*p1); }
+        { (*p2) = casinf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11397,21 +11398,21 @@ ca_monop_tan_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = ctan(*p1); }
+      { (*p2) = casinf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = ctan(*p1); }
+      { (*p2) = casinf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_tan_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_asin_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11427,7 +11428,7 @@ ca_monop_tan_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = ctan(*p1); }
+        { (*p2) = casin(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11442,21 +11443,21 @@ ca_monop_tan_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = ctan(*p1); }
+      { (*p2) = casin(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = ctan(*p1); }
+      { (*p2) = casin(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_tan_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_asin_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11475,9 +11476,9 @@ ca_monop_tan_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(tan(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(asin(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("tan"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("asin"), 0);
   }
 }
  }
@@ -11498,9 +11499,9 @@ ca_monop_tan_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(tan(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(asin(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("tan"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("asin"), 0);
   }
 }
  }
@@ -11513,9 +11514,9 @@ ca_monop_tan_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(tan(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(asin(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("tan"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("asin"), 0);
   }
 }
  }
@@ -11525,7 +11526,7 @@ ca_monop_tan_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_tan[CA_NTYPE] = {
+ca_monop_asin[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -11536,33 +11537,33 @@ ca_monop_tan[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_tan_float32_t,
-  ca_monop_tan_float64_t,
+  ca_monop_asin_float32_t,
+  ca_monop_asin_float64_t,
   ca_monop_not_implement,
-  ca_monop_tan_cmplx64_t,
-  ca_monop_tan_cmplx128_t,
+  ca_monop_asin_cmplx64_t,
+  ca_monop_asin_cmplx128_t,
   ca_monop_not_implement,
-  ca_monop_tan_VALUE,
+  ca_monop_asin_VALUE,
 };
 
-VALUE rb_ca_tan (VALUE self)
+VALUE rb_ca_asin (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_tan);
+  return rb_ca_call_monop(self, ca_monop_asin);
 }
 
-VALUE rb_ca_tan_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_tan); }
+VALUE rb_ca_asin_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_asin); }
 
-VALUE rb_cmath_tan (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("tan")); }
+VALUE rb_cmath_asin (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("asin")); }
 
-/* ===== monop asin ============================================ */
+/* ===== monop acos ============================================ */
 
 static void
-ca_monop_asin_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_acos_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11578,7 +11579,7 @@ ca_monop_asin_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = asin(*p1); }
+        { (*p2) = acosf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11593,21 +11594,21 @@ ca_monop_asin_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = asin(*p1); }
+      { (*p2) = acosf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = asin(*p1); }
+      { (*p2) = acosf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_asin_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_acos_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11623,7 +11624,7 @@ ca_monop_asin_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = asin(*p1); }
+        { (*p2) = acos(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11638,21 +11639,21 @@ ca_monop_asin_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = asin(*p1); }
+      { (*p2) = acos(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = asin(*p1); }
+      { (*p2) = acos(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_asin_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_acos_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11668,7 +11669,7 @@ ca_monop_asin_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = casin(*p1); }
+        { (*p2) = cacosf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11683,21 +11684,21 @@ ca_monop_asin_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = casin(*p1); }
+      { (*p2) = cacosf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = casin(*p1); }
+      { (*p2) = cacosf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_asin_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_acos_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11713,7 +11714,7 @@ ca_monop_asin_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = casin(*p1); }
+        { (*p2) = cacos(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11728,21 +11729,21 @@ ca_monop_asin_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = casin(*p1); }
+      { (*p2) = cacos(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = casin(*p1); }
+      { (*p2) = cacos(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_asin_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_acos_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11761,9 +11762,9 @@ ca_monop_asin_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(asin(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(acos(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("asin"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("acos"), 0);
   }
 }
  }
@@ -11784,9 +11785,9 @@ ca_monop_asin_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(asin(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(acos(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("asin"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("acos"), 0);
   }
 }
  }
@@ -11799,9 +11800,9 @@ ca_monop_asin_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(asin(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(acos(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("asin"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("acos"), 0);
   }
 }
  }
@@ -11811,7 +11812,7 @@ ca_monop_asin_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_asin[CA_NTYPE] = {
+ca_monop_acos[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -11822,33 +11823,33 @@ ca_monop_asin[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_asin_float32_t,
-  ca_monop_asin_float64_t,
+  ca_monop_acos_float32_t,
+  ca_monop_acos_float64_t,
   ca_monop_not_implement,
-  ca_monop_asin_cmplx64_t,
-  ca_monop_asin_cmplx128_t,
+  ca_monop_acos_cmplx64_t,
+  ca_monop_acos_cmplx128_t,
   ca_monop_not_implement,
-  ca_monop_asin_VALUE,
+  ca_monop_acos_VALUE,
 };
 
-VALUE rb_ca_asin (VALUE self)
+VALUE rb_ca_acos (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_asin);
+  return rb_ca_call_monop(self, ca_monop_acos);
 }
 
-VALUE rb_ca_asin_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_asin); }
+VALUE rb_ca_acos_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_acos); }
 
-VALUE rb_cmath_asin (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("asin")); }
+VALUE rb_cmath_acos (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("acos")); }
 
-/* ===== monop acos ============================================ */
+/* ===== monop atan ============================================ */
 
 static void
-ca_monop_acos_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_atan_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11864,7 +11865,7 @@ ca_monop_acos_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = acos(*p1); }
+        { (*p2) = atanf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11879,21 +11880,21 @@ ca_monop_acos_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = acos(*p1); }
+      { (*p2) = atanf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = acos(*p1); }
+      { (*p2) = atanf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_acos_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_atan_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11909,7 +11910,7 @@ ca_monop_acos_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = acos(*p1); }
+        { (*p2) = atan(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11924,21 +11925,21 @@ ca_monop_acos_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = acos(*p1); }
+      { (*p2) = atan(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = acos(*p1); }
+      { (*p2) = atan(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_acos_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_atan_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11954,7 +11955,7 @@ ca_monop_acos_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cacos(*p1); }
+        { (*p2) = catanf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -11969,21 +11970,21 @@ ca_monop_acos_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cacos(*p1); }
+      { (*p2) = catanf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cacos(*p1); }
+      { (*p2) = catanf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_acos_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_atan_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -11999,7 +12000,7 @@ ca_monop_acos_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cacos(*p1); }
+        { (*p2) = catan(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12014,21 +12015,21 @@ ca_monop_acos_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cacos(*p1); }
+      { (*p2) = catan(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cacos(*p1); }
+      { (*p2) = catan(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_acos_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_atan_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -12047,9 +12048,9 @@ ca_monop_acos_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(acos(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(atan(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("acos"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("atan"), 0);
   }
 }
  }
@@ -12070,9 +12071,9 @@ ca_monop_acos_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(acos(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(atan(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("acos"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("atan"), 0);
   }
 }
  }
@@ -12085,9 +12086,9 @@ ca_monop_acos_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(acos(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(atan(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("acos"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("atan"), 0);
   }
 }
  }
@@ -12097,7 +12098,7 @@ ca_monop_acos_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_acos[CA_NTYPE] = {
+ca_monop_atan[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -12108,33 +12109,33 @@ ca_monop_acos[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_acos_float32_t,
-  ca_monop_acos_float64_t,
+  ca_monop_atan_float32_t,
+  ca_monop_atan_float64_t,
   ca_monop_not_implement,
-  ca_monop_acos_cmplx64_t,
-  ca_monop_acos_cmplx128_t,
+  ca_monop_atan_cmplx64_t,
+  ca_monop_atan_cmplx128_t,
   ca_monop_not_implement,
-  ca_monop_acos_VALUE,
+  ca_monop_atan_VALUE,
 };
 
-VALUE rb_ca_acos (VALUE self)
+VALUE rb_ca_atan (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_acos);
+  return rb_ca_call_monop(self, ca_monop_atan);
 }
 
-VALUE rb_ca_acos_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_acos); }
+VALUE rb_ca_atan_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_atan); }
 
-VALUE rb_cmath_acos (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("acos")); }
+VALUE rb_cmath_atan (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("atan")); }
 
-/* ===== monop atan ============================================ */
+/* ===== monop log ============================================ */
 
 static void
-ca_monop_atan_float32_t (ca_size_t n, boolean8_t *m,
+ca_monop_log_float32_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -12150,7 +12151,7 @@ ca_monop_atan_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = atan(*p1); }
+        { (*p2) = logf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12165,21 +12166,21 @@ ca_monop_atan_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = atan(*p1); }
+      { (*p2) = logf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = atan(*p1); }
+      { (*p2) = logf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_atan_float64_t (ca_size_t n, boolean8_t *m,
+ca_monop_log_float64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -12195,7 +12196,7 @@ ca_monop_atan_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = atan(*p1); }
+        { (*p2) = log(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12210,21 +12211,21 @@ ca_monop_atan_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = atan(*p1); }
+      { (*p2) = log(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = atan(*p1); }
+      { (*p2) = log(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_atan_cmplx64_t (ca_size_t n, boolean8_t *m,
+ca_monop_log_cmplx64_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -12240,7 +12241,7 @@ ca_monop_atan_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = catan(*p1); }
+        { (*p2) = clog(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12255,21 +12256,21 @@ ca_monop_atan_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = catan(*p1); }
+      { (*p2) = clog(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = catan(*p1); }
+      { (*p2) = clog(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_atan_cmplx128_t (ca_size_t n, boolean8_t *m,
+ca_monop_log_cmplx128_t (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -12285,7 +12286,7 @@ ca_monop_atan_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = catan(*p1); }
+        { (*p2) = clog(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12300,21 +12301,21 @@ ca_monop_atan_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = catan(*p1); }
+      { (*p2) = clog(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = catan(*p1); }
+      { (*p2) = clog(*p1); }
       p1 += i1; p2 += i2;
     }
   }
 }
 
 static void
-ca_monop_atan_VALUE (ca_size_t n, boolean8_t *m,
+ca_monop_log_VALUE (ca_size_t n, boolean8_t *m,
                                  char * __restrict__ ptr1, ca_size_t i1,
                                  char * __restrict__ ptr2, ca_size_t i2)
 {
@@ -12333,9 +12334,9 @@ ca_monop_atan_VALUE (ca_size_t n, boolean8_t *m,
         { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(atan(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(log(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("atan"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("log"), 0);
   }
 }
  }
@@ -12356,9 +12357,9 @@ ca_monop_atan_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(atan(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(log(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("atan"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("log"), 0);
   }
 }
  }
@@ -12371,9 +12372,9 @@ ca_monop_atan_VALUE (ca_size_t n, boolean8_t *m,
       { {
   VALUE _obj_arg = (*p1);
   if (RB_FLOAT_TYPE_P(_obj_arg) || RB_INTEGER_TYPE_P(_obj_arg) || RB_TYPE_P(_obj_arg, T_RATIONAL)) {
-    (*p2) = DBL2NUM(atan(NUM2DBL(_obj_arg)));
+    (*p2) = DBL2NUM(log(NUM2DBL(_obj_arg)));
   } else {
-    (*p2) = rb_funcall(_obj_arg, rb_intern("atan"), 0);
+    (*p2) = rb_funcall(_obj_arg, rb_intern("log"), 0);
   }
 }
  }
@@ -12383,7 +12384,7 @@ ca_monop_atan_VALUE (ca_size_t n, boolean8_t *m,
 }
 
 ca_monop_func_t
-ca_monop_atan[CA_NTYPE] = {
+ca_monop_log[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
@@ -12394,28 +12395,28 @@ ca_monop_atan[CA_NTYPE] = {
   ca_monop_not_implement,
   ca_monop_not_implement,
   ca_monop_not_implement,
-  ca_monop_atan_float32_t,
-  ca_monop_atan_float64_t,
+  ca_monop_log_float32_t,
+  ca_monop_log_float64_t,
   ca_monop_not_implement,
-  ca_monop_atan_cmplx64_t,
-  ca_monop_atan_cmplx128_t,
+  ca_monop_log_cmplx64_t,
+  ca_monop_log_cmplx128_t,
   ca_monop_not_implement,
-  ca_monop_atan_VALUE,
+  ca_monop_log_VALUE,
 };
 
-VALUE rb_ca_atan (VALUE self)
+VALUE rb_ca_log (VALUE self)
 {
   if ( rb_ca_is_integer_type(self) ) {
     self = rb_ca_wrap_readonly(self, INT2NUM(CA_FLOAT64));
   }
-  return rb_ca_call_monop(self, ca_monop_atan);
+  return rb_ca_call_monop(self, ca_monop_log);
 }
 
-VALUE rb_ca_atan_bang (VALUE self)
-{ return rb_ca_call_monop_bang(self, ca_monop_atan); }
+VALUE rb_ca_log_bang (VALUE self)
+{ return rb_ca_call_monop_bang(self, ca_monop_log); }
 
-VALUE rb_cmath_atan (VALUE mod, VALUE arg)
-{ return ca_math_call(mod, arg, rb_intern("atan")); }
+VALUE rb_cmath_log (VALUE mod, VALUE arg)
+{ return ca_math_call(mod, arg, rb_intern("log")); }
 
 /* ===== monop exp2 ============================================ */
 
@@ -12436,7 +12437,7 @@ ca_monop_exp2_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = exp2(*p1); }
+        { (*p2) = exp2f(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12451,14 +12452,14 @@ ca_monop_exp2_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = exp2(*p1); }
+      { (*p2) = exp2f(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = exp2(*p1); }
+      { (*p2) = exp2f(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -12722,7 +12723,7 @@ ca_monop_log10_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = log10(*p1); }
+        { (*p2) = log10f(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12737,14 +12738,14 @@ ca_monop_log10_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = log10(*p1); }
+      { (*p2) = log10f(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = log10(*p1); }
+      { (*p2) = log10f(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -12918,7 +12919,7 @@ ca_monop_log2_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = log2(*p1); }
+        { (*p2) = log2f(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -12933,14 +12934,14 @@ ca_monop_log2_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = log2(*p1); }
+      { (*p2) = log2f(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = log2(*p1); }
+      { (*p2) = log2f(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13114,7 +13115,7 @@ ca_monop_logb_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = logb(*p1); }
+        { (*p2) = logbf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13129,14 +13130,14 @@ ca_monop_logb_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = logb(*p1); }
+      { (*p2) = logbf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = logb(*p1); }
+      { (*p2) = logbf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13310,7 +13311,7 @@ ca_monop_exp10_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = pow(10, (*p1)); }
+        { (*p2) = powf(10, (*p1)); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13325,14 +13326,14 @@ ca_monop_exp10_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = pow(10, (*p1)); }
+      { (*p2) = powf(10, (*p1)); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = pow(10, (*p1)); }
+      { (*p2) = powf(10, (*p1)); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13596,7 +13597,7 @@ ca_monop_sinh_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = sinh(*p1); }
+        { (*p2) = sinhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13611,14 +13612,14 @@ ca_monop_sinh_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = sinh(*p1); }
+      { (*p2) = sinhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = sinh(*p1); }
+      { (*p2) = sinhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13686,7 +13687,7 @@ ca_monop_sinh_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = sinh(*p1); }
+        { (*p2) = csinhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13701,14 +13702,14 @@ ca_monop_sinh_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = sinh(*p1); }
+      { (*p2) = csinhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = sinh(*p1); }
+      { (*p2) = csinhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13731,7 +13732,7 @@ ca_monop_sinh_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = sinh(*p1); }
+        { (*p2) = csinh(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13746,14 +13747,14 @@ ca_monop_sinh_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = sinh(*p1); }
+      { (*p2) = csinh(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = sinh(*p1); }
+      { (*p2) = csinh(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13882,7 +13883,7 @@ ca_monop_cosh_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cosh(*p1); }
+        { (*p2) = coshf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13897,14 +13898,14 @@ ca_monop_cosh_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cosh(*p1); }
+      { (*p2) = coshf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cosh(*p1); }
+      { (*p2) = coshf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -13972,7 +13973,7 @@ ca_monop_cosh_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cosh(*p1); }
+        { (*p2) = ccoshf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -13987,14 +13988,14 @@ ca_monop_cosh_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cosh(*p1); }
+      { (*p2) = ccoshf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cosh(*p1); }
+      { (*p2) = ccoshf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14017,7 +14018,7 @@ ca_monop_cosh_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = cosh(*p1); }
+        { (*p2) = ccosh(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14032,14 +14033,14 @@ ca_monop_cosh_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = cosh(*p1); }
+      { (*p2) = ccosh(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = cosh(*p1); }
+      { (*p2) = ccosh(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14168,7 +14169,7 @@ ca_monop_tanh_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = tanh(*p1); }
+        { (*p2) = tanhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14183,14 +14184,14 @@ ca_monop_tanh_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = tanh(*p1); }
+      { (*p2) = tanhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = tanh(*p1); }
+      { (*p2) = tanhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14258,7 +14259,7 @@ ca_monop_tanh_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = tanh(*p1); }
+        { (*p2) = ctanhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14273,14 +14274,14 @@ ca_monop_tanh_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = tanh(*p1); }
+      { (*p2) = ctanhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = tanh(*p1); }
+      { (*p2) = ctanhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14303,7 +14304,7 @@ ca_monop_tanh_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = tanh(*p1); }
+        { (*p2) = ctanh(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14318,14 +14319,14 @@ ca_monop_tanh_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = tanh(*p1); }
+      { (*p2) = ctanh(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = tanh(*p1); }
+      { (*p2) = ctanh(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14454,7 +14455,7 @@ ca_monop_asinh_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = asinh(*p1); }
+        { (*p2) = asinhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14469,14 +14470,14 @@ ca_monop_asinh_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = asinh(*p1); }
+      { (*p2) = asinhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = asinh(*p1); }
+      { (*p2) = asinhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14544,7 +14545,7 @@ ca_monop_asinh_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = asinh(*p1); }
+        { (*p2) = casinhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14559,14 +14560,14 @@ ca_monop_asinh_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = asinh(*p1); }
+      { (*p2) = casinhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = asinh(*p1); }
+      { (*p2) = casinhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14589,7 +14590,7 @@ ca_monop_asinh_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = asinh(*p1); }
+        { (*p2) = casinh(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14604,14 +14605,14 @@ ca_monop_asinh_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = asinh(*p1); }
+      { (*p2) = casinh(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = asinh(*p1); }
+      { (*p2) = casinh(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14740,7 +14741,7 @@ ca_monop_acosh_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = acosh(*p1); }
+        { (*p2) = acoshf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14755,14 +14756,14 @@ ca_monop_acosh_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = acosh(*p1); }
+      { (*p2) = acoshf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = acosh(*p1); }
+      { (*p2) = acoshf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14830,7 +14831,7 @@ ca_monop_acosh_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = acosh(*p1); }
+        { (*p2) = cacoshf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14845,14 +14846,14 @@ ca_monop_acosh_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = acosh(*p1); }
+      { (*p2) = cacoshf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = acosh(*p1); }
+      { (*p2) = cacoshf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -14875,7 +14876,7 @@ ca_monop_acosh_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = acosh(*p1); }
+        { (*p2) = cacosh(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -14890,14 +14891,14 @@ ca_monop_acosh_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = acosh(*p1); }
+      { (*p2) = cacosh(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = acosh(*p1); }
+      { (*p2) = cacosh(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15026,7 +15027,7 @@ ca_monop_atanh_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = atanh(*p1); }
+        { (*p2) = atanhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15041,14 +15042,14 @@ ca_monop_atanh_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = atanh(*p1); }
+      { (*p2) = atanhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = atanh(*p1); }
+      { (*p2) = atanhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15116,7 +15117,7 @@ ca_monop_atanh_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = atanh(*p1); }
+        { (*p2) = catanhf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15131,14 +15132,14 @@ ca_monop_atanh_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = atanh(*p1); }
+      { (*p2) = catanhf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = atanh(*p1); }
+      { (*p2) = catanhf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15161,7 +15162,7 @@ ca_monop_atanh_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = atanh(*p1); }
+        { (*p2) = catanh(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15176,14 +15177,14 @@ ca_monop_atanh_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = atanh(*p1); }
+      { (*p2) = catanh(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = atanh(*p1); }
+      { (*p2) = catanh(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15312,7 +15313,7 @@ ca_monop_expm1_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = expm1(*p1); }
+        { (*p2) = expm1f(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15327,14 +15328,14 @@ ca_monop_expm1_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = expm1(*p1); }
+      { (*p2) = expm1f(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = expm1(*p1); }
+      { (*p2) = expm1f(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15508,7 +15509,7 @@ ca_monop_log1p_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = log1p(*p1); }
+        { (*p2) = log1pf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15523,14 +15524,14 @@ ca_monop_log1p_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = log1p(*p1); }
+      { (*p2) = log1pf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = log1p(*p1); }
+      { (*p2) = log1pf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15704,7 +15705,7 @@ ca_monop_rsqrt_float32_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = 1.0 / sqrt(*p1); }
+        { (*p2) = ((float)1.0) / sqrtf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15719,14 +15720,14 @@ ca_monop_rsqrt_float32_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = 1.0 / sqrt(*p1); }
+      { (*p2) = ((float)1.0) / sqrtf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     float  *p1 = q1;
     float *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = 1.0 / sqrt(*p1); }
+      { (*p2) = ((float)1.0) / sqrtf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15749,7 +15750,7 @@ ca_monop_rsqrt_float64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = 1.0 / sqrt(*p1); }
+        { (*p2) = ((double)1.0) / sqrt(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15764,14 +15765,14 @@ ca_monop_rsqrt_float64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = 1.0 / sqrt(*p1); }
+      { (*p2) = ((double)1.0) / sqrt(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     double  *p1 = q1;
     double *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = 1.0 / sqrt(*p1); }
+      { (*p2) = ((double)1.0) / sqrt(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15794,7 +15795,7 @@ ca_monop_rsqrt_cmplx64_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = 1.0 / csqrt(*p1); }
+        { (*p2) = ((float)1.0) / csqrtf(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15809,14 +15810,14 @@ ca_monop_rsqrt_cmplx64_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = 1.0 / csqrt(*p1); }
+      { (*p2) = ((float)1.0) / csqrtf(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx64_t  *p1 = q1;
     cmplx64_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = 1.0 / csqrt(*p1); }
+      { (*p2) = ((float)1.0) / csqrtf(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -15839,7 +15840,7 @@ ca_monop_rsqrt_cmplx128_t (ca_size_t n, boolean8_t *m,
       if ( ! *pm ) {
         p1 = q1 + k*i1;
         p2 = q2 + k*i2;
-        { (*p2) = 1.0 / csqrt(*p1); }
+        { (*p2) = ((double)1.0) / csqrt(*p1); }
       }
     }
   } else if ( i1 == 1 && i2 == 1 ) {
@@ -15854,14 +15855,14 @@ ca_monop_rsqrt_cmplx128_t (ca_size_t n, boolean8_t *m,
     #pragma omp simd
     for (k=0; k<n; k++) {
       p1 = q1 + k; p2 = q2 + k;
-      { (*p2) = 1.0 / csqrt(*p1); }
+      { (*p2) = ((double)1.0) / csqrt(*p1); }
     }
   } else {
     /* strided branch: numo-shape `p += s` increment form. */
     cmplx128_t  *p1 = q1;
     cmplx128_t *p2 = q2;
     for (k=0; k<n; k++) {
-      { (*p2) = 1.0 / csqrt(*p1); }
+      { (*p2) = ((double)1.0) / csqrt(*p1); }
       p1 += i1; p2 += i2;
     }
   }
@@ -18228,9 +18229,6 @@ Init_carray_kernels_monop (void)
   rb_define_method(rb_cCArray, "exp", rb_ca_exp, 0);
   rb_define_method(rb_cCArray, "exp!", rb_ca_exp_bang, 0);
   rb_define_module_function(rb_mCAMath, "exp", rb_cmath_exp, 1);
-  rb_define_method(rb_cCArray, "log", rb_ca_log, 0);
-  rb_define_method(rb_cCArray, "log!", rb_ca_log_bang, 0);
-  rb_define_module_function(rb_mCAMath, "log", rb_cmath_log, 1);
   rb_define_method(rb_cCArray, "sin", rb_ca_sin, 0);
   rb_define_method(rb_cCArray, "sin!", rb_ca_sin_bang, 0);
   rb_define_module_function(rb_mCAMath, "sin", rb_cmath_sin, 1);
@@ -18249,6 +18247,9 @@ Init_carray_kernels_monop (void)
   rb_define_method(rb_cCArray, "atan", rb_ca_atan, 0);
   rb_define_method(rb_cCArray, "atan!", rb_ca_atan_bang, 0);
   rb_define_module_function(rb_mCAMath, "atan", rb_cmath_atan, 1);
+  rb_define_method(rb_cCArray, "log", rb_ca_log, 0);
+  rb_define_method(rb_cCArray, "log!", rb_ca_log_bang, 0);
+  rb_define_module_function(rb_mCAMath, "log", rb_cmath_log, 1);
   rb_define_method(rb_cCArray, "exp2", rb_ca_exp2, 0);
   rb_define_method(rb_cCArray, "exp2!", rb_ca_exp2_bang, 0);
   rb_define_module_function(rb_mCAMath, "exp2", rb_cmath_exp2, 1);