This is the mail archive of the
newlib@sourceware.org
mailing list for the newlib project.
[PATCH] libm/complex: add clog10, clog10f
- From: Yaakov Selkowitz <yselkowi at redhat dot com>
- To: newlib at sourceware dot org
- Date: Thu, 18 Feb 2016 10:53:09 -0600
- Subject: [PATCH] libm/complex: add clog10, clog10f
- Authentication-results: sourceware.org; auth=none
newlib/
* libc/include/complex.h (clog10, clog10f): Declare.
* libm/complex/Makefile.am (src): Add clog10.c.
(fsrc): Add clog10f.c.
(CHEWOUT_FILES): Add clog10.def.
* libm/complex/Makefile.in: Regenerate.
* libm/complex/clog10.c: New file.
* libm/complex/clog10f.c: New file.
* libm/complex/complex.tex: Add references to clog10 docs.
Signed-off-by: Yaakov Selkowitz <yselkowi@redhat.com>
---
newlib/libc/include/complex.h | 5 ++++
newlib/libm/complex/Makefile.am | 6 ++---
newlib/libm/complex/Makefile.in | 54 +++++++++++++++++++++++++----------------
newlib/libm/complex/clog10.c | 46 +++++++++++++++++++++++++++++++++++
newlib/libm/complex/clog10f.c | 15 ++++++++++++
newlib/libm/complex/complex.tex | 3 +++
6 files changed, 105 insertions(+), 24 deletions(-)
create mode 100644 newlib/libm/complex/clog10.c
create mode 100644 newlib/libm/complex/clog10f.c
diff --git a/newlib/libc/include/complex.h b/newlib/libc/include/complex.h
index c3ae3cf..3150f34 100644
--- a/newlib/libc/include/complex.h
+++ b/newlib/libc/include/complex.h
@@ -120,6 +120,11 @@ double creal(double complex);
float crealf(float complex);
long double creall(long double complex);
+#if __GNU_VISIBLE
+double complex clog10(double complex);
+float complex clog10f(float complex);
+#endif
+
__END_DECLS
#endif /* ! _COMPLEX_H */
diff --git a/newlib/libm/complex/Makefile.am b/newlib/libm/complex/Makefile.am
index d3db959..6ba49b6 100644
--- a/newlib/libm/complex/Makefile.am
+++ b/newlib/libm/complex/Makefile.am
@@ -6,14 +6,14 @@ INCLUDES = -I$(srcdir)/../common -I$(srcdir)/ $(NEWLIB_CFLAGS) $(CROSS_CFLAGS) $
src = cabs.c cacos.c cacosh.c carg.c casin.c casinh.c \
catan.c catanh.c ccos.c ccosh.c cephes_subr.c \
- cexp.c cimag.c clog.c conj.c \
+ cexp.c cimag.c clog.c clog10.c conj.c \
cpow.c cproj.c creal.c \
csin.c csinh.c csqrt.c ctan.c ctanh.c
lsrc = cabsl.c creall.c cimagl.c
fsrc = cabsf.c casinf.c ccosf.c cimagf.c cprojf.c \
- csqrtf.c cacosf.c casinhf.c ccoshf.c clogf.c \
+ csqrtf.c cacosf.c casinhf.c ccoshf.c clogf.c clog10f.c \
crealf.c ctanf.c cacoshf.c catanf.c \
cephes_subrf.c conjf.c csinf.c ctanhf.c \
cargf.c catanhf.c cexpf.c cpowf.c csinhf.c
@@ -36,7 +36,7 @@ include $(srcdir)/../../Makefile.shared
CHEWOUT_FILES = cabs.def cacos.def cacosh.def carg.def \
casin.def casinh.def catan.def catanh.def \
ccos.def ccosh.def cexp.def cimag.def clog.def \
- conj.def cpow.def cproj.def creal.def \
+ clog10.def conj.def cpow.def cproj.def creal.def \
csin.def csinh.def csqrt.def ctan.def ctanh.def
CHAPTERS = complex.tex
diff --git a/newlib/libm/complex/Makefile.in b/newlib/libm/complex/Makefile.in
index 49e35ac..0954736 100644
--- a/newlib/libm/complex/Makefile.in
+++ b/newlib/libm/complex/Makefile.in
@@ -79,23 +79,23 @@ am__objects_1 = lib_a-cabs.$(OBJEXT) lib_a-cacos.$(OBJEXT) \
lib_a-ccos.$(OBJEXT) lib_a-ccosh.$(OBJEXT) \
lib_a-cephes_subr.$(OBJEXT) lib_a-cexp.$(OBJEXT) \
lib_a-cimag.$(OBJEXT) lib_a-clog.$(OBJEXT) \
- lib_a-conj.$(OBJEXT) lib_a-cpow.$(OBJEXT) \
- lib_a-cproj.$(OBJEXT) lib_a-creal.$(OBJEXT) \
- lib_a-csin.$(OBJEXT) lib_a-csinh.$(OBJEXT) \
- lib_a-csqrt.$(OBJEXT) lib_a-ctan.$(OBJEXT) \
- lib_a-ctanh.$(OBJEXT)
+ lib_a-clog10.$(OBJEXT) lib_a-conj.$(OBJEXT) \
+ lib_a-cpow.$(OBJEXT) lib_a-cproj.$(OBJEXT) \
+ lib_a-creal.$(OBJEXT) lib_a-csin.$(OBJEXT) \
+ lib_a-csinh.$(OBJEXT) lib_a-csqrt.$(OBJEXT) \
+ lib_a-ctan.$(OBJEXT) lib_a-ctanh.$(OBJEXT)
am__objects_2 = lib_a-cabsf.$(OBJEXT) lib_a-casinf.$(OBJEXT) \
lib_a-ccosf.$(OBJEXT) lib_a-cimagf.$(OBJEXT) \
lib_a-cprojf.$(OBJEXT) lib_a-csqrtf.$(OBJEXT) \
lib_a-cacosf.$(OBJEXT) lib_a-casinhf.$(OBJEXT) \
lib_a-ccoshf.$(OBJEXT) lib_a-clogf.$(OBJEXT) \
- lib_a-crealf.$(OBJEXT) lib_a-ctanf.$(OBJEXT) \
- lib_a-cacoshf.$(OBJEXT) lib_a-catanf.$(OBJEXT) \
- lib_a-cephes_subrf.$(OBJEXT) lib_a-conjf.$(OBJEXT) \
- lib_a-csinf.$(OBJEXT) lib_a-ctanhf.$(OBJEXT) \
- lib_a-cargf.$(OBJEXT) lib_a-catanhf.$(OBJEXT) \
- lib_a-cexpf.$(OBJEXT) lib_a-cpowf.$(OBJEXT) \
- lib_a-csinhf.$(OBJEXT)
+ lib_a-clog10f.$(OBJEXT) lib_a-crealf.$(OBJEXT) \
+ lib_a-ctanf.$(OBJEXT) lib_a-cacoshf.$(OBJEXT) \
+ lib_a-catanf.$(OBJEXT) lib_a-cephes_subrf.$(OBJEXT) \
+ lib_a-conjf.$(OBJEXT) lib_a-csinf.$(OBJEXT) \
+ lib_a-ctanhf.$(OBJEXT) lib_a-cargf.$(OBJEXT) \
+ lib_a-catanhf.$(OBJEXT) lib_a-cexpf.$(OBJEXT) \
+ lib_a-cpowf.$(OBJEXT) lib_a-csinhf.$(OBJEXT)
am__objects_3 = lib_a-cabsl.$(OBJEXT) lib_a-creall.$(OBJEXT) \
lib_a-cimagl.$(OBJEXT)
@USE_LIBTOOL_FALSE@am_lib_a_OBJECTS = $(am__objects_1) \
@@ -105,13 +105,13 @@ LTLIBRARIES = $(noinst_LTLIBRARIES)
libcomplex_la_LIBADD =
am__objects_4 = cabs.lo cacos.lo cacosh.lo carg.lo casin.lo casinh.lo \
catan.lo catanh.lo ccos.lo ccosh.lo cephes_subr.lo cexp.lo \
- cimag.lo clog.lo conj.lo cpow.lo cproj.lo creal.lo csin.lo \
- csinh.lo csqrt.lo ctan.lo ctanh.lo
+ cimag.lo clog.lo clog10.lo conj.lo cpow.lo cproj.lo creal.lo \
+ csin.lo csinh.lo csqrt.lo ctan.lo ctanh.lo
am__objects_5 = cabsf.lo casinf.lo ccosf.lo cimagf.lo cprojf.lo \
- csqrtf.lo cacosf.lo casinhf.lo ccoshf.lo clogf.lo crealf.lo \
- ctanf.lo cacoshf.lo catanf.lo cephes_subrf.lo conjf.lo \
- csinf.lo ctanhf.lo cargf.lo catanhf.lo cexpf.lo cpowf.lo \
- csinhf.lo
+ csqrtf.lo cacosf.lo casinhf.lo ccoshf.lo clogf.lo clog10f.lo \
+ crealf.lo ctanf.lo cacoshf.lo catanf.lo cephes_subrf.lo \
+ conjf.lo csinf.lo ctanhf.lo cargf.lo catanhf.lo cexpf.lo \
+ cpowf.lo csinhf.lo
am__objects_6 = cabsl.lo creall.lo cimagl.lo
@USE_LIBTOOL_TRUE@am_libcomplex_la_OBJECTS = $(am__objects_4) \
@USE_LIBTOOL_TRUE@ $(am__objects_5) $(am__objects_6)
@@ -272,13 +272,13 @@ AUTOMAKE_OPTIONS = cygnus
INCLUDES = -I$(srcdir)/../common -I$(srcdir)/ $(NEWLIB_CFLAGS) $(CROSS_CFLAGS) $(TARGET_CFLAGS)
src = cabs.c cacos.c cacosh.c carg.c casin.c casinh.c \
catan.c catanh.c ccos.c ccosh.c cephes_subr.c \
- cexp.c cimag.c clog.c conj.c \
+ cexp.c cimag.c clog.c clog10.c conj.c \
cpow.c cproj.c creal.c \
csin.c csinh.c csqrt.c ctan.c ctanh.c
lsrc = cabsl.c creall.c cimagl.c
fsrc = cabsf.c casinf.c ccosf.c cimagf.c cprojf.c \
- csqrtf.c cacosf.c casinhf.c ccoshf.c clogf.c \
+ csqrtf.c cacosf.c casinhf.c ccoshf.c clogf.c clog10f.c \
crealf.c ctanf.c cacoshf.c catanf.c \
cephes_subrf.c conjf.c csinf.c ctanhf.c \
cargf.c catanhf.c cexpf.c cpowf.c csinhf.c
@@ -301,7 +301,7 @@ CLEANFILES = $(CHEWOUT_FILES) *.ref
CHEWOUT_FILES = cabs.def cacos.def cacosh.def carg.def \
casin.def casinh.def catan.def catanh.def \
ccos.def ccosh.def cexp.def cimag.def clog.def \
- conj.def cpow.def cproj.def creal.def \
+ clog10.def conj.def cpow.def cproj.def creal.def \
csin.def csinh.def csqrt.def ctan.def ctanh.def
CHAPTERS = complex.tex
@@ -458,6 +458,12 @@ lib_a-clog.o: clog.c
lib_a-clog.obj: clog.c
$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-clog.obj `if test -f 'clog.c'; then $(CYGPATH_W) 'clog.c'; else $(CYGPATH_W) '$(srcdir)/clog.c'; fi`
+lib_a-clog10.o: clog10.c
+ $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-clog10.o `test -f 'clog10.c' || echo '$(srcdir)/'`clog10.c
+
+lib_a-clog10.obj: clog10.c
+ $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-clog10.obj `if test -f 'clog10.c'; then $(CYGPATH_W) 'clog10.c'; else $(CYGPATH_W) '$(srcdir)/clog10.c'; fi`
+
lib_a-conj.o: conj.c
$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-conj.o `test -f 'conj.c' || echo '$(srcdir)/'`conj.c
@@ -572,6 +578,12 @@ lib_a-clogf.o: clogf.c
lib_a-clogf.obj: clogf.c
$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-clogf.obj `if test -f 'clogf.c'; then $(CYGPATH_W) 'clogf.c'; else $(CYGPATH_W) '$(srcdir)/clogf.c'; fi`
+lib_a-clog10f.o: clog10f.c
+ $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-clog10f.o `test -f 'clog10f.c' || echo '$(srcdir)/'`clog10f.c
+
+lib_a-clog10f.obj: clog10f.c
+ $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-clog10f.obj `if test -f 'clog10f.c'; then $(CYGPATH_W) 'clog10f.c'; else $(CYGPATH_W) '$(srcdir)/clog10f.c'; fi`
+
lib_a-crealf.o: crealf.c
$(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(lib_a_CFLAGS) $(CFLAGS) -c -o lib_a-crealf.o `test -f 'crealf.c' || echo '$(srcdir)/'`crealf.c
diff --git a/newlib/libm/complex/clog10.c b/newlib/libm/complex/clog10.c
new file mode 100644
index 0000000..c857a65
--- /dev/null
+++ b/newlib/libm/complex/clog10.c
@@ -0,0 +1,46 @@
+/*
+FUNCTION
+ <<clog10>>, <<clog10f>>---complex base-10 logarithm
+
+INDEX
+ clog10
+INDEX
+ clog10f
+
+ANSI_SYNOPSIS
+ #define _GNU_SOURCE
+ #include <complex.h>
+ double complex clog10(double complex <[z]>);
+ float complex clog10f(float complex <[z]>);
+
+
+DESCRIPTION
+ These functions compute the complex base-10 logarithm of <[z]>.
+ <<clog10>> is equivalent to <<clog>>(<[z]>)/<<log>>(10).
+
+ <<clog10f>> is identical to <<clog10>>, except that it performs
+ its calculations on <<floats complex>>.
+
+RETURNS
+ The clog10 functions return the complex base-10 logarithm value.
+
+PORTABILITY
+ <<clog10>> and <<clog10f>> are GNU extensions.
+
+*/
+
+#include <complex.h>
+#include <math.h>
+
+double complex
+clog10(double complex z)
+{
+ double complex w;
+ double p, rr;
+
+ rr = cabs(z);
+ p = log10(rr);
+ rr = atan2(cimag(z), creal(z)) * M_IVLN10;
+ w = p + rr * I;
+ return w;
+}
diff --git a/newlib/libm/complex/clog10f.c b/newlib/libm/complex/clog10f.c
new file mode 100644
index 0000000..98cf606
--- /dev/null
+++ b/newlib/libm/complex/clog10f.c
@@ -0,0 +1,15 @@
+#include <complex.h>
+#include <math.h>
+
+float complex
+clog10f(float complex z)
+{
+ float complex w;
+ float p, rr;
+
+ rr = cabsf(z);
+ p = log10f(rr);
+ rr = atan2f(cimagf(z), crealf(z)) * M_IVLN10;
+ w = p + rr * I;
+ return w;
+}
diff --git a/newlib/libm/complex/complex.tex b/newlib/libm/complex/complex.tex
index ec3258a..e91937a 100644
--- a/newlib/libm/complex/complex.tex
+++ b/newlib/libm/complex/complex.tex
@@ -19,6 +19,7 @@ Functions and documentations are taken from NetBSD.
* cexp:: Complex exponent
* cimag:: Imaginary part
* clog:: Complex natural (base e) logarithm
+* clog10:: Complex base 10 logarithm
* conj:: Complex conjugate
* cpow:: Complex power function
* cproj:: Compute a projection on the Riemann sphere
@@ -57,6 +58,8 @@ Functions and documentations are taken from NetBSD.
@page
@include complex/clog.def
@page
+@include complex/clog10.def
+@page
@include complex/conj.def
@page
@include complex/cpow.def
--
2.7.0