jn_yn.c File Reference

#include <math.h>
#include <float.h>
#include "ieee754/ieee754.h"

Include dependency graph for jn_yn.c:

Go to the source code of this file.

Typedefs
typedef int	fpclass_t

Functions
fpclass_t	_fpclass (double __d)
int *	_errno (void)
double	_jn (int n, double num)
double	_yn (int order, double num)

Typedef Documentation

fpclass_t

typedef int fpclass_t

Definition at line 5 of file jn_yn.c.

Function Documentation

_errno()

int * _errno

(

void

)

Definition at line 17 of file errno.c.

{
    return &(msvcrt_get_thread_data()->thread_errno);
}

Referenced by _get_errno(), _set_errno(), and _yn().

_fpclass()

fpclass_t _fpclass

(

double

__d

)

Definition at line 17 of file fpclass.c.

{
    union
    {
        double*   __d;
        double_s*   d;
    } d;
    d.__d = &__d;
 
 
    /* With 0x7ff, it can only be infinity or NaN */
    if (d.d->exponent == 0x7ff)
    {
        if (d.d->mantissah == 0 && d.d->mantissal == 0)
        {
            return (d.d->sign == 0) ? _FPCLASS_PINF : _FPCLASS_NINF;
        }
        /* Windows will never return Signaling NaN */
        else
        {
            return _FPCLASS_QNAN;
        }
    }
 
    /* With 0, it can only be zero or denormalized number */
    if (d.d->exponent == 0)
    {
        if (d.d->mantissah == 0 && d.d->mantissal == 0)
        {
            return (d.d->sign == 0) ? _FPCLASS_PZ : _FPCLASS_NZ;
        }
        else
        {
            return (d.d->sign == 0) ? _FPCLASS_PD : _FPCLASS_ND;
        }
    }
    /* Only remain normalized numbers */
    else
    {
        return (d.d->sign == 0) ? _FPCLASS_PN : _FPCLASS_NN;
    }
}

Referenced by _yn().

_jn()

double _jn	(	int	n,
		double	num
	)

Definition at line 12 of file jn_yn.c.

{
  /* FIXME: errno handling */
  return __ieee754_jn(n, num);
}

_yn()

double _yn	(	int	order,
		double	num
	)

Definition at line 21 of file jn_yn.c.

{
  double retval;
  if (!_finite(num)) *_errno() = EDOM;
  retval  = __ieee754_yn(order,num);
  if (_fpclass(retval) == _FPCLASS_NINF)
  {
    *_errno() = EDOM;
    retval = sqrt(-1);
  }
  return retval;
}