#ifndef INC__MATH_AS
#define INC__MATH_AS

// 定数
#define M_E                     2.7182818284590452354           ; e
#define M_LOG2E         1.4426950408889634074           ; log_2(e)
#define M_LOG10E        0.43429448190325182765          ; log_10(e)
#define M_LN2           0.69314718055994530942          ; log_e(2)
#define M_LN10          2.30258509299404568402          ; log_e(10)
#define M_PI            3.14159265358979323846          ; π
#define M_SQRTPI        1.77245385090551602792981       ; √π
#define M_SQRT2         1.41421356237309504880          ; √2
#define M_SQRT3         1.73205080756887719000          ; √3

// doubleの限度
#const  DBL_DIG                 15                                              ; 10進数で有効な桁数
#define DBL_EPSILON             0.00000000000000022204460492503131      ; 1.0とその次に大きい値との差
#const  DBL_MANT_DIG    53                                              ; 仮数部のbit数
#define DBL_MAX                 double("1.7976931348623158e+308")       ; 最大値
#const  DBL_MAX_10_EXP  308                                             ; 10進数での指数部の最大値
#const  DBL_MAX_EXP             1024                                    ; 2進数での指数部の最大値
#define DBL_MIN                 double("2.2250738585072014e-308")       ; 0を超える最小の値
#const  DBL_MIN_10_EXP  -307                                    ; 10進数での指数部の最小値
#const  DBL_MIN_EXP             -1021                                   ; 2進数での指数部の最小値

// intの限度
#const  INT_DIGIT               32                                              ; 2進数で有効な桁数
#const  INT_DIGIT10             9                                               ; 10進数で有効な桁数
#const  INT_MAX                 2147483647                              ; 最大値
#const  INT_MIN                 -INT_MAX-1                              ; 最小値

#define ctype rad2deg(%1)       (57.295779513082320877*(%1))    ; ラジアンを度に変換
#define ctype deg2rad(%1)       (0.017453292519943295769*(%1))  ; 度をラジアンに変換

#define ctype pow(%1,%2)        expf(logf(%1)*(%2))     ; x^y
#define ctype log10(%1)         (logf(%1)/logf(10))     ; log_10(x)
#define ctype log2(%1)          (logf(%1)/logf(2))      ; log_2(x)

#define ctype asin(%1)          atan((%1),sqrt(1.0-double(%1)*(%1)))    ; sin(y)=xとなるy
#define ctype acos(%1)          atan(sqrt(1.0-double(%1)*(%1)),(%1))    ; cos(y)=xとなるy

#define ctype sinh(%1)          ((expf(%1)-expf(-(%1)))/2)              ; sinh(x)
#define ctype cosh(%1)          ((expf(%1)+expf(-(%1)))/2)              ; cosh(x)
#define ctype tanh(%1)          (sinh(%1)/cosh(%1))                             ; tanh(x)

#define ctype asinh(%1)         logf(sqrt(double(%1)*(%1)+1.0)+(%1))    ; sinh(y)=xとなるy
#define ctype acosh(%1)         logf(sqrt(double(%1)*(%1)-1.0)+(%1))    ; cosh(y)=xとなるy
#define ctype atanh(%1)         (logf((1.0+(%1)) / (1.0-(%1)))/2.0)             ; tanh(y)=xとなるy

#define ctype isfinite(%1)      ((%1)*0=0)                              ; 有限／無限･非数の判定
#define ctype isnan(%1)         ((%1)!(%1))                             ; 非数の判定

#define ctype round(%1)         double(strf("%%0.0f", %1))              ; 四捨五入
#define ctype sgn(%1)           (((%1)>0) - ((%1)<0))                   ; 符号
#define ctype intf(%1)          round((-0.5)*sgn(%1)+(%1))              ; 0の方向へ丸め
#define ctype floor(%1)         ((intf(%1)-((%1)<intf(%1))))    ; 負の方向へ丸め
#define ctype ceil(%1)          (-floor(-(%1)))                                 ; 正の方向へ丸め
#define ctype fmod(%1,%2)       (double(%1)-(%2)*intf(double(%1)/(%2)))         ; モジュロ

#define ctype distance2(%1,%2)  sqrt((double(%1)-%2)*(double(%1)-%2)+(double(%1(1))-%2(1))*(double(%1(1))-%2(1)))       ; 2次元距離

#endif  ; !defined(INC__MATH_AS)