lib/sdk/crt/math/jn_yn.c

   1 #include <math.h>
   2
   3 typedef int fpclass_t;
   4 fpclass_t _fpclass(double __d);
   5 int *_errno(void);
   6
   7 /*
   8  * @unimplemented
   9  */
  10 double _jn(int n, double num)
  11 {
  12   /* FIXME: errno handling */
  13   return jn(n, num);
  14 }
  15
  16 /*
  17  * @implemented
  18  */
  19 double _yn(int order, double num)
  20 {
  21   double retval;
  22   if (!isfinite(num)) *_errno() = EDOM;
  23   retval  = yn(order,num);
  24   if (_fpclass(retval) == _FPCLASS_NINF)
  25   {
  26     *_errno() = EDOM;
  27     retval = sqrt(-1);
  28   }
  29   return retval;
  30 }