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
));