Feature #10474 ยป combine_some_macros_into_one_at_math.c.patch
| math.c | ||
|---|---|---|
|
VALUE rb_mMath;
|
||
|
VALUE rb_eMathDomainError;
|
||
|
#define Need_Float(x) do {if (!RB_TYPE_P(x, T_FLOAT)) {(x) = rb_to_float(x);}} while(0)
|
||
|
#define Need_Float2(x,y) do {\
|
||
|
Need_Float(x);\
|
||
|
Need_Float(y);\
|
||
|
} while (0)
|
||
|
#define Get_RFloat_Value(x) RFLOAT_VALUE(RB_TYPE_P(x, T_FLOAT) ? x : rb_to_float(x))
|
||
|
#define domain_error(msg) \
|
||
|
rb_raise(rb_eMathDomainError, "Numerical argument is out of domain - " #msg)
|
||
| ... | ... | |
|
# define M_PI 3.14159265358979323846
|
||
|
#endif
|
||
|
double dx, dy;
|
||
|
Need_Float2(y, x);
|
||
|
dx = RFLOAT_VALUE(x);
|
||
|
dy = RFLOAT_VALUE(y);
|
||
|
dx = Get_RFloat_Value(x);
|
||
|
dy = Get_RFloat_Value(y);
|
||
|
if (dx == 0.0 && dy == 0.0) {
|
||
|
if (!signbit(dx))
|
||
|
return DBL2NUM(dy);
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_cos(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(cos(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(cos(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_sin(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(sin(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(sin(Get_RFloat_Value(x)));
|
||
|
}
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_tan(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(tan(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(tan(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
{
|
||
|
double d0, d;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < -1.0 || 1.0 < d0) domain_error("acos");
|
||
|
d = acos(d0);
|
||
| ... | ... | |
|
{
|
||
|
double d0, d;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < -1.0 || 1.0 < d0) domain_error("asin");
|
||
|
d = asin(d0);
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_atan(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(atan(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(atan(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
#ifndef HAVE_COSH
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_cosh(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(cosh(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(cosh(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
#ifndef HAVE_SINH
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_sinh(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(sinh(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(sinh(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
#ifndef HAVE_TANH
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_tanh(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(tanh(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(tanh(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
{
|
||
|
double d0, d;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < 1.0) domain_error("acosh");
|
||
|
d = acosh(d0);
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_asinh(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(asinh(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(asinh(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
{
|
||
|
double d0, d;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < -1.0 || +1.0 < d0) domain_error("atanh");
|
||
|
/* check for pole error */
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_exp(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(exp(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(exp(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
#if defined __CYGWIN__
|
||
| ... | ... | |
|
numbits = 0;
|
||
|
}
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < 0.0) domain_error("log");
|
||
|
/* check for pole error */
|
||
| ... | ... | |
|
numbits = 0;
|
||
|
}
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < 0.0) domain_error("log2");
|
||
|
/* check for pole error */
|
||
| ... | ... | |
|
numbits = 0;
|
||
|
}
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < 0.0) domain_error("log10");
|
||
|
/* check for pole error */
|
||
| ... | ... | |
|
{
|
||
|
double d0, d;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (d0 < 0.0) domain_error("sqrt");
|
||
|
if (d0 == 0.0) return DBL2NUM(0.0);
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_cbrt(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(cbrt(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(cbrt(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
double d;
|
||
|
int exp;
|
||
|
Need_Float(x);
|
||
|
d = frexp(RFLOAT_VALUE(x), &exp);
|
||
|
d = frexp(Get_RFloat_Value(x), &exp);
|
||
|
return rb_assoc_new(DBL2NUM(d), INT2NUM(exp));
|
||
|
}
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_ldexp(VALUE obj, VALUE x, VALUE n)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(ldexp(RFLOAT_VALUE(x), NUM2INT(n)));
|
||
|
return DBL2NUM(ldexp(Get_RFloat_Value(x), NUM2INT(n)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_hypot(VALUE obj, VALUE x, VALUE y)
|
||
|
{
|
||
|
Need_Float2(x, y);
|
||
|
return DBL2NUM(hypot(RFLOAT_VALUE(x), RFLOAT_VALUE(y)));
|
||
|
return DBL2NUM(hypot(Get_RFloat_Value(x), Get_RFloat_Value(y)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_erf(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(erf(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(erf(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
static VALUE
|
||
|
math_erfc(VALUE obj, VALUE x)
|
||
|
{
|
||
|
Need_Float(x);
|
||
|
return DBL2NUM(erfc(RFLOAT_VALUE(x)));
|
||
|
return DBL2NUM(erfc(Get_RFloat_Value(x)));
|
||
|
}
|
||
|
/*
|
||
| ... | ... | |
|
};
|
||
|
double d0, d;
|
||
|
double intpart, fracpart;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (isinf(d0) && signbit(d0)) domain_error("gamma");
|
||
|
fracpart = modf(d0, &intpart);
|
||
| ... | ... | |
|
double d0, d;
|
||
|
int sign=1;
|
||
|
VALUE v;
|
||
|
Need_Float(x);
|
||
|
d0 = RFLOAT_VALUE(x);
|
||
|
d0 = Get_RFloat_Value(x);
|
||
|
/* check for domain error */
|
||
|
if (isinf(d0)) {
|
||
|
if (signbit(d0)) domain_error("lgamma");
|
||