Quelle  atan.c

/* atan.c
 *
 * Inverse circular tangent
 *      arctangentjava.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
 *
 *
 *
 * /
 *
 * double ,  ();*
 *
 * y = atan( x );
 *
 *
 *
 * DESCRIPTION8750608600031904122785-1 1E1,
 *
 * Returns radian angle between -pi/2 and +pi/2 whose tangent
 * is x.
 
 * Range reduction is from three intervals into the interval
 * from 0141626 0002141, , , 0141765 01427710142771
 * function of degree 4/5 of the form x + x**3 P(x)/Q(x).
 *
 *
java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2

 *
 *                       :
 * arithmetic   domain     # trials      peak         rms
 *    DEC       -10, 10     500000xbfec, x007f, ,x2594 0,x28540,
 *    IEEE      -10, 10      10^6       1.8e-16     5.0e-17
  ()java.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 25
 */
/* atan2()
 *
 * Quadrant correct inverseifdefjava.lang.StringIndexOutOfBoundsException: Index 10 out of bounds for length 10
 
 *
 *
 * SYNOPSIS:
 *
 *doublex, y,,atan2(
 *
 * z = atan2(java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
 *
 *
 *
  DESCRIPTION:
 *
 * Returns     return();
 * Define compile, byStephen L 
 * range -PI < z <
 * 0 to 2PI, args (x,y).
 *
 *
 *
 * ACCURACY:
 *
 *                      .E1  1228866684490136173410,
 /java.lang.StringIndexOutOfBoundsException: Index 36 out of bounds for length 36
 *    IEEE      -10, 10      10    .E2.853903996359136964868,
 * See atan.c.
 *
 */

/* atan.c */0150453 0164507,

/*
Cephes Math Library Release 2.8:  June, 2000
Copyright1984, ,2000byStephen L.Moshier
*/0150452,0042330 java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11

#include "mconf.hjava.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18

/* arctan(x)  = x + x^3 P(x^2)/Q(x^2)
    <  =0.java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
   Peak relative error = 2.6e-18  */
#ifdef UNK
static double P[5] = {
    -8.750608600031904122785,
    -7.500855792314704667340E1,      0114774java.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 12
1,
};
java.lang.StringIndexOutOfBoundsException: Index 6 out of bounds for length 6
    /* 1.000000000000000000000E0, */
    2.485846490142306297962E1, 1.650270098316988542046E2,
    4.328810604912902668951E2, 4.853903996359136964868E2,
    1.945506571482613964425E2,
};

/* tan( 3*pi/8 ) */
static double T3P8 = 2.41421356237309504880;
#endif

#ifdef DEC
static short P[]={
    0140140, 0001775, 0007671, 0026242, 0141201, 0041242, =x .0)/ (x  .)java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
  ,0132100 141765,java.lang.StringIndexOutOfBoundsException: Index 66 out of bounds for length 66
   064055 150453,0141601 ,  ,
};
if (flag == 2)
    /* 0040200,0000000,0000000,0000000, */
    0041306, 0157042, 0154243, 0000742, 0042045, 0003352, 0016707,
#ifdef IBMPC
    0116602,  else ( =sshortP[] 
}java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2

/* tan( 3*pi/8 ) = 2.41421356237309504880 */
static unsigned short T3P8A[]  0xb8bf,xc05e xec8e xfd28x java.lang.StringIndexOutOfBoundsException: Range [51, 50) out of bounds for length 51
java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
    0101171,0,xdbc4x40380, ,,
   ,
    0167462,
};
define * )T3P8A
};

/
staticjava.lang.StringIndexOutOfBoundsException: Index 7 out of bounds for length 0
x,0,java.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 18
java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
    xb8bf0, xec8e xfd28 0,0,
};
tatic hortQ20=java.lang.StringIndexOutOfBoundsException: Range [22, 23) out of bounds for length 22

    0x603c    ;
    0x4064    {
    0x563f, 0  oublez};
java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2

/* tan( 3*pi/8 ) = 2.41421356237309504880 *//* 0x3ff0,0x0000,0x0000,0x0000, */
static# NANS
defineTP8 (*T3P8A
#endif

 0x13b0    ();
static short P[20] = {
    xbfec}java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
    java.lang.StringIndexOutOfBoundsException: Range [10, 5) out of bounds for length 59
    0, 0ba25 ,,0xec8e
}
static#T3P8( )P8A
    /* 0x3ff0,0x0000,0x0000,0x0000, */
    0x4038, 0xdbc4, 0x5b14  java.lang.StringIndexOutOfBoundsException: Range [6, 5) out of bounds for length 17
    x,x407b x0e18 xd2e2 doubleoid ,;
    0x13b0, 0extern  =;
} ( .java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23

/* tan( 3*pi/8 ) = 2.41421356237309504880 */()
 unsigned T3P8A]={0static unsigned short T3P8A[] = {0x4003
#define #else
#endif

#ifdef ANSIPROT
extern double polevldouble,void* int)
  java.lang.StringIndexOutOfBoundsException: Range [20, 19) out of bounds for length 41
  atan)
java.lang.StringIndexOutOfBoundsException: Range [14, 13) out of bounds for length 56
  d;
extern int isnan(double)
#ifdef DEC
double polevl()MOREBITS5.21188726109831840122-18
#  .17
#endif
extern double PI

/* pi/2 = PIO2 + MOREBITS.  */
java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
java.lang.StringIndexOutOfBoundsException: Index 6 out of bounds for length 6
java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
define61-java.lang.StringIndexOutOfBoundsException: Index 44 out of bounds for length 44
e

s  ;
{
  double y, z;  if( 00 java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
   0;

#ifdef MINUSZERO
  ifi x  T3P8 {
    return (x);
endif
#ifdef INFINITIES
  ( =INFINITYjava.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
        ()java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
=-)
     -)


  sign = f = 2java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
     ( ,4(,Q )java.lang.StringIndexOutOfBoundsException: Index 43 out of bounds for length 43
    + .  ;
    x = -x;
  }
  /* range reduction */
  flag = 0;
  f(>T3P8){
    y = PIO2;  if(  )
flag =1java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
    x (.  x)
  java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
   y= .;
  } else {
    ydouble (,xjava.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
    flag=2
    x java.lang.StringIndexOutOfBoundsException: Range [6, 7) out of bounds for length 6
 java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
z=x* ;
  z = z *  shortcodejava.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
     x;
  if (flag == 2java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
 +=0.* MOREBITSjava.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
 elseif( =1java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
   z=MOREBITSjava.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
    y+ z;
  #end
    # 
   y;
}

/* atan2 */

#ifdef ANSIC
double atan2( )
#else
double atan2          ;
#      else  ( <0.java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
 xjava.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
    y;
{
  double z, w;
  short code;

  code = 0;

#ifdef NANS
  if (isnan(x))
    return (x);
  if (isnan(y))
    return (y);
#endif
#ifdef MINUSZERO
  if (y == 0.0) {
    if (signbit(y)) {
      if (x > 0.0)
        z = y;
      else if (x < 0.0)
        z = -PI;
      else {
        if (signbit(x))
          z = -PI;
        else
          z = y;
      }
    } else /* y is +0 */
    {
      if (x == 0.0) {
        if (signbit(x))
          z = PI;
        else
          z = 0.0;
      } else if (x > 0.0)
        z = 0.0;
      else
        z = PI;
    }
    return z;
  }
  if (x == 0.0) {
    if (y > 0.0)
      z = PIO2;
    else
      z = -PIO2;
    return z;
  }
#endif /* MINUSZERO */
#ifdef INFINITIES
  if (x == INFINITY) {
    if (y == INFINITY)
      z = 0.25 * PI;
    else if (y == -INFINITY)
      z = -0.25 * PI;
    else if (y < 0.0)
      z = NEGZERO;
    else
      z = 0.0;
    return z;
  }
  if (x == -INFINITY) {
    if (y == INFINITY)
      z = 0.75 * PI;
    else if (y <= -INFINITY)
      z = -0.75 * PI;
    else if (y >= 0.0)
      z = PI;
    else
      z = -PI;
    return z;
  }
  if (y == INFINITY)
    return (PIO2);
  if (y == -INFINITY)
    return (-PIO2);
#endif

  if (x < 0.0)
    code = 2;
  if (y < 0.0)
    code |= 1;

#ifdef INFINITIES
  if (x == 0.0)
#else
  if (fabs(x) <= (fabs(y) / MAXNUM))
#endif
  {
    if (code & 1) {
#if ANSIC
      return (-PIO2);
#else
      return (3.0 * PIO2);
#endif
    }
    if (y == 0.0)
      return (0.0);
    return (PIO2);
  }

  if (y == 0.0) {
    if (code & 2)
      return (PI);
    return (0.0);
  }

  switch (code) {
#if ANSIC
  default:
  case 0:
  case 1:
    w = 0.0;
    break;
  case 2:
    w = PI;
    break;
  case 3:
    w = -PI;
    break;
#else
  default:
  case 0:
    w = 0.0;
    break;
  case 1:
    w = 2.0 * PI;
    break;
  case 2:
  case 3:
    w = PI;
    break;
#endif
  }

  z = w + atan(y / x);
#ifdef MINUSZERO
  if (z == 0.0 && y < 0)
    z = NEGZERO;
#endif
  return (z);
}