7 #include "cl_complex.h"
15 const cl_N sqrt (const cl_N& x)
18 // x reell -> Für x>=0 klar, für x<0: sqrt(-x)*i.
20 // Bestimme r=abs(x)=sqrt(a*a+b*b).
21 // Falls a>=0: Setze c:=sqrt((r+a)/2), d:=(b/(2*c) falls c>0, c falls c=0).
22 // Falls a<0: Setze d:=sqrt((r-a)/2)*(1 falls b>=0, -1 falls b<0), c:=b/(2*d).
23 // Damit ist c>=0, 2*c*d=b, c*c=(r+a)/2, d*d=(r-a)/2, c*c-d*d=a, c*c+d*d=r,
24 // also c+di die gesuchte Wurzel.
30 return complex_C(0,sqrt(-x));
33 var const cl_R& a = realpart(x);
34 var const cl_R& b = imagpart(x);
35 var cl_R r = cl_hypot(a,b); // r = (abs x)
38 var cl_R c = sqrt((r+a)/2);
39 var cl_R d = (!zerop(c) ? b/(2*c) : c);
40 return complex_C(c,d);
42 var cl_R d = sqrt((r-a)/2);
46 return complex_C(c,d);