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 }