Mercurial > mxe-octave

annotate src/mingw-w64-complex-inverse-trig.patch @ 6354:9d1395334e68

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
mingw-w64: Handle negative zero input for cacosh (bug #49091). * src/mingw-w64-complex-inverse-trig.patch: Handle negative zero input for cacosh.
author Markus Mützel <markus.muetzel@gmx.de>
date Fri, 05 Aug 2022 19:32:03 +0200
parents 0825abaf61a7
children 6f3c099c0d38
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6335
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
1 From 59fd6e6ab7c93545175fa6ff46532b7a80dbdbfb Mon Sep 17 00:00:00 2001
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
2 From: =?UTF-8?q?Markus=20M=C3=BCtzel?= <markus.muetzel@gmx.de>
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
3 Date: Sat, 23 Jul 2022 15:30:46 +0200
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
4 Subject: [PATCH 1/2] casinh: Use approximation for large input.
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
5
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
6 ---
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
7 mingw-w64-crt/complex/casinh.def.h | 21 +++++++++++++++++++++
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
8 1 file changed, 21 insertions(+)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
9
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
10 diff --git a/mingw-w64-crt/complex/casinh.def.h b/mingw-w64-crt/complex/casinh.def.h
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
11 index 050d885a0..cce2d864e 100644
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
12 --- a/mingw-w64-crt/complex/casinh.def.h
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
13 +++ b/mingw-w64-crt/complex/casinh.def.h
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
14 @@ -87,6 +87,27 @@ __FLT_ABI(casinh) (__FLT_TYPE __complex__ z)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
15 if (r_class == FP_ZERO && i_class == FP_ZERO)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
16 return z;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
17
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
18 + /* casinh(z) = log(z + sqrt(z*z + 1)) */
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
19 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
20 + if (__FLT_ABI(fabs) (__real__ z) >= __FLT_CST(1.0)/__FLT_EPSILON
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
21 + || __FLT_ABI(fabs) (__imag__ z) >= __FLT_CST(1.0)/__FLT_EPSILON)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
22 + {
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
23 + /* For large z, z + sqrt(z*z + 1) is approximately 2*z.
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
24 + Use that approximation to avoid overflow when squaring.
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
25 + Additionally, use symmetries to perform the calculation in the first
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
26 + quadrant. */
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
27 + __real__ x = __FLT_ABI(fabs) (__real__ z);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
28 + __imag__ x = __FLT_ABI(fabs) (__imag__ z);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
29 + x = __FLT_ABI(clog) (x);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
30 + __real__ x += M_LN2;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
31 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
32 + /* adjust signs for input quadrant */
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
33 + __real__ ret = __FLT_ABI(copysign) (__real__ x, __real__ z);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
34 + __imag__ ret = __FLT_ABI(copysign) (__imag__ x, __imag__ z);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
35 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
36 + return ret;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
37 + }
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
38 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
39 __real__ x = (__real__ z - __imag__ z) * (__real__ z + __imag__ z) + __FLT_CST(1.0);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
40 __imag__ x = __FLT_CST(2.0) * __real__ z * __imag__ z;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
41
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
42 --
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
43 2.35.3.windows.1
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
44
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
45 From 4b5a229573ca98a622cc15d20bd561dbaec9bbfc Mon Sep 17 00:00:00 2001
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
46 From: =?UTF-8?q?Markus=20M=C3=BCtzel?= <markus.muetzel@gmx.de>
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
47 Date: Sat, 23 Jul 2022 18:01:29 +0200
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
48 Subject: [PATCH 2/2] cacosh: Use approximation for large input.
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
49
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
50 ---
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
51 mingw-w64-crt/complex/cacosh.def.h | 21 +++++++++++++++++++++
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
52 1 file changed, 21 insertions(+)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
53
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
54 diff --git a/mingw-w64-crt/complex/cacosh.def.h b/mingw-w64-crt/complex/cacosh.def.h
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
55 index f4ea2da07..a3edd5288 100644
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
56 --- a/mingw-w64-crt/complex/cacosh.def.h
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
57 +++ b/mingw-w64-crt/complex/cacosh.def.h
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
58 @@ -80,6 +80,27 @@ __FLT_ABI(cacosh) (__FLT_TYPE __complex__ z)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
59 return ret;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
60 }
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
61
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
62 + /* cacosh(z) = log(z + sqrt(z*z - 1)) */
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
63 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
64 + if (__FLT_ABI(fabs) (__real__ z) >= __FLT_CST(1.0)/__FLT_EPSILON
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
65 + || __FLT_ABI(fabs) (__imag__ z) >= __FLT_CST(1.0)/__FLT_EPSILON)
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
66 + {
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
67 + /* For large z, z + sqrt(z*z - 1) is approximately 2*z.
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
68 + Use that approximation to avoid overflow when squaring.
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
69 + Additionally, use symmetries to perform the calculation in the positive
6354
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
70 + half plane. */
6335
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
71 + __real__ x = __real__ z;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
72 + __imag__ x = __FLT_ABI(fabs) (__imag__ z);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
73 + x = __FLT_ABI(clog) (x);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
74 + __real__ x += M_LN2;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
75 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
76 + /* adjust signs for input */
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
77 + __real__ ret = __real__ x;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
78 + __imag__ ret = __FLT_ABI(copysign) (__imag__ x, __imag__ z);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
79 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
80 + return ret;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
81 + }
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
82 +
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
83 __real__ x = (__real__ z - __imag__ z) * (__real__ z + __imag__ z) - __FLT_CST(1.0);
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
84 __imag__ x = __FLT_CST(2.0) * __real__ z * __imag__ z;
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
85
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
86 --
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
87 2.35.3.windows.1
0825abaf61a7 mingw-w64: Special case for large input for complex inverse trigonometric functions (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents:
diff changeset
88
6354
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
89 From 9fcac91fc50120a5e6acfa94ba0a917c7faf2306 Mon Sep 17 00:00:00 2001
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
90 From: =?UTF-8?q?Markus=20M=C3=BCtzel?= <markus.muetzel@gmx.de>
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
91 Date: Fri, 5 Aug 2022 18:36:49 +0200
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
92 Subject: [PATCH] cacosh: Handle negative zero input
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
93
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
94 ---
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
95 mingw-w64-crt/complex/cacosh.def.h | 4 ++--
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
96 1 file changed, 2 insertions(+), 2 deletions(-)
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
97
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
98 diff --git a/mingw-w64-crt/complex/cacosh.def.h b/mingw-w64-crt/complex/cacosh.def.h
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
99 index 1f03befae..addf2958e 100644
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
100 --- a/mingw-w64-crt/complex/cacosh.def.h
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
101 +++ b/mingw-w64-crt/complex/cacosh.def.h
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
102 @@ -106,7 +106,7 @@ __FLT_ABI(cacosh) (__FLT_TYPE __complex__ z)
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
103
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
104 x = __FLT_ABI(csqrt) (x);
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
105
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
106 - if (__real__ z < __FLT_CST(0.0))
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
107 + if (__FLT_ABI(signbit) (__real__ z))
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
108 x = -x;
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
109
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
110 __real__ x += __real__ z;
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
111 @@ -114,7 +114,7 @@ __FLT_ABI(cacosh) (__FLT_TYPE __complex__ z)
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
112
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
113 ret = __FLT_ABI(clog) (x);
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
114
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
115 - if (__real__ ret < __FLT_CST(0.0))
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
116 + if (__FLT_ABI(signbit) (__real__ ret))
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
117 ret = -ret;
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
118
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
119 return ret;
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
120 --
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
121 2.35.3.windows.1
9d1395334e68 mingw-w64: Handle negative zero input for cacosh (bug #49091).
Markus Mützel <markus.muetzel@gmx.de>
parents: 6335
diff changeset
122