int within_distance_cutoff(double x1, double y1, double z1, double dim1,
                           double x2, double y2, double z2, double dim2,
                           double ext) {
  double d1 = dim1 + ext, d2 = -(dim2 + ext);
  if ((x2 - x1 >= d1) || (x2 - x1 <= d2) || (y2 - y1 >= d1) ||
      (y2 - y1 <= d2) || (z2 - z1 >= d1) || (z2 - z1 <= d2))
    return 0;
  d1 = 0.5 * (dim2 - dim1);
  d2 = (x2 - x1 + d1) * (x2 - x1 + d1) + (y2 - y1 + d1) * (y2 - y1 + d1) +
       (z2 - z1 + d1) * (z2 - z1 + d1);
  d1 = (0.86602540378443864676372317075294 * (dim1 + dim2) + ext);
  d1 *= d1;
  return (d2 < d1);
}


x86 assembly (clang, -O0):
<within_distance_cutoff>:
push   %rbp
mov    %rsp,%rbp
movsd  0x10(%rbp),%xmm8
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  -0x28(%rbp),%xmm0
addsd  0x10(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
ucomisd -0x50(%rbp),%xmm0
jae    11f1 <within_distance_cutoff+0xf1>
movsd  -0x30(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jae    11f1 <within_distance_cutoff+0xf1>
movsd  -0x38(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
ucomisd -0x50(%rbp),%xmm0
jae    11f1 <within_distance_cutoff+0xf1>
movsd  -0x38(%rbp),%xmm1
subsd  -0x18(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jae    11f1 <within_distance_cutoff+0xf1>
movsd  -0x40(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
ucomisd -0x50(%rbp),%xmm0
jae    11f1 <within_distance_cutoff+0xf1>
movsd  -0x40(%rbp),%xmm1
subsd  -0x20(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     11fd <within_distance_cutoff+0xfd>
movl   $0x0,-0x4(%rbp)
jmp    12d7 <within_distance_cutoff+0x1d7>
movsd  -0x48(%rbp),%xmm1
subsd  -0x28(%rbp),%xmm1
movsd  0xdf9(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
addsd  -0x50(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm2
subsd  -0x10(%rbp),%xmm2
addsd  -0x50(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm3
subsd  -0x18(%rbp),%xmm3
addsd  -0x50(%rbp),%xmm3
mulsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x40(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
subsd  -0x20(%rbp),%xmm2
addsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  0x10(%rbp),%xmm1
movsd  0xd5e(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<within_distance_cutoff>:
movsd  0x8(%rsp),%xmm8
movapd %xmm3,%xmm9
addsd  %xmm8,%xmm9
subsd  %xmm0,%xmm4
xor    %eax,%eax
ucomisd %xmm9,%xmm4
jae    11a9 <within_distance_cutoff+0xa9>
movapd %xmm7,%xmm0
addsd  %xmm8,%xmm0
xorpd  0xecd(%rip),%xmm0
ucomisd %xmm4,%xmm0
jae    11a9 <within_distance_cutoff+0xa9>
subsd  %xmm1,%xmm5
ucomisd %xmm9,%xmm5
jae    11a9 <within_distance_cutoff+0xa9>
ucomisd %xmm5,%xmm0
jae    11a9 <within_distance_cutoff+0xa9>
subsd  %xmm2,%xmm6
ucomisd %xmm9,%xmm6
jae    11a9 <within_distance_cutoff+0xa9>
ucomisd %xmm6,%xmm0
jae    11a9 <within_distance_cutoff+0xa9>
movapd %xmm7,%xmm0
subsd  %xmm3,%xmm0
mulsd  0xea5(%rip),%xmm0
addsd  %xmm0,%xmm4
addsd  %xmm0,%xmm5
mulsd  %xmm5,%xmm5
mulsd  %xmm4,%xmm4
addsd  %xmm5,%xmm4
addsd  %xmm0,%xmm6
mulsd  %xmm6,%xmm6
addsd  %xmm7,%xmm3
mulsd  0xe85(%rip),%xmm3
addsd  %xmm4,%xmm6
addsd  %xmm8,%xmm3
mulsd  %xmm3,%xmm3
xor    %eax,%eax
ucomisd %xmm6,%xmm3
seta   %al
ret

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

LLVM IR (-O2):
define i32 @within_distance_cutoff(double %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8) {
%10 = fadd double %3, %8
%11 = fadd double %7, %8
%12 = fneg double %11
%13 = fsub double %4, %0
%14 = fcmp ult double %13, %10
%15 = fcmp ugt double %13, %12
%16 = and i1 %14, %15
br i1 %16, label %17, label %41
17: ; preds = %9
%18 = fsub double %5, %1
%19 = fcmp ult double %18, %10
%20 = fcmp ugt double %18, %12
%21 = and i1 %19, %20
br i1 %21, label %22, label %41
22: ; preds = %17
%23 = fsub double %6, %2
%24 = fcmp ult double %23, %10
%25 = fcmp ugt double %23, %12
%26 = and i1 %24, %25
br i1 %26, label %27, label %41
27: ; preds = %22
%28 = fsub double %7, %3
%29 = fmul double %28, 5.000000e-01
%30 = fadd double %13, %29
%31 = fadd double %18, %29
%32 = fmul double %31, %31
%33 = tail call double @llvm.fmuladd.f64(double %30, double %30, double %32)
%34 = fadd double %23, %29
%35 = tail call double @llvm.fmuladd.f64(double %34, double %34, double %33)
%36 = fadd double %3, %7
%37 = tail call double @llvm.fmuladd.f64(double %36, double 0x3FEBB67AE8584CAA, double %8)
%38 = fmul double %37, %37
%39 = fcmp olt double %35, %38
%40 = zext i1 %39 to i32
br label %41
41: ; preds = %9, %17, %22, %27
%42 = phi i32 [ %40, %27 ], [ 0, %22 ], [ 0, %17 ], [ 0, %9 ]
ret i32 %42
}