summaryrefslogblamecommitdiff
path: root/src/complex/casin.c
blob: 79aff2789de96c0e8728fb0d3b6b25f9d114c800 (plain) (tree)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16















                                                   
#include "libm.h"

// FIXME

/* asin(z) = -i log(i z + sqrt(1 - z*z)) */

double complex casin(double complex z)
{
	double complex w;
	double x, y;

	x = creal(z);
	y = cimag(z);
	w = cpack(1.0 - (x - y)*(x + y), -2.0*x*y);
	return clog(cpack(-y, x) + csqrt(w));
}