/*
    Copyright (C) 2001 Nick Kurshev <nickols_k@mail.ru>
    Based on libm-test.c testsuite from glibc-2.1.3

    $Id: testm.c,v 1.2 2001/07/20 07:02:56 konst Exp $
*/
#ifndef _GNU_SOURCE
# define _GNU_SOURCE
#endif

# define __NO_MATH_INLINES	1
# define __USE_ISOC9X 1
# define __USE_MISC 1
# define __USE_XOPEN_EXTENDED 1
# define __USE_GNU 1
# define __USE_BSD 1
# define __USE_XOPEN 1
# define __USE_SVID 1
#include <stdio.h>
#include <stdlib.h>
#include <assert.h>

#define	M_LOG2E		1.4426950408889634074	/* log 2e */
#define	M_LOG10E	0.43429448190325182765	/* log 10e */
#define	M_LN2		0.69314718055994530942	/* log e2 */
#define	M_LN10		2.30258509299404568402	/* log e10 */
#define	M_E		2.7182818284590452354	/* e */
#define	M_PI		3.14159265358979323846	/* pi */
#if __FreeBSD_version < 501103
extern char __infinity[];
#define	HUGE_VAL	(*(double *) __infinity)
#else
#define	HUGE_VAL	__builtin_huge_val()
#endif

#ifndef NAN
#define NAN (0./0.)
#endif
#ifndef INF
#define INF HUGE_VAL
#endif

typedef double MATHTYPE;
#define FUNC(name) name
/*
typedef MATHTYPE long double;
#define FUNC(name) name##l

typedef MATHTYPE float;
#define FUNC(name) name##f
*/

#define EPS 10e-10

#define MATH_ASSERT(name,expr,result)\
  printf("%s - %s computed: %10.15e should be: %10.15e\n",name,\
  (expr >= result-EPS && expr <= result+EPS) ? "ok" : "fail",(double)expr,(double)result)
#define INT_ASSERT(name,expr,result)\
  printf("%s - %s computed: %i should be: %i\n",name,\
  (expr == result) ? "ok" : "fail",(int)expr,(int)result)
#define LLINT_ASSERT(name,expr,result)\
  printf("%s - %s computed: %lli should be: %lli\n",name,\
  (expr == result) ? "ok" : "fail",(long long int)expr,(long long int)result)
#define NAN_ASSERT(name,expr)\
  printf("%s - %s (NAN version)\n",name,\
  isnan(expr) ? "ok" : "fail")

int main ( void )
{
  MATHTYPE xtmp,ytmp;
  int	itmp;
  printf("*** It's test for new mathlib *** !!!\n\n");
  MATH_ASSERT("acos",acos(1.) , 0.);
  MATH_ASSERT("acos",acos(-1.) , M_PI);
  MATH_ASSERT("acos",acos(-0.5) , 2.*M_PI/3.);
  printf("\n");
  MATH_ASSERT("acosh",acosh(1.) , 0.);
  MATH_ASSERT("acosh",acosh(7.) , 2.6339157938496334172L);
  MATH_ASSERT("acosh",acosh(INF) , INF);
  printf("\n");
  MATH_ASSERT("asin",asin(0.) , 0.);
  MATH_ASSERT("asin",asin(1.) , M_PI/2.);
  MATH_ASSERT("asin",asin(0.5) , M_PI/6.);
  printf("\n");
  MATH_ASSERT("asinh",asinh(0.) , 0.);
  MATH_ASSERT("asinh",asinh(0.7) , 0.652666566082355786L);
  MATH_ASSERT("asinh",asinh(INF) , INF);
  printf("\n");
  MATH_ASSERT("atan",atan(0.) , 0.);
  MATH_ASSERT("atan",atan(1.) , M_PI/4);
  MATH_ASSERT("atan",atan(0.7) , 0.6107259643892086165L);
  printf("\n");
  MATH_ASSERT("atan2",atan2 (-0.,+0.) , -0.);
  MATH_ASSERT("atan2",atan2 (0.4,0.0003) , 1.5700463269355215718L);
  MATH_ASSERT("atan2",atan2 (3., INF) , 0.);
  MATH_ASSERT("atan2",atan2 (-3., INF) , -0.);
  MATH_ASSERT("atan2",atan2 (3., -INF) , M_PI);
  MATH_ASSERT("atan2",atan2 (-3., -INF) , -M_PI);
  MATH_ASSERT("atan2",atan2 (INF, 2.) , M_PI/2);
  MATH_ASSERT("atan2",atan2 (-INF, -3.) , -M_PI/2);
  MATH_ASSERT("atan2",atan2 (INF, INF) , M_PI/4);
  MATH_ASSERT("atan2",atan2 (-INF, INF) , -M_PI/4);
  MATH_ASSERT("atan2",atan2 (INF, -INF) , 3*M_PI/4);
  MATH_ASSERT("atan2",atan2 (-INF, -INF) , -3*M_PI/4);
  printf("\n");
  MATH_ASSERT("atanh",atanh(0.) , 0.);
  MATH_ASSERT("atanh",atanh(0.7) , 0.8673005276940531944L);
  printf("\n");
  MATH_ASSERT("cbrt",cbrt (-0.) , -0.);
  MATH_ASSERT("cbrt",cbrt (-0.001) , -0.1);
  MATH_ASSERT("cbrt",cbrt (27.) , 3.);
  MATH_ASSERT("cbrt",cbrt (0.7) , 0.8879040017426007084L);
  MATH_ASSERT("cbrt",cbrt (-8.) , -2.);
  MATH_ASSERT("cbrt",cbrt (INF) , INF);
  MATH_ASSERT("cbrt",cbrt (-INF) , -INF);
  NAN_ASSERT("cbrt",cbrt (NAN));
  printf("\n");
  MATH_ASSERT("ceil",ceil (-0.) , -0.);
  MATH_ASSERT("ceil",ceil (M_PI) , 4.);
  MATH_ASSERT("ceil",ceil (-M_PI) , -3.);
  MATH_ASSERT("ceil",ceil (INF) , INF);
  MATH_ASSERT("ceil",ceil (-INF) , -INF);
  printf("\n");
  MATH_ASSERT("copysign",copysign (0., -4.) , -0.);
  MATH_ASSERT("copysign",copysign (0., 4.) , 0.);
  MATH_ASSERT("copysign",copysign (INF, -4.) , -INF);
  MATH_ASSERT("copysign",copysign (INF, 4.) , INF);
  MATH_ASSERT("copysign",copysign (-INF, -4.) , -INF);
  MATH_ASSERT("copysign",copysign (-INF, 4.) , INF);
  NAN_ASSERT("copysign",copysign (NAN, -4.));
  NAN_ASSERT("copysign",copysign (NAN, 4.));
  NAN_ASSERT("copysign",copysign (-NAN, -4.));
  NAN_ASSERT("copysign",copysign (-NAN, 4.));
  printf("\n");
  MATH_ASSERT("cos",cos (+0.) , 1.);
  MATH_ASSERT("cos",cos (M_PI/2.) , 0.);
  MATH_ASSERT("cos",cos (2.*M_PI/3.) , -0.5);
  printf("\n");
  MATH_ASSERT("cosh",cosh (0.) , 1.);
  MATH_ASSERT("cosh",cosh (0.7) , 1.255169005630943018L);
  INT_ASSERT("cosh",isinf(cosh (INF)) , 1);
  INT_ASSERT("cosh",isinf(cosh (-INF)) , 1);
  printf("\n");
/* erf */
/* erfc */
  MATH_ASSERT("exp",exp (0.) , 1.);
  MATH_ASSERT("exp",exp (1.) , M_E);
  MATH_ASSERT("exp",exp (0.7) , 2.0137527074704765216L);
  INT_ASSERT("exp",isinf(exp (INF)) , 1);
  MATH_ASSERT("exp",exp (-INF) , 0.);
  printf("\n");
  MATH_ASSERT("exp10",exp10 (0.) , 1.);
  MATH_ASSERT("exp10",exp10 (3.) , 1000.);
  MATH_ASSERT("exp10",exp10 (0.7) , 5.0118723362727228500L);
  INT_ASSERT("exp10",isinf(exp10 (INF)) , 1);
  MATH_ASSERT("exp10",exp10 (-INF) , 0.);
  printf("\n");
  MATH_ASSERT("exp2",exp2 (0.) , 1.);
  MATH_ASSERT("exp2",exp2 (10.) , 1024.);
  MATH_ASSERT("exp2",exp2 (0.7) , 1.6245047927124710452L);
  INT_ASSERT("exp2",isinf(exp2 (INF)) , 1);
  MATH_ASSERT("exp2",exp2 (-INF) , 0.);
  printf("\n");
  MATH_ASSERT("expm1",expm1(0.), 0.);
  MATH_ASSERT("expm1",expm1(INF), INF);
  MATH_ASSERT("expm1",expm1(-INF), -1.);
  MATH_ASSERT("expm1",expm1(1.), M_E -1.);
  MATH_ASSERT("expm1",expm1(0.7), 1.0137527074704765216L);
  printf("\n");
  MATH_ASSERT("fabs",fabs(-0.) , 0.);
  MATH_ASSERT("fabs",fabs(-1.22) , 1.22);
  MATH_ASSERT("fabs",fabs(38.) , 38.);
  MATH_ASSERT("fabs",fabs(INF) , INF);
  MATH_ASSERT("fabs",fabs(-INF) , INF);
  printf("\n");
  MATH_ASSERT("fdim",fdim (0., 0.) , 0.);
  MATH_ASSERT("fdim",fdim (9., 0.) , 9.);
  MATH_ASSERT("fdim",fdim (0., 9.) , 9.);
  MATH_ASSERT("fdim",fdim (-9., 0.) , 9.);
  MATH_ASSERT("fdim",fdim (0., -9.) , 9.);
  printf("\n");
  MATH_ASSERT("floor",floor (-0.) , -0.);
  MATH_ASSERT("floor",floor (M_PI) , 3.);
  MATH_ASSERT("floor",floor (-M_PI) , -4.);
  MATH_ASSERT("floor",floor(INF) , INF);
  MATH_ASSERT("floor",floor(-INF) , -INF);
  printf("\n");
  MATH_ASSERT("fma",fma(1.0, 2.0, 3.0) , 5.0);
  printf("\n");
  MATH_ASSERT("fmax",fmax (0., 0.) , 0.);
  MATH_ASSERT("fmax",fmax (9., 0.) , 9.);
  MATH_ASSERT("fmax",fmax (0., -9.) , 0.);
  MATH_ASSERT("fmax",fmax (9., INF) , INF);
  MATH_ASSERT("fmax",fmax (0., -INF) , 0.);
  MATH_ASSERT("fmax",fmax (9., NAN) , 9.);
  MATH_ASSERT("fmax",fmax (-9., NAN) , -9.);
  MATH_ASSERT("fmax",fmax (INF, NAN) , INF);
  MATH_ASSERT("fmax",fmax (-INF, NAN) , -INF);
  printf("\n");
  MATH_ASSERT("fmin",fmin (0., 0.) , 0.);
  MATH_ASSERT("fmin",fmin (9., 0.) , 0.);
  MATH_ASSERT("fmin",fmin (0., -9.) , -9.);
  MATH_ASSERT("fmin",fmin (9., INF) , 9.);
  MATH_ASSERT("fmin",fmin (0., -INF) , -INF);
  MATH_ASSERT("fmin",fmin (9., NAN) , 9.);
  MATH_ASSERT("fmin",fmin (-9., NAN) , -9.);
  MATH_ASSERT("fmin",fmin (INF, NAN) , INF);
  MATH_ASSERT("fmin",fmin (-INF, NAN) , -INF);
  printf("\n");
  MATH_ASSERT("fmod",fmod (6.5, 2.3) , 1.9);
  MATH_ASSERT("fmod",fmod (-6.5, 2.3) , -1.9);
  MATH_ASSERT("fmod",fmod (6.5, -2.3) , 1.9);
  MATH_ASSERT("fmod",fmod (-6.5, -2.3) , -1.9);
  MATH_ASSERT("fmod",fmod (6.5, INF) , 6.5);
  MATH_ASSERT("fmod",fmod (-6.5, -INF) , -6.5);
  printf("\n");
/* fpclassify */
  MATH_ASSERT("frexp",frexp(-27.34, &itmp), -0.854375L);
  INT_ASSERT("frexp",itmp, 5);
  MATH_ASSERT("frexp",frexp(0., &itmp), 0.);
  INT_ASSERT("frexp",itmp, 0);
  MATH_ASSERT("frexp",frexp(12.8, &itmp), 0.8L);
  INT_ASSERT("frexp",itmp, 4);
  MATH_ASSERT("frexp",frexp(INF, &itmp), INF);
  MATH_ASSERT("frexp",frexp(-INF, &itmp), -INF);
  printf("\n");
/* gamma */
  MATH_ASSERT("hypot",hypot(3., 4.) , 5.);
  MATH_ASSERT("hypot",hypot(0.7, 1.2) , 1.3892443989449804508L);
  MATH_ASSERT("hypot",hypot(INF, NAN) , INF);
  MATH_ASSERT("hypot",hypot(-INF, NAN) , INF);
  NAN_ASSERT("hypot",hypot(NAN, NAN));
  printf("\n");
  INT_ASSERT("ilogb",ilogb (1.) , 0);
  INT_ASSERT("ilogb",ilogb (M_E) , 1);
  INT_ASSERT("ilogb",ilogb (1024) , 10);
  INT_ASSERT("ilogb",ilogb (-2000.) , 10);
  printf("\n");
  INT_ASSERT("isfinite", isfinite(0.), 1);
  INT_ASSERT("isfinite", isfinite(-0.), 1);
  INT_ASSERT("isfinite", isfinite(10.), 1);
  INT_ASSERT("isfinite", isfinite(INF), 0);
  INT_ASSERT("isfinite", isfinite(-INF), 0);
  INT_ASSERT("isfinite", isfinite(NAN), 0);
  printf("\n");
#if 0  
  INT_ASSERT("isnormal", isnormal(0.), 0);
  INT_ASSERT("isnormal", isnormal(-0.), 0);
  INT_ASSERT("isnormal", isnormal(10.), 1);
  INT_ASSERT("isnormal", isnormal(INF), 0);
  INT_ASSERT("isnormal", isnormal(-INF), 0);
  INT_ASSERT("isnormal", isnormal(NAN), 0);
  printf("\n");
#endif
  INT_ASSERT("isnan", isnan(0.), 0);
  INT_ASSERT("isnan", isnan(-0.), 0);
  INT_ASSERT("isnan", isnan(10.), 0);
  INT_ASSERT("isnan", isnan(INF), 0);
  INT_ASSERT("isnan", isnan(-INF), 0);
  INT_ASSERT("isnan", isnan(NAN), 1);
  printf("\n");
  INT_ASSERT("isinf", isinf(0.), 0);
  INT_ASSERT("isinf", isinf(-0.), 0);
  INT_ASSERT("isinf", isinf(10.), 0);
  INT_ASSERT("isinf", isinf(INF), 1);
  INT_ASSERT("isinf", isinf(-INF), -1);
  INT_ASSERT("isinf", isinf(NAN), 0);
  printf("\n");
/* lgamma */
  MATH_ASSERT("ldexp",ldexp (0., 0.) , 0.);
  MATH_ASSERT("ldexp",ldexp (0.8, 4.) , 12.8);
  MATH_ASSERT("ldexp",ldexp (-0.854375L, 5.) , -27.34L);
  MATH_ASSERT("ldexp",ldexp (INF, 1.) , INF);
  MATH_ASSERT("ldexp",ldexp (-INF, 1.) , -INF);
  NAN_ASSERT("ldexp",ldexp (NAN, 1.));
  printf("\n");
  LLINT_ASSERT("llrint", llrint(0.), 0);
  LLINT_ASSERT("llrint", llrint(-0.2), 0);
  LLINT_ASSERT("llrint", llrint(-1.4), -1);
  LLINT_ASSERT("llrint", llrint(72057594037927936.0), 72057594037927936LL);
  printf("\n");
  MATH_ASSERT("log",log (M_E) , 1.);
  MATH_ASSERT("log",log (1./M_E) , -1.);
  MATH_ASSERT("log",log (2.) , M_LN2);
  MATH_ASSERT("log",log (10.) , M_LN10);
  MATH_ASSERT("log",log (INF) , INF);
  printf("\n");
  MATH_ASSERT("log10",log10 (0.1) , -1.);
  MATH_ASSERT("log10",log10 (10.) , 1.);
  MATH_ASSERT("log10",log10 (M_E) , M_LOG10E);
  MATH_ASSERT("log10",log10 (10000.) , 4.);
  MATH_ASSERT("log10",log10 (INF) , INF);
  printf("\n");
  MATH_ASSERT("log2",log2 (1.) , 0.);
  MATH_ASSERT("log2",log2 (M_E) , M_LOG2E);
  MATH_ASSERT("log2",log2 (2.) , 1.);
  MATH_ASSERT("log2",log2 (16.) , 4.);
  MATH_ASSERT("log2",log2 (INF) , INF);
  printf("\n");
  MATH_ASSERT("log1p",log1p (0.) , 0.);
  MATH_ASSERT("log1p",log1p (M_E - 1.0) , 1.);
  MATH_ASSERT("log1p",log1p (-0.3) , -0.35667494393873237891L);
  MATH_ASSERT("log1p",log1p (INF) , INF);
  printf("\n");
  MATH_ASSERT("logb",logb (1.) , 0.);
  MATH_ASSERT("logb",logb (M_E) , 1.);
  MATH_ASSERT("logb",logb (1024.) , 10.);
  MATH_ASSERT("logb",logb (-2000.) , 10.);
  MATH_ASSERT("logb",logb (INF) , INF);
  printf("\n");
  MATH_ASSERT("lrint",lrint(0.) , 0.);
  MATH_ASSERT("lrint",lrint(0.4) , 0.);
  MATH_ASSERT("lrint",lrint(1.4) , 1.);
  MATH_ASSERT("lrint",lrint(-1.4) , -1.);
  printf("\n");
/* modf */
#if 0
  /* Sometime those functions raise exceptions that is wrong */
  MATH_ASSERT("nearbyint",nearbyint(copysign(0., -1.)), copysign(0., -1.));
  MATH_ASSERT("nearbyint",nearbyint(+0.), +0.);
  MATH_ASSERT("nearbyint",nearbyint(INF), INF);
  MATH_ASSERT("nearbyint",nearbyint(-INF), -INF);
  MATH_ASSERT("nearbyint",nearbyint(+5.1), +5.0);
  MATH_ASSERT("nearbyint",nearbyint(+4.9), +5.0);
  printf("\n");
#endif
/* nextafter */
  MATH_ASSERT("pow",pow(0., 1.25) , 0.);
  MATH_ASSERT("pow",pow(3.16, 0.) , 1.);
  MATH_ASSERT("pow",pow(5., 3.) , 125.);
  MATH_ASSERT("pow",pow(36., 0.5) , 6.);
  MATH_ASSERT("pow",pow(10., 3.) , 1000.);
  MATH_ASSERT("pow",pow(2., 3.) , 8.);
  NAN_ASSERT("pow",pow(NAN, 1.));
  NAN_ASSERT("pow",pow(NAN, NAN));
  NAN_ASSERT("pow",pow(0., NAN));
  MATH_ASSERT("pow",pow(INF, 1.) , INF);
  MATH_ASSERT("pow",pow(INF, INF) , INF);
  MATH_ASSERT("pow",pow(1., INF) , 1.);
  MATH_ASSERT("pow",pow(0., INF) , 0.);
  printf("\n");
  MATH_ASSERT("pow10",pow10(0.) , 1.);
  MATH_ASSERT("pow10",pow10(1.) , 10.);
  MATH_ASSERT("pow10",pow10(3.) , 1000.);
  printf("\n");
  MATH_ASSERT("remainder",remainder(1.625, 1.0) , -0.375);
  MATH_ASSERT("remainder",remainder(-1.625, 1.0) , 0.375);
  MATH_ASSERT("remainder",remainder(1.625, -1.0) , -0.375);
  MATH_ASSERT("remainder",remainder(-1.625, -1.0) , 0.375);
  printf("\n");
  MATH_ASSERT("rint",rint(-0.) , -0.);
  printf("\n");
/* remquo */
  MATH_ASSERT("scalb",scalb(0.8, 4.), 12.8);
  MATH_ASSERT("scalb",scalb(-0.854375L, 5), -27.34L);
  MATH_ASSERT("scalb",scalb(0., 2.), 0.);
  NAN_ASSERT("scalb",scalb(3., -2.5));
  NAN_ASSERT("scalb",scalb(2., 0.5));
  MATH_ASSERT("scalb",scalb(INF, 0.), INF);
  MATH_ASSERT("scalb",scalb(-INF, 0.), -INF);
  NAN_ASSERT("scalb",scalb(INF, -INF));
  NAN_ASSERT("scalb",scalb(-INF, -INF));
  MATH_ASSERT("scalb",scalb(3.,  -INF), 0.);
  MATH_ASSERT("scalb",scalb(-2., -INF), -0.);
  printf("\n");
  MATH_ASSERT("scalbn",scalbn (0., 0.) , 0.);
  MATH_ASSERT("scalbn",scalbn (0.8, 4.) , 12.8);
  MATH_ASSERT("scalbn",scalbn (-0.854375L, 5.) , -27.34L);
  MATH_ASSERT("scalbn",scalbn (INF, 1.) , INF);
  MATH_ASSERT("scalbn",scalbn (-INF, 1.) , -INF);
  NAN_ASSERT("scalbn",scalbn (NAN, 1.));
  printf("\n");
/*
  MATH_ASSERT("significand",significand (0.) , 0.);
  Later
  printf("\n");
*/
/* signbit */
  MATH_ASSERT("sin",sin (-0.) , -0.);
  MATH_ASSERT("sin",sin (M_PI/2.) , 1.);
  MATH_ASSERT("sin",sin (-M_PI/2.) , -1.);
  MATH_ASSERT("sin",sin (0.7) , 0.64421768723769105367L);
  printf("\n");
  sincos(0.,&xtmp,&ytmp);
  MATH_ASSERT("sincos",xtmp, 0.);
  MATH_ASSERT("sincos",ytmp, 1.);
  sincos(M_PI/2.,&xtmp,&ytmp);
  MATH_ASSERT("sincos",xtmp, 1.);
  MATH_ASSERT("sincos",ytmp, 0.);
  sincos(0.7,&xtmp,&ytmp);
  MATH_ASSERT("sincos",xtmp, 0.64421768723769105367L);
  MATH_ASSERT("sincos",ytmp, 0.76484218728448842626L);
  printf("\n");
  MATH_ASSERT("sinh",sinh (0.) , 0.);
  MATH_ASSERT("sinh",sinh (0.7) , 0.75858370183953350346L);
  MATH_ASSERT("sinh",sinh (INF) , INF);
  MATH_ASSERT("sinh",sinh (-INF) , -INF);
  printf("\n");
  MATH_ASSERT("sqrt",sqrt(9.) , 3.);
  MATH_ASSERT("sqrt",sqrt(0.25) , 0.5);
  MATH_ASSERT("sqrt",sqrt(0.7) , 0.83666002653407554798L);
  NAN_ASSERT("sqrt",sqrt(NAN));
  MATH_ASSERT("sqrt",sqrt(INF) , INF);
  printf("\n");
  MATH_ASSERT("tan",tan (-0.) , -0.);
  MATH_ASSERT("tan",tan (M_PI/4.) , 1.);
  MATH_ASSERT("tan",tan (0.7) , 0.84228838046307944813L);
  printf("\n");
  MATH_ASSERT("tanh",tanh (0.) , 0.);
  MATH_ASSERT("tanh",tanh (0.7) , 0.60436777711716349631L);
  MATH_ASSERT("tanh",tanh (INF) , 1.);
  MATH_ASSERT("tanh",tanh (-INF) , -1.);
  printf("\n");
/* tgamma */
  MATH_ASSERT("trunc",trunc(0.) , 0.);
  MATH_ASSERT("trunc",trunc(0.625) , 0.);
  MATH_ASSERT("trunc",trunc(1048580.625L) , 1048580L);
  MATH_ASSERT("trunc",trunc(INF) , INF);
  MATH_ASSERT("trunc",trunc(-INF) , -INF);
  NAN_ASSERT("trunc",trunc(NAN));
  printf("\n");
  exit(0); /* no errors */
  return 0;
}