diff options
| author | Vladimir Azarov <avm@intermediate-node.net> | 2024-10-01 15:47:05 +0200 |
|---|---|---|
| committer | Vladimir Azarov <avm@intermediate-node.net> | 2024-10-01 15:47:05 +0200 |
| commit | 4abab5ad6c8465a7528ccdd5f49367da05f78bbd (patch) | |
| tree | ebf009bf1376a5a223a915bc27cbbd791a1316bc /src/complex/casinl.c | |
Initial version
Diffstat (limited to 'src/complex/casinl.c')
| -rw-r--r-- | src/complex/casinl.c | 21 |
1 files changed, 21 insertions, 0 deletions
diff --git a/src/complex/casinl.c b/src/complex/casinl.c new file mode 100644 index 0000000..072adc4 --- /dev/null +++ b/src/complex/casinl.c @@ -0,0 +1,21 @@ +#include "complex_impl.h" + +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 +long double complex casinl(long double complex z) +{ + return casin(z); +} +#else +// FIXME +long double complex casinl(long double complex z) +{ + long double complex w; + long double x, y; + + x = creall(z); + y = cimagl(z); + w = CMPLXL(1.0 - (x - y)*(x + y), -2.0*x*y); + long double complex r = clogl(CMPLXL(-y, x) + csqrtl(w)); + return CMPLXL(cimagl(r), -creall(r)); +} +#endif |