Function line_line_intersect

code/geometry/primitives3d.cpp:64–71

    return
      abs((A - *this) * (B - *this) % (C - *this)) < EPS; } };
int line_line_intersect(L(A, B), L(C, D), point3d &O){
  if (abs((B - A) * (C - A) % (D - A)) > EPS) return 0;
  if (((A - B) * (C - D)).length() < EPS)
    return A.isOnLine(C, D) ? 2 : 0;
  point3d normal = ((A - B) * (C - B)).normalize();
  double s1 = (C - A) * (D - A) % normal;
  O = A + ((B - A) / (s1 + ((D - B) * (C - B) % normal))) * s1;
  return 1; }
int line_plane_intersect(L(A, B), PL(C, D, E), point3d & O) {
  double V1 = (C - A) * (D - A) % (E - A);
  double V2 = (D - B) * (C - B) % (E - B);