mirror of
https://gitlab.alpinelinux.org/alpine/aports.git
synced 2026-01-07 01:32:26 +01:00
102 lines
3.2 KiB
Diff
102 lines
3.2 KiB
Diff
From 8733cd43293bfa04b4334db7873c4de613f5cf75 Mon Sep 17 00:00:00 2001
|
|
From: William Pitcock <nenolod@dereferenced.org>
|
|
Date: Tue, 1 Nov 2011 18:41:30 -0500
|
|
Subject: [PATCH 1/2] libm: implement generic cexp(), cexpf() and cexpl().
|
|
|
|
The cexp*() family of functions is defined by the C99 math standard as
|
|
implementing exponential functions for complex types.
|
|
|
|
Signed-off-by: William Pitcock <nenolod@dereferenced.org>
|
|
---
|
|
libm/Makefile.in | 3 +-
|
|
libm/cexp.c | 63 ++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
2 files changed, 65 insertions(+), 1 deletions(-)
|
|
create mode 100644 libm/cexp.c
|
|
|
|
diff --git a/libm/Makefile.in b/libm/Makefile.in
|
|
index f0aaa2a..af949e8 100644
|
|
--- a/libm/Makefile.in
|
|
+++ b/libm/Makefile.in
|
|
@@ -72,7 +72,8 @@ libm_CSRC := \
|
|
s_fpclassify.c s_fpclassifyf.c s_signbit.c s_signbitf.c \
|
|
s_isnan.c s_isnanf.c s_isinf.c s_isinff.c s_finitef.c \
|
|
s_fdim.c s_fma.c s_fmax.c s_fmin.c \
|
|
- s_remquo.c w_exp2.c
|
|
+ s_remquo.c w_exp2.c \
|
|
+ cexp.c
|
|
|
|
# Not implemented [yet?], see comment in float_wrappers.c:
|
|
# fdimf.o fmaf.o fmaxf.o fminf.o
|
|
diff --git a/libm/cexp.c b/libm/cexp.c
|
|
new file mode 100644
|
|
index 0000000..87512b7
|
|
--- /dev/null
|
|
+++ b/libm/cexp.c
|
|
@@ -0,0 +1,63 @@
|
|
+/*
|
|
+ * Copyright (c) 2011 William Pitcock <nenolod@dereferenced.org>
|
|
+ *
|
|
+ * Permission to use, copy, modify, and/or distribute this software for any
|
|
+ * purpose with or without fee is hereby granted, provided that the above
|
|
+ * copyright notice and this permission notice appear in all copies.
|
|
+ *
|
|
+ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
|
|
+ * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
|
|
+ * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
|
|
+ * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR 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.
|
|
+ */
|
|
+
|
|
+#include <features.h>
|
|
+#include <math.h>
|
|
+#include <complex.h>
|
|
+
|
|
+__complex__ double cexp(__complex__ double z)
|
|
+{
|
|
+ __complex__ double ret;
|
|
+ double r_exponent = exp(__real__ z);
|
|
+
|
|
+ __real__ ret = r_exponent * cos(__imag__ z);
|
|
+ __imag__ ret = r_exponent * sin(__imag__ z);
|
|
+
|
|
+ return ret;
|
|
+}
|
|
+libm_hidden_def(cexp)
|
|
+
|
|
+libm_hidden_proto(cexpf)
|
|
+__complex__ float cexpf(__complex__ float z)
|
|
+{
|
|
+ __complex__ float ret;
|
|
+ double r_exponent = exp(__real__ z);
|
|
+
|
|
+ __real__ ret = r_exponent * cosf(__imag__ z);
|
|
+ __imag__ ret = r_exponent * sinf(__imag__ z);
|
|
+
|
|
+ return ret;
|
|
+}
|
|
+libm_hidden_def(cexpf)
|
|
+
|
|
+#if defined __UCLIBC_HAS_LONG_DOUBLE_MATH__ && !defined __NO_LONG_DOUBLE_MATH
|
|
+libm_hidden_proto(cexpl)
|
|
+__complex__ long double cexpl(__complex__ long double z)
|
|
+{
|
|
+ __complex__ long double ret;
|
|
+ long double r_exponent = expl(__real__ z);
|
|
+
|
|
+ __real__ ret = r_exponent * cosl(__imag__ z);
|
|
+ __imag__ ret = r_exponent * sinl(__imag__ z);
|
|
+
|
|
+ return ret;
|
|
+}
|
|
+libm_hidden_def(cexpl)
|
|
+#endif
|
|
--
|
|
1.7.7.1
|
|
|