1 

     2 /*

     3  * Copyright (c) 1998, 2001, Oracle and/or its affiliates. All rights reserved.

     4  * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.

     5  *

     6  * This code is free software; you can redistribute it and/or modify it

     7  * under the terms of the GNU General Public License version 2 only, as

     8  * published by the Free Software Foundation.  Oracle designates this

     9  * particular file as subject to the "Classpath" exception as provided

    10  * by Oracle in the LICENSE file that accompanied this code.

    11  *

    12  * This code is distributed in the hope that it will be useful, but WITHOUT

    13  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

    14  * FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License

    15  * version 2 for more details (a copy is included in the LICENSE file that

    16  * accompanied this code).

    17  *

    18  * You should have received a copy of the GNU General Public License version

    19  * 2 along with this work; if not, write to the Free Software Foundation,

    20  * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.

    21  *

    22  * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA

    23  * or visit www.oracle.com if you need additional information or have any

    24  * questions.

    25  */

    26 

    27 /*

    28  * wrapper of j1,y1

    29  */

    30 

    31 #include "fdlibm.h"

    32 

    33 #ifdef __STDC__

    34         double j1(double x)             /* wrapper j1 */

    35 #else

    36         double j1(x)                    /* wrapper j1 */

    37         double x;

    38 #endif

    39 {

    40 #ifdef _IEEE_LIBM

    41         return __ieee754_j1(x);

    42 #else

    43         double z;

    44         z = __ieee754_j1(x);

    45         if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;

    46         if(fabs(x)>X_TLOSS) {

    47                 return __kernel_standard(x,x,36); /* j1(|x|>X_TLOSS) */

    48         } else

    49             return z;

    50 #endif

    51 }

    52 

    53 #ifdef __STDC__

    54         double y1(double x)             /* wrapper y1 */

    55 #else

    56         double y1(x)                    /* wrapper y1 */

    57         double x;

    58 #endif

    59 {

    60 #ifdef _IEEE_LIBM

    61         return __ieee754_y1(x);

    62 #else

    63         double z;

    64         z = __ieee754_y1(x);

    65         if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;

    66         if(x <= 0.0){

    67                 if(x==0.0)

    68                     /* d= -one/(x-x); */

    69                     return __kernel_standard(x,x,10);

    70                 else

    71                     /* d = zero/(x-x); */

    72                     return __kernel_standard(x,x,11);

    73         }

    74         if(x>X_TLOSS) {

    75                 return __kernel_standard(x,x,37); /* y1(x>X_TLOSS) */

    76         } else

    77             return z;

    78 #endif

    79 }