Home | History | Annotate | Line # | Download | only in libm 
      1 #!/usr/bin/lua
      2 
      3 local lm = require("libm")
      4 
      5 local fmtf = "%-24s%+2.13f"
      6 local fmti = "%-24s%d"
      7 
      8 local function test(s, x)
      9    print(string.format(fmtf, s, x))
     10 end
     11 
     12 local function testb(s, x)
     13    print(string.format(fmti, s, x and 1 or 0))
     14 end
     15 
     16 test("M_E", lm.M_E);
     17 test("M_LOG2E", lm.M_LOG2E)
     18 test("M_LOG10E", lm.M_LOG10E)
     19 test("M_LN2", lm.M_LN2)
     20 test("M_LN10", lm.M_LN10)
     21 test("M_PI", lm.M_PI)
     22 test("M_PI_2", lm.M_PI_2)
     23 test("M_PI_4", lm.M_PI_4)
     24 test("M_1_PI", lm.M_1_PI)
     25 test("M_2_PI", lm.M_2_PI)
     26 test("M_2_SQRTPI", lm.M_2_SQRTPI)
     27 test("M_SQRT2", lm.M_SQRT2)
     28 test("M_SQRT1_2", lm.M_SQRT1_2)
     29 
     30 test("cos(M_PI_2)", lm.cos(lm.M_PI_2))
     31 test("acos(cos(M_PI_2))", lm.acos(lm.cos(lm.M_PI_2)))
     32 test("cosh(M_SQRT1_2)", lm.cosh(lm.M_SQRT1_2))
     33 test("acosh(cosh(M_SQRT1_2))", lm.acosh(lm.cosh(lm.M_SQRT1_2)))
     34 test("sin(M_PI_2)", lm.sin(lm.M_PI_2))
     35 test("asin(sin(M_PI_2))", lm.asin(lm.sin(lm.M_PI_2)))
     36 test("sinh(M_PI_2)", lm.sinh(lm.M_PI_2))
     37 test("asinh(sinh(M_PI_2))", lm.asinh(lm.sinh(lm.M_PI_2)))
     38 test("tan(M_SQRT2)", lm.tan(lm.M_SQRT2))
     39 test("atan(tan(M_SQRT2))", lm.atan(lm.tan(lm.M_SQRT2)))
     40 test("tanh(M_SQRT2)", lm.tanh(lm.M_SQRT2))
     41 test("atanh(tanh(M_SQRT2))", lm.atanh(lm.tanh(lm.M_SQRT2)))
     42 test("atan2(M_SQRT2, M_SQRT2)", lm.atan2(lm.M_SQRT2, lm.M_SQRT2))
     43 test("cbrt(8.0)", lm.cbrt(8.0))
     44 test("ceil(M_PI)", lm.ceil(lm.M_PI))
     45 test("ceil(M_E)", lm.ceil(lm.M_E))
     46 test("copysign(-2.0, 10.0)", lm.copysign(-2.0, 10.0))
     47 test("copysign(2.0, -10.0)", lm.copysign(2.0, -10.0))
     48 test("erf(M_SQRT1_2)", lm.erf(lm.M_SQRT1_2))
     49 test("erfc(M_SQRT1_2)", lm.erfc(lm.M_SQRT1_2))
     50 test("exp(M_PI_4)", lm.exp(lm.M_PI_4))
     51 test("exp2(3.0)", lm.exp2(3.0))
     52 test("expm1(M_PI_4)", lm.expm1(lm.M_PI_4))
     53 test("fabs(-M_SQRT2)", lm.fabs(-lm.M_SQRT2))
     54 test("fabs(M_SQRT2)", lm.fabs(lm.M_SQRT2))
     55 test("fdim(M_PI, M_PI_2)", lm.fdim(lm.M_PI, lm.M_PI_2))
     56 test("fdim(M_PI, -M_PI_2)", lm.fdim(lm.M_PI, -lm.M_PI_2))
     57 testb("finite(1.0 / 0.0)", lm.finite(1.0 / 0.0))
     58 test("floor(M_PI)", lm.floor(lm.M_PI))
     59 test("floor(M_E)", lm.floor(lm.M_E))
     60 test("fma(M_PI, M_E, M_SQRT2)", lm.fma(lm.M_PI, lm.M_E, lm.M_SQRT2))
     61 test("fmax(M_PI, M_E)", lm.fmax(lm.M_PI, lm.M_E))
     62 test("fmin(M_PI, M_E)", lm.fmin(lm.M_PI, lm.M_E))
     63 test("fmod(100.5, 10.0)", lm.fmod(100.5, 10.0))
     64 test("gamma(5.0)", lm.gamma(5.0))
     65 test("hypot(3.0, 4.0)", lm.hypot(3.0, 4.0))
     66 print(string.format(fmti, "ilogb(10.0)", lm.ilogb(10.0)))
     67 testb("isinf(1.0 / 0.0)", lm.isinf(1.0 / 0.0))
     68 testb("isnan(log(-1.0))", lm.isnan(lm.log(-1.0)))
     69 test("j0(M_PI)", lm.j0(lm.M_PI))
     70 test("j1(M_PI)", lm.j1(lm.M_PI))
     71 test("jn(3, M_PI)", lm.jn(3, lm.M_PI))
     72 test("lgamma(M_PI)", lm.lgamma(lm.M_PI))
     73 test("log(M_E)", lm.log(lm.M_E))
     74 test("log10(100.0)", lm.log10(100.0))
     75 test("log1p(M_PI)", lm.log1p(lm.M_PI))
     76 test("nan(\"\")", lm.nan(""))
     77 -- XXX vax test("nextafter(1.0e-14, 1.0)", lm.nextafter(1.0e-14, 1.0))
     78 test("pow(M_SQRT2, 2.0)", lm.pow(lm.M_SQRT2, 2.0))
     79 test("remainder(M_PI, M_E)", lm.remainder(lm.M_PI, lm.M_E))
     80 test("rint(M_PI)", lm.rint(lm.M_PI))
     81 test("rint(M_E)", lm.rint(lm.M_E))
     82 print(string.format(fmtf, "scalbn(1.0, 2)", lm.scalbn(1.0, 2)))
     83 test("sin(M_PI_4)", lm.sin(lm.M_PI_4))
     84 test("sinh(M_PI_4)", lm.sinh(lm.M_PI_4))
     85 test("sqrt(9.0)", lm.sqrt(9.0))
     86 test("tan(M_PI_4)", lm.tan(lm.M_PI_4))
     87 test("tanh(M_PI_4)", lm.tanh(lm.M_PI_4))
     88 test("trunc(M_PI)", lm.trunc(lm.M_PI))
     89 test("trunc(M_E)", lm.trunc(lm.M_E))
     90 test("y0(M_PI)", lm.y0(lm.M_PI))
     91 test("y1(M_PI)", lm.y1(lm.M_PI))
     92 test("yn(3, M_PI)", lm.yn(3, lm.M_PI))
     93