double cbrt(double x) {
  if (x == 0) {
    return 0;
  }
  double b = 1;
  double last_b_1 = 0;
  double last_b_2 = 0;
  while (last_b_1 != b && last_b_2 != b) {
    last_b_1 = b;
    b = (b + x / b / b) / 2;
    last_b_2 = b;
    b = (b + x / b / b) / 2;
  }
  return b;
}


x86 assembly (clang, -O0):
<cbrt>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    112e <cbrt+0x2e>
jp     112e <cbrt+0x2e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1205 <cbrt+0x105>
movsd  0xeca(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
xor    %eax,%eax
ucomisd -0x18(%rbp),%xmm0
mov    %al,-0x29(%rbp)
jne    116b <cbrt+0x6b>
jp     116b <cbrt+0x6b>
jmp    1180 <cbrt+0x80>
movsd  -0x28(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
setne  %al
setp   %cl
or     %cl,%al
mov    %al,-0x29(%rbp)
mov    -0x29(%rbp),%al
test   $0x1,%al
jne    1190 <cbrt+0x90>
jmp    11fb <cbrt+0xfb>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  0xe4e(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  0xe1b(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    114b <cbrt+0x4b>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<cbrt>:
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    110c <cbrt+0xc>
jnp    115c <cbrt+0x5c>
movsd  0xeec(%rip),%xmm1
movsd  0xeec(%rip),%xmm2
nopl   0x0(%rax)
movapd %xmm1,%xmm4
movapd %xmm0,%xmm3
divsd  %xmm1,%xmm3
divsd  %xmm1,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm0,%xmm1
divsd  %xmm3,%xmm1
divsd  %xmm3,%xmm1
addsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm4
jne    1154 <cbrt+0x54>
jnp    115c <cbrt+0x5c>
ucomisd %xmm1,%xmm3
jne    1120 <cbrt+0x20>
jp     1120 <cbrt+0x20>
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @cbrt(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
store double %0, double* %3
%7 = load double, double* %3
%8 = fcmp oeq double %7, 0.000000e+00
br i1 %8, label %9, label %10
9: ; preds = %1
store double 0.000000e+00, double* %2
br label %42
10: ; preds = %1
store double 1.000000e+00, double* %4
store double 0.000000e+00, double* %5
store double 0.000000e+00, double* %6
br label %11
11: ; preds = %21, %10
%12 = load double, double* %5
%13 = load double, double* %4
%14 = fcmp une double %12, %13
br i1 %14, label %15, label %19
15: ; preds = %11
%16 = load double, double* %6
%17 = load double, double* %4
%18 = fcmp une double %16, %17
br label %19
19: ; preds = %15, %11
%20 = phi i1 [ false, %11 ], [ %18, %15 ]
br i1 %20, label %21, label %40
21: ; preds = %19
%22 = load double, double* %4
store double %22, double* %5
%23 = load double, double* %4
%24 = load double, double* %3
%25 = load double, double* %4
%26 = fdiv double %24, %25
%27 = load double, double* %4
%28 = fdiv double %26, %27
%29 = fadd double %23, %28
%30 = fdiv double %29, 2.000000e+00
store double %30, double* %4
%31 = load double, double* %4
store double %31, double* %6
%32 = load double, double* %4
%33 = load double, double* %3
%34 = load double, double* %4
%35 = fdiv double %33, %34
%36 = load double, double* %4
%37 = fdiv double %35, %36
%38 = fadd double %32, %37
%39 = fdiv double %38, 2.000000e+00
store double %39, double* %4
br label %11
40: ; preds = %19
%41 = load double, double* %4
store double %41, double* %2
br label %42
42: ; preds = %40, %9
%43 = load double, double* %2
ret double %43
}

LLVM IR (-O2):
define double @cbrt(double %0) {
%2 = fcmp oeq double %0, 0.000000e+00
br i1 %2, label %16, label %3
3: ; preds = %1, %3
%4 = phi double [ %12, %3 ], [ 1.000000e+00, %1 ]
%5 = fdiv double %0, %4
%6 = fdiv double %5, %4
%7 = fadd double %4, %6
%8 = fmul double %7, 5.000000e-01
%9 = fdiv double %0, %8
%10 = fdiv double %9, %8
%11 = fadd double %8, %10
%12 = fmul double %11, 5.000000e-01
%13 = fcmp une double %4, %12
%14 = fcmp une double %8, %12
%15 = select i1 %13, i1 %14, i1 false
br i1 %15, label %3, label %16
16: ; preds = %3, %1
%17 = phi double [ 0.000000e+00, %1 ], [ %12, %3 ]
ret double %17
}