summaryrefslogtreecommitdiff
path: root/lib/mlibc/options/ansi/generic/complex
diff options
context:
space:
mode:
authorIan Moffett <ian@osmora.org>2024-03-07 17:28:52 -0500
committerIan Moffett <ian@osmora.org>2024-03-07 18:24:51 -0500
commitf5e48e94a2f4d4bbd6e5628c7f2afafc6dbcc459 (patch)
tree93b156621dc0303816b37f60ba88051b702d92f6 /lib/mlibc/options/ansi/generic/complex
parentbd5969fc876a10b18613302db7087ef3c40f18e1 (diff)
build: Build mlibc + add distclean target
Signed-off-by: Ian Moffett <ian@osmora.org>
Diffstat (limited to 'lib/mlibc/options/ansi/generic/complex')
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cabs.c53
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cabsf.c19
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cacos.c99
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cacosf.c46
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cacosh.c93
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cacoshf.c48
-rw-r--r--lib/mlibc/options/ansi/generic/complex/carg.c59
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cargf.c19
-rw-r--r--lib/mlibc/options/ansi/generic/complex/casin.c165
-rw-r--r--lib/mlibc/options/ansi/generic/complex/casinf.c122
-rw-r--r--lib/mlibc/options/ansi/generic/complex/casinh.c97
-rw-r--r--lib/mlibc/options/ansi/generic/complex/casinhf.c44
-rw-r--r--lib/mlibc/options/ansi/generic/complex/catan.c130
-rw-r--r--lib/mlibc/options/ansi/generic/complex/catanf.c79
-rw-r--r--lib/mlibc/options/ansi/generic/complex/catanh.c90
-rw-r--r--lib/mlibc/options/ansi/generic/complex/catanhf.c44
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ccos.c81
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ccosf.c48
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ccosh.c81
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ccoshf.c48
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cephes_subr.c126
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cephes_subr.h9
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cephes_subrf.c125
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cephes_subrf.h9
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cexp.c82
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cexpf.c49
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cimag.c54
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cimagf.c21
-rw-r--r--lib/mlibc/options/ansi/generic/complex/clog.c91
-rw-r--r--lib/mlibc/options/ansi/generic/complex/clogf.c49
-rw-r--r--lib/mlibc/options/ansi/generic/complex/conj.c56
-rw-r--r--lib/mlibc/options/ansi/generic/complex/conjf.c23
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cpow.c101
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cpowf.c59
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cproj.c105
-rw-r--r--lib/mlibc/options/ansi/generic/complex/cprojf.c67
-rw-r--r--lib/mlibc/options/ansi/generic/complex/creal.c54
-rw-r--r--lib/mlibc/options/ansi/generic/complex/crealf.c21
-rw-r--r--lib/mlibc/options/ansi/generic/complex/csin.c81
-rw-r--r--lib/mlibc/options/ansi/generic/complex/csinf.c48
-rw-r--r--lib/mlibc/options/ansi/generic/complex/csinh.c80
-rw-r--r--lib/mlibc/options/ansi/generic/complex/csinhf.c48
-rw-r--r--lib/mlibc/options/ansi/generic/complex/csqrt.c137
-rw-r--r--lib/mlibc/options/ansi/generic/complex/csqrtf.c102
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ctan.c91
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ctanf.c58
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ctanh.c83
-rw-r--r--lib/mlibc/options/ansi/generic/complex/ctanhf.c50
-rw-r--r--lib/mlibc/options/ansi/generic/complex/fdlibm.h17
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