#include "dk3all.h"
#include "dk3md.h"


DK3_MAIN
{
  long double res;
  int		ec = 0;
  $!trace-init test-md.deb
  res = dk3md_add_ok(3.0L, 5.0L, &ec);
  printf("%Lg\n", res);
  res = dk3md_sub_ok(3.0L, 5.0L, &ec);
  printf("%Lg\n", res);
  res = dk3md_mul_ok(3.0L, 5.0L, &ec);
  printf("%Lg\n", res);
  res = dk3md_div_ok(3.0L, 5.0L, &ec);
  printf("%Lg\n", res);
  res = dk3md_sin(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_cos(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_tan(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_asin(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_acos(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_atan(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_atan2(5.0L, 3.0L);
  printf("%Lg\n", res);
  res = dk3md_rint(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_floor(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_ceil(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_exp(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_log(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_sqrt(0.5 * M_PI_2);
  printf("%Lg\n", res);

  res = dk3md_sinh(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_cosh(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_tanh(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_asinh(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_acosh(0.5 * M_PI_2);
  printf("%Lg\n", res);
  res = dk3md_atanh(0.5 * M_PI_2);
  printf("%Lg\n", res);

  $!trace-end
}