#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 */
double pow(int, double);
double series(int m, int id) {
  int k;
  double ak, eps, p, s, t;
  double expm(double x, double y);
  s = 0.;
  for (k = 0; k < id; k++) {
    ak = 8 * k + m;
    p = id - k;
    t = expm(p, ak);
    s = s + t / ak;
    s = s - (int)s;
  }
  for (k = id; k <= id + 100; k++) {
    ak = 8 * k + m;
    t = pow(16., (double)(id - k)) / ak;
    if (t < 1e-17)
      break;
    s = s + t;
    s = s - (int)s;
  }
  return s;
}


x86 assembly (clang, -O0):
<series>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11bb <series+0x9b>
mov    -0xc(%rbp),%eax
shl    $0x3,%eax
add    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
call   1040 <expm@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
cvttsd2si -0x30(%rbp),%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    113d <series+0x1d>
mov    -0x8(%rbp),%eax
mov    %eax,-0xc(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x8(%rbp),%ecx
add    $0x64,%ecx
cmp    %ecx,%eax
jg     124e <series+0x12e>
mov    -0xc(%rbp),%eax
shl    $0x3,%eax
add    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%eax
sub    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    $0x10,%edi
call   1030 <pow@plt>
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xdf6(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    121a <series+0xfa>
jmp    124e <series+0x12e>
movsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
cvttsd2si -0x30(%rbp),%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    11c1 <series+0xa1>
movsd  -0x30(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<series>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %esi,%r12d
mov    %edi,%r14d
test   %esi,%esi
jle    1189 <series+0x69>
xorpd  %xmm0,%xmm0
mov    %r12d,%ebp
mov    %r14d,%ebx
movsd  %xmm0,0x8(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
movsd  %xmm1,0x10(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
call   1040 <expm@plt>
movapd %xmm0,%xmm1
divsd  0x10(%rsp),%xmm1
addsd  0x8(%rsp),%xmm1
cvttpd2dq %xmm1,%xmm0
cvtdq2pd %xmm0,%xmm0
subsd  %xmm0,%xmm1
add    $0x8,%ebx
movapd %xmm1,%xmm0
add    $0xffffffff,%ebp
jne    1140 <series+0x20>
jmp    118d <series+0x6d>
xorpd  %xmm1,%xmm1
lea    0x64(%r12),%r15d
cmp    %r15d,%r12d
cmovg  %r12d,%r15d
not    %r15d
add    %r12d,%r15d
lea    (%r14,%r12,8),%ebx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm1,0x10(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,0x8(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
mov    $0x10,%edi
call   1030 <pow@plt>
divsd  0x8(%rsp),%xmm0
movsd  0xe1e(%rip),%xmm1
ucomisd %xmm0,%xmm1
ja     120b <series+0xeb>
movapd 0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
cvttpd2dq %xmm1,%xmm0
cvtdq2pd %xmm0,%xmm0
subsd  %xmm0,%xmm1
add    $0xffffffff,%ebp
add    $0x8,%ebx
cmp    %ebp,%r15d
jne    11b0 <series+0x90>
jmp    1211 <series+0xf1>
movapd 0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @series(i32 %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %8, label %4
4: ; preds = %8, %2
%5 = phi double [ 0.000000e+00, %2 ], [ %21, %8 ]
%6 = add i32 %1, 100
%7 = call i32 @llvm.smax.i32(i32 %1, i32 %6)
br label %24
8: ; preds = %2, %8
%9 = phi double [ %21, %8 ], [ 0.000000e+00, %2 ]
%10 = phi i32 [ %22, %8 ], [ 0, %2 ]
%11 = shl nsw i32 %10, 3
%12 = add nsw i32 %11, %0
%13 = sitofp i32 %12 to double
%14 = sub nsw i32 %1, %10
%15 = sitofp i32 %14 to double
%16 = tail call double @expm(double %15, double %13)
%17 = fdiv double %16, %13
%18 = fadd double %9, %17
%19 = fptosi double %18 to i32
%20 = sitofp i32 %19 to double
%21 = fsub double %18, %20
%22 = add nuw nsw i32 %10, 1
%23 = icmp eq i32 %22, %1
br i1 %23, label %4, label %8
24: ; preds = %4, %35
%25 = phi double [ %5, %4 ], [ %39, %35 ]
%26 = phi i32 [ %1, %4 ], [ %40, %35 ]
%27 = shl nsw i32 %26, 3
%28 = add nsw i32 %27, %0
%29 = sitofp i32 %28 to double
%30 = sub nsw i32 %1, %26
%31 = sitofp i32 %30 to double
%32 = tail call double @pow(i32 16, double %31)
%33 = fdiv double %32, %29
%34 = fcmp olt double %33, 1.000000e-17
br i1 %34, label %42, label %35
35: ; preds = %24
%36 = fadd double %25, %33
%37 = fptosi double %36 to i32
%38 = sitofp i32 %37 to double
%39 = fsub double %36, %38
%40 = add i32 %26, 1
%41 = icmp eq i32 %26, %7
br i1 %41, label %42, label %24
42: ; preds = %24, %35
%43 = phi double [ %25, %24 ], [ %39, %35 ]
ret double %43
}