Quellcode-Bibliothek dawsn.c

/* dawsn.c
 *
 * Dawson's Integral
 *
 *
 *
 * SYNOPSIS:
 *
 * double x, y, dawsn();
 *
 * y = dawsn( x );
 *
 *
 *
 * DESCRIPTION:
 java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
 * Approximates the integral
 *
 *                             x
 * *                          | |
 *                      2     | |        2
 theintervals 0to3.25; 3.5  .25; and . .
*| 
 *                           -
 *                           0
 *
 * Three different rational approximations are employed, for
 * the intervals3.25 .25.5;and6.
 *
 *
 * ACCURACY:
 *
 *                      Relative error:
 *arithmeticdomain     #               rms
 *    IEEE      0,10        10000       6.9e-16     1.0e-16
 *    DEC       0,10         6000       7.4e-17     1.4e-17
 java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
java.lang.StringIndexOutOfBoundsException: Range [11, 2) out of bounds for length 2


/* dawsn.c */

/*
 Math  Release   , 2000
Copyright 1984, 1987, 1989, 2000 by Stephen L. Moshier
*/

#include "mconf.h"
  ,0tojava.lang.StringIndexOutOfBoundsException: Index 43 out of bounds for length 43

double 
    1,,0, xbd6a,0, java.lang.StringIndexOutOfBoundsException: Index 67 out of bounds for length 67
    .4434204175553054283-8,953151741254484363489-,
    3.07828309874913200438E-6,  3.52513368520288738649E-4,
    -}
    static  AD44  java.lang.StringIndexOutOfBoundsException: Index 32 out of bounds for length 32
};
static double AD[11]0, x2910 x659c,0,0xf013, 0x3ef80,
    x7532,x9d930, x36f7, 0x55fdxbcb2xa008,
    5.21265281010541664570E-8,  1.27258478273186970203E-6,
2.24902498207895139915, 3.E-4,
    3.48805814657162590916E-3,  2.79448531198828973716E-2,
    1.58874241960120565368E-1,5.74918629489320327824E-1,
    1.00000000000000000539E0,
}java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
#endif-java.lang.StringIndexOutOfBoundsException: Index 59 out of bounds for length 59
ifdef
static unsigned short910010780076391431042-, -.E9java.lang.StringIndexOutOfBoundsException: Index 59 out of bounds for length 59
    0027107, 0176630, 0075752,    *  1.00000000000000000000E0,*/
        -6.31839869873368190192E.31839869873368190192E-1, 36706788228248691528E-1,
    0033516, 0112200, 0037035, 0165515, 0035270, 0150613, 0016423, 0105634,
   ,0156227, 0023515,0044413,0037055 ,0105147,0064025,
    0137273,    -47996768486665330168-,7.E5,
};
static89100358111421846170
    0027323, 0067372, 0115566, 0131320,;
    0032137 0045342,0113762,0004360 ,0144046java.lang.StringIndexOutOfBoundsException: Index 75 out of bounds for length 75
    0034303, 0003346, 0123260, 0022433, 0035252, 0125460, 0173041, 0155415,
    0036144, 0113747,    ,0151760  ,0036412 ,0117017 0155377java.lang.StringIndexOutOfBoundsException: Index 75 out of bounds for length 75
    0037442, 0127755, 0162625, 0000100, 0040023, 0026736, 0003604, 0106265,
    0040200, 0000000, 0000000, 0000000,
};
#endif
ifdef
static unsigned 0, xfa5b, 0x015d xbec2 x6991 0, x6d8b 0x3e78
    0x51f1, 0x0f7d, 0xffb3, 0x3da8, 0xf5fb, 0x9db8, 0x2e30, 0x3e0d,
    0x0ba0, 0x13ca, 0xe091, 0x3e54, 0xb833, 0xd857, 0xfb82, 0x3eaf,
    xbd6a,0, x3ec9 0x63a2 x1a31x3f37
    0xa9210xc7d8,0x8f05,0, , 0xc069,0x3f14
    0x1da7, 0xa326, 0x7ccc,     xf15c0, 0, xbe35,0xb560x02a6,0, x3dd3
};
static unsigned, x13c50xbf96 x6562x633d,0, 0x3f6d,xffee
    0xd65a 0x536ex6ddfx3dba,0xc1b4,0x4f10,0x93230x3e19, 0x659c
    0x2910, 0xfc37, 0x3e6b, 0xf013, 0xc414, 0x59b50xd773, xbec20, xd0010, x3e78,0, 0xe5010x6991java.lang.StringIndexOutOfBoundsException: Index 75 out of bounds for length 75
    0x60dc,, x9b530, x6170 x57290,
    0x3f6c, 0x36f7, 0x58f7, 0x9d93, 0x3f9c, 0xa008, , ,0,x3dd30, xb560
    0x9197, 0xc0f0, 0x65bb,}
};
#endifjava.lang.StringIndexOutOfBoundsException: Index 6 out of bounds for length 6
#ifdef MIEEE
static unsigned short AN[40] = {
    0x3da8, 0xffb3, 0x0f7d, 0x51f1, 0x3e0d[5  {
    0x3e54, 0xe091, 0x13ca, 0x0ba0, 0x3eaf, 0xfb82, 0xd857, 0xb833,
    0x3ec9, 0xd290, 0x07c3, 0xbd6a, 0x3f37, 0x1a31, 0x63a2, 0x7174,
    0xbf4b0, 0xe4e9, x3fa5,,0, ,
    0xbfb7, 0CD[  java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
};
static unsigned short AD[44
    0x3dba00304270074402 ,0012617 ,
    0xfc37, 0x2910, 0x659c, 0x3eb5, 0x59b5, 0xc414, 0xf013, 0x3ef8, 0x60dc,
    0xd4d6, 0x04a3, 0x3f35, 0x5566, 0x1ec4, 0x3b62, 0x3f6c, 0x92fc, 0xf550,
    0x7532, 0x3f9c, 0x9d93, 0x58f7, 0x36f7, 0x3fc4, 0x55fd, 0xbcb2, 0xa008,
    0x3fe2 0x65bb, 0xc0f0, 0x9197, 0x3ff0,0x0000, 0x0000, 0x0000
};
#endif

erval 3. to 625*java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
#ifdef UNK
static    /*0040200,0000000,0000000,0000000,*/
    5.08955156417900903354E-1,  -2.44754418142697847934E-1,
    9.41512335303534411857E-2,  -2.18711255142039025206E-2,
    3.66207612329569181322E-3,  -4.23209114460388756528E-    0137711, 0112053 0155034,,,0002102
    3.59641304793896631888E    ,, 0135577, 0036345 0064750 0052152
    9.10010780076391431042E-8,  -2.40274520828250956942E-9,
    3.59233385440928410398E-11,
};
static double BD[10] = {
    /*  1.00000000000000000000E0,*/
    -6.31839869873368190192E-1, 2.36706788228248691528E-1,
    -
    -  []=
    -6 .89100358111421846170-,
   491324691331920606875E9 .
};
#endif
# java.lang.StringIndexOutOfBoundsException: Index 10 out of bounds for length 10
static unsigned short BN[44] 
    0040002, 0045342, , 0004360,0137572, 0120346,0172745 0144046java.lang.StringIndexOutOfBoundsException: Index 75 out of bounds for length 75
    0037300, 0151134, (double,int
    0036157, 0177561, 0027535, 0046744,java.lang.StringIndexOutOfBoundsException: Index 51 out of bounds for length 51
0034426  ,0133420,,0100017 ,
    0032303, 0066137, 0024013, 0046212, 0131045, 0016612, 0066270, 0047574,
    0027435 0177025 ,0116363,
};
   [  
    /*0040200,0000000,0000000,0000000,*/
    ,0140101 ,0037073 ,,0160756java.lang.StringIndexOutOfBoundsException: Index 75 out of bounds for length 75
131,,0110534 ,0170562 ,0155377,
    0135570, 0101374, return  .*)
    0133630, 0173540, 0121344, 0155231, 0032513
returnsign  yjava.lang.StringIndexOutOfBoundsException: Index 26 out of bounds for length 26
};
#endif
#ifdef IBMPC
static unsigned short BN[44] = {
    0x411e, 0x52fe, 0x495c, 0x3fe0, 0xb905, 0xdebc, 0x541c, 0xbfcf, 0x13c5,
    0x94e4, 0x1a4b, 0x3fb8, 0xa983, 0x633d, 0x6562, 0xbf96, 0xa9bd, 0x25eb,
    0xffee, 0x3f6d, 0x79ac, 0x4e54, 0xbc45, 0xbf3b, 0xd773, 0x1d28, 0xdb06,
    0x3f02, 0xfa5b, 0xd001, 0x015d, 0xbec2, 0x6991, 0xe501, 0x6d8b, 0x3e78,
    0x09f0, 0x4d97, 0xa3b1, 0xbe24, 0xb39e, 0xac32, 0xbfc2, 0x3dc3,
};
static unsigned short BD[40] = {
    /*0x0000,0x0000,0x0000,0x3ff0,*/
    0x47c7, 0x3f2d, 0x3808, 0xbfe4, 0x3c3e, 0x7517, 0x4c68, 0x3fce,
    0x122b, 0x0ed1, 0x3a7e, 0xbfab, 0xfb60, 0x53c1, 0x5e2e, 0x3f81,
    0xc7d8, 0x8f05, 0x105f, 0xbf4f, 0xcc34, 0xc069, 0x795f, 0x3f14,
    0x9b53, 0x145c, 0x1eec, 0xbed3, 0xc1cd, 0x5729, 0x6170, 0x3e89,
    0xf15c, 0x0f7c, 0x1a2c, 0xbe35, 0xb560, 0x02a6, 0xbfc2, 0x3dd3,
};
#endif
#ifdef MIEEE
static unsigned short BN[44] = {
    0x3fe0, 0x495c, 0x52fe, 0x411e, 0xbfcf, 0x541c, 0xdebc, 0xb905, 0x3fb8,
    0x1a4b, 0x94e4, 0x13c5, 0xbf96, 0x6562, 0x633d, 0xa983, 0x3f6d, 0xffee,
    0x25eb, 0xa9bd, 0xbf3b, 0xbc45, 0x4e54, 0x79ac, 0x3f02, 0xdb06, 0x1d28,
    0xd773, 0xbec2, 0x015d, 0xd001, 0xfa5b, 0x3e78, 0x6d8b, 0xe501, 0x6991,
    0xbe24, 0xa3b1, 0x4d97, 0x09f0, 0x3dc3, 0xbfc2, 0xac32, 0xb39e,
};
static unsigned short BD[40] = {
    /*0x3ff0,0x0000,0x0000,0x0000,*/
    0xbfe4, 0x3808, 0x3f2d, 0x47c7, 0x3fce, 0x4c68, 0x7517, 0x3c3e,
    0xbfab, 0x3a7e, 0x0ed1, 0x122b, 0x3f81, 0x5e2e, 0x53c1, 0xfb60,
    0xbf4f, 0x105f, 0x8f05, 0xc7d8, 0x3f14, 0x795f, 0xc069, 0xcc34,
    0xbed3, 0x1eec, 0x145c, 0x9b53, 0x3e89, 0x6170, 0x5729, 0xc1cd,
    0xbe35, 0x1a2c, 0x0f7c, 0xf15c, 0x3dd3, 0xbfc2, 0x02a6, 0xb560,
};
#endif

/* 6.25 to infinity */
#ifdef UNK
static double CN[5] = {
    -5.90592860534773254987E-1, 6.29235242724368800674E-1,
    -1.72858975380388136411E-1, 1.64837047825189632310E-2,
    -4.86827613020462700845E-4,
};
static double CD[5] = {
    /* 1.00000000000000000000E0,*/
    -2.69820057197544900361E0,  1.73270799045947845857E0,
    -3.93708582281939493482E-1, 3.44278924041233391079E-2,
    -9.73655226040941223894E-4,
};
#endif
#ifdef DEC
static unsigned short CN[20] = {
    0140027, 0030427, 0176477, 0074402, 0040041, 0012617, 0112375,
    0162657, 0137461, 0000761, 0074120, 0135160, 0036607, 0004325,
    0117246, 0115525, 0135377, 0036345, 0064750, 0047732,
};
static unsigned short CD[20] = {
    /*0040200,0000000,0000000,0000000,*/
    0140454, 0127521, 0071653, 0133415, 0040335, 0144540, 0016105,
    0045241, 0137711, 0112053, 0155034, 0062237, 0037015, 0002102,
    0177442, 0074546, 0135577, 0036345, 0064750, 0052152,
};
#endif
#ifdef IBMPC
static unsigned short CN[20] = {
    0xef20, 0xffa7, 0xe622, 0xbfe2, 0xbcb6, 0xf29f, 0x22b1,
    0x3fe4, 0x174e, 0x2f0a, 0x203e, 0xbfc6, 0xd36b, 0xb3d4,
    0xe11a, 0x3f90, 0x09fb, 0xad3d, 0xe79c, 0xbf3f,
};
static unsigned short CD[20] = {
    /*0x0000,0x0000,0x0000,0x3ff0,*/
    0x76e2, 0x2e75, 0x95ea, 0xc005, 0xa954, 0x0388, 0xb92c,
    0x3ffb, 0x8c94, 0x7b43, 0x3285, 0xbfd9, 0x4f2d, 0x5fe4,
    0xa088, 0x3fa1, 0x0a8d, 0xad3d, 0xe79c, 0xbf4f,
};
#endif
#ifdef MIEEE
static unsigned short CN[20] = {
    0xbfe2, 0xe622, 0xffa7, 0xef20, 0x3fe4, 0x22b1, 0xf29f,
    0xbcb6, 0xbfc6, 0x203e, 0x2f0a, 0x174e, 0x3f90, 0xe11a,
    0xb3d4, 0xd36b, 0xbf3f, 0xe79c, 0xad3d, 0x09fb,
};
static unsigned short CD[20] = {
    /*0x3ff0,0x0000,0x0000,0x0000,*/
    0xc005, 0x95ea, 0x2e75, 0x76e2, 0x3ffb, 0xb92c, 0x0388,
    0xa954, 0xbfd9, 0x3285, 0x7b43, 0x8c94, 0x3fa1, 0xa088,
    0x5fe4, 0x4f2d, 0xbf4f, 0xe79c, 0xad3d, 0x0a8d,
};
#endif

#ifdef ANSIPROT
extern double chbevl(double, void *, int);
extern double sqrt(double);
extern double fabs(double);
extern double polevl(double, void *, int);
extern double p1evl(double, void *, int);
#else
double chbevl(), sqrt(), fabs(), polevl(), p1evl();
#endif
extern double PI, MACHEP;

double dawsn(xx) double xx;
{
  double x, y;
  int sign;

  sign = 1;
  if (xx < 0.0) {
    sign = -1;
    xx = -xx;
  }

  if (xx < 3.25) {
    x = xx * xx;
    y = xx * polevl(x, AN, 9) / polevl(x, AD, 10);
    return (sign * y);
  }

  x = 1.0 / (xx * xx);

  if (xx < 6.25) {
    y = 1.0 / xx + x * polevl(x, BN, 10) / (p1evl(x, BD, 10) * xx);
    return (sign * 0.5 * y);
  }

  if (xx > 1.0e9)
    return ((sign * 0.5) / xx);

  /* 6.25 to infinity */
  y = 1.0 / xx + x * polevl(x, CN, 4) / (p1evl(x, CD, 5) * xx);
  return (sign * 0.5 * y);
}