10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  19  * CDDL HEADER END
  20  */
  21 
  22 /*
  23  * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
  24  */
  25 /*
  26  * Copyright 2006 Sun Microsystems, Inc.  All rights reserved.
  27  * Use is subject to license terms.
  28  */
  29 
  30 #pragma weak tgammaf = __tgammaf
  31 
  32 /*
  33  * True gamma function
  34  *
  35  * float tgammaf(float x)
  36  *
  37  * Algorithm: see tgamma.c
  38  *
  39  * Maximum error observed: 0.87ulp (both positive and negative arguments)
  40  */
  41 
  42 #include "libm.h"
  43 #include "libm_synonyms.h"
  44 #include <math.h>
  45 #if defined(__SUNPRO_C)
  46 #include <sunmath.h>
  47 #endif
  48 #include <sys/isa_defs.h>
  49 
  50 #if defined(_BIG_ENDIAN)
  51 #define HIWORD  0
  52 #define LOWORD  1
  53 #else
  54 #define HIWORD  1
  55 #define LOWORD  0
  56 #endif
  57 #define __HI(x) ((int *) &x)[HIWORD]
  58 #define __LO(x) ((unsigned *) &x)[LOWORD]
  59 
  60 /* Coefficients for primary intervals GTi() */
  61 static const double cr[] = {
  62         /* p1 */
  63         +7.09087253435088360271451613398019280077561279443e-0001,

  10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  19  * CDDL HEADER END
  20  */
  21 
  22 /*
  23  * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
  24  */
  25 /*
  26  * Copyright 2006 Sun Microsystems, Inc.  All rights reserved.
  27  * Use is subject to license terms.
  28  */
  29 
  30 #pragma weak __tgammaf = tgammaf
  31 
  32 /*
  33  * True gamma function
  34  *
  35  * float tgammaf(float x)
  36  *
  37  * Algorithm: see tgamma.c
  38  *
  39  * Maximum error observed: 0.87ulp (both positive and negative arguments)
  40  */
  41 
  42 #include "libm.h"

  43 #include <math.h>
  44 #if defined(__SUNPRO_C)
  45 #include <sunmath.h>
  46 #endif
  47 #include <sys/isa_defs.h>
  48 
  49 #if defined(_BIG_ENDIAN)
  50 #define HIWORD  0
  51 #define LOWORD  1
  52 #else
  53 #define HIWORD  1
  54 #define LOWORD  0
  55 #endif
  56 #define __HI(x) ((int *) &x)[HIWORD]
  57 #define __LO(x) ((unsigned *) &x)[LOWORD]
  58 
  59 /* Coefficients for primary intervals GTi() */
  60 static const double cr[] = {
  61         /* p1 */
  62         +7.09087253435088360271451613398019280077561279443e-0001,