1 extern double fabs (double);
2 extern double sqrt (double);
4 typedef struct complexm {
9 setCom (double r, double i)
12 ct.re=fabs(r)<1E-300?0.0:r;
13 ct.im=fabs(i)<1E-300?0.0:i;
18 csqrt_crash (double x)
20 return (x>=0) ? setCom(sqrt(x),0) : setCom(0,sqrt(-x));