diff options
author | Ian Moffett <ian@osmora.org> | 2024-03-07 17:28:52 -0500 |
---|---|---|
committer | Ian Moffett <ian@osmora.org> | 2024-03-07 18:24:51 -0500 |
commit | f5e48e94a2f4d4bbd6e5628c7f2afafc6dbcc459 (patch) | |
tree | 93b156621dc0303816b37f60ba88051b702d92f6 /lib/mlibc/options/ansi/generic/complex | |
parent | bd5969fc876a10b18613302db7087ef3c40f18e1 (diff) |
build: Build mlibc + add distclean target
Signed-off-by: Ian Moffett <ian@osmora.org>
Diffstat (limited to 'lib/mlibc/options/ansi/generic/complex')
49 files changed, 0 insertions, 3361 deletions
diff --git a/lib/mlibc/options/ansi/generic/complex/cabs.c b/lib/mlibc/options/ansi/generic/complex/cabs.c deleted file mode 100644 index 2750fab..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cabs.c +++ /dev/null @@ -1,53 +0,0 @@ -/* $NetBSD: cabs.c,v 1.1 2007/08/20 16:01:30 drochner Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cabs>>, <<cabsf>>---complex absolute-value - -INDEX - cabs -INDEX - cabsf - -ANSI_SYNOPSIS - #include <complex.h> - double cabs(double complex <[z]>); - float cabsf(float complex <[z]>); - - -DESCRIPTION - These functions compute compute the complex absolute value - (also called norm, modulus, or magnitude) of <[z]>. - - <<cabsf>> is identical to <<cabs>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The cabs functions return the complex absolute value. - -PORTABILITY - <<cabs>> and <<cabsf>> are ISO C99 - -QUICKREF - <<cabs>> and <<cabsf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -double -cabs(double complex z) -{ - - return hypot( creal(z), cimag(z) ); -} diff --git a/lib/mlibc/options/ansi/generic/complex/cabsf.c b/lib/mlibc/options/ansi/generic/complex/cabsf.c deleted file mode 100644 index 635e23e..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cabsf.c +++ /dev/null @@ -1,19 +0,0 @@ -/* $NetBSD: cabsf.c,v 1.1 2007/08/20 16:01:30 drochner Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float -cabsf(float complex z) -{ - - return hypotf( crealf(z), cimagf(z) ); -} diff --git a/lib/mlibc/options/ansi/generic/complex/cacos.c b/lib/mlibc/options/ansi/generic/complex/cacos.c deleted file mode 100644 index 86e1198..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cacos.c +++ /dev/null @@ -1,99 +0,0 @@ -/* $NetBSD: cacos.c,v 1.1 2007/08/20 16:01:30 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cacos>>, <<cacosf>>---complex arc cosine - -INDEX - cacos -INDEX - cacosf - -ANSI_SYNOPSIS - #include <complex.h> - double complex cacos(double complex <[z]>); - float complex cacosf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex arc cosine of <[z]>, - with branch cuts outside the interval [-1, +1] along the real axis. - - <<cacosf>> is identical to <<cacos>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - These functions return the complex arc cosine value, in the range - of a strip mathematically unbounded along the imaginary axis - and in the interval [0, pi] along the real axis. - @end ifnottex - @tex - These functions return the complex arc cosine value, in the range - of a strip mathematically unbounded along the imaginary axis - and in the interval [<<0>>, $\pi$] along the real axis. - @end tex - -PORTABILITY - <<cacos>> and <<cacosf>> are ISO C99 - -QUICKREF - <<cacos>> and <<cacosf>> are ISO C99 - -*/ - -#include <complex.h> -#include <math.h> - -double complex -cacos(double complex z) -{ - double complex w; - - /* FIXME: The original NetBSD code results in an ICE when trying to - build this function on ARM/Thumb using gcc 4.5.1. For now we use - a hopefully temporary workaround. */ -#if 0 - w = casin(z); - w = (M_PI_2 - creal(w)) - cimag(w) * I; -#else - double complex tmp0, tmp1; - - tmp0 = casin(z); - tmp1 = M_PI_2 - creal(tmp0); - w = tmp1 - (cimag(tmp0) * I); -#endif - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cacosf.c b/lib/mlibc/options/ansi/generic/complex/cacosf.c deleted file mode 100644 index 3874dd5..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cacosf.c +++ /dev/null @@ -1,46 +0,0 @@ -/* $NetBSD: cacosf.c,v 1.1 2007/08/20 16:01:30 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -cacosf(float complex z) -{ - float complex w; - - w = casinf(z); - w = ((float)M_PI_2 - crealf(w)) - cimagf(w) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cacosh.c b/lib/mlibc/options/ansi/generic/complex/cacosh.c deleted file mode 100644 index 3d42c40..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cacosh.c +++ /dev/null @@ -1,93 +0,0 @@ -/* $NetBSD: cacosh.c,v 1.2 2009/08/03 19:41:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cacosh>>, <<cacoshf>>---complex arc hyperbolic cosine - -INDEX - cacosh -INDEX - cacoshf - -ANSI_SYNOPSIS - #include <complex.h> - double complex cacosh(double complex <[z]>); - float complex cacoshf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex arc hyperbolic cosine of <[z]>, - with a branch cut at values less than 1 along the real axis. - - <<cacoshf>> is identical to <<cacosh>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - These functions return the complex arc hyperbolic cosine value, - in the range of a half-strip of non-negative values along the - real axis and in the interval [-i * pi, +i * pi] along the - imaginary axis. - @end ifnottex - @tex - These functions return the complex arc hyperbolic cosine value, - in the range of a half-strip of non-negative values along the - real axis and in the interval [$-i\pi$, $+i\pi$] along the - imaginary axis. - @end tex - -PORTABILITY - <<cacosh>> and <<cacoshf>> are ISO C99 - -QUICKREF - <<cacosh>> and <<cacoshf>> are ISO C99 - -*/ - - -#include <complex.h> - -double complex -cacosh(double complex z) -{ - double complex w; - -#if 0 /* does not give the principal value */ - w = I * cacos(z); -#else - w = clog(z + csqrt(z + 1) * csqrt(z - 1)); -#endif - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cacoshf.c b/lib/mlibc/options/ansi/generic/complex/cacoshf.c deleted file mode 100644 index 41a557a..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cacoshf.c +++ /dev/null @@ -1,48 +0,0 @@ -/* $NetBSD: cacoshf.c,v 1.2 2009/08/03 19:41:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> - -float complex -cacoshf(float complex z) -{ - float complex w; - -#if 0 /* does not give the principal value */ - w = I * cacosf(z); -#else - w = clogf(z + csqrtf(z + 1) * csqrtf(z - 1)); -#endif - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/carg.c b/lib/mlibc/options/ansi/generic/complex/carg.c deleted file mode 100644 index 0447420..0000000 --- a/lib/mlibc/options/ansi/generic/complex/carg.c +++ /dev/null @@ -1,59 +0,0 @@ -/* $NetBSD: carg.c,v 1.1 2007/08/20 16:01:31 drochner Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<carg>>, <<cargf>>---argument (phase angle) - -INDEX - carg -INDEX - cargf - -ANSI_SYNOPSIS - #include <complex.h> - double carg(double complex <[z]>); - float cargf(float complex <[z]>); - - -DESCRIPTION - These functions compute the argument (also called phase angle) - of <[z]>, with a branch cut along the negative real axis. - - <<cargf>> is identical to <<carg>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - The carg functions return the value of the argument in the - interval [-pi, +pi] - @end ifnottex - @tex - The carg functions return the value of the argument in the - interval [$-\pi$, $+\pi$] - @end tex - -PORTABILITY - <<carg>> and <<cargf>> are ISO C99 - -QUICKREF - <<carg>> and <<cargf>> are ISO C99 - -*/ - -#include <complex.h> -#include <math.h> - -double -carg(double complex z) -{ - - return atan2( cimag(z) , creal(z) ); -} diff --git a/lib/mlibc/options/ansi/generic/complex/cargf.c b/lib/mlibc/options/ansi/generic/complex/cargf.c deleted file mode 100644 index 1683d21..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cargf.c +++ /dev/null @@ -1,19 +0,0 @@ -/* $NetBSD: cargf.c,v 1.1 2007/08/20 16:01:31 drochner Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float -cargf(float complex z) -{ - - return atan2f( cimagf(z), crealf(z) ); -} diff --git a/lib/mlibc/options/ansi/generic/complex/casin.c b/lib/mlibc/options/ansi/generic/complex/casin.c deleted file mode 100644 index 5019fd8..0000000 --- a/lib/mlibc/options/ansi/generic/complex/casin.c +++ /dev/null @@ -1,165 +0,0 @@ -/* $NetBSD: casin.c,v 1.1 2007/08/20 16:01:31 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<casin>>, <<casinf>>---complex arc sine - -INDEX - casin -INDEX - casinf - -ANSI_SYNOPSIS - #include <complex.h> - double complex casin(double complex <[z]>); - float complex casinf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex arc sine of <[z]>, - with branch cuts outside the interval [-1, +1] along the real axis. - - <<casinf>> is identical to <<casin>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - These functions return the complex arc sine value, in the range - of a strip mathematically unbounded along the imaginary axis - and in the interval [-pi/2, +pi/2] along the real axis. - @end ifnottex - @tex - These functions return the complex arc sine value, in the range - of a strip mathematically unbounded along the imaginary axis - and in the interval [$-\pi/2$, $+\pi/2$] along the real axis. - @end tex - -PORTABILITY - <<casin>> and <<casinf>> are ISO C99 - -QUICKREF - <<casin>> and <<casinf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -#ifdef __weak_alias -__weak_alias(casin, _casin) -#endif - -double complex -casin(double complex z) -{ - double complex w; - double complex ca, ct, zz, z2; - double x, y; - - x = creal(z); - y = cimag(z); - -#if 0 /* MD: test is incorrect, casin(>1) is defined */ - if (y == 0.0) { - if (fabs(x) > 1.0) { - w = M_PI_2 + 0.0 * I; -#if 0 - mtherr ("casin", DOMAIN); -#endif - } else { - w = asin(x) + 0.0 * I; - } - return w; - } -#endif - -/* Power series expansion */ -/* -b = cabs(z); -if( b < 0.125 ) -{ -z2.r = (x - y) * (x + y); -z2.i = 2.0 * x * y; - -cn = 1.0; -n = 1.0; -ca.r = x; -ca.i = y; -sum.r = x; -sum.i = y; -do - { - ct.r = z2.r * ca.r - z2.i * ca.i; - ct.i = z2.r * ca.i + z2.i * ca.r; - ca.r = ct.r; - ca.i = ct.i; - - cn *= n; - n += 1.0; - cn /= n; - n += 1.0; - b = cn/n; - - ct.r *= b; - ct.i *= b; - sum.r += ct.r; - sum.i += ct.i; - b = fabs(ct.r) + fabs(ct.i); - } -while( b > MACHEP ); -w->r = sum.r; -w->i = sum.i; -return; -} -*/ - - - ca = x + y * I; - ct = ca * I; - /* sqrt( 1 - z*z) */ - /* cmul( &ca, &ca, &zz ) */ - /*x * x - y * y */ - zz = (x - y) * (x + y) + (2.0 * x * y) * I; - - zz = 1.0 - creal(zz) - cimag(zz) * I; - z2 = csqrt(zz); - - zz = ct + z2; - zz = clog(zz); - /* multiply by 1/i = -i */ - w = zz * (-1.0 * I); - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/casinf.c b/lib/mlibc/options/ansi/generic/complex/casinf.c deleted file mode 100644 index 9a9f759..0000000 --- a/lib/mlibc/options/ansi/generic/complex/casinf.c +++ /dev/null @@ -1,122 +0,0 @@ -/* $NetBSD: casinf.c,v 1.1 2007/08/20 16:01:31 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -#ifdef __weak_alias -__weak_alias(casinf, _casinf) -#endif - -float complex -casinf(float complex z) -{ - float complex w; - float complex ca, ct, zz, z2; - float x, y; - - x = crealf(z); - y = cimagf(z); - -#if 0 /* MD: test is incorrect, casin(>1) is defined */ - if (y == 0.0f) { - if (fabsf(x) > 1.0) { - w = M_PI_2 + 0.0f * I; -#if 0 - mtherr ("casin", DOMAIN); -#endif - } else { - w = asinf(x) + 0.0f * I; - } - return w; - } -#endif - -/* Power series expansion */ -/* -b = cabsf(z); -if( b < 0.125 ) -{ -z2.r = (x - y) * (x + y); -z2.i = 2.0 * x * y; - -cn = 1.0; -n = 1.0; -ca.r = x; -ca.i = y; -sum.r = x; -sum.i = y; -do - { - ct.r = z2.r * ca.r - z2.i * ca.i; - ct.i = z2.r * ca.i + z2.i * ca.r; - ca.r = ct.r; - ca.i = ct.i; - - cn *= n; - n += 1.0; - cn /= n; - n += 1.0; - b = cn/n; - - ct.r *= b; - ct.i *= b; - sum.r += ct.r; - sum.i += ct.i; - b = fabsf(ct.r) + fabsf(ct.i); - } -while( b > MACHEP ); -w->r = sum.r; -w->i = sum.i; -return; -} -*/ - - - ca = x + y * I; - ct = ca * I; - /* sqrt( 1 - z*z) */ - /* cmul( &ca, &ca, &zz ) */ - /*x * x - y * y */ - zz = (x - y) * (x + y) + (2.0f * x * y) * I; - - zz = 1.0f - crealf(zz) - cimagf(zz) * I; - z2 = csqrtf(zz); - - zz = ct + z2; - zz = clogf(zz); - /* multiply by 1/i = -i */ - w = zz * (-1.0f * I); - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/casinh.c b/lib/mlibc/options/ansi/generic/complex/casinh.c deleted file mode 100644 index 16238a6..0000000 --- a/lib/mlibc/options/ansi/generic/complex/casinh.c +++ /dev/null @@ -1,97 +0,0 @@ -/* $NetBSD: casinh.c,v 1.1 2007/08/20 16:01:31 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<casinh>>, <<casinhf>>---complex arc hyperbolic sine - -INDEX - casinh -INDEX - casinhf - -ANSI_SYNOPSIS - #include <complex.h> - double complex casinh(double complex <[z]>); - float complex casinhf(float complex <[z]>); - - -DESCRIPTION - @ifnottex - These functions compute the complex arc hyperbolic sine of <[z]>, - with branch cuts outside the interval [-i, +i] along the - imaginary axis. - @end ifnottex - @tex - These functions compute the complex arc hyperbolic sine of <[z]>, - with branch cuts outside the interval [$-i$, $+i$] along the - imaginary axis. - @end tex - - <<casinhf>> is identical to <<casinh>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - These functions return the complex arc hyperbolic sine value, - in the range of a strip mathematically unbounded along the - real axis and in the interval [-i*p/2, +i*p/2] along the - imaginary axis. - @end ifnottex - @tex - These functions return the complex arc hyperbolic sine value, - in the range of a strip mathematically unbounded along the - real axis and in the interval [$-i\pi/2$, $+i\pi/2$] along the - imaginary axis. - @end tex - -PORTABILITY - <<casinh>> and <<casinhf>> are ISO C99 - -QUICKREF - <<casinh>> and <<casinhf>> are ISO C99 - -*/ - - -#include <complex.h> - -double complex -casinh(double complex z) -{ - double complex w; - - w = -1.0 * I * casin(z * I); - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/casinhf.c b/lib/mlibc/options/ansi/generic/complex/casinhf.c deleted file mode 100644 index 0db55a0..0000000 --- a/lib/mlibc/options/ansi/generic/complex/casinhf.c +++ /dev/null @@ -1,44 +0,0 @@ -/* $NetBSD: casinhf.c,v 1.1 2007/08/20 16:01:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> - -float complex -casinhf(float complex z) -{ - float complex w; - - w = -1.0f * I * casinf(z * I); - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/catan.c b/lib/mlibc/options/ansi/generic/complex/catan.c deleted file mode 100644 index 0cf4739..0000000 --- a/lib/mlibc/options/ansi/generic/complex/catan.c +++ /dev/null @@ -1,130 +0,0 @@ -/* $NetBSD: catan.c,v 1.1 2007/08/20 16:01:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<catan>>, <<catanf>>---complex arc tangent - -INDEX - catan -INDEX - catanf - -ANSI_SYNOPSIS - #include <complex.h> - double complex catan(double complex <[z]>); - float complex catanf(float complex <[z]>); - - -DESCRIPTION - @ifnottex - These functions compute the complex arc tangent of <[z]>, - with branch cuts outside the interval [-i, +i] along the - imaginary axis. - @end ifnottex - @tex - These functions compute the complex arc tangent of <[z]>, - with branch cuts outside the interval [$-i$, $+i$] along the - imaginary axis. - @end tex - - <<catanf>> is identical to <<catan>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - These functions return the complex arc tangent value, in the range - of a strip mathematically unbounded along the imaginary axis - and in the interval [-pi/2, +pi/2] along the real axis. - @end ifnottex - @tex - These functions return the complex arc tangent, in the range - of a strip mathematically unbounded along the imaginary axis - and in the interval [$-\pi/2$, $+\pi/2$] along the real axis. - @end tex - -PORTABILITY - <<catan>> and <<catanf>> are ISO C99 - -QUICKREF - <<catan>> and <<catanf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> -#include "cephes_subr.h" - -#ifdef __weak_alias -__weak_alias(catan, _catan) -#endif - -double complex -catan(double complex z) -{ - double complex w; - double a, t, x, x2, y; - - x = creal(z); - y = cimag(z); - - if ((x == 0.0) && (y > 1.0)) - goto ovrf; - - x2 = x * x; - a = 1.0 - x2 - (y * y); - if (a == 0.0) - goto ovrf; - - t = 0.5 * atan2(2.0 * x, a); - w = __mlibc_redupi(t); - - t = y - 1.0; - a = x2 + (t * t); - if (a == 0.0) - goto ovrf; - - t = y + 1.0; - a = (x2 + (t * t))/a; - w = w + (0.25 * log(a)) * I; - return w; - -ovrf: -#if 0 - mtherr ("catan", OVERFLOW); -#endif - w = HUGE_VAL + HUGE_VAL * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/catanf.c b/lib/mlibc/options/ansi/generic/complex/catanf.c deleted file mode 100644 index 33c47df..0000000 --- a/lib/mlibc/options/ansi/generic/complex/catanf.c +++ /dev/null @@ -1,79 +0,0 @@ -/* $NetBSD: catanf.c,v 1.1 2007/08/20 16:01:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> -#include "cephes_subrf.h" - -#ifdef __weak_alias -__weak_alias(catanf, _catanf) -#endif - -float complex -catanf(float complex z) -{ - float complex w; - float a, t, x, x2, y; - - x = crealf(z); - y = cimagf(z); - - if ((x == 0.0f) && (y > 1.0f)) - goto ovrf; - - x2 = x * x; - a = 1.0f - x2 - (y * y); - if (a == 0.0f) - goto ovrf; - - t = 0.5f * atan2f(2.0f * x, a); - w = __mlibc_redupif(t); - - t = y - 1.0f; - a = x2 + (t * t); - if (a == 0.0f) - goto ovrf; - - t = y + 1.0f; - a = (x2 + (t * t))/a; - w = w + (0.25f * logf(a)) * I; - return w; - -ovrf: -#if 0 - mtherr ("catan", OVERFLOW); -#endif - w = HUGE_VALF + HUGE_VALF * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/catanh.c b/lib/mlibc/options/ansi/generic/complex/catanh.c deleted file mode 100644 index 2b9ef9e..0000000 --- a/lib/mlibc/options/ansi/generic/complex/catanh.c +++ /dev/null @@ -1,90 +0,0 @@ -/* $NetBSD: catanh.c,v 1.1 2007/08/20 16:01:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<catanh>>, <<catanhf>>---complex arc hyperbolic tangent - -INDEX - catanh -INDEX - catanhf - -ANSI_SYNOPSIS - #include <complex.h> - double complex catanh(double complex <[z]>); - float complex catanhf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex arc hyperbolic tan of <[z]>, - with branch cuts outside the interval [-1, +1] along the - real axis. - - <<catanhf>> is identical to <<catanh>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - These functions return the complex arc hyperbolic tangent value, - in the range of a strip mathematically unbounded along the - real axis and in the interval [-i*p/2, +i*p/2] along the - imaginary axis. - @end ifnottex - @tex - These functions return the complex arc hyperbolic tangent value, - in the range of a strip mathematically unbounded along the - real axis and in the interval [$-i\pi/2$, $+i\pi/2$] along the - imaginary axis. - @end tex - -PORTABILITY - <<catanh>> and <<catanhf>> are ISO C99 - -QUICKREF - <<catanh>> and <<catanhf>> are ISO C99 - -*/ - - -#include <complex.h> - -double complex -catanh(double complex z) -{ - double complex w; - - w = -1.0 * I * catan(z * I); - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/catanhf.c b/lib/mlibc/options/ansi/generic/complex/catanhf.c deleted file mode 100644 index fe6127a..0000000 --- a/lib/mlibc/options/ansi/generic/complex/catanhf.c +++ /dev/null @@ -1,44 +0,0 @@ -/* $NetBSD: catanhf.c,v 1.1 2007/08/20 16:01:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> - -float complex -catanhf(float complex z) -{ - float complex w; - - w = -1.0f * I * catanf(z * I); - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ccos.c b/lib/mlibc/options/ansi/generic/complex/ccos.c deleted file mode 100644 index ebb52bf..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ccos.c +++ /dev/null @@ -1,81 +0,0 @@ -/* $NetBSD: ccos.c,v 1.1 2007/08/20 16:01:32 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<ccos>>, <<ccosf>>---complex cosine - -INDEX - ccos -INDEX - ccosf - -ANSI_SYNOPSIS - #include <complex.h> - double complex ccos(double complex <[z]>); - float complex ccosf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex cosine of <[z]>. - - <<ccosf>> is identical to <<ccos>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - These functions return the complex cosine value. - -PORTABILITY - <<ccos>> and <<ccosf>> are ISO C99 - -QUICKREF - <<ccos>> and <<ccosf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> -#include "cephes_subr.h" - -double complex -ccos(double complex z) -{ - double complex w; - double ch, sh; - - __mlibc_cchsh(cimag(z), &ch, &sh); - w = cos(creal(z)) * ch - (sin(creal(z)) * sh) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ccosf.c b/lib/mlibc/options/ansi/generic/complex/ccosf.c deleted file mode 100644 index db7fab3..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ccosf.c +++ /dev/null @@ -1,48 +0,0 @@ -/* $NetBSD: ccosf.c,v 1.1 2007/08/20 16:01:33 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> -#include "cephes_subrf.h" - -float complex -ccosf(float complex z) -{ - float complex w; - float ch, sh; - - __mlibc_cchshf(cimagf(z), &ch, &sh); - w = cosf(crealf(z)) * ch - (sinf(crealf(z)) * sh) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ccosh.c b/lib/mlibc/options/ansi/generic/complex/ccosh.c deleted file mode 100644 index 223a5ed..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ccosh.c +++ /dev/null @@ -1,81 +0,0 @@ -/* $NetBSD: ccosh.c,v 1.1 2007/08/20 16:01:33 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<ccosh>>, <<ccoshf>>---complex hyperbolic cosine - -INDEX - ccosh -INDEX - ccoshf - -ANSI_SYNOPSIS - #include <complex.h> - double complex ccosh(double complex <[z]>); - float complex ccoshf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex hyperbolic cosine of <[z]>. - - <<ccoshf>> is identical to <<ccosh>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - These functions return the complex hyperbolic cosine value. - -PORTABILITY - <<ccosh>> and <<ccoshf>> are ISO C99 - -QUICKREF - <<ccosh>> and <<ccoshf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -double complex -ccosh(double complex z) -{ - double complex w; - double x, y; - - x = creal(z); - y = cimag(z); - w = cosh(x) * cos(y) + (sinh(x) * sin(y)) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ccoshf.c b/lib/mlibc/options/ansi/generic/complex/ccoshf.c deleted file mode 100644 index af11353..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ccoshf.c +++ /dev/null @@ -1,48 +0,0 @@ -/* $NetBSD: ccoshf.c,v 1.1 2007/08/20 16:01:33 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -ccoshf(float complex z) -{ - float complex w; - float x, y; - - x = crealf(z); - y = cimagf(z); - w = coshf(x) * cosf(y) + (sinhf(x) * sinf(y)) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cephes_subr.c b/lib/mlibc/options/ansi/generic/complex/cephes_subr.c deleted file mode 100644 index fe08b42..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cephes_subr.c +++ /dev/null @@ -1,126 +0,0 @@ -/* $NetBSD: cephes_subr.c,v 1.1 2007/08/20 16:01:33 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> -#include "cephes_subr.h" - -/* calculate cosh and sinh */ - -void -__mlibc_cchsh(double x, double *c, double *s) -{ - double e, ei; - - if (fabs(x) <= 0.5) { - *c = cosh(x); - *s = sinh(x); - } else { - e = exp(x); - ei = 0.5 / e; - e = 0.5 * e; - *s = e - ei; - *c = e + ei; - } -} - -/* Program to subtract nearest integer multiple of PI */ - -/* extended precision value of PI: */ -static const double DP1 = 3.14159265160560607910E0; -static const double DP2 = 1.98418714791870343106E-9; -static const double DP3 = 1.14423774522196636802E-17; -#define MACHEP 1.1e-16 - -double -__mlibc_redupi(double x) -{ - double t; - long i; - - t = x / M_PI; - if (t >= 0.0) - t += 0.5; - else - t -= 0.5; - - i = t; /* the multiple */ - t = i; - t = ((x - t * DP1) - t * DP2) - t * DP3; - return t; -} - -/* Taylor series expansion for cosh(2y) - cos(2x) */ - -double -__mlibc_ctans(double complex z) -{ - double f, x, x2, y, y2, rn, t; - double d; - - x = fabs(2.0 * creal(z)); - y = fabs(2.0 * cimag(z)); - - x = __mlibc_redupi(x); - - x = x * x; - y = y * y; - x2 = 1.0; - y2 = 1.0; - f = 1.0; - rn = 0.0; - d = 0.0; - do { - rn += 1.0; - f *= rn; - rn += 1.0; - f *= rn; - x2 *= x; - y2 *= y; - t = y2 + x2; - t /= f; - d += t; - - rn += 1.0; - f *= rn; - rn += 1.0; - f *= rn; - x2 *= x; - y2 *= y; - t = y2 - x2; - t /= f; - d += t; - } while (fabs(t/d) > MACHEP); - return d; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cephes_subr.h b/lib/mlibc/options/ansi/generic/complex/cephes_subr.h deleted file mode 100644 index 719075e..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cephes_subr.h +++ /dev/null @@ -1,9 +0,0 @@ -/* $NetBSD: cephes_subr.h,v 1.1 2007/08/20 16:01:33 drochner Exp $ */ - -#ifndef __MLIBC_ABI_ONLY - -void __mlibc_cchsh(double, double *, double *); -double __mlibc_redupi(double); -double __mlibc_ctans(double complex); - -#endif /* !__MLIBC_ABI_ONLY */ diff --git a/lib/mlibc/options/ansi/generic/complex/cephes_subrf.c b/lib/mlibc/options/ansi/generic/complex/cephes_subrf.c deleted file mode 100644 index 1ce18e5..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cephes_subrf.c +++ /dev/null @@ -1,125 +0,0 @@ -/* $NetBSD: cephes_subrf.c,v 1.1 2007/08/20 16:01:34 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> -#include "cephes_subrf.h" - -/* calculate cosh and sinh */ - -void -__mlibc_cchshf(float x, float *c, float *s) -{ - float e, ei; - - if (fabsf(x) <= 0.5f) { - *c = coshf(x); - *s = sinhf(x); - } else { - e = expf(x); - ei = 0.5f / e; - e = 0.5f * e; - *s = e - ei; - *c = e + ei; - } -} - -/* Program to subtract nearest integer multiple of PI */ - -/* extended precision value of PI: */ -static const double DP1 = 3.140625; -static const double DP2 = 9.67502593994140625E-4; -static const double DP3 = 1.509957990978376432E-7; -#define MACHEPF 3.0e-8 - -float -__mlibc_redupif(float x) -{ - float t; - long i; - - t = x / (float)M_PI; - if (t >= 0.0f) - t += 0.5f; - else - t -= 0.5f; - - i = t; /* the multiple */ - t = i; - t = ((x - t * DP1) - t * DP2) - t * DP3; - return t; -} - -/* Taylor series expansion for cosh(2y) - cos(2x) */ - -float -__mlibc_ctansf(float complex z) -{ - float f, x, x2, y, y2, rn, t, d; - - x = fabsf(2.0f * crealf(z)); - y = fabsf(2.0f * cimagf(z)); - - x = __mlibc_redupif(x); - - x = x * x; - y = y * y; - x2 = 1.0f; - y2 = 1.0f; - f = 1.0f; - rn = 0.0f; - d = 0.0f; - do { - rn += 1.0f; - f *= rn; - rn += 1.0f; - f *= rn; - x2 *= x; - y2 *= y; - t = y2 + x2; - t /= f; - d += t; - - rn += 1.0f; - f *= rn; - rn += 1.0f; - f *= rn; - x2 *= x; - y2 *= y; - t = y2 - x2; - t /= f; - d += t; - } while (fabsf(t/d) > MACHEPF); - return d; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cephes_subrf.h b/lib/mlibc/options/ansi/generic/complex/cephes_subrf.h deleted file mode 100644 index 84cdd82..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cephes_subrf.h +++ /dev/null @@ -1,9 +0,0 @@ -/* $NetBSD: cephes_subrf.h,v 1.1 2007/08/20 16:01:34 drochner Exp $ */ - -#ifndef __MLIBC_ABI_ONLY - -void __mlibc_cchshf(float, float *, float *); -float __mlibc_redupif(float); -float __mlibc_ctansf(float complex); - -#endif /* !__MLIBC_ABI_ONLY */ diff --git a/lib/mlibc/options/ansi/generic/complex/cexp.c b/lib/mlibc/options/ansi/generic/complex/cexp.c deleted file mode 100644 index b9a3fd0..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cexp.c +++ /dev/null @@ -1,82 +0,0 @@ -/* $NetBSD: cexp.c,v 1.1 2007/08/20 16:01:34 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cexp>>, <<cexpf>>---complex base-e exponential - -INDEX - cexp -INDEX - cexpf - -ANSI_SYNOPSIS - #include <complex.h> - double complex cexp(double complex <[z]>); - float complex cexpf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex base-<[e]> exponential of <[z]>. - - <<cexpf>> is identical to <<cexp>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The cexp functions return the complex base-<[e]> exponential value. - -PORTABILITY - <<cexp>> and <<cexpf>> are ISO C99 - -QUICKREF - <<cexp>> and <<cexpf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -double complex -cexp(double complex z) -{ - double complex w; - double r, x, y; - - x = creal(z); - y = cimag(z); - r = exp(x); - w = r * cos(y) + r * sin(y) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cexpf.c b/lib/mlibc/options/ansi/generic/complex/cexpf.c deleted file mode 100644 index 07fab1f..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cexpf.c +++ /dev/null @@ -1,49 +0,0 @@ -/* $NetBSD: cexpf.c,v 1.1 2007/08/20 16:01:34 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -cexpf(float complex z) -{ - float complex w; - float r, x, y; - - x = crealf(z); - y = cimagf(z); - r = expf(x); - w = r * cosf(y) + r * sinf(y) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cimag.c b/lib/mlibc/options/ansi/generic/complex/cimag.c deleted file mode 100644 index 24619f0..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cimag.c +++ /dev/null @@ -1,54 +0,0 @@ -/* $NetBSD: cimag.c,v 1.2 2010/09/15 16:11:29 christos Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cimag>>, <<cimagf>>---imaginary part - -INDEX - cimag -INDEX - cimagf - -ANSI_SYNOPSIS - #include <complex.h> - double cimag(double complex <[z]>); - float cimagf(float complex <[z]>); - - -DESCRIPTION - These functions compute the imaginary part of <[z]>. - - <<cimagf>> is identical to <<cimag>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The cimag functions return the imaginary part value (as a real). - -PORTABILITY - <<cimag>> and <<cimagf>> are ISO C99 - -QUICKREF - <<cimag>> and <<cimagf>> are ISO C99 - -*/ - - -#include <complex.h> - -#include "fdlibm.h" - -double -cimag(double complex z) -{ - double_complex w = { .z = z }; - - return (IMAG_PART(w)); -} diff --git a/lib/mlibc/options/ansi/generic/complex/cimagf.c b/lib/mlibc/options/ansi/generic/complex/cimagf.c deleted file mode 100644 index 28ed81c..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cimagf.c +++ /dev/null @@ -1,21 +0,0 @@ -/* $NetBSD: cimagf.c,v 1.2 2010/09/15 16:11:29 christos Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> - -#include "fdlibm.h" - -float -cimagf(float complex z) -{ - float_complex w = { .z = z }; - - return (IMAG_PART(w)); -} diff --git a/lib/mlibc/options/ansi/generic/complex/clog.c b/lib/mlibc/options/ansi/generic/complex/clog.c deleted file mode 100644 index f7ad3d2..0000000 --- a/lib/mlibc/options/ansi/generic/complex/clog.c +++ /dev/null @@ -1,91 +0,0 @@ -/* $NetBSD: clog.c,v 1.1 2007/08/20 16:01:35 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<clog>>, <<clogf>>---complex base-e logarithm - -INDEX - clog -INDEX - clogf - -ANSI_SYNOPSIS - #include <complex.h> - double complex clog(double complex <[z]>); - float complex clogf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex natural (base-<[e]>) logarithm - of <[z]>, with a branch cut along the negative real axis. - - <<clogf>> is identical to <<clog>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - @ifnottex - The clog functions return the complex natural logarithm value, in - the range of a strip mathematically unbounded along the real axis - and in the interval [-i*pi , +i*pi] along the imaginary axis. - @end ifnottex - @tex - The clog functions return the complex natural logarithm value, in - the range of a strip mathematically unbounded along the real axis - and in the interval [$-i\pi$, $+i\pi$] along the imaginary axis. - @end tex - -PORTABILITY - <<clog>> and <<clogf>> are ISO C99 - -QUICKREF - <<clog>> and <<clogf>> are ISO C99 - -*/ - -#include <complex.h> -#include <math.h> - -double complex -clog(double complex z) -{ - double complex w; - double p, rr; - - rr = cabs(z); - p = log(rr); - rr = atan2(cimag(z), creal(z)); - w = p + rr * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/clogf.c b/lib/mlibc/options/ansi/generic/complex/clogf.c deleted file mode 100644 index 078cea5..0000000 --- a/lib/mlibc/options/ansi/generic/complex/clogf.c +++ /dev/null @@ -1,49 +0,0 @@ -/* $NetBSD: clogf.c,v 1.1 2007/08/20 16:01:35 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -clogf(float complex z) -{ - float complex w; - float p, rr; - - rr = cabsf(z); - p = logf(rr); - rr = atan2f(cimagf(z), crealf(z)); - w = p + rr * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/conj.c b/lib/mlibc/options/ansi/generic/complex/conj.c deleted file mode 100644 index a761b5a..0000000 --- a/lib/mlibc/options/ansi/generic/complex/conj.c +++ /dev/null @@ -1,56 +0,0 @@ -/* $NetBSD: conj.c,v 1.2 2010/09/15 16:11:29 christos Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<conj>>, <<conjf>>---complex conjugate - -INDEX - conj -INDEX - conjf - -ANSI_SYNOPSIS - #include <complex.h> - double complex conj(double complex <[z]>); - float complex conjf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex conjugate of <[z]>, - by reversing the sign of its imaginary part. - - <<conjf>> is identical to <<conj>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The conj functions return the complex conjugate value. - -PORTABILITY - <<conj>> and <<conjf>> are ISO C99 - -QUICKREF - <<conj>> and <<conjf>> are ISO C99 - -*/ - -#include <complex.h> - -#include "fdlibm.h" - -double complex -conj(double complex z) -{ - double_complex w = { .z = z }; - - IMAG_PART(w) = -IMAG_PART(w); - - return (w.z); -} diff --git a/lib/mlibc/options/ansi/generic/complex/conjf.c b/lib/mlibc/options/ansi/generic/complex/conjf.c deleted file mode 100644 index 0ca71ef..0000000 --- a/lib/mlibc/options/ansi/generic/complex/conjf.c +++ /dev/null @@ -1,23 +0,0 @@ -/* $NetBSD: conjf.c,v 1.2 2010/09/15 16:11:29 christos Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> - -#include "fdlibm.h" - -float complex -conjf(float complex z) -{ - float_complex w = { .z = z }; - - IMAG_PART(w) = -IMAG_PART(w); - - return (w.z); -} diff --git a/lib/mlibc/options/ansi/generic/complex/cpow.c b/lib/mlibc/options/ansi/generic/complex/cpow.c deleted file mode 100644 index b60f7be..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cpow.c +++ /dev/null @@ -1,101 +0,0 @@ -/* $NetBSD: cpow.c,v 1.1 2007/08/20 16:01:35 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cpow>>, <<cpowf>>---complex power - -INDEX - cpow -INDEX - cpowf - -ANSI_SYNOPSIS - #include <complex.h> - double complex cpow(double complex <[x]>, double complex <[y]>); - float complex cpowf(float complex <[x]>, float complex <[y]>); - - -DESCRIPTION - @ifnottex - The cpow functions compute the complex power function x^y - power, with a branch cut for the first parameter along the - negative real axis. - @end ifnottex - @tex - The cpow functions compute the complex power function $x^y$ - power, with a branch cut for the first parameter along the - negative real axis. - @end tex - - <<cpowf>> is identical to <<cpow>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The cpow functions return the complex power function value. - -PORTABILITY - <<cpow>> and <<cpowf>> are ISO C99 - -QUICKREF - <<cpow>> and <<cpowf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -double complex -cpow(double complex a, double complex z) -{ - double complex w; - double x, y, r, theta, absa, arga; - - x = creal(z); - y = cimag(z); - absa = cabs(a); - if (absa == 0.0) { - return (0.0 + 0.0 * I); - } - arga = carg(a); - r = pow(absa, x); - theta = x * arga; - if (y != 0.0) { - r = r * exp(-y * arga); - theta = theta + y * log(absa); - } - w = r * cos(theta) + (r * sin(theta)) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cpowf.c b/lib/mlibc/options/ansi/generic/complex/cpowf.c deleted file mode 100644 index 1e736af..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cpowf.c +++ /dev/null @@ -1,59 +0,0 @@ -/* $NetBSD: cpowf.c,v 1.1 2007/08/20 16:01:36 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -cpowf(float complex a, float complex z) -{ - float complex w; - float x, y, r, theta, absa, arga; - - x = crealf(z); - y = cimagf(z); - absa = cabsf(a); - if (absa == 0.0f) { - return (0.0f + 0.0f * I); - } - arga = cargf(a); - r = powf(absa, x); - theta = x * arga; - if (y != 0.0f) { - r = r * expf(-y * arga); - theta = theta + y * logf(absa); - } - w = r * cosf(theta) + (r * sinf(theta)) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/cproj.c b/lib/mlibc/options/ansi/generic/complex/cproj.c deleted file mode 100644 index 0ed50f2..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cproj.c +++ /dev/null @@ -1,105 +0,0 @@ -/* $NetBSD: cproj.c,v 1.3 2010/09/20 17:51:38 christos Exp $ */ - -/*- - * Copyright (c) 2010 The NetBSD Foundation, Inc. - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<cproj>>, <<cprojf>>--- Riemann sphere projection - -INDEX - cproj -INDEX - cprojf - -ANSI_SYNOPSIS - #include <complex.h> - double complex cproj(double complex <[z]>); - float complex cprojf(float complex <[z]>); - - -DESCRIPTION - These functions compute a projection of <[z]> onto the Riemann - sphere: <[z]> projects to <[z]> except that all complex infinities - (even those with one infinite part and one NaN part) project - to positive infinity on the real axis. If <[z]> has an infinite part, - then <<cproj>>(<[z]>) is equivalent to - - INFINITY + I * copysign(0.0, cimag(z)) - - <<cprojf>> is identical to <<cproj>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The cproj functions return the value of the projection onto - the Riemann sphere. - -PORTABILITY - <<cproj>> and <<cprojf>> are ISO C99 - -QUICKREF - <<cproj>> and <<cprojf>> are ISO C99 - -*/ - -/*__RCSID("$NetBSD: cproj.c,v 1.3 2010/09/20 17:51:38 christos Exp $"); */ - -#include <complex.h> -#include <math.h> - -#include "fdlibm.h" - -/* - * cproj(double complex z) - * - * These functions return the value of the projection (not stereographic!) - * onto the Riemann sphere. - * - * z projects to z, except that all complex infinities (even those with one - * infinite part and one NaN part) project to positive infinity on the real axis. - * If z has an infinite part, then cproj(z) shall be equivalent to: - * - * INFINITY + I * copysign(0.0, cimag(z)) - */ -double complex -cproj(double complex z) -{ - double_complex w = { .z = z }; - - if (isinf(creal(z)) || isinf(cimag(z))) { -#ifdef __INFINITY - REAL_PART(w) = __INFINITY; -#else - REAL_PART(w) = INFINITY; -#endif - IMAG_PART(w) = copysign(0.0, cimag(z)); - } - - return (w.z); -} diff --git a/lib/mlibc/options/ansi/generic/complex/cprojf.c b/lib/mlibc/options/ansi/generic/complex/cprojf.c deleted file mode 100644 index 76c3d8a..0000000 --- a/lib/mlibc/options/ansi/generic/complex/cprojf.c +++ /dev/null @@ -1,67 +0,0 @@ -/* $NetBSD: cprojf.c,v 1.3 2010/09/20 17:51:38 christos Exp $ */ - -/*- - * Copyright (c) 2010 The NetBSD Foundation, Inc. - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/*__RCSID("$NetBSD: cprojf.c,v 1.3 2010/09/20 17:51:38 christos Exp $"); */ - -#include <complex.h> -#include <math.h> - -#include "fdlibm.h" - -/* - * cprojf(float complex z) - * - * These functions return the value of the projection (not stereographic!) - * onto the Riemann sphere. - * - * z projects to z, except that all complex infinities (even those with one - * infinite part and one NaN part) project to positive infinity on the real axis. - * If z has an infinite part, then cproj(z) shall be equivalent to: - * - * INFINITY + I * copysign(0.0, cimag(z)) - */ - -float complex -cprojf(float complex z) -{ - float_complex w = { .z = z }; - - if (isinf(crealf(z)) || isinf(cimagf(z))) { -#ifdef __INFINITY - REAL_PART(w) = __INFINITY; -#else - REAL_PART(w) = INFINITY; -#endif - IMAG_PART(w) = copysignf(0.0, cimagf(z)); - } - - return (w.z); -} diff --git a/lib/mlibc/options/ansi/generic/complex/creal.c b/lib/mlibc/options/ansi/generic/complex/creal.c deleted file mode 100644 index 07bf96f..0000000 --- a/lib/mlibc/options/ansi/generic/complex/creal.c +++ /dev/null @@ -1,54 +0,0 @@ -/* $NetBSD: creal.c,v 1.2 2010/09/15 16:11:29 christos Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<creal>>, <<crealf>>---real part - -INDEX - creal -INDEX - crealf - -ANSI_SYNOPSIS - #include <complex.h> - double creal(double complex <[z]>); - float crealf(float complex <[z]>); - - -DESCRIPTION - These functions compute the real part of <[z]>. - - <<crealf>> is identical to <<creal>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The creal functions return the real part value. - -PORTABILITY - <<creal>> and <<crealf>> are ISO C99 - -QUICKREF - <<creal>> and <<crealf>> are ISO C99 - -*/ - - -#include <complex.h> - -#include "fdlibm.h" - -double -creal(double complex z) -{ - double_complex w = { .z = z }; - - return (REAL_PART(w)); -} diff --git a/lib/mlibc/options/ansi/generic/complex/crealf.c b/lib/mlibc/options/ansi/generic/complex/crealf.c deleted file mode 100644 index 245986d..0000000 --- a/lib/mlibc/options/ansi/generic/complex/crealf.c +++ /dev/null @@ -1,21 +0,0 @@ -/* $NetBSD: crealf.c,v 1.2 2010/09/15 16:11:29 christos Exp $ */ - -/* - * Written by Matthias Drochner <drochner@NetBSD.org>. - * Public domain. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> - -#include "fdlibm.h" - -float -crealf(float complex z) -{ - float_complex w = { .z = z }; - - return (REAL_PART(w)); -} diff --git a/lib/mlibc/options/ansi/generic/complex/csin.c b/lib/mlibc/options/ansi/generic/complex/csin.c deleted file mode 100644 index b32d057..0000000 --- a/lib/mlibc/options/ansi/generic/complex/csin.c +++ /dev/null @@ -1,81 +0,0 @@ -/* $NetBSD: csin.c,v 1.1 2007/08/20 16:01:36 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<csin>>, <<csinf>>---complex sine - -INDEX - csin -INDEX - csinf - -ANSI_SYNOPSIS - #include <complex.h> - double complex csin(double complex <[z]>); - float complex csinf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex sine of <[z]>. - - <<csinf>> is identical to <<csin>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - These functions return the complex sine value. - -PORTABILITY - <<csin>> and <<csinf>> are ISO C99 - -QUICKREF - <<csin>> and <<csinf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> -#include "cephes_subr.h" - -double complex -csin(double complex z) -{ - double complex w; - double ch, sh; - - __mlibc_cchsh(cimag(z), &ch, &sh); - w = sin(creal(z)) * ch + (cos(creal(z)) * sh) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/csinf.c b/lib/mlibc/options/ansi/generic/complex/csinf.c deleted file mode 100644 index 0d81d41..0000000 --- a/lib/mlibc/options/ansi/generic/complex/csinf.c +++ /dev/null @@ -1,48 +0,0 @@ -/* $NetBSD: csinf.c,v 1.1 2007/08/20 16:01:36 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> -#include "cephes_subrf.h" - -float complex -csinf(float complex z) -{ - float complex w; - float ch, sh; - - __mlibc_cchshf(cimagf(z), &ch, &sh); - w = sinf(crealf(z)) * ch + (cosf(crealf(z)) * sh) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/csinh.c b/lib/mlibc/options/ansi/generic/complex/csinh.c deleted file mode 100644 index f117162..0000000 --- a/lib/mlibc/options/ansi/generic/complex/csinh.c +++ /dev/null @@ -1,80 +0,0 @@ -/* $NetBSD: csinh.c,v 1.1 2007/08/20 16:01:36 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<csinh>>, <<csinhf>>---complex hyperbolic sine - -INDEX - csinh -INDEX - csinhf - -ANSI_SYNOPSIS - #include <complex.h> - double complex csinh(double complex <[z]>); - float complex csinhf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex hyperbolic sine of <[z]>. - - <<ccoshf>> is identical to <<ccosh>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - These functions return the complex hyperbolic sine value. - -PORTABILITY - <<csinh>> and <<csinhf>> are ISO C99 - -QUICKREF - <<csinh>> and <<csinhf>> are ISO C99 - -*/ - -#include <complex.h> -#include <math.h> - -double complex -csinh(double complex z) -{ - double complex w; - double x, y; - - x = creal(z); - y = cimag(z); - w = sinh(x) * cos(y) + (cosh(x) * sin(y)) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/csinhf.c b/lib/mlibc/options/ansi/generic/complex/csinhf.c deleted file mode 100644 index 3cd6ba7..0000000 --- a/lib/mlibc/options/ansi/generic/complex/csinhf.c +++ /dev/null @@ -1,48 +0,0 @@ -/* $NetBSD: csinhf.c,v 1.1 2007/08/20 16:01:37 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -csinhf(float complex z) -{ - float complex w; - float x, y; - - x = crealf(z); - y = cimagf(z); - w = sinhf(x) * cosf(y) + (coshf(x) * sinf(y)) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/csqrt.c b/lib/mlibc/options/ansi/generic/complex/csqrt.c deleted file mode 100644 index b144b7c..0000000 --- a/lib/mlibc/options/ansi/generic/complex/csqrt.c +++ /dev/null @@ -1,137 +0,0 @@ -/* $NetBSD: csqrt.c,v 1.1 2007/08/20 16:01:37 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<csqrt>>, <<csqrtf>>---complex square root - -INDEX - csqrt -INDEX - csqrtf - -ANSI_SYNOPSIS - #include <complex.h> - double complex csqrt(double complex <[z]>); - float complex csqrtf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex square root of <[z]>, with - a branch cut along the negative real axis. - - <<csqrtf>> is identical to <<csqrt>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - The csqrt functions return the complex square root value, in - the range of the right halfplane (including the imaginary axis). - -PORTABILITY - <<csqrt>> and <<csqrtf>> are ISO C99 - -QUICKREF - <<csqrt>> and <<csqrtf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -double complex -csqrt(double complex z) -{ - double complex w; - double x, y, r, t, scale; - - x = creal (z); - y = cimag (z); - - if (y == 0.0) { - if (x == 0.0) { - w = 0.0 + y * I; - } else { - r = fabs(x); - r = sqrt(r); - if (x < 0.0) { - w = 0.0 + r * I; - } else { - w = r + y * I; - } - } - return w; - } - if (x == 0.0) { - r = fabs(y); - r = sqrt(0.5 * r); - if (y > 0) - w = r + r * I; - else - w = r - r * I; - return w; - } - /* Rescale to avoid internal overflow or underflow. */ - if ((fabs(x) > 4.0) || (fabs(y) > 4.0)) { - x *= 0.25; - y *= 0.25; - scale = 2.0; - } else { -#if 1 - x *= 1.8014398509481984e16; /* 2^54 */ - y *= 1.8014398509481984e16; - scale = 7.450580596923828125e-9; /* 2^-27 */ -#else - x *= 4.0; - y *= 4.0; - scale = 0.5; -#endif - } - w = x + y * I; - r = cabs(w); - if (x > 0) { - t = sqrt(0.5 * r + 0.5 * x); - r = scale * fabs((0.5 * y) / t ); - t *= scale; - } else { - r = sqrt(0.5 * r - 0.5 * x); - t = scale * fabs((0.5 * y) / r); - r *= scale; - } - if (y < 0) - w = t - r * I; - else - w = t + r * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/csqrtf.c b/lib/mlibc/options/ansi/generic/complex/csqrtf.c deleted file mode 100644 index 13451fa..0000000 --- a/lib/mlibc/options/ansi/generic/complex/csqrtf.c +++ /dev/null @@ -1,102 +0,0 @@ -/* $NetBSD: csqrtf.c,v 1.1 2007/08/20 16:01:37 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -csqrtf(float complex z) -{ - float complex w; - float x, y, r, t, scale; - - x = crealf (z); - y = cimagf (z); - - if (y == 0.0f) { - if (x < 0.0f) { - w = 0.0f + sqrtf(-x) * I; - return w; - } else if (x == 0.0f) { - return (0.0f + y * I); - } else { - w = sqrtf(x) + y * I; - return w; - } - } - - if (x == 0.0f) { - r = fabsf(y); - r = sqrtf(0.5f * r); - if (y > 0) - w = r + r * I; - else - w = r - r * I; - return w; - } - - /* Rescale to avoid internal overflow or underflow. */ - if ((fabsf(x) > 4.0f) || (fabsf(y) > 4.0f)) { - x *= 0.25f; - y *= 0.25f; - scale = 2.0f; - } else { -#if 1 - x *= 6.7108864e7f; /* 2^26 */ - y *= 6.7108864e7f; - scale = 1.220703125e-4f; /* 2^-13 */ -#else - x *= 4.0f; - y *= 4.0f; - scale = 0.5f; -#endif - } - w = x + y * I; - r = cabsf(w); - if( x > 0 ) { - t = sqrtf(0.5f * r + 0.5f * x); - r = scale * fabsf((0.5f * y) / t); - t *= scale; - } else { - r = sqrtf(0.5f * r - 0.5f * x); - t = scale * fabsf((0.5f * y) / r); - r *= scale; - } - - if (y < 0) - w = t - r * I; - else - w = t + r * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ctan.c b/lib/mlibc/options/ansi/generic/complex/ctan.c deleted file mode 100644 index 600989d..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ctan.c +++ /dev/null @@ -1,91 +0,0 @@ -/* $NetBSD: ctan.c,v 1.1 2007/08/20 16:01:37 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<ctan>>, <<ctanf>>---complex tangent - -INDEX - ctan -INDEX - ctanf - -ANSI_SYNOPSIS - #include <complex.h> - double complex ctan(double complex <[z]>); - float complex ctanf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex tangent of <[z]>. - - <<ctanf>> is identical to <<ctan>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - These functions return the complex tangent value. - -PORTABILITY - <<ctan>> and <<ctanf>> are ISO C99 - -QUICKREF - <<ctan>> and <<ctanf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> -#include "cephes_subr.h" - -double complex -ctan(double complex z) -{ - double complex w; - double d; - - d = cos(2.0 * creal(z)) + cosh(2.0 * cimag(z)); - - if (fabs(d) < 0.25) - d = __mlibc_ctans(z); - - if (d == 0.0) { - /* mtherr ("ctan", OVERFLOW); */ - w = HUGE_VAL + HUGE_VAL * I; - return w; - } - - w = sin(2.0 * creal(z)) / d + (sinh(2.0 * cimag(z)) / d) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ctanf.c b/lib/mlibc/options/ansi/generic/complex/ctanf.c deleted file mode 100644 index 52360e0..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ctanf.c +++ /dev/null @@ -1,58 +0,0 @@ -/* $NetBSD: ctanf.c,v 1.1 2007/08/20 16:01:38 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> -#include "cephes_subrf.h" - -float complex -ctanf(float complex z) -{ - float complex w; - float d; - - d = cosf(2.0f * crealf(z)) + coshf(2.0f * cimagf(z)); - - if (fabsf(d) < 0.25f) - d = __mlibc_ctansf(z); - - if (d == 0.0f) { - /* mtherr ("ctan", OVERFLOW); */ - w = HUGE_VALF + HUGE_VALF * I; - return w; - } - - w = sinf(2.0f * crealf(z)) / d + (sinhf(2.0f * cimagf(z)) / d) * I; - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ctanh.c b/lib/mlibc/options/ansi/generic/complex/ctanh.c deleted file mode 100644 index db27e5b..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ctanh.c +++ /dev/null @@ -1,83 +0,0 @@ -/* $NetBSD: ctanh.c,v 1.1 2007/08/20 16:01:38 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -/* -FUNCTION - <<ctanh>>, <<ctanf>>---complex hyperbolic tangent - -INDEX - ctanh -INDEX - ctanhf - -ANSI_SYNOPSIS - #include <complex.h> - double complex ctanh(double complex <[z]>); - float complex ctanhf(float complex <[z]>); - - -DESCRIPTION - These functions compute the complex hyperbolic tangent of <[z]>. - - <<ctanhf>> is identical to <<ctanh>>, except that it performs - its calculations on <<floats complex>>. - -RETURNS - These functions return the complex hyperbolic tangent value. - -PORTABILITY - <<ctanh>> and <<ctanhf>> are ISO C99 - -QUICKREF - <<ctanh>> and <<ctanhf>> are ISO C99 - -*/ - - -#include <complex.h> -#include <math.h> - -double complex -ctanh(double complex z) -{ - double complex w; - double x, y, d; - - x = creal(z); - y = cimag(z); - d = cosh(2.0 * x) + cos(2.0 * y); - w = sinh(2.0 * x) / d + (sin(2.0 * y) / d) * I; - - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/ctanhf.c b/lib/mlibc/options/ansi/generic/complex/ctanhf.c deleted file mode 100644 index 6aaf20f..0000000 --- a/lib/mlibc/options/ansi/generic/complex/ctanhf.c +++ /dev/null @@ -1,50 +0,0 @@ -/* $NetBSD: ctanhf.c,v 1.1 2007/08/20 16:01:38 drochner Exp $ */ - -/*- - * Copyright (c) 2007 The NetBSD Foundation, Inc. - * All rights reserved. - * - * This code is derived from software written by Stephen L. Moshier. - * It is redistributed by the NetBSD Foundation by permission of the author. - * - * 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 NETBSD FOUNDATION, INC. 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 FOUNDATION 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. - * - * imported and modified include for newlib 2010/10/03 - * Marco Atzeri <marco_atzeri@yahoo.it> - */ - -#include <complex.h> -#include <math.h> - -float complex -ctanhf(float complex z) -{ - float complex w; - float x, y, d; - - x = crealf(z); - y = cimagf(z); - d = coshf(2.0f * x) + cosf(2.0f * y); - w = sinhf(2.0f * x) / d + (sinf(2.0f * y) / d) * I; - - return w; -} diff --git a/lib/mlibc/options/ansi/generic/complex/fdlibm.h b/lib/mlibc/options/ansi/generic/complex/fdlibm.h deleted file mode 100644 index 75cdd2a..0000000 --- a/lib/mlibc/options/ansi/generic/complex/fdlibm.h +++ /dev/null @@ -1,17 +0,0 @@ -#ifndef __MLIBC_FDLIBM_H -#define __MLIBC_FDLIBM_H - -#define REAL_PART(z) ((z).parts[0]) -#define IMAG_PART(z) ((z).parts[1]) - -typedef union { - float complex z; - float parts[2]; -} float_complex; - -typedef union { - double complex z; - double parts[2]; -} double_complex; - -#endif |