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 tgamma = __tgamma
31
32 /* INDENT OFF */
33 /*
34 * True gamma function
35 * double tgamma(double x)
36 *
37 * Error:
38 * ------
39 * Less that one ulp for both positive and negative arguments.
40 *
41 * Algorithm:
42 * ---------
43 * A: For negative argument
44 * (1) gamma(-n or -inf) is NaN
45 * (2) Underflow Threshold
46 * (3) Reduction to gamma(1+x)
47 * B: For x between 1 and 2
48 * C: For x between 0 and 1
49 * D: For x between 2 and 8
50 * E: Overflow thresold {see over.c}
|
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 __tgamma = tgamma
31
32 /* INDENT OFF */
33 /*
34 * True gamma function
35 * double tgamma(double x)
36 *
37 * Error:
38 * ------
39 * Less that one ulp for both positive and negative arguments.
40 *
41 * Algorithm:
42 * ---------
43 * A: For negative argument
44 * (1) gamma(-n or -inf) is NaN
45 * (2) Underflow Threshold
46 * (3) Reduction to gamma(1+x)
47 * B: For x between 1 and 2
48 * C: For x between 0 and 1
49 * D: For x between 2 and 8
50 * E: Overflow thresold {see over.c}
|