unsigned mandel_double(double cr, double ci) {
  double zr = cr, zi = ci;
  double tmp;
  unsigned i;
  for (i = 0; i < (1 << 14); i++) {
    tmp = zr * zr - zi * zi + cr;
    zi *= 2 * zr;
    zi += ci;
    zr = tmp;
    if (zr * zr + zi * zi > 4.0)
      break;
  }
  return i;
}


x86 assembly (clang, -O0):
<mandel_double>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
cmpl   $0x4000,-0x2c(%rbp)
jae    11e9 <mandel_double+0xe9>
movsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe8d(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe39(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11d6 <mandel_double+0xd6>
jmp    11e9 <mandel_double+0xe9>
jmp    11db <mandel_double+0xdb>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1129 <mandel_double+0x29>
mov    -0x2c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<mandel_double>:
xor    %eax,%eax
movsd  0xef6(%rip),%xmm2
movapd %xmm0,%xmm3
movapd %xmm1,%xmm4
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm3,%xmm5
movapd %xmm4,%xmm6
mulsd  %xmm4,%xmm6
mulsd  %xmm3,%xmm3
subsd  %xmm6,%xmm3
addsd  %xmm0,%xmm3
addsd  %xmm5,%xmm5
mulsd  %xmm5,%xmm4
addsd  %xmm1,%xmm4
movapd %xmm4,%xmm5
mulsd  %xmm4,%xmm5
movapd %xmm3,%xmm6
mulsd  %xmm3,%xmm6
addsd  %xmm5,%xmm6
ucomisd %xmm2,%xmm6
ja     116d <mandel_double+0x6d>
add    $0x1,%eax
cmp    $0x4000,%eax
jne    1120 <mandel_double+0x20>
mov    $0x4000,%eax
ret

LLVM IR (-O0):
define i32 @mandel_double(double %0, double %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca i32
store double %0, double* %3
store double %1, double* %4
%9 = load double, double* %3
store double %9, double* %5
%10 = load double, double* %4
store double %10, double* %6
store i32 0, i32* %8
br label %11
11: ; preds = %41, %2
%12 = load i32, i32* %8
%13 = icmp ult i32 %12, 16384
br i1 %13, label %14, label %44
14: ; preds = %11
%15 = load double, double* %5
%16 = load double, double* %5
%17 = load double, double* %6
%18 = load double, double* %6
%19 = fmul double %17, %18
%20 = fneg double %19
%21 = call double @llvm.fmuladd.f64(double %15, double %16, double %20)
%22 = load double, double* %3
%23 = fadd double %21, %22
store double %23, double* %7
%24 = load double, double* %5
%25 = fmul double 2.000000e+00, %24
%26 = load double, double* %6
%27 = fmul double %26, %25
store double %27, double* %6
%28 = load double, double* %4
%29 = load double, double* %6
%30 = fadd double %29, %28
store double %30, double* %6
%31 = load double, double* %7
store double %31, double* %5
%32 = load double, double* %5
%33 = load double, double* %5
%34 = load double, double* %6
%35 = load double, double* %6
%36 = fmul double %34, %35
%37 = call double @llvm.fmuladd.f64(double %32, double %33, double %36)
%38 = fcmp ogt double %37, 4.000000e+00
br i1 %38, label %39, label %40
39: ; preds = %14
br label %44
40: ; preds = %14
br label %41
41: ; preds = %40
%42 = load i32, i32* %8
%43 = add i32 %42, 1
store i32 %43, i32* %8
br label %11
44: ; preds = %39, %11
%45 = load i32, i32* %8
ret i32 %45
}

LLVM IR (-O2):
define i32 @mandel_double(double %0, double %1) {
br label %3
3: ; preds = %2, %17
%4 = phi i32 [ 0, %2 ], [ %18, %17 ]
%5 = phi double [ %0, %2 ], [ %10, %17 ]
%6 = phi double [ %1, %2 ], [ %13, %17 ]
%7 = fneg double %6
%8 = fmul double %6, %7
%9 = tail call double @llvm.fmuladd.f64(double %5, double %5, double %8)
%10 = fadd double %9, %0
%11 = fmul double %5, 2.000000e+00
%12 = fmul double %6, %11
%13 = fadd double %12, %1
%14 = fmul double %13, %13
%15 = tail call double @llvm.fmuladd.f64(double %10, double %10, double %14)
%16 = fcmp ogt double %15, 4.000000e+00
br i1 %16, label %20, label %17
17: ; preds = %3
%18 = add nuw nsw i32 %4, 1
%19 = icmp eq i32 %18, 16384
br i1 %19, label %20, label %3
20: ; preds = %3, %17
%21 = phi i32 [ %4, %3 ], [ 16384, %17 ]
ret i32 %21
}