#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ dblcmp(double, double);
double nearbyint(double);
double round_to_3digits(double val) {
  int i, exp = 0, polarity = 1;
  if (!dblcmp(val, 0.0)) {
    return 0;
  }
  if (val < 0) {
    polarity = -1;
    val *= -1;
  }
  while (val < 99.999) {
    val *= 10;
    exp--;
  }
  while (val > 999.999) {
    val /= 10;
    exp++;
  }
  val = nearbyint(val);
  for (i = 0; i < exp; i++) {
    val *= 10;
  }
  for (i = 0; i > exp; i--) {
    val /= 10;
  }
  return val * polarity;
}


x86 assembly (clang, -O0):
<round_to_3digits>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movsd  %xmm0,-0x10(%rbp)
movl   $0x0,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
call   1030 <dblcmp@plt>
cmp    $0x0,%eax
jne    115e <round_to_3digits+0x3e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1289 <round_to_3digits+0x169>
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    1185 <round_to_3digits+0x65>
movl   $0xffffffff,-0x1c(%rbp)
movsd  0xe85(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
jmp    118a <round_to_3digits+0x6a>
movsd  0xe76(%rip),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    11bd <round_to_3digits+0x9d>
movsd  0xe73(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x18(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x18(%rbp)
jmp    118a <round_to_3digits+0x6a>
jmp    11c2 <round_to_3digits+0xa2>
movsd  -0x10(%rbp),%xmm0
movsd  0xe41(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11fd <round_to_3digits+0xdd>
movsd  -0x10(%rbp),%xmm0
movsd  0xe32(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11c2 <round_to_3digits+0xa2>
movsd  -0x10(%rbp),%xmm0
call   1040 <nearbyint@plt>
movsd  %xmm0,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    123f <round_to_3digits+0x11f>
movsd  0xdf1(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1213 <round_to_3digits+0xf3>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
jle    1276 <round_to_3digits+0x156>
movsd  -0x10(%rbp),%xmm0
movsd  0xdb9(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    1246 <round_to_3digits+0x126>
movsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<round_to_3digits>:
push   %rbx
sub    $0x10,%rsp
xorpd  %xmm1,%xmm1
movapd %xmm0,(%rsp)
call   1030 <dblcmp@plt>
xorpd  %xmm0,%xmm0
test   %eax,%eax
je     12de <round_to_3digits+0x1be>
movapd 0xeb9(%rip),%xmm1
movapd (%rsp),%xmm3
xorpd  %xmm3,%xmm1
xorpd  %xmm2,%xmm2
movapd %xmm3,%xmm0
cmpltsd %xmm2,%xmm0
andpd  %xmm0,%xmm1
andnpd %xmm3,%xmm0
orpd   %xmm1,%xmm0
xor    %ebx,%ebx
movsd  0xe9d(%rip),%xmm1
ucomisd %xmm0,%xmm1
jbe    119d <round_to_3digits+0x7d>
xor    %ebx,%ebx
movsd  0xe95(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mulsd  %xmm2,%xmm0
add    $0xffffffff,%ebx
ucomisd %xmm0,%xmm1
ja     1190 <round_to_3digits+0x70>
ucomisd 0xe7b(%rip),%xmm0
jbe    11cd <round_to_3digits+0xad>
movsd  0xe69(%rip),%xmm1
movsd  0xe69(%rip),%xmm2
nopw   0x0(%rax,%rax,1)
divsd  %xmm1,%xmm0
add    $0x1,%ebx
ucomisd %xmm2,%xmm0
ja     11c0 <round_to_3digits+0xa0>
call   1040 <nearbyint@plt>
test   %ebx,%ebx
jle    11ea <round_to_3digits+0xca>
lea    -0x1(%rbx),%ecx
mov    %ebx,%eax
and    $0x7,%eax
cmp    $0x7,%ecx
jae    11f1 <round_to_3digits+0xd1>
movapd (%rsp),%xmm3
jmp    1235 <round_to_3digits+0x115>
movapd (%rsp),%xmm3
jmp    1259 <round_to_3digits+0x139>
mov    %ebx,%ecx
and    $0xfffffff8,%ecx
movsd  0xe1a(%rip),%xmm1
movapd (%rsp),%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   (%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,%ecx
jne    1210 <round_to_3digits+0xf0>
test   %eax,%eax
je     1259 <round_to_3digits+0x139>
movsd  0xdd7(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mulsd  %xmm1,%xmm0
add    $0xffffffff,%eax
jne    1250 <round_to_3digits+0x130>
test   %ebx,%ebx
jns    12b9 <round_to_3digits+0x199>
mov    %ebx,%ecx
neg    %ecx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0xfffffffc,%ebx
ja     1295 <round_to_3digits+0x175>
and    $0xfffffffc,%ecx
movsd  0xda2(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
add    $0xfffffffc,%ecx
jne    1280 <round_to_3digits+0x160>
test   %eax,%eax
je     12b9 <round_to_3digits+0x199>
movsd  0xd77(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
divsd  %xmm1,%xmm0
add    $0xffffffff,%eax
jne    12b0 <round_to_3digits+0x190>
movapd 0xd3f(%rip),%xmm1
xorpd  %xmm0,%xmm1
xorpd  %xmm2,%xmm2
cmpltsd %xmm2,%xmm3
andpd  %xmm3,%xmm1
andnpd %xmm0,%xmm3
orpd   %xmm1,%xmm3
movapd %xmm3,%xmm0
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define double @round_to_3digits(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
store double %0, double* %3
store i32 0, i32* %5
store i32 1, i32* %6
%7 = load double, double* %3
%8 = call i32 @dblcmp(double %7, double 0.000000e+00)
%9 = icmp ne i32 %8, 0
br i1 %9, label %11, label %10
10: ; preds = %1
store double 0.000000e+00, double* %2
br label %64
11: ; preds = %1
%12 = load double, double* %3
%13 = fcmp olt double %12, 0.000000e+00
br i1 %13, label %14, label %17
14: ; preds = %11
store i32 -1, i32* %6
%15 = load double, double* %3
%16 = fmul double %15, -1.000000e+00
store double %16, double* %3
br label %17
17: ; preds = %14, %11
br label %18
18: ; preds = %21, %17
%19 = load double, double* %3
%20 = fcmp olt double %19, 0x4058FFEF9DB22D0E
br i1 %20, label %21, label %26
21: ; preds = %18
%22 = load double, double* %3
%23 = fmul double %22, 1.000000e+01
store double %23, double* %3
%24 = load i32, i32* %5
%25 = add nsw i32 %24, -1
store i32 %25, i32* %5
br label %18
26: ; preds = %18
br label %27
27: ; preds = %30, %26
%28 = load double, double* %3
%29 = fcmp ogt double %28, 9.999990e+02
br i1 %29, label %30, label %35
30: ; preds = %27
%31 = load double, double* %3
%32 = fdiv double %31, 1.000000e+01
store double %32, double* %3
%33 = load i32, i32* %5
%34 = add nsw i32 %33, 1
store i32 %34, i32* %5
br label %27
35: ; preds = %27
%36 = load double, double* %3
%37 = call double @llvm.nearbyint.f64(double %36)
store double %37, double* %3
store i32 0, i32* %4
br label %38
38: ; preds = %45, %35
%39 = load i32, i32* %4
%40 = load i32, i32* %5
%41 = icmp slt i32 %39, %40
br i1 %41, label %42, label %48
42: ; preds = %38
%43 = load double, double* %3
%44 = fmul double %43, 1.000000e+01
store double %44, double* %3
br label %45
45: ; preds = %42
%46 = load i32, i32* %4
%47 = add nsw i32 %46, 1
store i32 %47, i32* %4
br label %38
48: ; preds = %38
store i32 0, i32* %4
br label %49
49: ; preds = %56, %48
%50 = load i32, i32* %4
%51 = load i32, i32* %5
%52 = icmp sgt i32 %50, %51
br i1 %52, label %53, label %59
53: ; preds = %49
%54 = load double, double* %3
%55 = fdiv double %54, 1.000000e+01
store double %55, double* %3
br label %56
56: ; preds = %53
%57 = load i32, i32* %4
%58 = add nsw i32 %57, -1
store i32 %58, i32* %4
br label %49
59: ; preds = %49
%60 = load double, double* %3
%61 = load i32, i32* %6
%62 = sitofp i32 %61 to double
%63 = fmul double %60, %62
store double %63, double* %2
br label %64
64: ; preds = %59, %10
%65 = load double, double* %2
ret double %65
}

LLVM IR (-O2):
define double @round_to_3digits(double %0) {
%2 = tail call i32 @dblcmp(double %0, double 0.000000e+00)
%3 = icmp eq i32 %2, 0
br i1 %3, label %91, label %4
4: ; preds = %1
%5 = fcmp olt double %0, 0.000000e+00
%6 = fneg double %0
%7 = select i1 %5, double %6, double %0
%8 = fcmp olt double %7, 0x4058FFEF9DB22D0E
br i1 %8, label %13, label %9
9: ; preds = %13, %4
%10 = phi double [ %7, %4 ], [ %16, %13 ]
%11 = phi i32 [ 0, %4 ], [ %17, %13 ]
%12 = fcmp ogt double %10, 9.999990e+02
br i1 %12, label %19, label %25
13: ; preds = %4, %13
%14 = phi i32 [ %17, %13 ], [ 0, %4 ]
%15 = phi double [ %16, %13 ], [ %7, %4 ]
%16 = fmul double %15, 1.000000e+01
%17 = add nsw i32 %14, -1
%18 = fcmp olt double %16, 0x4058FFEF9DB22D0E
br i1 %18, label %13, label %9
19: ; preds = %9, %19
%20 = phi i32 [ %23, %19 ], [ %11, %9 ]
%21 = phi double [ %22, %19 ], [ %10, %9 ]
%22 = fdiv double %21, 1.000000e+01
%23 = add nsw i32 %20, 1
%24 = fcmp ogt double %22, 9.999990e+02
br i1 %24, label %19, label %25
25: ; preds = %19, %9
%26 = phi double [ %10, %9 ], [ %22, %19 ]
%27 = phi i32 [ %11, %9 ], [ %23, %19 ]
%28 = tail call double @llvm.nearbyint.f64(double %26)
%29 = icmp sgt i32 %27, 0
br i1 %29, label %30, label %46
30: ; preds = %25
%31 = add i32 %27, -1
%32 = and i32 %27, 7
%33 = icmp ult i32 %31, 7
br i1 %33, label %36, label %34
34: ; preds = %30
%35 = and i32 %27, -8
br label %55
36: ; preds = %55, %30
%37 = phi double [ undef, %30 ], [ %65, %55 ]
%38 = phi double [ %28, %30 ], [ %65, %55 ]
%39 = icmp eq i32 %32, 0
br i1 %39, label %46, label %40
40: ; preds = %36, %40
%41 = phi double [ %43, %40 ], [ %38, %36 ]
%42 = phi i32 [ %44, %40 ], [ 0, %36 ]
%43 = fmul double %41, 1.000000e+01
%44 = add i32 %42, 1
%45 = icmp eq i32 %44, %32
br i1 %45, label %46, label %40
46: ; preds = %36, %40, %25
%47 = phi double [ %28, %25 ], [ %37, %36 ], [ %43, %40 ]
%48 = icmp slt i32 %27, 0
br i1 %48, label %49, label %87
49: ; preds = %46
%50 = sub i32 0, %27
%51 = and i32 %50, 3
%52 = icmp ugt i32 %27, -4
br i1 %52, label %77, label %53
53: ; preds = %49
%54 = and i32 %50, -4
br label %68
55: ; preds = %55, %34
%56 = phi double [ %28, %34 ], [ %65, %55 ]
%57 = phi i32 [ 0, %34 ], [ %66, %55 ]
%58 = fmul double %56, 1.000000e+01
%59 = fmul double %58, 1.000000e+01
%60 = fmul double %59, 1.000000e+01
%61 = fmul double %60, 1.000000e+01
%62 = fmul double %61, 1.000000e+01
%63 = fmul double %62, 1.000000e+01
%64 = fmul double %63, 1.000000e+01
%65 = fmul double %64, 1.000000e+01
%66 = add i32 %57, 8
%67 = icmp eq i32 %66, %35
br i1 %67, label %36, label %55
68: ; preds = %68, %53
%69 = phi double [ %47, %53 ], [ %74, %68 ]
%70 = phi i32 [ 0, %53 ], [ %75, %68 ]
%71 = fdiv double %69, 1.000000e+01
%72 = fdiv double %71, 1.000000e+01
%73 = fdiv double %72, 1.000000e+01
%74 = fdiv double %73, 1.000000e+01
%75 = add i32 %70, 4
%76 = icmp eq i32 %75, %54
br i1 %76, label %77, label %68
77: ; preds = %68, %49
%78 = phi double [ undef, %49 ], [ %74, %68 ]
%79 = phi double [ %47, %49 ], [ %74, %68 ]
%80 = icmp eq i32 %51, 0
br i1 %80, label %87, label %81
81: ; preds = %77, %81
%82 = phi double [ %84, %81 ], [ %79, %77 ]
%83 = phi i32 [ %85, %81 ], [ 0, %77 ]
%84 = fdiv double %82, 1.000000e+01
%85 = add i32 %83, 1
%86 = icmp eq i32 %85, %51
br i1 %86, label %87, label %81
87: ; preds = %77, %81, %46
%88 = phi double [ %47, %46 ], [ %78, %77 ], [ %84, %81 ]
%89 = fneg double %88
%90 = select i1 %5, double %89, double %88
br label %91
91: ; preds = %1, %87
%92 = phi double [ %90, %87 ], [ 0.000000e+00, %1 ]
ret double %92
}