11 double foo (S *, S *);
12 void bar (S *, S *, S *, double &, double &, double &);
17 return __builtin_sqrt ((a1->x - a2->x) * (a1->x - a2->x)
18 + (a1->y - a2->y) * (a1->y - a2->y)
19 + (a1->z - a2->z) * (a1->z - a2->z));
23 bar (S *p, S *q, S *r, double &x, double &y, double &z)
25 if (foo (p, q) == 0.0)
32 if (foo (p, r) == 0.0)
39 if (foo (q, r) == 0.0)
47 double a1, b1, c1, d1, e1;
48 double dx, dy, dz, dw, dv;
53 e1 = __builtin_sqrt (a1 * a1 + b1 * b1 + c1 * c1);
60 dw = dx * dx + dy * dy + dz * dz;
61 dv = 2.0 * dx * a1 + 2.0 * dy * b1 + 2.0 * dz * c1;
72 S a, b, c, d, *p, *q, *r;
86 bar (p, q, r, d.x, d.y, d.z);