Sha256: 7ecda3560c35d900e738cb8421335129d5f0cdc7c0d7e83c8e0c4d12cbaa20d2
Contents?: true
Size: 1.16 KB
Versions: 19
Compression:
Stored size: 1.16 KB
Contents
#include "xint_macro.h"
#define m_sign(x) (((x)==0) ? 0 : (((x)>0) ? 1 : -1))
static inline dtype m_abs(dtype x) {
if (x==DATA_MIN) {
rb_raise(cumo_na_eValueError, "cannot convert the minimum integer");
}
return (x<0)?-x:x;
}
static inline dtype int_reciprocal(dtype x) {
switch (x) {
case 1:
return 1;
case -1:
return -1;
case 0:
rb_raise(rb_eZeroDivError, "divided by 0");
default:
return 0;
}
}
/*
static dtype pow_int(dtype x, int p)
{
dtype r = m_one;
switch(p) {
case 0: return 1;
case 1: return x;
case 2: return x*x;
case 3: return x*x*x;
}
if (p<0) return 0;
while (p) {
if (p&1) r *= x;
x *= x;
p >>= 1;
}
return r;
}
*/
static inline int64_t f_sum(size_t n, char *p, ssize_t stride)
{
int64_t x,y=0;
size_t i=n;
for (; i--;) {
x = *(dtype*)p;
y += x;
p += stride;
}
return y;
}
static inline int64_t f_prod(size_t n, char *p, ssize_t stride)
{
int64_t x,y=1;
size_t i=n;
for (; i--;) {
x = *(dtype*)p;
y *= x;
p += stride;
}
return y;
}
Version data entries
19 entries across 19 versions & 1 rubygems