diff options
author | Danny Smith <dannysmith@users.sourceforge.net> | 2005-01-10 08:03:55 +0300 |
---|---|---|
committer | Danny Smith <dannysmith@users.sourceforge.net> | 2005-01-10 08:03:55 +0300 |
commit | 4073a59300a4411813a4a1c4de5e0c4e3d77ee09 (patch) | |
tree | be52f78e7dc02863f8ccfa8141421a8b5b1b437a /winsup/mingw/mingwex/complex/casinl.c | |
parent | 76eb61923900be1b6fd4f83bb62f48670fe7293c (diff) |
* mingwex/complex/(cabsl.c cacosl.c cacoshl.c cargl.c casinl.c
casinhl.c catanl.c catanhl.c ccosl.c ccoshl.c cexpl.c cimagl.c
clogl.c cpowl.c cprojl.c creall.c csinl.c csinhl.c csqrtl.c
ctanl.c ctanhl.c): New files.
* mingwex/Makefile.in (COMPLEX_DISTFILES): Adjust.
(COMPLEX_OBJS): Adjust.
* include/complex.h (cabsl, cacosl, cacoshl, cargl, casinf.
casinhl, catanl, catanhl, ccosl, ccoshl, cexpl, cimagl, clogl,
cpowl, cprojl, creall, csinl, csinhl, csqrtl, ctanl, ctanhl):
Declare.
Diffstat (limited to 'winsup/mingw/mingwex/complex/casinl.c')
-rwxr-xr-x | winsup/mingw/mingwex/complex/casinl.c | 48 |
1 files changed, 48 insertions, 0 deletions
diff --git a/winsup/mingw/mingwex/complex/casinl.c b/winsup/mingw/mingwex/complex/casinl.c new file mode 100755 index 000000000..3093fa542 --- /dev/null +++ b/winsup/mingw/mingwex/complex/casinl.c @@ -0,0 +1,48 @@ +/* + casinl.c + Contributed by Danny Smith + 2005-01-04 +*/ + +#include <math.h> +#include <complex.h> + +/* casin (Z ) = -I * clog(I * Z + csqrt (1.0 - Z * Z))) */ + +long double complex casinl (long double complex Z) +{ + long double complex Res; + long double x, y; + + x = __real__ Z; + y = __imag__ Z; + + if (y == 0.0L) + { + __real__ Res = asinl (x); + __imag__ Res = 0.0L; + } + else /* -I * clog(I * Z + csqrt(1.0 - Z * Z))) */ + { + long double complex ZZ; + + /* Z * Z = ((x - y) * (x + y)) + (2.0 * x * y) * I */ + /* calculate 1 - Z * Z */ + __real__ ZZ = 1.0L - (x - y) * (x + y); + __imag__ ZZ = -2.0L * x * y; + ZZ = csqrtl (ZZ); + + + /* add I * Z to ZZ */ + + __real__ ZZ -= y; + __imag__ ZZ += x; + + ZZ = clogl (ZZ); + + /* mult by -I */ + __real__ Res = __imag__ ZZ; + __imag__ Res = - __real__ ZZ; + } + return (Res); +} |