improve pixie.math scope (adapted from work of @edw found here:

https://gist.github.com/edw/87a4e234a33b5a1e0377)

Author: mpenet

pixie/math.pxi

@@ -4,52 +4,93 @@
 (i/with-config {:library "m"
                 :cxx-flags ["-lm"]
                 :includes ["math.h"]}
-  (i/defcfn sin)
-  (i/defcfn cos)
-  (i/defcfn tan)
+  (i/defconst M_E)
+  (i/defconst M_LOG2E)
+  (i/defconst M_LOG10E)
+  (i/defconst M_LN2)
+  (i/defconst M_LN10)
+  (i/defconst M_PI)
+  (i/defconst M_PI_2)
+  (i/defconst M_PI_4)
+  (i/defconst M_1_PI)
+  (i/defconst M_2_PI)
+  (i/defconst M_2_SQRTPI)
+  (i/defconst M_SQRT2)
+  (i/defconst M_SQRT1_2)
 
-  (i/defcfn asin)
-  (i/defcfn acos)
-  (i/defcfn atan)
-  (i/defcfn atan2) ; Arc tangent function of two variables.
+  (i/defcfn nan)
+  (i/defcfn ceil)
+  (i/defcfn floor)
+  (i/defcfn nearbyint)
+  (i/defcfn rint)
+  (i/defcfn lround)
+  (i/defcfn llrint)
+  (i/defcfn llround)
+  (i/defcfn trunc)
 
-  (i/defcfn sinh)
-  (i/defcfn cosh)
-  (i/defcfn tanh)
+  (i/defcfn fmod)
+  (i/defcfn remainder)
+  (i/defcfn remquo)
 
-  (i/defcfn asinh)
-  (i/defcfn acosh)
-  (i/defcfn atanh)
+  (i/defcfn fdim)
+  (i/defcfn fmax)
+  (i/defcfn fmin)
+
+  (i/defcfn fma)
+
+  (i/defcfn fabs)
+  (i/defcfn sqrt)
+  (i/defcfn cbrt)
+  (i/defcfn hypot)
 
   (i/defcfn exp)
-  (i/defcfn ldexp)
+  (i/defcfn exp2)
+  (i/defcfn exp10)
+  (i/defcfn expm1)
 
   (i/defcfn log)
   (i/defcfn log2)
   (i/defcfn log10)
   (i/defcfn log1p)
+
   (i/defcfn logb)
   (i/defcfn ilogb)
 
-  ;; (i/defcfn modf) ;; Needs ffi support
+  (i/defcfn modf)
+  (i/defcfn frexp)
+
+  (i/defcfn ldexp)
+  (i/defcfn scalbn)
+  (i/defcfn scalbln)
+
   (i/defcfn pow)
-  (i/defcfn sqrt)
-  (i/defcfn ceil)
-  (i/defcfn fabs)
-  (i/defcfn floor)
-  (i/defcfn fmod)
 
-  (i/defconst M_E)             ; base of natural logarithm, e
-  (i/defconst M_LOG2E)         ; log2(e)
-  (i/defconst M_LOG10E)        ; log10(e)
-  (i/defconst M_LN2)           ; ln(2)
-  (i/defconst M_LN10)          ; ln(10)
-  (i/defconst M_PI)            ; pi
-  (i/defconst M_PI_2)          ; pi / 2
-  (i/defconst M_PI_4)          ; pi / 4
-  (i/defconst M_1_PI)          ; 1 / pi
-  (i/defconst M_2_PI)          ; 2 / pi
-  (i/defconst M_2_SQRTPI)      ; 2 / sqrt(pi)
-  (i/defconst M_SQRT2)         ; sqrt(2)
-  (i/defconst M_SQRT1_2))      ; sqrt(1/2)
+  (i/defcfn cos)
+  (i/defcfn sin)
+  (i/defcfn tan)
+
+  (i/defcfn cosh)
+  (i/defcfn sinh)
+  (i/defcfn tanh)
+
+  (i/defcfn acos)
+  (i/defcfn asin)
+  (i/defcfn atan)
+  (i/defcfn atan2)
+
+  (i/defcfn acosh)
+  (i/defcfn asinh)
+  (i/defcfn atanh)
+
+  (i/defcfn tgamma)
+  (i/defcfn lgamma)
+
+  (i/defcfn j0)
+  (i/defcfn j1)
+  (i/defcfn jn)
+  (i/defcfn y0)
+  (i/defcfn y1)
+  (i/defcfn yn)
 
+  (i/defcfn erf)
+  (i/defcfn erfc))