ccmath 0.2.7 → 0.2.71

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 83d2fef732aa63e2e110ef2dacbfdad43e9367f0
-  data.tar.gz: f3b1d96ddfa548493c1f6172e26334d43af65ee6
+  metadata.gz: c2fda0ca3c87789fe9fbba2ab15250a2e4cf7c6d
+  data.tar.gz: 73c23f5a3218c774c21035689ca9ded2d6d30b14
 SHA512:
-  metadata.gz: 5d97f22e2547750e60adf92c264683fa0733f67ed207bab79b7494b22a94a83fdf99ed89d307eadd50632e8f72a7e22ca346b8392f3488e1d0f628f3a06f52cf
-  data.tar.gz: 6babd532fb5e80db788d875ca051c08ebb33ca2247982594d3451db7a64ff3f0c05cca0d4b3e36723a6e7d007d1dd09ed379d1f075b87cd095b6d3772ebc17db
+  metadata.gz: eae6b0686f5e9bcc96a12d557ae39830badc4c3863c9ac80baf8d1a4385bd3ea8cb2819ccca506d60ada66b5520cbaac87ab9a2eadb7ddcbf583627070fd5576
+  data.tar.gz: 2792755a0be3d4849243eca410c013eb33b302d99ce0022191abc10f235b8bd002f1fb7ff49264e4b5f497a755dcc03637e98923b6616f74a4c1ee8f21a29914

data/CHANGELOG.md

@@ -0,0 +1,9 @@
+## CHANGELOG
+
+### v2.71
+
+* `CCMath.gamma` can receive `Complex` object.
+
+### v2.7
+
+* Initial release.

data/README.md

@@ -1,28 +1,19 @@
-# CMath implemented with C
+# Hopefully, CCMath is better than CMath... (・ω・;)
 
 [![Build Status](https://travis-ci.org/gogotanaka/ccmath.svg?branch=master)](https://travis-ci.org/gogotanaka/ccmath)
 
-Ruby has two official libraries related with mathematics.
-
-* **Math** ...Implemented with C. Only for Real number.
-* **CMath**...Implemented with Ruby. Accept Complex number.
-
-And now,
-
-* **CCMath**...Implemented with C. Accept Complex number.
-
-You can declare codomain.
+CCMath is CMath implemented with C-lang.
 
+## Demo
 ```ruby
-CCMath.set = 'R' # equivalent to Math
-CCMath.set = 'C' # equivalent to CMath
-```
+CCMath.cos(1i)
+#=> (1.5430806348152437+0.0i)
 
-## Demo
-```
 CCMath.cos(1)
-#=> 0.5403023058681398
+#=> (0.5403023058681398+0.0i)
+```
 
+## Bench
 
 ```
 Calculating -------------------------------------
@@ -93,3 +84,17 @@ AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
+
+# Contributor Code of Conduct
+
+As contributors and maintainers of this project, we pledge to respect all people who contribute through reporting issues, posting feature requests, updating documentation, submitting pull requests or patches, and other activities.
+
+We are committed to making participation in this project a harassment-free experience for everyone, regardless of level of experience, gender, gender identity and expression, sexual orientation, disability, personal appearance, body size, race, ethnicity, age, or religion.
+
+Examples of unacceptable behavior by participants include the use of sexual language or imagery, derogatory comments or personal attacks, trolling, public or private harassment, insults, or other unprofessional conduct.
+
+Project maintainers have the right and responsibility to remove, edit, or reject comments, commits, code, wiki edits, issues, and other contributions that are not aligned to this Code of Conduct. Project maintainers who do not follow the Code of Conduct may be removed from the project team.
+
+Instances of abusive, harassing, or otherwise unacceptable behavior may be reported by opening an issue or contacting one or more of the project maintainers.
+
+This Code of Conduct is adapted from the [Contributor Covenant](http://contributor-covenant.org), version 1.0.0, available at [http://contributor-covenant.org/version/1/0/0/](http://contributor-covenant.org/version/1/0/0/)

data/ext/ccmath/ccmath.c

@@ -3,7 +3,7 @@
 VALUE rb_mCCMath;
 VALUE rb_eMathDomainError;
 
-static ID id_real_p, id_codomain;
+static ID id_real_p, id_set, id_power;
 
 #define RB_BIGNUM_TYPE_P(x) RB_TYPE_P((x), T_BIGNUM)
 #define BIGNUM_POSITIVE_P(b) (FIX2LONG(rb_big_cmp((b), INT2FIX(0))) >= 0)
@@ -28,12 +28,9 @@ num2dbl_without_to_f(VALUE num)
     }
     else {
         switch (BUILTIN_TYPE(num)) {
-          case T_FLOAT:
-              return RFLOAT_VALUE(num);
-          case T_BIGNUM:
-              return big2dbl_without_to_f(num);
-          case T_RATIONAL:
-              return rat2dbl_without_to_f(num);
+          case T_FLOAT:    return RFLOAT_VALUE(num);
+          case T_BIGNUM:   return big2dbl_without_to_f(num);
+          case T_RATIONAL: return rat2dbl_without_to_f(num);
         }
     }
     return RFLOAT_VALUE(rb_to_float(num));
@@ -47,6 +44,19 @@ f_real_p(VALUE x)
     return rb_funcall(x, id_real_p, 0);
 }
 
+#define binop(n,op) \
+inline static VALUE \
+f_##n(VALUE x, VALUE y)\
+{\
+    return rb_funcall(x, (op), 1, y);\
+}
+
+binop(add, '+')
+binop(div, '/')
+binop(sub, '-')
+binop(mul, '*')
+binop(pow, id_power)
+
 inline static VALUE
 DBLS2COMP(double real, double imag)
 {
@@ -71,73 +81,43 @@ DBLS2COMP(double real, double imag)
 
 #define EXTRACT_DBLS(z);              \
     double z##_real, z##_imag;        \
-    do{                               \
+    do {                              \
         if (f_real_p(z)) {            \
             EXTRACT_DBLS_FROM_REAL(z) \
         } else {                      \
             EXTRACT_DBLS_FROM_COMP(z) \
         }                             \
-    }while(0)
+    } while(0)
 
 #define domain_error(msg) \
     rb_raise(rb_eMathDomainError, "Numerical argument is out of domain - " #msg)
 
-static inline int
-codomain_check(VALUE z) {
-    if (!rb_respond_to(z, id_real_p)) rb_raise(rb_eTypeError, "Numeric Number required");
-
-    char codomain = *(RSTRING_PTR(rb_ivar_get(rb_mCCMath, id_codomain)));
-    if (codomain == 'R') {
-        if (!f_real_p(z)) rb_raise(rb_eMathDomainError, "Real Number required");
-        return 1;
-    }
-    else if (codomain == 'C') {
-        return 0;
-    }
-    else {
-        domain_error("wow");
-    }
-}
-
 static VALUE
 ccmath_sqrt(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d = NUM2DBL_F(z);
-        if (d < 0.0) domain_error("sqrt");
-        if (d == 0.0) return DBL2NUM(0.0);
-        return DBL2NUM(sqrt(d));
-    }
-    else {
-        EXTRACT_DBLS(z);
-        if (z_imag == 0.0) {
-            if (z_real < 0.0) {
-                return DBLS2COMP(0.0, sqrt(fabs(z_real)));
-            }
-            else {
-                return DBLS2COMP(sqrt(z_real), 0.0);
-            }
+    EXTRACT_DBLS(z);
+    if (z_imag == 0.0) {
+        if (z_real < 0.0) {
+            return DBLS2COMP(0.0, sqrt(fabs(z_real)));
         }
         else {
-            double s = sqrt((hypot(z_real, z_imag) + z_real) / 2.0);
-            return DBLS2COMP(s, z_imag / (2 * s));
+            return DBLS2COMP(sqrt(z_real), 0.0);
         }
     }
+    else {
+        double s = sqrt((hypot(z_real, z_imag) + z_real) / 2.0);
+        return DBLS2COMP(s, z_imag / (2 * s));
+    }
 }
 
 static VALUE
 ccmath_exp(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(exp(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);
-        if (z_imag == 0.0) return DBLS2COMP(exp(z_real), 0.0);
+    EXTRACT_DBLS(z);
+    if (z_imag == 0.0) return DBLS2COMP(exp(z_real), 0.0);
 
-        double ere = exp(z_real);
-        return DBLS2COMP(ere * cos(z_imag), ere * sin(z_imag));
-    }
+    double ere = exp(z_real);
+    return DBLS2COMP(ere * cos(z_imag), ere * sin(z_imag));
 }
 
 void
@@ -173,291 +153,229 @@ ccmath_log(int argc, const VALUE* argv, VALUE obj)
     VALUE z, base;
     rb_scan_args(argc, argv, "11", &z, &base);
 
-    if (codomain_check(z)) {
-        double d = internal_log(z);
-        if (argc == 2) {
-            d /= internal_log(base);
-        }
-        return DBL2NUM(d);
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        float r = hypot(z_real, z_imag);
-        if (argc == 2) {
-            if (!rb_respond_to(base, id_real_p)) rb_raise(rb_eTypeError, "Numeric Number required");
-            EXTRACT_DBLS(base);
-            if (base_imag != 0.0) domain_error("log");
-            if (base_real > 0.0) {
-                double ln_base = log(base_real);
-                return DBLS2COMP(log(r) / ln_base, m_atan2(z_imag, z_real) / ln_base);
-            }
-            else {
-                VALUE ln_base = DBLS2COMP(log(fabs(base_real)), M_PI);
-                return rb_funcall(DBLS2COMP(log(r), m_atan2(z_imag, z_real)), '/', 1, ln_base);
-            }
+    EXTRACT_DBLS(z);
+    float r = hypot(z_real, z_imag);
+    if (argc == 2) {
+        if (!rb_respond_to(base, id_real_p)) rb_raise(rb_eTypeError, "Numeric Number required");
+        EXTRACT_DBLS(base);
+        if (base_imag != 0.0) domain_error("log");
+        if (base_real > 0.0) {
+            double ln_base = log(base_real);
+            return DBLS2COMP(log(r) / ln_base, m_atan2(z_imag, z_real) / ln_base);
         }
         else {
-            return DBLS2COMP(log(r), m_atan2(z_imag, z_real));
+            VALUE ln_base = DBLS2COMP(log(fabs(base_real)), M_PI);
+            return f_div(DBLS2COMP(log(r), m_atan2(z_imag, z_real)), ln_base);
         }
     }
+    else {
+        return DBLS2COMP(log(r), m_atan2(z_imag, z_real));
+    }
 }
 
 static VALUE
 ccmath_log2(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d;
-        size_t numbits;
-
-        assing_numbits(z, &d, &numbits);
-
-        if (d < 0.0) domain_error("log2");
-        if (d == 0.0) return DBL2NUM(-INFINITY);
-
-        return DBL2NUM(log2(d) + numbits);
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        float r = hypot(z_real, z_imag);
-        return DBLS2COMP(log(r) / M_LN2, m_atan2(z_imag, z_real) / M_LN2);
-    }
+    EXTRACT_DBLS(z);
+    float r = hypot(z_real, z_imag);
+    return DBLS2COMP(log(r) / M_LN2, m_atan2(z_imag, z_real) / M_LN2);
 }
 
 static VALUE
 ccmath_log10(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d;
-        size_t numbits;
-
-        assing_numbits(z, &d, &numbits);
-
-        if (d < 0.0) domain_error("log10");
-        if (d == 0.0) return DBL2NUM(-INFINITY);
-
-        return DBL2NUM(log10(d) + numbits * log10(2));
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        float r = hypot(z_real, z_imag);
-        return DBLS2COMP(log(r) / M_LN10, m_atan2(z_imag, z_real) / M_LN10);
-    }
+    EXTRACT_DBLS(z);
+    float r = hypot(z_real, z_imag);
+    return DBLS2COMP(log(r) / M_LN10, m_atan2(z_imag, z_real) / M_LN10);
 }
 
 static VALUE
 ccmath_cos(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(cos(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        if (z_real == 0.0) return DBLS2COMP(cosh(z_imag), 0.0);
-        if (z_imag == 0.0) return DBLS2COMP(cos(z_real), 0.0);
-        return DBLS2COMP(cos(z_real) * cosh(z_imag), -sin(z_real) * sinh(z_imag));
-    }
+    EXTRACT_DBLS(z);
+    if (z_real == 0.0) return DBLS2COMP(cosh(z_imag), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(cos(z_real), 0.0);
+    return DBLS2COMP(cos(z_real) * cosh(z_imag), -sin(z_real) * sinh(z_imag));
 }
 
 static VALUE
 ccmath_sin(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(sin(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        if (z_real == 0.0) return DBLS2COMP(0.0, sinh(z_imag));
-        if (z_imag == 0.0) return DBLS2COMP(sin(z_real), 0.0);
-        return DBLS2COMP(sin(z_real) * cosh(z_imag), cos(z_real) * sinh(z_imag));
-    }
-
+    EXTRACT_DBLS(z);
+    if (z_real == 0.0) return DBLS2COMP(0.0, sinh(z_imag));
+    if (z_imag == 0.0) return DBLS2COMP(sin(z_real), 0.0);
+    return DBLS2COMP(sin(z_real) * cosh(z_imag), cos(z_real) * sinh(z_imag));
 }
 
 static VALUE
 ccmath_tan(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(tan(NUM2DBL_F(z)));
-    }
-    else {
-        return rb_funcall(ccmath_sin(obj, z), '/', 1, ccmath_cos(obj, z));
-    }
+    return f_div(ccmath_sin(obj, z), ccmath_cos(obj, z));
 }
 
 static VALUE
 ccmath_cosh(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(cosh(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        if (z_real == 0.0) return DBLS2COMP(cos(z_imag), 0.0);
-        if (z_imag == 0.0) return DBLS2COMP(cosh(z_real), 0.0);
-        return DBLS2COMP(cosh(z_real) * cos(z_imag), sinh(z_real) * sin(z_imag));
-    }
+    EXTRACT_DBLS(z);
+    if (z_real == 0.0) return DBLS2COMP(cos(z_imag), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(cosh(z_real), 0.0);
+    return DBLS2COMP(cosh(z_real) * cos(z_imag), sinh(z_real) * sin(z_imag));
 }
 
 static VALUE
 ccmath_sinh(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(sinh(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);;
-        if (z_real == 0.0) return DBLS2COMP(0.0, sin(z_imag));
-        if (z_imag == 0.0) return DBLS2COMP(sinh(z_real), 0.0);
-        return DBLS2COMP(sinh(z_real) * cos(z_imag), cosh(z_real) * sin(z_imag));
-    }
+    EXTRACT_DBLS(z);
+    if (z_real == 0.0) return DBLS2COMP(0.0, sin(z_imag));
+    if (z_imag == 0.0) return DBLS2COMP(sinh(z_real), 0.0);
+    return DBLS2COMP(sinh(z_real) * cos(z_imag), cosh(z_real) * sin(z_imag));
 }
 
 static VALUE
 ccmath_tanh(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(tanh(NUM2DBL_F(z)));
-    }
-    else {
-        return rb_funcall(ccmath_sinh(obj, z), '/', 1, ccmath_cosh(obj, z));
-    }
+    return f_div(ccmath_sinh(obj, z), ccmath_cosh(obj, z));
 }
 
 static VALUE
 ccmath_asinh(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(asinh(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);;
+    EXTRACT_DBLS(z);
 
-        if (z_imag == 0.0) return DBLS2COMP(asinh(z_real), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(asinh(z_real), 0.0);
 
-        VALUE s1 = DBLS2COMP(1.0 + z_imag, -z_real);
-        VALUE s2 = DBLS2COMP(1.0 - z_imag, z_real);
+    VALUE s1 = DBLS2COMP(1.0 + z_imag, -z_real);
+    VALUE s2 = DBLS2COMP(1.0 - z_imag, z_real);
 
-        s1 = ccmath_sqrt(obj, s1);
-        s2 = ccmath_sqrt(obj, s2);
-        EXTRACT_DBLS(s1);
-        EXTRACT_DBLS(s2);
+    s1 = ccmath_sqrt(obj, s1);
+    s2 = ccmath_sqrt(obj, s2);
+    EXTRACT_DBLS(s1);
+    EXTRACT_DBLS(s2);
 
-        return DBLS2COMP(asinh(s1_real * s2_imag - s2_real * s1_imag), atan2(z_imag, s1_real * s2_real - s1_imag * s2_imag));
-    }
+    return DBLS2COMP(asinh(s1_real * s2_imag - s2_real * s1_imag), atan2(z_imag, s1_real * s2_real - s1_imag * s2_imag));
 }
 
 static VALUE
 ccmath_asin(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d = NUM2DBL_F(z);
-        if (d < -1.0 || 1.0 < d) domain_error("asin");
-        return DBL2NUM(asin(d));
-    }
-    else {
-        EXTRACT_DBLS(z);;
+    EXTRACT_DBLS(z);
 
-        if (z_imag == 0.0) return DBLS2COMP(asin(z_real), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(asin(z_real), 0.0);
 
-        VALUE s = ccmath_asinh(obj, DBLS2COMP(-z_imag, z_real));
-        EXTRACT_DBLS(s);
-        return DBLS2COMP(s_imag, -s_real);
-    }
+    VALUE s = ccmath_asinh(obj, DBLS2COMP(-z_imag, z_real));
+    EXTRACT_DBLS(s);
+    return DBLS2COMP(s_imag, -s_real);
 }
 
 static VALUE
 ccmath_acosh(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d = NUM2DBL_F(z);
-        if (d < 1.0) domain_error("acosh");
-        return DBL2NUM(acosh(d));
-    }
-    else {
-        EXTRACT_DBLS(z);;
+    EXTRACT_DBLS(z);
 
-        if (z_imag == 0.0) return DBLS2COMP(acosh(z_real), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(acosh(z_real), 0.0);
 
-        VALUE s1 = DBLS2COMP(z_real - 1.0, z_imag);
-        VALUE s2 = DBLS2COMP(z_real + 1.0, z_imag);
+    VALUE s1 = DBLS2COMP(z_real - 1.0, z_imag);
+    VALUE s2 = DBLS2COMP(z_real + 1.0, z_imag);
 
-        s1 = ccmath_sqrt(obj, s1);
-        s2 = ccmath_sqrt(obj, s2);
-        EXTRACT_DBLS(s1);
-        EXTRACT_DBLS(s2);
+    s1 = ccmath_sqrt(obj, s1);
+    s2 = ccmath_sqrt(obj, s2);
+    EXTRACT_DBLS(s1);
+    EXTRACT_DBLS(s2);
 
-        return DBLS2COMP(asinh(s1_real * s2_real + s1_imag * s2_imag), 2.0 * atan2(s1_imag, s2_real));
-    }
+    return DBLS2COMP(asinh(s1_real * s2_real + s1_imag * s2_imag), 2.0 * atan2(s1_imag, s2_real));
 }
 
 static VALUE
 ccmath_acos(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d = NUM2DBL_F(z);
-        if (d < -1.0 || 1.0 < d) domain_error("acos");
-        return DBL2NUM(acos(d));
-    }
-    else {
-        EXTRACT_DBLS(z);;
+    EXTRACT_DBLS(z);
 
-        if (z_imag == 0.0) return DBLS2COMP(acos(z_real), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(acos(z_real), 0.0);
 
-        return rb_funcall(DBL2NUM(M_PI / 2.0), '-', 1, ccmath_asin(obj, z));
-    }
+    return f_sub(DBL2NUM(M_PI / 2.0), ccmath_asin(obj, z));
 }
 
 static VALUE
 ccmath_atanh(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        double d = NUM2DBL_F(z);
+    EXTRACT_DBLS(z);
 
-        if (d <  -1.0 || +1.0 <  d) domain_error("atanh");
-        if (d == -1.0) return DBL2NUM(-INFINITY);
-        if (d == +1.0) return DBL2NUM(+INFINITY);
-        return DBL2NUM(atanh(d));
-    }
-    else {
-        EXTRACT_DBLS(z);;
+    if (z_imag == 0.0) return DBLS2COMP(atanh(z_real), 0.0);
 
-        if (z_imag == 0.0) return DBLS2COMP(atanh(z_real), 0.0);
+    double sq_imag = z_imag * z_imag;
 
-        double sq_imag = z_imag * z_imag;
-
-        return DBLS2COMP(m_log1p(4.0 * z_real / ((1.0 - z_real) * (1.0 - z_real) + sq_imag)) / 4.0,
-                        -m_atan2(-2.0 * z_imag, (1.0 - z_real) * (1.0 + z_real) - sq_imag) / 2.0);
-    }
+    return DBLS2COMP(m_log1p(4.0 * z_real / ((1.0 - z_real) * (1.0 - z_real) + sq_imag)) / 4.0,
+                    -m_atan2(-2.0 * z_imag, (1.0 - z_real) * (1.0 + z_real) - sq_imag) / 2.0);
 }
 
 static VALUE
 ccmath_atan(VALUE obj, VALUE z)
 {
-    if (codomain_check(z)) {
-        return DBL2NUM(atan(NUM2DBL_F(z)));
-    }
-    else {
-        EXTRACT_DBLS(z);;
+    EXTRACT_DBLS(z);
 
-        if (z_imag == 0.0) return DBLS2COMP(atan(z_real), 0.0);
+    if (z_imag == 0.0) return DBLS2COMP(atan(z_real), 0.0);
 
-        VALUE s = ccmath_atanh(obj, DBLS2COMP(-z_imag, z_real));
-        EXTRACT_DBLS(s);
-        return DBLS2COMP(s_imag, -s_real);
-    }
+    VALUE s = ccmath_atanh(obj, DBLS2COMP(-z_imag, z_real));
+    EXTRACT_DBLS(s);
+    return DBLS2COMP(s_imag, -s_real);
 }
 
 static VALUE
-ccmath_define_set(VALUE obj, VALUE str)
+ccmath_gamma(VALUE obj, VALUE z)
 {
-    rb_ivar_set(obj, id_codomain, str);
-    return obj;
+    EXTRACT_DBLS(z);
+    if (z_real < 0.5) {
+        VALUE s1, s2;
+        s1 = ccmath_gamma(obj, f_sub(DBL2NUM(1), z));
+        s2 = ccmath_sin(obj, f_mul(DBL2NUM(M_PI), z));
+        return f_div(DBL2NUM(M_PI), f_mul(s1, s2));
+    }
+    else {
+        static const double lanczos_coef[] = {
+            0.99999999999980993,
+            676.5203681218851,
+            -1259.1392167224028,
+            771.32342877765313,
+            -176.61502916214059,
+            12.507343278686905,
+            -0.13857109526572012,
+            9.9843695780195716e-6,
+            1.5056327351493116e-7
+        };
+        int i;
+        double g, s, x_real, x_imag;
+
+        g = 7.0;
+        z_real -= 1;
+        x_real = lanczos_coef[0];
+        x_imag = 0.0;
+        for(i=1; i<g+2; i++) {
+            s = lanczos_coef[i] / ((z_real+i) * (z_real+i) + (z_imag * z_imag));
+            x_real += s * (z_real+i);
+            x_imag -= s * z_imag;
+        }
+        VALUE x = DBLS2COMP(x_real, x_imag);
+        VALUE sqrt_2_pi = DBL2NUM(sqrt(2 * M_PI));
+        VALUE t = DBLS2COMP(z_real + g + 0.5, z_imag);
+        return f_mul(sqrt_2_pi,
+               f_mul(f_pow(t, DBLS2COMP(z_real+0.5, z_imag)),
+               f_mul(f_div(DBL2NUM(1), ccmath_exp(obj,t)), x)));
+
+    }
 }
 
+// static VALUE
+// ccmath_define_set(VALUE obj, VALUE str)
+// {
+//     rb_ivar_set(obj, id_set, str);
+//     return obj;
+// }
+
 void Init_ccmath(void)
 {
     id_real_p = rb_intern("real?");
-    id_codomain = rb_intern("codomain");
+    id_power = rb_intern("**");
+    id_set = rb_intern("set");
     rb_mCCMath = rb_define_module("CCMath");
 
     rb_eMathDomainError = rb_define_class_under(rb_mCCMath, "DomainError", rb_eStandardError);
@@ -490,7 +408,8 @@ void Init_ccmath(void)
     rb_define_module_function(rb_mCCMath, "acos", ccmath_acos, 1);
     rb_define_module_function(rb_mCCMath, "atanh", ccmath_atanh, 1);
     rb_define_module_function(rb_mCCMath, "atan", ccmath_atan, 1);
-    rb_define_module_function(rb_mCCMath, "set=", ccmath_define_set, 1);
+    rb_define_module_function(rb_mCCMath, "gamma", ccmath_gamma, 1);
+    // rb_define_module_function(rb_mCCMath, "set=", ccmath_define_set, 1);
 
-    rb_ivar_set(rb_mCCMath, id_codomain, rb_str_new2("R"));
+    // rb_ivar_set(rb_mCCMath, id_set, rb_str_new2("R"));
 }

data/lib/ccmath.rb

@@ -2,21 +2,37 @@ require "ccmath/version"
 require "ccmath/ccmath"
 
 module CCMath
-  class << self
-    %w[
-      atan2
-      cbrt
-      frexp
-      ldexp
-      hypot
-      erf
-      erfc
-      gamma
-      lgamma
-    ].each do |meth|
-      define_method(meth) do |*args, &blk|
-        Math.send(meth, *args, &blk)
-      end
+  def cbrt(z)
+    z ** (1.0/3)
+  end
+
+  def atan2(y,x)
+    if y.real? and x.real?
+      atan2!(y,x)
+    else
+      (-1.0).i * log((x + 1.0.i * y) / sqrt(x * x + y * y))
+    end
+  end
+
+  %w[
+    frexp
+    ldexp
+    hypot
+    erf
+    erfc
+    lgamma
+  ].each do |meth|
+    define_method(meth) do |*args, &blk|
+      Math.send(meth, *args, &blk)
     end
   end
+
+  module_function :cbrt
+  module_function :atan2
+  module_function :frexp
+  module_function :ldexp
+  module_function :hypot
+  module_function :erf
+  module_function :erfc
+  module_function :lgamma
 end

data/lib/ccmath/version.rb

@@ -1,3 +1,3 @@
 module Ccmath
-  VERSION = "0.2.7"
+  VERSION = "0.2.71"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: ccmath
 version: !ruby/object:Gem::Version
-  version: 0.2.7
+  version: 0.2.71
 platform: ruby
 authors:
 - gogotanaka
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2015-07-05 00:00:00.000000000 Z
+date: 2015-08-07 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -76,7 +76,7 @@ extra_rdoc_files: []
 files:
 - ".gitignore"
 - ".travis.yml"
-- CODE_OF_CONDUCT.md
+- CHANGELOG.md
 - Gemfile
 - README.md
 - Rakefile

data/CODE_OF_CONDUCT.md

@@ -1,13 +0,0 @@
-# Contributor Code of Conduct
-
-As contributors and maintainers of this project, we pledge to respect all people who contribute through reporting issues, posting feature requests, updating documentation, submitting pull requests or patches, and other activities.
-
-We are committed to making participation in this project a harassment-free experience for everyone, regardless of level of experience, gender, gender identity and expression, sexual orientation, disability, personal appearance, body size, race, ethnicity, age, or religion.
-
-Examples of unacceptable behavior by participants include the use of sexual language or imagery, derogatory comments or personal attacks, trolling, public or private harassment, insults, or other unprofessional conduct.
-
-Project maintainers have the right and responsibility to remove, edit, or reject comments, commits, code, wiki edits, issues, and other contributions that are not aligned to this Code of Conduct. Project maintainers who do not follow the Code of Conduct may be removed from the project team.
-
-Instances of abusive, harassing, or otherwise unacceptable behavior may be reported by opening an issue or contacting one or more of the project maintainers.
-
-This Code of Conduct is adapted from the [Contributor Covenant](http://contributor-covenant.org), version 1.0.0, available at [http://contributor-covenant.org/version/1/0/0/](http://contributor-covenant.org/version/1/0/0/)