void binlensd1(double x1, double x2, double d, double m1, double *y11,
               double *y12, double *y22) {
  double m2;
  double del11, del12, x2sq, n1, n2, n1sq, n2sq;
  double sum, sum11, sum12, sum22;
  m2 = 1.0 - m1;
  if (m2 < 0.0) {
    m2 = 0.0;
    m1 = 1.0;
  }
  del11 = x1 + m2 * d;
  del12 = x1 - m1 * d;
  x2sq = x2 * x2;
  n1 = del11 * del11 + x2sq;
  n2 = del12 * del12 + x2sq;
  n1sq = n1 * n1;
  n2sq = n2 * n2;
  sum = m1 / n1 + m2 / n2;
  sum11 = 2.0 * m1 * del11 * del11 / n1sq + 2.0 * m2 * del12 * del12 / n2sq;
  sum12 = 2.0 * m1 * del11 * x2 / n1sq + 2.0 * m2 * del12 * x2 / n2sq;
  sum22 = 2.0 * m1 * x2sq / n1sq + 2.0 * m2 * x2sq / n2sq;
  *y11 = 1.0 - sum + sum11;
  *y12 = sum12;
  *y22 = 1.0 - sum + sum22;
}


x86 assembly (clang, -O0):
<binlensd1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x18,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
mov    %rdi,-0x28(%rbp)
mov    %rsi,-0x30(%rbp)
mov    %rdx,-0x38(%rbp)
movsd  0xed0(%rip),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    115d <binlensd1+0x5d>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xea8(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x8(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x8(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x60(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x20(%rbp),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
divsd  -0x68(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  0xdd2(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  0xdb6(%rip),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
divsd  -0x78(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xd8e(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  0xd72(%rip),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
divsd  -0x78(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0xd4a(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  0xd33(%rip),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  -0x58(%rbp),%xmm1
divsd  -0x78(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xd08(%rip),%xmm0
subsd  -0x80(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  0xcdb(%rip),%xmm0
subsd  -0x80(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x18,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<binlensd1>:
movsd  0xef7(%rip),%xmm8
movapd %xmm8,%xmm6
subsd  %xmm3,%xmm6
xorpd  %xmm5,%xmm5
xorpd  %xmm11,%xmm11
maxsd  %xmm6,%xmm11
cmpltsd %xmm5,%xmm6
movapd %xmm6,%xmm5
andpd  %xmm8,%xmm5
andnpd %xmm3,%xmm6
orpd   %xmm5,%xmm6
movapd %xmm11,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm6,%xmm2
subsd  %xmm2,%xmm0
movapd %xmm3,%xmm2
movapd %xmm0,%xmm7
movapd %xmm6,%xmm9
movapd %xmm11,%xmm10
addsd  %xmm6,%xmm6
movapd %xmm6,%xmm4
mulsd  %xmm3,%xmm4
mulsd  %xmm4,%xmm3
addsd  %xmm11,%xmm11
movapd %xmm11,%xmm5
mulsd  %xmm0,%xmm5
mulsd  %xmm5,%xmm0
mulsd  %xmm1,%xmm4
mulsd  %xmm1,%xmm5
mulsd  %xmm1,%xmm1
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm7,%xmm7
addsd  %xmm1,%xmm7
divsd  %xmm2,%xmm9
mulsd  %xmm2,%xmm2
divsd  %xmm7,%xmm10
mulsd  %xmm7,%xmm7
addsd  %xmm9,%xmm10
divsd  %xmm2,%xmm3
divsd  %xmm7,%xmm0
addsd  %xmm3,%xmm0
divsd  %xmm2,%xmm4
divsd  %xmm7,%xmm5
addsd  %xmm4,%xmm5
mulsd  %xmm1,%xmm6
divsd  %xmm2,%xmm6
mulsd  %xmm1,%xmm11
divsd  %xmm7,%xmm11
addsd  %xmm6,%xmm11
subsd  %xmm10,%xmm8
addsd  %xmm8,%xmm0
movsd  %xmm0,(%rdi)
movsd  %xmm5,(%rsi)
addsd  %xmm11,%xmm8
movsd  %xmm8,(%rdx)
ret

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

LLVM IR (-O2):
define void @binlensd1(double %0, double %1, double %2, double %3, double* nocapture writeonly %4, double* nocapture writeonly %5, double* nocapture writeonly %6) {
%8 = fsub double 1.000000e+00, %3
%9 = fcmp olt double %8, 0.000000e+00
%10 = select i1 %9, double 0.000000e+00, double %8
%11 = select i1 %9, double 1.000000e+00, double %3
%12 = tail call double @llvm.fmuladd.f64(double %10, double %2, double %0)
%13 = fneg double %11
%14 = tail call double @llvm.fmuladd.f64(double %13, double %2, double %0)
%15 = fmul double %1, %1
%16 = tail call double @llvm.fmuladd.f64(double %12, double %12, double %15)
%17 = tail call double @llvm.fmuladd.f64(double %14, double %14, double %15)
%18 = fmul double %16, %16
%19 = fmul double %17, %17
%20 = fdiv double %11, %16
%21 = fdiv double %10, %17
%22 = fadd double %20, %21
%23 = fmul double %11, 2.000000e+00
%24 = fmul double %23, %12
%25 = fmul double %12, %24
%26 = fdiv double %25, %18
%27 = fmul double %10, 2.000000e+00
%28 = fmul double %27, %14
%29 = fmul double %14, %28
%30 = fdiv double %29, %19
%31 = fadd double %26, %30
%32 = fmul double %24, %1
%33 = fdiv double %32, %18
%34 = fmul double %28, %1
%35 = fdiv double %34, %19
%36 = fadd double %33, %35
%37 = fmul double %15, %23
%38 = fdiv double %37, %18
%39 = fmul double %15, %27
%40 = fdiv double %39, %19
%41 = fadd double %38, %40
%42 = fsub double 1.000000e+00, %22
%43 = fadd double %42, %31
store double %43, double* %4
store double %36, double* %5
%44 = fadd double %42, %41
store double %44, double* %6
ret void
}