overby.me
this repo has no description
1/*
2 * Copyright (c) 2002 Apple Computer, Inc. All rights reserved.
3 *
4 * @APPLE_LICENSE_HEADER_START@
5 *
6 * The contents of this file constitute Original Code as defined in and
7 * are subject to the Apple Public Source License Version 1.1 (the
8 * "License"). You may not use this file except in compliance with the
9 * License. Please obtain a copy of the License at
10 * http://www.apple.com/publicsource and read it before using this file.
11 *
12 * This Original Code and all software distributed under the License are
13 * distributed on an "AS IS" basis, WITHOUT WARRANTY OF ANY KIND, EITHER
14 * EXPRESS OR IMPLIED, AND APPLE HEREBY DISCLAIMS ALL SUCH WARRANTIES,
15 * INCLUDING WITHOUT LIMITATION, ANY WARRANTIES OF MERCHANTABILITY,
16 * FITNESS FOR A PARTICULAR PURPOSE OR NON-INFRINGEMENT. Please see the
17 * License for the specific language governing rights and limitations
18 * under the License.
19 *
20 * @APPLE_LICENSE_HEADER_END@
21 */
22/*******************************************************************************
23* *
24* File: gamma9.c *
25* *
26* Contains: Legacy gamma(). *
27* *
28* Copyright � 2001 Apple Computer, Inc. All rights reserved. *
29* *
30* Written by Stephen C. Peters, started in November 2001. *
31* *
32* A MathLib v5 file. *
33* *
34* Change History (most recent first): *
35* *
36* 12 Dec 01 scp First created. *
37* *
38* W A R N I N G: *
39* These routines require a 64-bit double precision IEEE-754 model. *
40* *
41* A version of gcc higher than 932 is required. *
42* *
43* GCC compiler options: *
44* optimization level 3 (-O3) *
45* -fschedule-insns -finline-functions -funroll-all-loops *
46* *
47*******************************************************************************/
48#ifndef _REENTRANT
49 #define _REENTRANT 1
50#endif
51
52#include "math.h"
53
54int signgam;
55
56double gamma ( double x )
57{
58 double g = tgamma ( x ); // return *True* gamma a la MacOS9 Mathlib
59
60 signgam = (g < 0.0 ? -1 : 1); // set signgam as a courtesy.
61
62 return g;
63}
64
65double gamma_r ( double x, int *psigngam )
66{
67 double g = tgamma ( x ); // return *True* gamma a la MacOS9 Mathlib
68
69 *psigngam = (g < 0.0 ? -1 : 1);
70
71 return g;
72}