#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 expm(double, double);
double pow(int, double);
double series(int m, int id) {
  int k;
  double ak, p, s, t;
  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    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%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,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
call   1040 <expm@plt>
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
divsd  -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
cvttsd2si -0x28(%rbp),%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%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,-0x30(%rbp)
movsd  0xdf6(%rip),%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    121a <series+0xfa>
jmp    124e <series+0x12e>
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
cvttsd2si -0x28(%rbp),%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    11c1 <series+0xa1>
movsd  -0x28(%rbp),%xmm0
add    $0x30,%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
store i32 %0, i32* %3
store i32 %1, i32* %4
store double 0.000000e+00, double* %8
store i32 0, i32* %5
br label %10
10: ; preds = %37, %2
%11 = load i32, i32* %5
%12 = load i32, i32* %4
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %40
14: ; preds = %10
%15 = load i32, i32* %5
%16 = mul nsw i32 8, %15
%17 = load i32, i32* %3
%18 = add nsw i32 %16, %17
%19 = sitofp i32 %18 to double
store double %19, double* %6
%20 = load i32, i32* %4
%21 = load i32, i32* %5
%22 = sub nsw i32 %20, %21
%23 = sitofp i32 %22 to double
store double %23, double* %7
%24 = load double, double* %7
%25 = load double, double* %6
%26 = call double @expm(double %24, double %25)
store double %26, double* %9
%27 = load double, double* %8
%28 = load double, double* %9
%29 = load double, double* %6
%30 = fdiv double %28, %29
%31 = fadd double %27, %30
store double %31, double* %8
%32 = load double, double* %8
%33 = load double, double* %8
%34 = fptosi double %33 to i32
%35 = sitofp i32 %34 to double
%36 = fsub double %32, %35
store double %36, double* %8
br label %37
37: ; preds = %14
%38 = load i32, i32* %5
%39 = add nsw i32 %38, 1
store i32 %39, i32* %5
br label %10
40: ; preds = %10
%41 = load i32, i32* %4
store i32 %41, i32* %5
br label %42
42: ; preds = %72, %40
%43 = load i32, i32* %5
%44 = load i32, i32* %4
%45 = add nsw i32 %44, 100
%46 = icmp sle i32 %43, %45
br i1 %46, label %47, label %75
47: ; preds = %42
%48 = load i32, i32* %5
%49 = mul nsw i32 8, %48
%50 = load i32, i32* %3
%51 = add nsw i32 %49, %50
%52 = sitofp i32 %51 to double
store double %52, double* %6
%53 = load i32, i32* %4
%54 = load i32, i32* %5
%55 = sub nsw i32 %53, %54
%56 = sitofp i32 %55 to double
%57 = call double @pow(i32 16, double %56)
%58 = load double, double* %6
%59 = fdiv double %57, %58
store double %59, double* %9
%60 = load double, double* %9
%61 = fcmp olt double %60, 1.000000e-17
br i1 %61, label %62, label %63
62: ; preds = %47
br label %75
63: ; preds = %47
%64 = load double, double* %8
%65 = load double, double* %9
%66 = fadd double %64, %65
store double %66, double* %8
%67 = load double, double* %8
%68 = load double, double* %8
%69 = fptosi double %68 to i32
%70 = sitofp i32 %69 to double
%71 = fsub double %67, %70
store double %71, double* %8
br label %72
72: ; preds = %63
%73 = load i32, i32* %5
%74 = add nsw i32 %73, 1
store i32 %74, i32* %5
br label %42
75: ; preds = %62, %42
%76 = load double, double* %8
ret double %76
}

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
}