double RoundUp(val)
double val;
{
  int exponent, index;
  exponent = (int)floor((val == 0.0 ? 0.0 : log10(val)));
  if (exponent < 0) {
    for (index = exponent; index < 0; index++) {
      val *= 10.0;
    }
  } else {
    for (index = 0; index < exponent; index++) {
      val /= 10.0;
    }
  }
  if (val > 5.0)
    val = 10.0;
  else if (val > 2.0)
    val = 5.0;
  else if (val > 1.0)
    val = 2.0;
  else
    val = 1.0;
  if (exponent < 0)
    for (index = exponent; index < 0; index++)
      val /= 10.0;
  else
    for (index = 0; index < exponent; index++)
      val *= 10.0;
  return (val);
}


x86 assembly (clang, -O0):
<RoundUp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1152 <RoundUp+0x32>
jp     1152 <RoundUp+0x32>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    1161 <RoundUp+0x41>
movsd  -0x8(%rbp),%xmm0
call   1030 <log10@plt>
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1040 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0xc(%rbp)
cmpl   $0x0,-0xc(%rbp)
jge    11b1 <RoundUp+0x91>
mov    -0xc(%rbp),%eax
mov    %eax,-0x10(%rbp)
cmpl   $0x0,-0x10(%rbp)
jge    11ac <RoundUp+0x8c>
movsd  0xe6c(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1182 <RoundUp+0x62>
jmp    11ed <RoundUp+0xcd>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11e8 <RoundUp+0xc8>
movsd  -0x8(%rbp),%xmm0
movsd  0xe2f(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    11b8 <RoundUp+0x98>
jmp    11ed <RoundUp+0xcd>
movsd  -0x8(%rbp),%xmm0
movsd  0xe0e(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1216 <RoundUp+0xf6>
movsd  0xdf4(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    127f <RoundUp+0x15f>
movsd  -0x8(%rbp),%xmm0
movsd  0xded(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    123f <RoundUp+0x11f>
movsd  0xdd3(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    127a <RoundUp+0x15a>
movsd  -0x8(%rbp),%xmm0
movsd  0xdcc(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1268 <RoundUp+0x148>
movsd  0xdb2(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1275 <RoundUp+0x155>
movsd  0xda8(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    127a <RoundUp+0x15a>
jmp    127f <RoundUp+0x15f>
cmpl   $0x0,-0xc(%rbp)
jge    12c2 <RoundUp+0x1a2>
mov    -0xc(%rbp),%eax
mov    %eax,-0x10(%rbp)
cmpl   $0x0,-0x10(%rbp)
jge    12bd <RoundUp+0x19d>
movsd  -0x8(%rbp),%xmm0
movsd  0xd5a(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    128f <RoundUp+0x16f>
jmp    12fa <RoundUp+0x1da>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12f5 <RoundUp+0x1d5>
movsd  0xd23(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    12c9 <RoundUp+0x1a9>
jmp    12fa <RoundUp+0x1da>
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<RoundUp>:
push   %rax
movapd %xmm0,%xmm1
xorpd  %xmm0,%xmm0
movsd  %xmm1,(%rsp)
ucomisd %xmm0,%xmm1
jne    1136 <RoundUp+0x16>
jnp    1140 <RoundUp+0x20>
movsd  (%rsp),%xmm0
call   1030 <log10@plt>
call   1040 <floor@plt>
cvttsd2si %xmm0,%eax
test   %eax,%eax
js     11be <RoundUp+0x9e>
je     11f6 <RoundUp+0xd6>
lea    -0x1(%rax),%edx
mov    %eax,%ecx
and    $0x3,%ecx
cmp    $0x3,%edx
movsd  (%rsp),%xmm1
jb     1195 <RoundUp+0x75>
mov    %eax,%edx
and    $0xfffffffc,%edx
movsd  0xe8e(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
divsd  %xmm0,%xmm1
divsd  %xmm0,%xmm1
divsd  %xmm0,%xmm1
divsd  %xmm0,%xmm1
add    $0xfffffffc,%edx
jne    1180 <RoundUp+0x60>
test   %ecx,%ecx
je     1245 <RoundUp+0x125>
movsd  0xe5b(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
divsd  %xmm0,%xmm1
add    $0xffffffff,%ecx
jne    11b0 <RoundUp+0x90>
jmp    1245 <RoundUp+0x125>
mov    %eax,%ecx
neg    %ecx
and    $0x7,%ecx
je     11fd <RoundUp+0xdd>
neg    %ecx
xor    %edx,%edx
movsd  0xe2d(%rip),%xmm0
movsd  (%rsp),%xmm1
nopl   0x0(%rax,%rax,1)
mulsd  %xmm0,%xmm1
add    $0xffffffff,%edx
cmp    %edx,%ecx
jne    11e0 <RoundUp+0xc0>
mov    %eax,%ecx
sub    %edx,%ecx
cmp    $0xfffffff8,%eax
jbe    1209 <RoundUp+0xe9>
jmp    1245 <RoundUp+0x125>
movsd  (%rsp),%xmm1
jmp    1245 <RoundUp+0x125>
mov    %eax,%ecx
movsd  (%rsp),%xmm1
cmp    $0xfffffff8,%eax
ja     1245 <RoundUp+0x125>
movsd  0xdef(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
add    $0x8,%ecx
jne    1220 <RoundUp+0x100>
ucomisd 0xdbb(%rip),%xmm1
jbe    1260 <RoundUp+0x140>
movsd  0xda9(%rip),%xmm0
test   %eax,%eax
jns    127e <RoundUp+0x15e>
jmp    130f <RoundUp+0x1ef>
ucomisd 0xda8(%rip),%xmm1
jbe    12ee <RoundUp+0x1ce>
movsd  0xd92(%rip),%xmm0
test   %eax,%eax
js     130f <RoundUp+0x1ef>
je     1375 <RoundUp+0x255>
lea    -0x1(%rax),%edx
mov    %eax,%ecx
and    $0x7,%ecx
cmp    $0x7,%edx
jb     12c5 <RoundUp+0x1a5>
and    $0xfffffff8,%eax
movsd  0xd64(%rip),%xmm1
nopl   0x0(%rax)
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
add    $0xfffffff8,%eax
jne    12a0 <RoundUp+0x180>
test   %ecx,%ecx
je     1375 <RoundUp+0x255>
movsd  0xd2b(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nop
mulsd  %xmm1,%xmm0
add    $0xffffffff,%ecx
jne    12e0 <RoundUp+0x1c0>
jmp    1375 <RoundUp+0x255>
xor    %ecx,%ecx
ucomisd 0xd20(%rip),%xmm1
seta   %cl
lea    0xd1e(%rip),%rdx
movsd  (%rdx,%rcx,8),%xmm0
test   %eax,%eax
jns    127e <RoundUp+0x15e>
mov    %eax,%ecx
neg    %ecx
and    $0x3,%ecx
je     1346 <RoundUp+0x226>
neg    %ecx
xor    %edx,%edx
movsd  0xcdc(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
divsd  %xmm1,%xmm0
add    $0xffffffff,%edx
cmp    %edx,%ecx
jne    1330 <RoundUp+0x210>
mov    %eax,%ecx
sub    %edx,%ecx
cmp    $0xfffffffc,%eax
jbe    134d <RoundUp+0x22d>
jmp    1375 <RoundUp+0x255>
mov    %eax,%ecx
cmp    $0xfffffffc,%eax
ja     1375 <RoundUp+0x255>
movsd  0xcab(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nop
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
add    $0x4,%ecx
jne    1360 <RoundUp+0x240>
pop    %rax
ret

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

LLVM IR (-O2):
define double @RoundUp(double %0) {
%2 = fcmp oeq double %0, 0.000000e+00
br i1 %2, label %5, label %3
3: ; preds = %1
%4 = tail call double @log10(double %0)
br label %5
5: ; preds = %1, %3
%6 = phi double [ %4, %3 ], [ 0.000000e+00, %1 ]
%7 = tail call double @llvm.floor.f64(double %6)
%8 = fptosi double %7 to i32
%9 = icmp slt i32 %8, 0
br i1 %9, label %10, label %27
10: ; preds = %5
%11 = sub i32 0, %8
%12 = and i32 %11, 7
%13 = icmp eq i32 %12, 0
br i1 %13, label %22, label %14
14: ; preds = %10, %14
%15 = phi i32 [ %19, %14 ], [ %8, %10 ]
%16 = phi double [ %18, %14 ], [ %0, %10 ]
%17 = phi i32 [ %20, %14 ], [ 0, %10 ]
%18 = fmul double %16, 1.000000e+01
%19 = add i32 %15, 1
%20 = add i32 %17, 1
%21 = icmp eq i32 %20, %12
br i1 %21, label %22, label %14
22: ; preds = %14, %10
%23 = phi double [ undef, %10 ], [ %18, %14 ]
%24 = phi i32 [ %8, %10 ], [ %19, %14 ]
%25 = phi double [ %0, %10 ], [ %18, %14 ]
%26 = icmp ugt i32 %8, -8
br i1 %26, label %67, label %35
27: ; preds = %5
%28 = icmp eq i32 %8, 0
br i1 %28, label %67, label %29
29: ; preds = %27
%30 = add i32 %8, -1
%31 = and i32 %8, 3
%32 = icmp ult i32 %30, 3
br i1 %32, label %57, label %33
33: ; preds = %29
%34 = and i32 %8, -4
br label %48
35: ; preds = %22, %35
%36 = phi i32 [ %46, %35 ], [ %24, %22 ]
%37 = phi double [ %45, %35 ], [ %25, %22 ]
%38 = fmul double %37, 1.000000e+01
%39 = fmul double %38, 1.000000e+01
%40 = fmul double %39, 1.000000e+01
%41 = fmul double %40, 1.000000e+01
%42 = fmul double %41, 1.000000e+01
%43 = fmul double %42, 1.000000e+01
%44 = fmul double %43, 1.000000e+01
%45 = fmul double %44, 1.000000e+01
%46 = add i32 %36, 8
%47 = icmp eq i32 %46, 0
br i1 %47, label %67, label %35
48: ; preds = %48, %33
%49 = phi double [ %0, %33 ], [ %54, %48 ]
%50 = phi i32 [ 0, %33 ], [ %55, %48 ]
%51 = fdiv double %49, 1.000000e+01
%52 = fdiv double %51, 1.000000e+01
%53 = fdiv double %52, 1.000000e+01
%54 = fdiv double %53, 1.000000e+01
%55 = add i32 %50, 4
%56 = icmp eq i32 %55, %34
br i1 %56, label %57, label %48
57: ; preds = %48, %29
%58 = phi double [ undef, %29 ], [ %54, %48 ]
%59 = phi double [ %0, %29 ], [ %54, %48 ]
%60 = icmp eq i32 %31, 0
br i1 %60, label %67, label %61
61: ; preds = %57, %61
%62 = phi double [ %64, %61 ], [ %59, %57 ]
%63 = phi i32 [ %65, %61 ], [ 0, %57 ]
%64 = fdiv double %62, 1.000000e+01
%65 = add i32 %63, 1
%66 = icmp eq i32 %65, %31
br i1 %66, label %67, label %61
67: ; preds = %57, %61, %22, %35, %27
%68 = phi double [ %0, %27 ], [ %23, %22 ], [ %45, %35 ], [ %58, %57 ], [ %64, %61 ]
%69 = fcmp ogt double %68, 5.000000e+00
br i1 %69, label %75, label %70
70: ; preds = %67
%71 = fcmp ogt double %68, 2.000000e+00
br i1 %71, label %75, label %72
72: ; preds = %70
%73 = fcmp ogt double %68, 1.000000e+00
%74 = select i1 %73, double 2.000000e+00, double 1.000000e+00
br label %75
75: ; preds = %72, %70, %67
%76 = phi double [ 1.000000e+01, %67 ], [ 5.000000e+00, %70 ], [ %74, %72 ]
br i1 %9, label %77, label %94
77: ; preds = %75
%78 = sub i32 0, %8
%79 = and i32 %78, 3
%80 = icmp eq i32 %79, 0
br i1 %80, label %89, label %81
81: ; preds = %77, %81
%82 = phi i32 [ %86, %81 ], [ %8, %77 ]
%83 = phi double [ %85, %81 ], [ %76, %77 ]
%84 = phi i32 [ %87, %81 ], [ 0, %77 ]
%85 = fdiv double %83, 1.000000e+01
%86 = add i32 %82, 1
%87 = add i32 %84, 1
%88 = icmp eq i32 %87, %79
br i1 %88, label %89, label %81
89: ; preds = %81, %77
%90 = phi double [ undef, %77 ], [ %85, %81 ]
%91 = phi i32 [ %8, %77 ], [ %86, %81 ]
%92 = phi double [ %76, %77 ], [ %85, %81 ]
%93 = icmp ugt i32 %8, -4
br i1 %93, label %134, label %102
94: ; preds = %75
%95 = icmp eq i32 %8, 0
br i1 %95, label %134, label %96
96: ; preds = %94
%97 = add i32 %8, -1
%98 = and i32 %8, 7
%99 = icmp ult i32 %97, 7
br i1 %99, label %124, label %100
100: ; preds = %96
%101 = and i32 %8, -8
br label %111
102: ; preds = %89, %102
%103 = phi i32 [ %109, %102 ], [ %91, %89 ]
%104 = phi double [ %108, %102 ], [ %92, %89 ]
%105 = fdiv double %104, 1.000000e+01
%106 = fdiv double %105, 1.000000e+01
%107 = fdiv double %106, 1.000000e+01
%108 = fdiv double %107, 1.000000e+01
%109 = add i32 %103, 4
%110 = icmp eq i32 %109, 0
br i1 %110, label %134, label %102
111: ; preds = %111, %100
%112 = phi double [ %76, %100 ], [ %121, %111 ]
%113 = phi i32 [ 0, %100 ], [ %122, %111 ]
%114 = fmul double %112, 1.000000e+01
%115 = fmul double %114, 1.000000e+01
%116 = fmul double %115, 1.000000e+01
%117 = fmul double %116, 1.000000e+01
%118 = fmul double %117, 1.000000e+01
%119 = fmul double %118, 1.000000e+01
%120 = fmul double %119, 1.000000e+01
%121 = fmul double %120, 1.000000e+01
%122 = add i32 %113, 8
%123 = icmp eq i32 %122, %101
br i1 %123, label %124, label %111
124: ; preds = %111, %96
%125 = phi double [ undef, %96 ], [ %121, %111 ]
%126 = phi double [ %76, %96 ], [ %121, %111 ]
%127 = icmp eq i32 %98, 0
br i1 %127, label %134, label %128
128: ; preds = %124, %128
%129 = phi double [ %131, %128 ], [ %126, %124 ]
%130 = phi i32 [ %132, %128 ], [ 0, %124 ]
%131 = fmul double %129, 1.000000e+01
%132 = add i32 %130, 1
%133 = icmp eq i32 %132, %98
br i1 %133, label %134, label %128
134: ; preds = %124, %128, %89, %102, %94
%135 = phi double [ %76, %94 ], [ %90, %89 ], [ %108, %102 ], [ %125, %124 ], [ %131, %128 ]
ret double %135
}