#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__ fabs(float);
void eulsum(float *sum, float term, int jterm, float wksp[]) {
  int j;
  static int nterm;
  float tmp, dum;
  if (jterm == 1) {
    nterm = 1;
    *sum = 0.5 * (wksp[1] = term);
  } else {
    tmp = wksp[1];
    wksp[1] = term;
    for (j = 1; j <= nterm - 1; j++) {
      dum = wksp[j + 1];
      wksp[j + 1] = 0.5 * (wksp[j] + tmp);
      tmp = dum;
    }
    wksp[nterm + 1] = 0.5 * (wksp[nterm] + tmp);
    if (fabs(wksp[nterm + 1]) <= fabs(wksp[nterm]))
      *sum += (0.5 * wksp[++nterm]);
    else
      *sum += wksp[nterm + 1];
  }
}


x86 assembly (clang, -O0):
<eulsum>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
cmpl   $0x1,-0x10(%rbp)
jne    116b <eulsum+0x5b>
movl   $0x1,0x2ef0(%rip)
movss  -0xc(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,0x4(%rax)
cvtss2sd %xmm0,%xmm1
movsd  0xeaa(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,(%rax)
jmp    12f8 <eulsum+0x1e8>
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0xc(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,0x4(%rax)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    0x2e95(%rip),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1209 <eulsum+0xf9>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe29(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    118e <eulsum+0x7e>
mov    -0x18(%rbp),%rax
movslq 0x2e18(%rip),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdd6(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
mov    0x2df0(%rip),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
mov    0x2ddb(%rip),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
call   1030 <fabs@plt>
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq 0x2dbc(%rip),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1030 <fabs@plt>
mov    %rax,%rcx
mov    -0x30(%rbp),%rax
cmp    %rcx,%rax
jg     12d2 <eulsum+0x1c2>
mov    -0x18(%rbp),%rax
mov    0x2d98(%rip),%ecx
add    $0x1,%ecx
mov    %ecx,0x2d8f(%rip)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xd43(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax)
jmp    12f3 <eulsum+0x1e3>
mov    -0x18(%rbp),%rax
mov    0x2d50(%rip),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
addss  (%rax),%xmm0
movss  %xmm0,(%rax)
jmp    12f8 <eulsum+0x1e8>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<eulsum>:
push   %r15
push   %r14
push   %rbx
mov    %rdx,%rbx
mov    %rdi,%r14
cmp    $0x1,%esi
jne    113c <eulsum+0x2c>
movl   $0x1,0x2f02(%rip)
movss  %xmm0,0x4(%rbx)
mulss  0xec9(%rip),%xmm0
jmp    12a3 <eulsum+0x193>
movss  0x4(%rbx),%xmm1
movss  %xmm0,0x4(%rbx)
movslq 0x2edf(%rip),%rcx
mov    %ecx,%eax
cmp    $0x2,%rcx
jl     122b <eulsum+0x11b>
lea    -0x1(%rax),%rdi
lea    -0x2(%rax),%rsi
mov    %edi,%edx
and    $0x3,%edx
cmp    $0x3,%rsi
jae    117b <eulsum+0x6b>
mov    $0x1,%esi
test   %rdx,%rdx
jne    11f2 <eulsum+0xe2>
jmp    122b <eulsum+0x11b>
and    $0xfffffffffffffffc,%rdi
xor    %esi,%esi
movss  0xe77(%rip),%xmm2
nopl   0x0(%rax)
addss  %xmm1,%xmm0
mulss  %xmm2,%xmm0
movss  0x8(%rbx,%rsi,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm0,0x8(%rbx,%rsi,4)
mulss  %xmm2,%xmm1
movss  0xc(%rbx,%rsi,4),%xmm3
addss  %xmm1,%xmm3
movss  %xmm1,0xc(%rbx,%rsi,4)
mulss  %xmm2,%xmm3
movss  0x10(%rbx,%rsi,4),%xmm0
addss  %xmm3,%xmm0
movss  %xmm3,0x10(%rbx,%rsi,4)
movss  0x14(%rbx,%rsi,4),%xmm1
mulss  %xmm2,%xmm0
movss  %xmm0,0x14(%rbx,%rsi,4)
add    $0x4,%rsi
cmp    %rsi,%rdi
jne    1190 <eulsum+0x80>
add    $0x1,%rsi
test   %rdx,%rdx
je     122b <eulsum+0x11b>
lea    (%rbx,%rsi,4),%rsi
add    $0x4,%rsi
xor    %edi,%edi
movss  0xdfc(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addss  %xmm1,%xmm0
movss  (%rsi,%rdi,4),%xmm1
mulss  %xmm2,%xmm0
movss  %xmm0,(%rsi,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rdx
jne    1210 <eulsum+0x100>
addss  (%rbx,%rcx,4),%xmm1
mulss  0xdc8(%rip),%xmm1
cltq
movss  %xmm1,0x4(%rbx,%rax,4)
movaps %xmm1,%xmm0
call   1030 <fabs@plt>
mov    %rax,%r15
movslq 0x2dda(%rip),%rax
movss  (%rbx,%rax,4),%xmm0
call   1030 <fabs@plt>
movslq 0x2dc9(%rip),%rcx
lea    0x1(%rcx),%edx
cmp    %rax,%r15
jle    127a <eulsum+0x16a>
movslq %edx,%rax
movss  (%rbx,%rax,4),%xmm0
addss  (%r14),%xmm0
jmp    12a3 <eulsum+0x193>
mov    %edx,0x2dac(%rip)
movss  0x4(%rbx,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movss  (%r14),%xmm1
mulsd  0xd71(%rip),%xmm0
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r14)
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
@eulsum.nterm = internal global i32 0
define void @eulsum(float* %0, float %1, i32 %2, float* %3) {
%5 = alloca float*
%6 = alloca float
%7 = alloca i32
%8 = alloca float*
%9 = alloca i32
%10 = alloca float
%11 = alloca float
store float* %0, float** %5
store float %1, float* %6
store i32 %2, i32* %7
store float* %3, float** %8
%12 = load i32, i32* %7
%13 = icmp eq i32 %12, 1
br i1 %13, label %14, label %22
14: ; preds = %4
store i32 1, i32* @eulsum.nterm
%15 = load float, float* %6
%16 = load float*, float** %8
%17 = getelementptr inbounds float, float* %16, i64 1
store float %15, float* %17
%18 = fpext float %15 to double
%19 = fmul double 5.000000e-01, %18
%20 = fptrunc double %19 to float
%21 = load float*, float** %5
store float %20, float* %21
br label %114
22: ; preds = %4
%23 = load float*, float** %8
%24 = getelementptr inbounds float, float* %23, i64 1
%25 = load float, float* %24
store float %25, float* %10
%26 = load float, float* %6
%27 = load float*, float** %8
%28 = getelementptr inbounds float, float* %27, i64 1
store float %26, float* %28
store i32 1, i32* %9
br label %29
29: ; preds = %57, %22
%30 = load i32, i32* %9
%31 = load i32, i32* @eulsum.nterm
%32 = sub nsw i32 %31, 1
%33 = icmp sle i32 %30, %32
br i1 %33, label %34, label %60
34: ; preds = %29
%35 = load float*, float** %8
%36 = load i32, i32* %9
%37 = add nsw i32 %36, 1
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds float, float* %35, i64 %38
%40 = load float, float* %39
store float %40, float* %11
%41 = load float*, float** %8
%42 = load i32, i32* %9
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
%45 = load float, float* %44
%46 = load float, float* %10
%47 = fadd float %45, %46
%48 = fpext float %47 to double
%49 = fmul double 5.000000e-01, %48
%50 = fptrunc double %49 to float
%51 = load float*, float** %8
%52 = load i32, i32* %9
%53 = add nsw i32 %52, 1
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %51, i64 %54
store float %50, float* %55
%56 = load float, float* %11
store float %56, float* %10
br label %57
57: ; preds = %34
%58 = load i32, i32* %9
%59 = add nsw i32 %58, 1
store i32 %59, i32* %9
br label %29
60: ; preds = %29
%61 = load float*, float** %8
%62 = load i32, i32* @eulsum.nterm
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %61, i64 %63
%65 = load float, float* %64
%66 = load float, float* %10
%67 = fadd float %65, %66
%68 = fpext float %67 to double
%69 = fmul double 5.000000e-01, %68
%70 = fptrunc double %69 to float
%71 = load float*, float** %8
%72 = load i32, i32* @eulsum.nterm
%73 = add nsw i32 %72, 1
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds float, float* %71, i64 %74
store float %70, float* %75
%76 = load float*, float** %8
%77 = load i32, i32* @eulsum.nterm
%78 = add nsw i32 %77, 1
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %76, i64 %79
%81 = load float, float* %80
%82 = call i64 @fabs(float %81)
%83 = load float*, float** %8
%84 = load i32, i32* @eulsum.nterm
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds float, float* %83, i64 %85
%87 = load float, float* %86
%88 = call i64 @fabs(float %87)
%89 = icmp sle i64 %82, %88
br i1 %89, label %90, label %103
90: ; preds = %60
%91 = load float*, float** %8
%92 = load i32, i32* @eulsum.nterm
%93 = add nsw i32 %92, 1
store i32 %93, i32* @eulsum.nterm
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds float, float* %91, i64 %94
%96 = load float, float* %95
%97 = fpext float %96 to double
%98 = load float*, float** %5
%99 = load float, float* %98
%100 = fpext float %99 to double
%101 = call double @llvm.fmuladd.f64(double 5.000000e-01, double %97, double %100)
%102 = fptrunc double %101 to float
store float %102, float* %98
br label %113
103: ; preds = %60
%104 = load float*, float** %8
%105 = load i32, i32* @eulsum.nterm
%106 = add nsw i32 %105, 1
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds float, float* %104, i64 %107
%109 = load float, float* %108
%110 = load float*, float** %5
%111 = load float, float* %110
%112 = fadd float %111, %109
store float %112, float* %110
br label %113
113: ; preds = %103, %90
br label %114
114: ; preds = %113, %14
ret void
}

LLVM IR (-O2):
@eulsum.nterm = internal unnamed_addr global i32 0
define void @eulsum(float* nocapture %0, float %1, i32 %2, float* nocapture %3) {
%5 = icmp eq i32 %2, 1
br i1 %5, label %6, label %9
6: ; preds = %4
store i32 1, i32* @eulsum.nterm
%7 = getelementptr inbounds float, float* %3, i64 1
store float %1, float* %7
%8 = fmul float %1, 5.000000e-01
br label %101
9: ; preds = %4
%10 = getelementptr inbounds float, float* %3, i64 1
%11 = load float, float* %10
store float %1, float* %10
%12 = load i32, i32* @eulsum.nterm
%13 = icmp sgt i32 %12, 1
br i1 %13, label %14, label %67
14: ; preds = %9
%15 = zext i32 %12 to i64
%16 = add nsw i64 %15, -1
%17 = add nsw i64 %15, -2
%18 = and i64 %16, 3
%19 = icmp ult i64 %17, 3
br i1 %19, label %49, label %20
20: ; preds = %14
%21 = and i64 %16, -4
br label %22
22: ; preds = %22, %20
%23 = phi float [ %1, %20 ], [ %46, %22 ]
%24 = phi i64 [ 1, %20 ], [ %42, %22 ]
%25 = phi float [ %11, %20 ], [ %44, %22 ]
%26 = phi i64 [ 0, %20 ], [ %47, %22 ]
%27 = add nuw nsw i64 %24, 1
%28 = getelementptr inbounds float, float* %3, i64 %27
%29 = load float, float* %28
%30 = fadd float %25, %23
%31 = fmul float %30, 5.000000e-01
store float %31, float* %28
%32 = add nuw nsw i64 %24, 2
%33 = getelementptr inbounds float, float* %3, i64 %32
%34 = load float, float* %33
%35 = fadd float %29, %31
%36 = fmul float %35, 5.000000e-01
store float %36, float* %33
%37 = add nuw nsw i64 %24, 3
%38 = getelementptr inbounds float, float* %3, i64 %37
%39 = load float, float* %38
%40 = fadd float %34, %36
%41 = fmul float %40, 5.000000e-01
store float %41, float* %38
%42 = add nuw nsw i64 %24, 4
%43 = getelementptr inbounds float, float* %3, i64 %42
%44 = load float, float* %43
%45 = fadd float %39, %41
%46 = fmul float %45, 5.000000e-01
store float %46, float* %43
%47 = add i64 %26, 4
%48 = icmp eq i64 %47, %21
br i1 %48, label %49, label %22
49: ; preds = %22, %14
%50 = phi float [ undef, %14 ], [ %44, %22 ]
%51 = phi float [ %1, %14 ], [ %46, %22 ]
%52 = phi i64 [ 1, %14 ], [ %42, %22 ]
%53 = phi float [ %11, %14 ], [ %44, %22 ]
%54 = icmp eq i64 %18, 0
br i1 %54, label %67, label %55
55: ; preds = %49, %55
%56 = phi float [ %64, %55 ], [ %51, %49 ]
%57 = phi i64 [ %60, %55 ], [ %52, %49 ]
%58 = phi float [ %62, %55 ], [ %53, %49 ]
%59 = phi i64 [ %65, %55 ], [ 0, %49 ]
%60 = add nuw nsw i64 %57, 1
%61 = getelementptr inbounds float, float* %3, i64 %60
%62 = load float, float* %61
%63 = fadd float %58, %56
%64 = fmul float %63, 5.000000e-01
store float %64, float* %61
%65 = add i64 %59, 1
%66 = icmp eq i64 %65, %18
br i1 %66, label %67, label %55
67: ; preds = %49, %55, %9
%68 = phi float [ %11, %9 ], [ %50, %49 ], [ %62, %55 ]
%69 = sext i32 %12 to i64
%70 = getelementptr inbounds float, float* %3, i64 %69
%71 = load float, float* %70
%72 = fadd float %68, %71
%73 = fmul float %72, 5.000000e-01
%74 = add nsw i32 %12, 1
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds float, float* %3, i64 %75
store float %73, float* %76
%77 = tail call i64 @fabs(float %73)
%78 = load i32, i32* @eulsum.nterm
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %3, i64 %79
%81 = load float, float* %80
%82 = tail call i64 @fabs(float %81)
%83 = icmp sgt i64 %77, %82
%84 = load i32, i32* @eulsum.nterm
%85 = add nsw i32 %84, 1
br i1 %83, label %95, label %86
86: ; preds = %67
store i32 %85, i32* @eulsum.nterm
%87 = sext i32 %85 to i64
%88 = getelementptr inbounds float, float* %3, i64 %87
%89 = load float, float* %88
%90 = fpext float %89 to double
%91 = load float, float* %0
%92 = fpext float %91 to double
%93 = tail call double @llvm.fmuladd.f64(double %90, double 5.000000e-01, double %92)
%94 = fptrunc double %93 to float
br label %101
95: ; preds = %67
%96 = sext i32 %85 to i64
%97 = getelementptr inbounds float, float* %3, i64 %96
%98 = load float, float* %97
%99 = load float, float* %0
%100 = fadd float %98, %99
br label %101
101: ; preds = %86, %95, %6
%102 = phi float [ %94, %86 ], [ %100, %95 ], [ %8, %6 ]
store float %102, float* %0
ret void
}