diff options
author | Corinna Vinschen <corinna@vinschen.de> | 2010-10-08 10:35:14 +0000 |
---|---|---|
committer | Corinna Vinschen <corinna@vinschen.de> | 2010-10-08 10:35:14 +0000 |
commit | 5af2a793bc8cc1a41783515ded2aee98d7e50644 (patch) | |
tree | d07f6b0863fe8607d1290980d9bd5ae0a75f2cff /newlib/libm/complex/ctan.c | |
parent | 1a91b8c0ac168f498576bd8833e8d4a7197fce57 (diff) | |
download | cygnal-5af2a793bc8cc1a41783515ded2aee98d7e50644.tar.gz cygnal-5af2a793bc8cc1a41783515ded2aee98d7e50644.tar.bz2 cygnal-5af2a793bc8cc1a41783515ded2aee98d7e50644.zip |
* libc/include/complex.h: New complex header.
* libm/common/fdlibm.h: Added ifdef _COMPLEX_H.
* libm/complex/*: New complex functions imported from NetBSD.
* libm/Makefile.am: Added complex subdir.
* libm/Makefile.in: Regenerate.
* libm/configure.in: Added complex subdir.
* libm/configure: Regenerate.
Diffstat (limited to 'newlib/libm/complex/ctan.c')
-rw-r--r-- | newlib/libm/complex/ctan.c | 60 |
1 files changed, 60 insertions, 0 deletions
diff --git a/newlib/libm/complex/ctan.c b/newlib/libm/complex/ctan.c new file mode 100644 index 000000000..33b39b497 --- /dev/null +++ b/newlib/libm/complex/ctan.c @@ -0,0 +1,60 @@ +/* $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> + */ + +#include <complex.h> +#include <math.h> +#include "cephes_subr.h" + +#define MAXNUM 1.0e308 + +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 = _ctans(z); + + if (d == 0.0) { + /* mtherr ("ctan", OVERFLOW); */ + w = MAXNUM + MAXNUM * I; + return w; + } + + w = sin(2.0 * creal(z)) / d + (sinh(2.0 * cimag(z)) / d) * I; + return w; +} |