mirror of
https://github.com/freebsd/freebsd-src.git
synced 2026-10-30 20:35:03 +01:00
Dimitry Andric
3d1716c561
cacos(3): correct spelling of 'I' In some cases we had 'i' instead of 'I'. PR: 195517 Submitted by: stephen MFC r329259 (by eadler): msun: signed overflow in atan2 As a component of atan2(y, x), the case of x == 1.0 is farmed out to atan(y). The current implementation of this comparison is vulnerable to signed integer underflow (that is, undefined behavior), and it's performed in a somewhat more complicated way than it need be. Change it to not be quite so cute, rather directly comparing the high/low bits of x to the specific IEEE-754 bit pattern that encodes 1.0. Note that while there are three different e_atan* files in the relevant directory, only this one needs fixing. e_atan2f.c already compares against the full bit pattern encoding 1.0f, while e_atan2l.cuses bitwise-ands/ors/nots and so doesn't require a change. Closes #130 Submitted by: Jeff Walden (@jswalden github PR #130) Reviewed by: bde MFC r334721 (by cem): clog.3, complex.3: Fix typos and igor style issues PR: 228783 Reported by: Karsten <freebsd-bugzilla AT kkoenig.net> MFC r336299 (by mmacy): msun: add ld80/ld128 powl, cpow, cpowf, cpowl from openbsd This corresponds to the latest status (hasn't changed in 9+ years) from openbsd of ld80/ld128 powl, and source cpowf, cpow, cpowl (the complex power functions for float complex, double complex, and long double complex) which are required for C99 compliance and were missing from FreeBSD. Also required for some numerical codes using complex numbered Hamiltonians. Thanks to jhb for tracking down the issue with making weak_reference compile on powerpc. When asked to review, bde said "I don't like it" - but provided no actionable feedback or superior implementations. Discussed with: jhb Submitted by: jmd Differential Revision: https://reviews.freebsd.org/D15919 MFC r336563: Recommit r336497: Fix powl, cpow, cpowf, and cpowl imports from OpenBSD This is a follow-up to r336299. * lib/msun/Makefile: . Remove polevll.c * lib/msun/ld80/e_powl.c: . Copy contents of polevll.c to here. This is the only consumer of these functions. Make functions 'static inline'. . Make reducl a 'static inline' function. * lib/msun/man/exp.3: . Remove BUGS section that no longer applies. * lib/msun/src/math_private.h: . Remove prototypes of __p1evll() and __polevll() * lib/msun/src/s_cpow.c: * lib/msun/src/s_cpowf.c: * lib/msun/src/s_cpowl.c . Include math_private.h. . Use the CMPLX macro from either C99 or math_private.h (depends on compiler support) instead of the problematic use of complex I. Submitted by: Steve Kargl <[email protected]> PR: 229876
125 lines
4.3 KiB
C
125 lines
4.3 KiB
C
/*-
|
|
* Copyright (c) 2001-2011 The FreeBSD Project.
|
|
* All rights reserved.
|
|
*
|
|
* Redistribution and use in source and binary forms, with or without
|
|
* modification, are permitted provided that the following conditions
|
|
* are met:
|
|
* 1. Redistributions of source code must retain the above copyright
|
|
* notice, this list of conditions and the following disclaimer.
|
|
* 2. Redistributions in binary form must reproduce the above copyright
|
|
* notice, this list of conditions and the following disclaimer in the
|
|
* documentation and/or other materials provided with the distribution.
|
|
*
|
|
* THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
|
|
* ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
* ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
|
|
* FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
|
|
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
|
|
* OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
|
|
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
|
|
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
|
|
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
|
|
* SUCH DAMAGE.
|
|
*
|
|
* $FreeBSD$
|
|
*/
|
|
|
|
#ifndef _COMPLEX_H
|
|
#define _COMPLEX_H
|
|
|
|
#include <sys/cdefs.h>
|
|
|
|
#ifdef __GNUC__
|
|
#if __STDC_VERSION__ < 199901
|
|
#define _Complex __complex__
|
|
#endif
|
|
#define _Complex_I ((float _Complex)1.0i)
|
|
#endif
|
|
|
|
#ifdef __generic
|
|
_Static_assert(__generic(_Complex_I, float _Complex, 1, 0),
|
|
"_Complex_I must be of type float _Complex");
|
|
#endif
|
|
|
|
#define complex _Complex
|
|
#define I _Complex_I
|
|
|
|
#if __ISO_C_VISIBLE >= 2011
|
|
#ifdef __clang__
|
|
#define CMPLX(x, y) ((double complex){ x, y })
|
|
#define CMPLXF(x, y) ((float complex){ x, y })
|
|
#define CMPLXL(x, y) ((long double complex){ x, y })
|
|
#elif __GNUC_PREREQ__(4, 7)
|
|
#define CMPLX(x, y) __builtin_complex((double)(x), (double)(y))
|
|
#define CMPLXF(x, y) __builtin_complex((float)(x), (float)(y))
|
|
#define CMPLXL(x, y) __builtin_complex((long double)(x), (long double)(y))
|
|
#endif
|
|
#endif /* __ISO_C_VISIBLE >= 2011 */
|
|
|
|
__BEGIN_DECLS
|
|
|
|
double cabs(double complex);
|
|
float cabsf(float complex);
|
|
long double cabsl(long double complex);
|
|
double complex cacos(double complex);
|
|
float complex cacosf(float complex);
|
|
double complex cacosh(double complex);
|
|
float complex cacoshf(float complex);
|
|
double carg(double complex);
|
|
float cargf(float complex);
|
|
long double cargl(long double complex);
|
|
double complex casin(double complex);
|
|
float complex casinf(float complex);
|
|
double complex casinh(double complex);
|
|
float complex casinhf(float complex);
|
|
double complex catan(double complex);
|
|
float complex catanf(float complex);
|
|
double complex catanh(double complex);
|
|
float complex catanhf(float complex);
|
|
double complex ccos(double complex);
|
|
float complex ccosf(float complex);
|
|
double complex ccosh(double complex);
|
|
float complex ccoshf(float complex);
|
|
double complex cexp(double complex);
|
|
float complex cexpf(float complex);
|
|
double cimag(double complex) __pure2;
|
|
float cimagf(float complex) __pure2;
|
|
long double cimagl(long double complex) __pure2;
|
|
double complex clog(double complex);
|
|
float complex clogf(float complex);
|
|
long double complex
|
|
clogl(long double complex);
|
|
double complex conj(double complex) __pure2;
|
|
float complex conjf(float complex) __pure2;
|
|
long double complex
|
|
conjl(long double complex) __pure2;
|
|
float complex cpowf(float complex, float complex);
|
|
double complex cpow(double complex, double complex);
|
|
long double complex
|
|
cpowl(long double complex, long double complex);
|
|
float complex cprojf(float complex) __pure2;
|
|
double complex cproj(double complex) __pure2;
|
|
long double complex
|
|
cprojl(long double complex) __pure2;
|
|
double creal(double complex) __pure2;
|
|
float crealf(float complex) __pure2;
|
|
long double creall(long double complex) __pure2;
|
|
double complex csin(double complex);
|
|
float complex csinf(float complex);
|
|
double complex csinh(double complex);
|
|
float complex csinhf(float complex);
|
|
double complex csqrt(double complex);
|
|
float complex csqrtf(float complex);
|
|
long double complex
|
|
csqrtl(long double complex);
|
|
double complex ctan(double complex);
|
|
float complex ctanf(float complex);
|
|
double complex ctanh(double complex);
|
|
float complex ctanhf(float complex);
|
|
|
|
__END_DECLS
|
|
|
|
#endif /* _COMPLEX_H */
|