#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 */
scalar_t__ ceil(float);
scalar_t__ floor(float);
int calc_label_digits(float start, float stop, float labels_ppm,
                      float steps_ppm) {
  float a, b, c;
  int j, k, ia, ib, init = 0, digits = 0;
  if (labels_ppm == 0 || steps_ppm == 0)
    return 1;
  k = ((int)(((steps_ppm / labels_ppm) < 0.0)
                 ? (ceil((steps_ppm / labels_ppm) - 0.5))
                 : (floor((steps_ppm / labels_ppm) + 0.5))));
  if (k < 1)
    k = 1;
  c = 1.0;
  for (j = start; j >= stop; j--) {
    if (!(j % k)) {
      b = ((float)j) / steps_ppm;
      if (!init)
        init = 1;
      else {
        do {
          ia = (int)(c * a);
          ib = (int)(c * b);
          if (ia == ib) {
            c *= 10;
            digits++;
          }
        } while (ia == ib);
      }
      a = b;
    }
  }
  return digits;
}


x86 assembly (clang, -O0):
<calc_label_digits>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
movss  %xmm2,-0x10(%rbp)
movss  %xmm3,-0x14(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
movss  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1166 <calc_label_digits+0x46>
jp     1166 <calc_label_digits+0x46>
jmp    117d <calc_label_digits+0x5d>
movss  -0x14(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1189 <calc_label_digits+0x69>
jp     1189 <calc_label_digits+0x69>
movl   $0x1,-0x4(%rbp)
jmp    12ef <calc_label_digits+0x1cf>
movss  -0x14(%rbp),%xmm0
divss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    11d0 <calc_label_digits+0xb0>
movss  -0x14(%rbp),%xmm0
divss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe46(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1030 <ceil@plt>
mov    %rax,-0x40(%rbp)
jmp    11f7 <calc_label_digits+0xd7>
movss  -0x14(%rbp),%xmm0
divss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe1a(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <floor@plt>
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    %eax,-0x28(%rbp)
cmpl   $0x1,-0x28(%rbp)
jge    120f <calc_label_digits+0xef>
movl   $0x1,-0x28(%rbp)
movss  0xdf1(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
cvttss2si -0x8(%rbp),%eax
mov    %eax,-0x24(%rbp)
cvtsi2ssl -0x24(%rbp),%xmm0
ucomiss -0xc(%rbp),%xmm0
jb     12e9 <calc_label_digits+0x1c9>
mov    -0x24(%rbp),%eax
cltd
idivl  -0x28(%rbp)
cmp    $0x0,%edx
jne    12d6 <calc_label_digits+0x1b6>
cvtsi2ssl -0x24(%rbp),%xmm0
divss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
cmpl   $0x0,-0x34(%rbp)
jne    1268 <calc_label_digits+0x148>
movl   $0x1,-0x34(%rbp)
jmp    12cc <calc_label_digits+0x1ac>
jmp    126d <calc_label_digits+0x14d>
movss  -0x20(%rbp),%xmm0
mulss  -0x18(%rbp),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x2c(%rbp)
movss  -0x20(%rbp),%xmm0
mulss  -0x1c(%rbp),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jne    12b6 <calc_label_digits+0x196>
movss  0xd69(%rip),%xmm0
mulss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12bb <calc_label_digits+0x19b>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
je     126d <calc_label_digits+0x14d>
jmp    12cc <calc_label_digits+0x1ac>
movss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
jmp    12db <calc_label_digits+0x1bb>
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1224 <calc_label_digits+0x104>
mov    -0x38(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_label_digits>:
sub    $0x28,%rsp
xorps  %xmm4,%xmm4
ucomiss %xmm4,%xmm2
mov    $0x1,%ecx
jne    1133 <calc_label_digits+0x13>
jnp    11a0 <calc_label_digits+0x80>
ucomiss %xmm4,%xmm3
jne    113a <calc_label_digits+0x1a>
jnp    11a0 <calc_label_digits+0x80>
movaps %xmm0,0x10(%rsp)
movss  %xmm1,0x8(%rsp)
movss  %xmm3,0xc(%rsp)
movaps %xmm3,%xmm0
divss  %xmm2,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm0,%xmm1
jbe    1169 <calc_label_digits+0x49>
addss  0xea2(%rip),%xmm0
call   1030 <ceil@plt>
jmp    1176 <calc_label_digits+0x56>
addss  0xe8f(%rip),%xmm0
call   1040 <floor@plt>
test   %eax,%eax
mov    $0x1,%edi
cmovg  %eax,%edi
movaps 0x10(%rsp),%xmm1
cvttps2dq %xmm1,%xmm0
cvtdq2ps %xmm0,%xmm0
xor    %ecx,%ecx
movss  0x8(%rsp),%xmm7
ucomiss %xmm7,%xmm0
movss  0xc(%rsp),%xmm8
jae    11a7 <calc_label_digits+0x87>
mov    %ecx,%eax
add    $0x28,%rsp
ret
cvttss2si %xmm1,%esi
movss  0xe55(%rip),%xmm1
xor    %ecx,%ecx
movss  0xe4f(%rip),%xmm2
xor    %r8d,%r8d
jmp    11e5 <calc_label_digits+0xc5>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movaps %xmm3,%xmm0
movaps %xmm0,%xmm3
add    $0xffffffff,%esi
xorps  %xmm0,%xmm0
cvtsi2ss %esi,%xmm0
ucomiss %xmm7,%xmm0
jb     11a0 <calc_label_digits+0x80>
mov    %esi,%eax
cltd
idiv   %edi
test   %edx,%edx
jne    11d0 <calc_label_digits+0xb0>
divss  %xmm8,%xmm0
test   %r8d,%r8d
je     122e <calc_label_digits+0x10e>
unpcklps %xmm0,%xmm3
add    $0xffffffff,%ecx
xchg   %ax,%ax
movaps %xmm1,%xmm4
movaps %xmm1,%xmm5
shufps $0x0,%xmm1,%xmm5
mulps  %xmm3,%xmm5
cvttps2dq %xmm5,%xmm5
pshufd $0x55,%xmm5,%xmm6
pcmpeqd %xmm5,%xmm6
movd   %xmm6,%eax
mulss  %xmm2,%xmm1
add    $0x1,%ecx
test   $0x1,%al
jne    1200 <calc_label_digits+0xe0>
movaps %xmm4,%xmm1
jmp    11d3 <calc_label_digits+0xb3>
mov    $0x1,%r8d
jmp    11d3 <calc_label_digits+0xb3>

LLVM IR (-O0):
define i32 @calc_label_digits(float %0, float %1, float %2, float %3) {
%5 = alloca i32
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
store float %0, float* %6
store float %1, float* %7
store float %2, float* %8
store float %3, float* %9
store i32 0, i32* %17
store i32 0, i32* %18
%19 = load float, float* %8
%20 = fcmp oeq float %19, 0.000000e+00
br i1 %20, label %24, label %21
21: ; preds = %4
%22 = load float, float* %9
%23 = fcmp oeq float %22, 0.000000e+00
br i1 %23, label %24, label %25
24: ; preds = %21, %4
store i32 1, i32* %5
br label %106
25: ; preds = %21
%26 = load float, float* %9
%27 = load float, float* %8
%28 = fdiv float %26, %27
%29 = fpext float %28 to double
%30 = fcmp olt double %29, 0.000000e+00
br i1 %30, label %31, label %39
31: ; preds = %25
%32 = load float, float* %9
%33 = load float, float* %8
%34 = fdiv float %32, %33
%35 = fpext float %34 to double
%36 = fsub double %35, 5.000000e-01
%37 = fptrunc double %36 to float
%38 = call i64 @ceil(float %37)
br label %47
39: ; preds = %25
%40 = load float, float* %9
%41 = load float, float* %8
%42 = fdiv float %40, %41
%43 = fpext float %42 to double
%44 = fadd double %43, 5.000000e-01
%45 = fptrunc double %44 to float
%46 = call i64 @floor(float %45)
br label %47
47: ; preds = %39, %31
%48 = phi i64 [ %38, %31 ], [ %46, %39 ]
%49 = trunc i64 %48 to i32
store i32 %49, i32* %14
%50 = load i32, i32* %14
%51 = icmp slt i32 %50, 1
br i1 %51, label %52, label %53
52: ; preds = %47
store i32 1, i32* %14
br label %53
53: ; preds = %52, %47
store float 1.000000e+00, float* %12
%54 = load float, float* %6
%55 = fptosi float %54 to i32
store i32 %55, i32* %13
br label %56
56: ; preds = %101, %53
%57 = load i32, i32* %13
%58 = sitofp i32 %57 to float
%59 = load float, float* %7
%60 = fcmp oge float %58, %59
br i1 %60, label %61, label %104
61: ; preds = %56
%62 = load i32, i32* %13
%63 = load i32, i32* %14
%64 = srem i32 %62, %63
%65 = icmp ne i32 %64, 0
br i1 %65, label %100, label %66
66: ; preds = %61
%67 = load i32, i32* %13
%68 = sitofp i32 %67 to float
%69 = load float, float* %9
%70 = fdiv float %68, %69
store float %70, float* %11
%71 = load i32, i32* %17
%72 = icmp ne i32 %71, 0
br i1 %72, label %74, label %73
73: ; preds = %66
store i32 1, i32* %17
br label %98
74: ; preds = %66
br label %75
75: ; preds = %93, %74
%76 = load float, float* %12
%77 = load float, float* %10
%78 = fmul float %76, %77
%79 = fptosi float %78 to i32
store i32 %79, i32* %15
%80 = load float, float* %12
%81 = load float, float* %11
%82 = fmul float %80, %81
%83 = fptosi float %82 to i32
store i32 %83, i32* %16
%84 = load i32, i32* %15
%85 = load i32, i32* %16
%86 = icmp eq i32 %84, %85
br i1 %86, label %87, label %92
87: ; preds = %75
%88 = load float, float* %12
%89 = fmul float %88, 1.000000e+01
store float %89, float* %12
%90 = load i32, i32* %18
%91 = add nsw i32 %90, 1
store i32 %91, i32* %18
br label %92
92: ; preds = %87, %75
br label %93
93: ; preds = %92
%94 = load i32, i32* %15
%95 = load i32, i32* %16
%96 = icmp eq i32 %94, %95
br i1 %96, label %75, label %97
97: ; preds = %93
br label %98
98: ; preds = %97, %73
%99 = load float, float* %11
store float %99, float* %10
br label %100
100: ; preds = %98, %61
br label %101
101: ; preds = %100
%102 = load i32, i32* %13
%103 = add nsw i32 %102, -1
store i32 %103, i32* %13
br label %56
104: ; preds = %56
%105 = load i32, i32* %18
store i32 %105, i32* %5
br label %106
106: ; preds = %104, %24
%107 = load i32, i32* %5
ret i32 %107
}

LLVM IR (-O2):
define i32 @calc_label_digits(float %0, float %1, float %2, float %3) {
%5 = fcmp oeq float %2, 0.000000e+00
%6 = fcmp oeq float %3, 0.000000e+00
%7 = or i1 %5, %6
br i1 %7, label %60, label %8
8: ; preds = %4
%9 = fdiv float %3, %2
%10 = fcmp olt float %9, 0.000000e+00
br i1 %10, label %11, label %14
11: ; preds = %8
%12 = fadd float %9, -5.000000e-01
%13 = tail call i64 @ceil(float %12)
br label %17
14: ; preds = %8
%15 = fadd float %9, 5.000000e-01
%16 = tail call i64 @floor(float %15)
br label %17
17: ; preds = %14, %11
%18 = phi i64 [ %13, %11 ], [ %16, %14 ]
%19 = trunc i64 %18 to i32
%20 = icmp sgt i32 %19, 1
%21 = select i1 %20, i32 %19, i32 1
%22 = fptosi float %0 to i32
%23 = sitofp i32 %22 to float
%24 = fcmp ult float %23, %1
br i1 %24, label %60, label %25
25: ; preds = %17, %52
%26 = phi float [ %58, %52 ], [ %23, %17 ]
%27 = phi i32 [ %56, %52 ], [ 0, %17 ]
%28 = phi i32 [ %55, %52 ], [ 0, %17 ]
%29 = phi float [ %54, %52 ], [ undef, %17 ]
%30 = phi i32 [ %57, %52 ], [ %22, %17 ]
%31 = phi float [ %53, %52 ], [ 1.000000e+00, %17 ]
%32 = srem i32 %30, %21
%33 = icmp eq i32 %32, 0
br i1 %33, label %34, label %52
34: ; preds = %25
%35 = fdiv float %26, %3
%36 = icmp eq i32 %28, 0
br i1 %36, label %52, label %37
37: ; preds = %34
%38 = insertelement <2 x float> poison, float %29, i64 0
%39 = insertelement <2 x float> %38, float %35, i64 1
br label %40
40: ; preds = %37, %40
%41 = phi float [ %50, %40 ], [ %31, %37 ]
%42 = phi i32 [ %51, %40 ], [ %27, %37 ]
%43 = insertelement <2 x float> poison, float %41, i64 0
%44 = shufflevector <2 x float> %43, <2 x float> poison, <2 x i32> zeroinitializer
%45 = fmul <2 x float> %39, %44
%46 = fptosi <2 x float> %45 to <2 x i32>
%47 = shufflevector <2 x i32> %46, <2 x i32> poison, <2 x i32> <i32 1, i32 undef>
%48 = icmp eq <2 x i32> %47, %46
%49 = extractelement <2 x i1> %48, i64 0
%50 = fmul float %41, 1.000000e+01
%51 = add nsw i32 %42, 1
br i1 %49, label %40, label %52
52: ; preds = %40, %34, %25
%53 = phi float [ %31, %25 ], [ %31, %34 ], [ %41, %40 ]
%54 = phi float [ %29, %25 ], [ %35, %34 ], [ %35, %40 ]
%55 = phi i32 [ %28, %25 ], [ 1, %34 ], [ %28, %40 ]
%56 = phi i32 [ %27, %25 ], [ %27, %34 ], [ %42, %40 ]
%57 = add nsw i32 %30, -1
%58 = sitofp i32 %57 to float
%59 = fcmp ult float %58, %1
br i1 %59, label %60, label %25
60: ; preds = %52, %17, %4
%61 = phi i32 [ 1, %4 ], [ 0, %17 ], [ %56, %52 ]
ret i32 %61
}