| 1 | /*- |
| 2 | * Copyright (c) 2008 David Schultz <das@FreeBSD.ORG> |
| 3 | * All rights reserved. |
| 4 | * |
| 5 | * Redistribution and use in source and binary forms, with or without |
| 6 | * modification, are permitted provided that the following conditions |
| 7 | * are met: |
| 8 | * 1. Redistributions of source code must retain the above copyright |
| 9 | * notice, this list of conditions and the following disclaimer. |
| 10 | * 2. Redistributions in binary form must reproduce the above copyright |
| 11 | * notice, this list of conditions and the following disclaimer in the |
| 12 | * documentation and/or other materials provided with the distribution. |
| 13 | * |
| 14 | * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND |
| 15 | * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE |
| 16 | * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE |
| 17 | * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE |
| 18 | * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL |
| 19 | * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS |
| 20 | * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) |
| 21 | * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT |
| 22 | * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY |
| 23 | * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF |
| 24 | * SUCH DAMAGE. |
| 25 | * |
| 26 | * $FreeBSD: head/lib/msun/ld80/invtrig.c 181074 2008-07-31 22:41:26Z das $ |
| 27 | */ |
| 28 | |
| 29 | #include "invtrig.h" |
| 30 | |
| 31 | /* |
| 32 | * asinl() and acosl() |
| 33 | */ |
| 34 | const long double |
| 35 | pS0 = 1.66666666666666666631e-01L, |
| 36 | pS1 = -4.16313987993683104320e-01L, |
| 37 | pS2 = 3.69068046323246813704e-01L, |
| 38 | pS3 = -1.36213932016738603108e-01L, |
| 39 | pS4 = 1.78324189708471965733e-02L, |
| 40 | pS5 = -2.19216428382605211588e-04L, |
| 41 | pS6 = -7.10526623669075243183e-06L, |
| 42 | qS1 = -2.94788392796209867269e+00L, |
| 43 | qS2 = 3.27309890266528636716e+00L, |
| 44 | qS3 = -1.68285799854822427013e+00L, |
| 45 | qS4 = 3.90699412641738801874e-01L, |
| 46 | qS5 = -3.14365703596053263322e-02L; |
| 47 | |
| 48 | /* |
| 49 | * atanl() |
| 50 | */ |
| 51 | const long double atanhi[] = { |
| 52 | 4.63647609000806116202e-01L, |
| 53 | 7.85398163397448309628e-01L, |
| 54 | 9.82793723247329067960e-01L, |
| 55 | 1.57079632679489661926e+00L, |
| 56 | }; |
| 57 | |
| 58 | const long double atanlo[] = { |
| 59 | 1.18469937025062860669e-20L, |
| 60 | -1.25413940316708300586e-20L, |
| 61 | 2.55232234165405176172e-20L, |
| 62 | -2.50827880633416601173e-20L, |
| 63 | }; |
| 64 | |
| 65 | const long double aT[] = { |
| 66 | 3.33333333333333333017e-01L, |
| 67 | -1.99999999999999632011e-01L, |
| 68 | 1.42857142857046531280e-01L, |
| 69 | -1.11111111100562372733e-01L, |
| 70 | 9.09090902935647302252e-02L, |
| 71 | -7.69230552476207730353e-02L, |
| 72 | 6.66661718042406260546e-02L, |
| 73 | -5.88158892835030888692e-02L, |
| 74 | 5.25499891539726639379e-02L, |
| 75 | -4.70119845393155721494e-02L, |
| 76 | 4.03539201366454414072e-02L, |
| 77 | -2.91303858419364158725e-02L, |
| 78 | 1.24822046299269234080e-02L, |
| 79 | }; |
| 80 | |
| 81 | const long double pi_lo = -5.01655761266833202345e-20L; |