int lines_intersect(double p1x, double p1y, double p2x, double p2y, double q1x,
                    double q1y, double q2x, double q2y) {
  double tmp;
  if (p1x == p2x)
    p1x -= 0.0001;
  else if (p1x > p2x) {
    tmp = p1x;
    p1x = p2x;
    p2x = tmp;
    tmp = p1y;
    p1y = p2y;
    p2y = tmp;
  }
  if (q1x == q2x)
    q2x += 0.0001;
  else if (q1x > q2x) {
    tmp = q1x;
    q1x = q2x;
    q2x = tmp;
    tmp = q1y;
    q1y = q2y;
    q2y = tmp;
  }
  double pa = (p2y - p1y) / (p2x - p1x), pc = p1y - pa * p1x;
  double qa = (q2y - q1y) / (q2x - q1x), qc = q1y - qa * q1x;
  if (pa == qa)
    return 0;
  double x = (qc - pc) / (pa - qa);
  return x + 0.0001 > p1x && x - 0.0001 < p2x && x + 0.0001 > q1x &&
         x - 0.0001 < q2x;
}


x86 assembly (clang, -O0):
<lines_intersect>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
movsd  %xmm4,-0x30(%rbp)
movsd  %xmm5,-0x38(%rbp)
movsd  %xmm6,-0x40(%rbp)
movsd  %xmm7,-0x48(%rbp)
movsd  -0x10(%rbp),%xmm0
ucomisd -0x20(%rbp),%xmm0
jne    115d <lines_intersect+0x5d>
jp     115d <lines_intersect+0x5d>
movsd  -0x10(%rbp),%xmm0
movsd  0xeb1(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x10(%rbp)
jmp    11ae <lines_intersect+0xae>
movsd  -0x10(%rbp),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    11a9 <lines_intersect+0xa9>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11ae <lines_intersect+0xae>
movsd  -0x30(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jne    11db <lines_intersect+0xdb>
jp     11db <lines_intersect+0xdb>
movsd  0xe34(%rip),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    122c <lines_intersect+0x12c>
movsd  -0x30(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1227 <lines_intersect+0x127>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    122c <lines_intersect+0x12c>
movsd  -0x28(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
subsd  -0x30(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x58(%rbp),%xmm0
ucomisd -0x68(%rbp),%xmm0
jne    12ee <lines_intersect+0x1ee>
jp     12ee <lines_intersect+0x1ee>
movl   $0x0,-0x4(%rbp)
jmp    1395 <lines_intersect+0x295>
movsd  -0x70(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
subsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xced(%rip),%xmm0
addsd  -0x78(%rbp),%xmm0
xor    %eax,%eax
ucomisd -0x10(%rbp),%xmm0
mov    %al,-0x79(%rbp)
jbe    138a <lines_intersect+0x28a>
movsd  -0x78(%rbp),%xmm1
movsd  0xccb(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
xor    %eax,%eax
ucomisd %xmm1,%xmm0
mov    %al,-0x79(%rbp)
jbe    138a <lines_intersect+0x28a>
movsd  0xcab(%rip),%xmm0
addsd  -0x78(%rbp),%xmm0
xor    %eax,%eax
ucomisd -0x30(%rbp),%xmm0
mov    %al,-0x79(%rbp)
jbe    138a <lines_intersect+0x28a>
movsd  -0x78(%rbp),%xmm1
movsd  0xc89(%rip),%xmm0
subsd  %xmm0,%xmm1
movsd  -0x40(%rbp),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x79(%rbp)
mov    -0x79(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<lines_intersect>:
ucomisd %xmm2,%xmm0
jne    1112 <lines_intersect+0x12>
jp     1112 <lines_intersect+0x12>
addsd  0xef0(%rip),%xmm0
jmp    1124 <lines_intersect+0x24>
jbe    1124 <lines_intersect+0x24>
movapd %xmm1,%xmm9
movapd %xmm0,%xmm8
movapd %xmm3,%xmm1
jmp    1132 <lines_intersect+0x32>
movapd %xmm3,%xmm9
movapd %xmm2,%xmm8
movapd %xmm0,%xmm2
ucomisd %xmm6,%xmm4
jne    1144 <lines_intersect+0x44>
jp     1144 <lines_intersect+0x44>
addsd  0xec6(%rip),%xmm6
jmp    1154 <lines_intersect+0x54>
jbe    1154 <lines_intersect+0x54>
movapd %xmm6,%xmm0
movapd %xmm7,%xmm3
movapd %xmm5,%xmm7
jmp    1160 <lines_intersect+0x60>
movapd %xmm4,%xmm0
movapd %xmm5,%xmm3
movapd %xmm6,%xmm4
unpcklpd %xmm1,%xmm3
unpcklpd %xmm9,%xmm7
subpd  %xmm3,%xmm7
movapd %xmm4,%xmm5
unpcklpd %xmm8,%xmm5
movapd %xmm0,%xmm1
unpcklpd %xmm2,%xmm1
subpd  %xmm1,%xmm5
divpd  %xmm5,%xmm7
movapd %xmm7,%xmm5
unpckhpd %xmm7,%xmm5
xor    %eax,%eax
ucomisd %xmm7,%xmm5
jne    1199 <lines_intersect+0x99>
jp     1199 <lines_intersect+0x99>
ret
mulpd  %xmm7,%xmm1
subpd  %xmm1,%xmm3
movapd %xmm7,%xmm1
movsd  %xmm3,%xmm1
shufpd $0x1,%xmm7,%xmm3
subpd  %xmm3,%xmm1
movapd %xmm1,%xmm3
unpckhpd %xmm1,%xmm3
divsd  %xmm3,%xmm1
movsd  0xe42(%rip),%xmm3
addsd  %xmm1,%xmm3
ucomisd %xmm2,%xmm3
jbe    1198 <lines_intersect+0x98>
addsd  0xe28(%rip),%xmm1
ucomisd %xmm1,%xmm8
jbe    1198 <lines_intersect+0x98>
ucomisd %xmm0,%xmm3
jbe    1198 <lines_intersect+0x98>
xor    %eax,%eax
ucomisd %xmm1,%xmm4
seta   %al
ret

LLVM IR (-O0):
define i32 @lines_intersect(double %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7) {
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
store double %0, double* %10
store double %1, double* %11
store double %2, double* %12
store double %3, double* %13
store double %4, double* %14
store double %5, double* %15
store double %6, double* %16
store double %7, double* %17
%24 = load double, double* %10
%25 = load double, double* %12
%26 = fcmp oeq double %24, %25
br i1 %26, label %27, label %30
27: ; preds = %8
%28 = load double, double* %10
%29 = fsub double %28, 1.000000e-04
store double %29, double* %10
br label %42
30: ; preds = %8
%31 = load double, double* %10
%32 = load double, double* %12
%33 = fcmp ogt double %31, %32
br i1 %33, label %34, label %41
34: ; preds = %30
%35 = load double, double* %10
store double %35, double* %18
%36 = load double, double* %12
store double %36, double* %10
%37 = load double, double* %18
store double %37, double* %12
%38 = load double, double* %11
store double %38, double* %18
%39 = load double, double* %13
store double %39, double* %11
%40 = load double, double* %18
store double %40, double* %13
br label %41
41: ; preds = %34, %30
br label %42
42: ; preds = %41, %27
%43 = load double, double* %14
%44 = load double, double* %16
%45 = fcmp oeq double %43, %44
br i1 %45, label %46, label %49
46: ; preds = %42
%47 = load double, double* %16
%48 = fadd double %47, 1.000000e-04
store double %48, double* %16
br label %61
49: ; preds = %42
%50 = load double, double* %14
%51 = load double, double* %16
%52 = fcmp ogt double %50, %51
br i1 %52, label %53, label %60
53: ; preds = %49
%54 = load double, double* %14
store double %54, double* %18
%55 = load double, double* %16
store double %55, double* %14
%56 = load double, double* %18
store double %56, double* %16
%57 = load double, double* %15
store double %57, double* %18
%58 = load double, double* %17
store double %58, double* %15
%59 = load double, double* %18
store double %59, double* %17
br label %60
60: ; preds = %53, %49
br label %61
61: ; preds = %60, %46
%62 = load double, double* %13
%63 = load double, double* %11
%64 = fsub double %62, %63
%65 = load double, double* %12
%66 = load double, double* %10
%67 = fsub double %65, %66
%68 = fdiv double %64, %67
store double %68, double* %19
%69 = load double, double* %11
%70 = load double, double* %19
%71 = load double, double* %10
%72 = fneg double %70
%73 = call double @llvm.fmuladd.f64(double %72, double %71, double %69)
store double %73, double* %20
%74 = load double, double* %17
%75 = load double, double* %15
%76 = fsub double %74, %75
%77 = load double, double* %16
%78 = load double, double* %14
%79 = fsub double %77, %78
%80 = fdiv double %76, %79
store double %80, double* %21
%81 = load double, double* %15
%82 = load double, double* %21
%83 = load double, double* %14
%84 = fneg double %82
%85 = call double @llvm.fmuladd.f64(double %84, double %83, double %81)
store double %85, double* %22
%86 = load double, double* %19
%87 = load double, double* %21
%88 = fcmp oeq double %86, %87
br i1 %88, label %89, label %90
89: ; preds = %61
store i32 0, i32* %9
br label %120
90: ; preds = %61
%91 = load double, double* %22
%92 = load double, double* %20
%93 = fsub double %91, %92
%94 = load double, double* %19
%95 = load double, double* %21
%96 = fsub double %94, %95
%97 = fdiv double %93, %96
store double %97, double* %23
%98 = load double, double* %23
%99 = fadd double %98, 1.000000e-04
%100 = load double, double* %10
%101 = fcmp ogt double %99, %100
br i1 %101, label %102, label %117
102: ; preds = %90
%103 = load double, double* %23
%104 = fsub double %103, 1.000000e-04
%105 = load double, double* %12
%106 = fcmp olt double %104, %105
br i1 %106, label %107, label %117
107: ; preds = %102
%108 = load double, double* %23
%109 = fadd double %108, 1.000000e-04
%110 = load double, double* %14
%111 = fcmp ogt double %109, %110
br i1 %111, label %112, label %117
112: ; preds = %107
%113 = load double, double* %23
%114 = fsub double %113, 1.000000e-04
%115 = load double, double* %16
%116 = fcmp olt double %114, %115
br label %117
117: ; preds = %112, %107, %102, %90
%118 = phi i1 [ false, %107 ], [ false, %102 ], [ false, %90 ], [ %116, %112 ]
%119 = zext i1 %118 to i32
store i32 %119, i32* %9
br label %120
120: ; preds = %117, %89
%121 = load i32, i32* %9
ret i32 %121
}

LLVM IR (-O2):
define i32 @lines_intersect(double %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7) {
%9 = fcmp oeq double %0, %2
br i1 %9, label %10, label %12
10: ; preds = %8
%11 = fadd double %0, -1.000000e-04
br label %15
12: ; preds = %8
%13 = fcmp ogt double %0, %2
br i1 %13, label %14, label %15
14: ; preds = %12
br label %15
15: ; preds = %12, %14, %10
%16 = phi double [ %3, %10 ], [ %1, %14 ], [ %3, %12 ]
%17 = phi double [ %2, %10 ], [ %0, %14 ], [ %2, %12 ]
%18 = phi double [ %1, %10 ], [ %3, %14 ], [ %1, %12 ]
%19 = phi double [ %11, %10 ], [ %2, %14 ], [ %0, %12 ]
%20 = fcmp oeq double %4, %6
br i1 %20, label %21, label %23
21: ; preds = %15
%22 = fadd double %6, 1.000000e-04
br label %26
23: ; preds = %15
%24 = fcmp ogt double %4, %6
br i1 %24, label %25, label %26
25: ; preds = %23
br label %26
26: ; preds = %23, %25, %21
%27 = phi double [ %4, %21 ], [ %6, %25 ], [ %4, %23 ]
%28 = phi double [ %5, %21 ], [ %7, %25 ], [ %5, %23 ]
%29 = phi double [ %22, %21 ], [ %4, %25 ], [ %6, %23 ]
%30 = phi double [ %7, %21 ], [ %5, %25 ], [ %7, %23 ]
%31 = insertelement <2 x double> poison, double %16, i64 1
%32 = insertelement <2 x double> %31, double %30, i64 0
%33 = insertelement <2 x double> poison, double %18, i64 1
%34 = insertelement <2 x double> %33, double %28, i64 0
%35 = fsub <2 x double> %32, %34
%36 = insertelement <2 x double> poison, double %17, i64 1
%37 = insertelement <2 x double> %36, double %29, i64 0
%38 = insertelement <2 x double> poison, double %19, i64 1
%39 = insertelement <2 x double> %38, double %27, i64 0
%40 = fsub <2 x double> %37, %39
%41 = fdiv <2 x double> %35, %40
%42 = extractelement <2 x double> %41, i64 0
%43 = extractelement <2 x double> %41, i64 1
%44 = fcmp oeq double %43, %42
br i1 %44, label %64, label %45
45: ; preds = %26
%46 = fneg <2 x double> %41
%47 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %46, <2 x double> %39, <2 x double> %34)
%48 = shufflevector <2 x double> %47, <2 x double> %41, <2 x i32> <i32 0, i32 3>
%49 = shufflevector <2 x double> %47, <2 x double> %41, <2 x i32> <i32 1, i32 2>
%50 = fsub <2 x double> %48, %49
%51 = extractelement <2 x double> %50, i64 0
%52 = extractelement <2 x double> %50, i64 1
%53 = fdiv double %51, %52
%54 = fadd double %53, 1.000000e-04
%55 = fcmp ogt double %54, %19
br i1 %55, label %56, label %64
56: ; preds = %45
%57 = fadd double %53, -1.000000e-04
%58 = fcmp olt double %57, %17
%59 = fcmp ogt double %54, %27
%60 = and i1 %58, %59
br i1 %60, label %61, label %64
61: ; preds = %56
%62 = fcmp olt double %57, %29
%63 = zext i1 %62 to i32
br label %64
64: ; preds = %45, %56, %61, %26
%65 = phi i32 [ 0, %26 ], [ 0, %56 ], [ 0, %45 ], [ %63, %61 ]
ret i32 %65
}