double get_limit(double value) {
  int i;
  double limit;
  if (value >= 1.0) {
    limit = 1.0;
    for (i = 0; i <= 20; i++) {
      if (limit >= value)
        return (limit);
      if (2.0 * limit >= value)
        return (2.0 * limit);
      if (5.0 * limit >= value)
        return (5.0 * limit);
      limit *= 10.0;
    }
    return (value);
  }
  limit = 1.0;
  for (i = 0; i <= 20; i++) {
    if (0.5 < value * limit)
      return (1.0 / limit);
    limit *= 10.0;
    if (2.0 < value * limit)
      return (5.0 / limit);
    if (1.0 < value * limit)
      return (2.0 / limit);
  }
  return (value);
}


x86 assembly (clang, -O0):
<get_limit>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  0xeea(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     11ea <get_limit+0xea>
movsd  0xed8(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x14,-0x14(%rbp)
jg     11db <get_limit+0xdb>
movsd  -0x20(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jb     115d <get_limit+0x5d>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
movsd  0xeab(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jb     118c <get_limit+0x8c>
movsd  0xe93(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
movsd  0xe8c(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jb     11bb <get_limit+0xbb>
movsd  0xe74(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
movsd  0xe55(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1134 <get_limit+0x34>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
movsd  0xe0e(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x14,-0x14(%rbp)
jg     12c6 <get_limit+0x1c6>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xdee(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    123b <get_limit+0x13b>
movsd  0xdd4(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
movsd  0xdd5(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xdb1(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1280 <get_limit+0x180>
movsd  0xdaf(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xd6e(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12b3 <get_limit+0x1b3>
movsd  0xd6c(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12d0 <get_limit+0x1d0>
jmp    12b8 <get_limit+0x1b8>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11fe <get_limit+0xfe>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_limit>:
ucomisd 0xef8(%rip),%xmm0
jae    11a9 <get_limit+0xa9>
movsd  0xeea(%rip),%xmm3
mov    $0x14,%eax
movsd  0xef5(%rip),%xmm4
movsd  0xee5(%rip),%xmm5
movsd  0xeed(%rip),%xmm1
movapd %xmm3,%xmm2
nopw   0x0(%rax,%rax,1)
movapd %xmm2,%xmm6
mulsd  %xmm0,%xmm6
ucomisd %xmm4,%xmm6
ja     1228 <get_limit+0x128>
mulsd  %xmm5,%xmm2
movapd %xmm2,%xmm6
mulsd  %xmm0,%xmm6
ucomisd %xmm1,%xmm6
ja     1232 <get_limit+0x132>
ucomisd %xmm3,%xmm6
ja     123a <get_limit+0x13a>
test   %eax,%eax
je     1242 <get_limit+0x142>
ucomisd %xmm4,%xmm6
ja     1228 <get_limit+0x128>
mulsd  %xmm5,%xmm2
movapd %xmm2,%xmm6
mulsd  %xmm0,%xmm6
ucomisd %xmm1,%xmm6
ja     1232 <get_limit+0x132>
ucomisd %xmm3,%xmm6
ja     123a <get_limit+0x13a>
add    $0xfffffffe,%eax
jmp    1140 <get_limit+0x40>
movsd  0xe4f(%rip),%xmm2
mov    $0x14,%eax
movsd  0xe4a(%rip),%xmm3
movsd  0xe4a(%rip),%xmm4
ucomisd %xmm0,%xmm2
jae    1223 <get_limit+0x123>
nopl   0x0(%rax)
movapd %xmm2,%xmm1
addsd  %xmm2,%xmm1
ucomisd %xmm0,%xmm1
jae    123e <get_limit+0x13e>
movapd %xmm2,%xmm1
mulsd  %xmm3,%xmm1
ucomisd %xmm0,%xmm1
jae    123e <get_limit+0x13e>
mulsd  %xmm4,%xmm2
test   %eax,%eax
je     1242 <get_limit+0x142>
ucomisd %xmm0,%xmm2
jae    1223 <get_limit+0x123>
movapd %xmm2,%xmm1
addsd  %xmm2,%xmm1
ucomisd %xmm0,%xmm1
jae    123e <get_limit+0x13e>
movapd %xmm2,%xmm1
mulsd  %xmm3,%xmm1
ucomisd %xmm0,%xmm1
jae    123e <get_limit+0x13e>
mulsd  %xmm4,%xmm2
add    $0xfffffffe,%eax
ucomisd %xmm0,%xmm2
jb     11d0 <get_limit+0xd0>
movapd %xmm2,%xmm0
ret
movsd  0xdd0(%rip),%xmm1
jmp    123a <get_limit+0x13a>
movsd  0xdce(%rip),%xmm1
divsd  %xmm2,%xmm1
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @get_limit(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca i32
%5 = alloca double
store double %0, double* %3
%6 = load double, double* %3
%7 = fcmp oge double %6, 1.000000e+00
br i1 %7, label %8, label %42
8: ; preds = %1
store double 1.000000e+00, double* %5
store i32 0, i32* %4
br label %9
9: ; preds = %37, %8
%10 = load i32, i32* %4
%11 = icmp sle i32 %10, 20
br i1 %11, label %12, label %40
12: ; preds = %9
%13 = load double, double* %5
%14 = load double, double* %3
%15 = fcmp oge double %13, %14
br i1 %15, label %16, label %18
16: ; preds = %12
%17 = load double, double* %5
store double %17, double* %2
br label %78
18: ; preds = %12
%19 = load double, double* %5
%20 = fmul double 2.000000e+00, %19
%21 = load double, double* %3
%22 = fcmp oge double %20, %21
br i1 %22, label %23, label %26
23: ; preds = %18
%24 = load double, double* %5
%25 = fmul double 2.000000e+00, %24
store double %25, double* %2
br label %78
26: ; preds = %18
%27 = load double, double* %5
%28 = fmul double 5.000000e+00, %27
%29 = load double, double* %3
%30 = fcmp oge double %28, %29
br i1 %30, label %31, label %34
31: ; preds = %26
%32 = load double, double* %5
%33 = fmul double 5.000000e+00, %32
store double %33, double* %2
br label %78
34: ; preds = %26
%35 = load double, double* %5
%36 = fmul double %35, 1.000000e+01
store double %36, double* %5
br label %37
37: ; preds = %34
%38 = load i32, i32* %4
%39 = add nsw i32 %38, 1
store i32 %39, i32* %4
br label %9
40: ; preds = %9
%41 = load double, double* %3
store double %41, double* %2
br label %78
42: ; preds = %1
store double 1.000000e+00, double* %5
store i32 0, i32* %4
br label %43
43: ; preds = %73, %42
%44 = load i32, i32* %4
%45 = icmp sle i32 %44, 20
br i1 %45, label %46, label %76
46: ; preds = %43
%47 = load double, double* %3
%48 = load double, double* %5
%49 = fmul double %47, %48
%50 = fcmp olt double 5.000000e-01, %49
br i1 %50, label %51, label %54
51: ; preds = %46
%52 = load double, double* %5
%53 = fdiv double 1.000000e+00, %52
store double %53, double* %2
br label %78
54: ; preds = %46
%55 = load double, double* %5
%56 = fmul double %55, 1.000000e+01
store double %56, double* %5
%57 = load double, double* %3
%58 = load double, double* %5
%59 = fmul double %57, %58
%60 = fcmp olt double 2.000000e+00, %59
br i1 %60, label %61, label %64
61: ; preds = %54
%62 = load double, double* %5
%63 = fdiv double 5.000000e+00, %62
store double %63, double* %2
br label %78
64: ; preds = %54
%65 = load double, double* %3
%66 = load double, double* %5
%67 = fmul double %65, %66
%68 = fcmp olt double 1.000000e+00, %67
br i1 %68, label %69, label %72
69: ; preds = %64
%70 = load double, double* %5
%71 = fdiv double 2.000000e+00, %70
store double %71, double* %2
br label %78
72: ; preds = %64
br label %73
73: ; preds = %72
%74 = load i32, i32* %4
%75 = add nsw i32 %74, 1
store i32 %75, i32* %4
br label %43
76: ; preds = %43
%77 = load double, double* %3
store double %77, double* %2
br label %78
78: ; preds = %76, %69, %61, %51, %40, %31, %23, %16
%79 = load double, double* %2
ret double %79
}

LLVM IR (-O2):
define double @get_limit(double %0) {
%2 = fcmp ult double %0, 1.000000e+00
br i1 %2, label %40, label %3
3: ; preds = %1, %24
%4 = phi double [ %25, %24 ], [ 1.000000e+00, %1 ]
%5 = phi i32 [ %26, %24 ], [ 0, %1 ]
%6 = fcmp ult double %4, %0
br i1 %6, label %7, label %60
7: ; preds = %3
%8 = fmul double %4, 2.000000e+00
%9 = fcmp ult double %8, %0
br i1 %9, label %10, label %60
10: ; preds = %7
%11 = fmul double %4, 5.000000e+00
%12 = fcmp ult double %11, %0
br i1 %12, label %13, label %60
13: ; preds = %10
%14 = fmul double %4, 1.000000e+01
%15 = icmp eq i32 %5, 20
br i1 %15, label %60, label %16
16: ; preds = %13
%17 = fcmp ult double %14, %0
br i1 %17, label %18, label %60
18: ; preds = %16
%19 = fmul double %14, 2.000000e+00
%20 = fcmp ult double %19, %0
br i1 %20, label %21, label %60
21: ; preds = %18
%22 = fmul double %14, 5.000000e+00
%23 = fcmp ult double %22, %0
br i1 %23, label %24, label %60
24: ; preds = %21
%25 = fmul double %14, 1.000000e+01
%26 = add nuw nsw i32 %5, 2
br label %3
27: ; preds = %55
%28 = icmp eq i32 %42, 20
br i1 %28, label %60, label %29
29: ; preds = %27
%30 = fmul double %49, %0
%31 = fcmp ogt double %30, 5.000000e-01
br i1 %31, label %45, label %32
32: ; preds = %29
%33 = fmul double %49, 1.000000e+01
%34 = fmul double %33, %0
%35 = fcmp ogt double %34, 2.000000e+00
br i1 %35, label %52, label %36
36: ; preds = %32
%37 = fcmp ogt double %34, 1.000000e+00
br i1 %37, label %57, label %38
38: ; preds = %36
%39 = add nuw nsw i32 %42, 2
br label %40
40: ; preds = %1, %38
%41 = phi double [ %33, %38 ], [ 1.000000e+00, %1 ]
%42 = phi i32 [ %39, %38 ], [ 0, %1 ]
%43 = fmul double %41, %0
%44 = fcmp ogt double %43, 5.000000e-01
br i1 %44, label %45, label %48
45: ; preds = %29, %40
%46 = phi double [ %41, %40 ], [ %49, %29 ]
%47 = fdiv double 1.000000e+00, %46
br label %60
48: ; preds = %40
%49 = fmul double %41, 1.000000e+01
%50 = fmul double %49, %0
%51 = fcmp ogt double %50, 2.000000e+00
br i1 %51, label %52, label %55
52: ; preds = %32, %48
%53 = phi double [ %49, %48 ], [ %33, %32 ]
%54 = fdiv double 5.000000e+00, %53
br label %60
55: ; preds = %48
%56 = fcmp ogt double %50, 1.000000e+00
br i1 %56, label %57, label %27
57: ; preds = %36, %55
%58 = phi double [ %49, %55 ], [ %33, %36 ]
%59 = fdiv double 2.000000e+00, %58
br label %60
60: ; preds = %13, %10, %7, %3, %16, %18, %21, %27, %57, %52, %45
%61 = phi double [ %47, %45 ], [ %54, %52 ], [ %59, %57 ], [ %0, %27 ], [ %0, %13 ], [ %11, %10 ], [ %8, %7 ], [ %4, %3 ], [ %14, %16 ], [ %19, %18 ], [ %22, %21 ]
ret double %61
}