#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 */
float *zeros_f(int);
float *cumtrapz(float *buf, int n, float dh) {
  float *s;
  int i;
  s = zeros_f(n);
  for (i = 1; i < n; i++) {
    s[i] = s[i - 1] + 0.5 * dh * (buf[i - 1] + buf[i]);
  }
  return s;
}


x86 assembly (clang, -O0):
<cumtrapz>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
mov    -0xc(%rbp),%edi
call   1030 <zeros_f@plt>
mov    %rax,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b8 <cumtrapz+0xa8>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm2
movsd  0xe96(%rip),%xmm0
mulsd  %xmm2,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1137 <cumtrapz+0x27>
mov    -0x18(%rbp),%rax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cumtrapz>:
push   %rbp
push   %rbx
push   %rax
movss  %xmm0,0x4(%rsp)
mov    %esi,%ebp
mov    %rdi,%rbx
mov    %esi,%edi
call   1030 <zeros_f@plt>
cmp    $0x2,%ebp
jl     11fa <cumtrapz+0xea>
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  0xec0(%rip),%xmm0
mov    %ebp,%edx
movss  (%rax),%xmm1
add    $0xffffffffffffffff,%rdx
mov    $0x1,%ecx
cmp    $0x2,%ebp
je     11d1 <cumtrapz+0xc1>
mov    %rdx,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
mov    $0x1,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  -0x4(%rbx,%rcx,4),%xmm2
addss  (%rbx,%rcx,4),%xmm2
cvtss2sd %xmm1,%xmm1
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
movss  %xmm1,(%rax,%rcx,4)
cvtss2sd %xmm1,%xmm1
movss  (%rbx,%rcx,4),%xmm2
addss  0x4(%rbx,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
movss  %xmm1,0x4(%rax,%rcx,4)
lea    (%rsi,%rcx,1),%rdi
add    $0x2,%rdi
add    $0x2,%rcx
cmp    $0x1,%rdi
jne    1170 <cumtrapz+0x60>
test   $0x1,%dl
je     11fa <cumtrapz+0xea>
cvtss2sd %xmm1,%xmm1
movss  -0x4(%rbx,%rcx,4),%xmm2
addss  (%rbx,%rcx,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
add    $0x8,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define float* @cumtrapz(float* %0, i32 %1, float %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca float
%7 = alloca float*
%8 = alloca i32
store float* %0, float** %4
store i32 %1, i32* %5
store float %2, float* %6
%9 = load i32, i32* %5
%10 = call float* @zeros_f(i32 %9)
store float* %10, float** %7
store i32 1, i32* %8
br label %11
11: ; preds = %45, %3
%12 = load i32, i32* %8
%13 = load i32, i32* %5
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %48
15: ; preds = %11
%16 = load float*, float** %7
%17 = load i32, i32* %8
%18 = sub nsw i32 %17, 1
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds float, float* %16, i64 %19
%21 = load float, float* %20
%22 = fpext float %21 to double
%23 = load float, float* %6
%24 = fpext float %23 to double
%25 = fmul double 5.000000e-01, %24
%26 = load float*, float** %4
%27 = load i32, i32* %8
%28 = sub nsw i32 %27, 1
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %26, i64 %29
%31 = load float, float* %30
%32 = load float*, float** %4
%33 = load i32, i32* %8
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %32, i64 %34
%36 = load float, float* %35
%37 = fadd float %31, %36
%38 = fpext float %37 to double
%39 = call double @llvm.fmuladd.f64(double %25, double %38, double %22)
%40 = fptrunc double %39 to float
%41 = load float*, float** %7
%42 = load i32, i32* %8
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
store float %40, float* %44
br label %45
45: ; preds = %15
%46 = load i32, i32* %8
%47 = add nsw i32 %46, 1
store i32 %47, i32* %8
br label %11
48: ; preds = %11
%49 = load float*, float** %7
ret float* %49
}

LLVM IR (-O2):
define float* @cumtrapz(float* nocapture %0, i32 %1, float %2) {
%4 = tail call float* @zeros_f(i32 %1)
%5 = icmp sgt i32 %1, 1
br i1 %5, label %6, label %61
6: ; preds = %3
%7 = fpext float %2 to double
%8 = fmul double %7, 5.000000e-01
%9 = zext i32 %1 to i64
%10 = load float, float* %4
%11 = add nsw i64 %9, -1
%12 = and i64 %11, 1
%13 = icmp eq i32 %1, 2
br i1 %13, label %45, label %14
14: ; preds = %6
%15 = and i64 %11, -2
br label %16
16: ; preds = %16, %14
%17 = phi float [ %10, %14 ], [ %40, %16 ]
%18 = phi i64 [ 1, %14 ], [ %42, %16 ]
%19 = phi i64 [ 0, %14 ], [ %43, %16 ]
%20 = add nsw i64 %18, -1
%21 = fpext float %17 to double
%22 = getelementptr inbounds float, float* %0, i64 %20
%23 = load float, float* %22
%24 = getelementptr inbounds float, float* %0, i64 %18
%25 = load float, float* %24
%26 = fadd float %23, %25
%27 = fpext float %26 to double
%28 = tail call double @llvm.fmuladd.f64(double %8, double %27, double %21)
%29 = fptrunc double %28 to float
%30 = getelementptr inbounds float, float* %4, i64 %18
store float %29, float* %30
%31 = add nuw nsw i64 %18, 1
%32 = fpext float %29 to double
%33 = getelementptr inbounds float, float* %0, i64 %18
%34 = load float, float* %33
%35 = getelementptr inbounds float, float* %0, i64 %31
%36 = load float, float* %35
%37 = fadd float %34, %36
%38 = fpext float %37 to double
%39 = tail call double @llvm.fmuladd.f64(double %8, double %38, double %32)
%40 = fptrunc double %39 to float
%41 = getelementptr inbounds float, float* %4, i64 %31
store float %40, float* %41
%42 = add nuw nsw i64 %18, 2
%43 = add i64 %19, 2
%44 = icmp eq i64 %43, %15
br i1 %44, label %45, label %16
45: ; preds = %16, %6
%46 = phi float [ %10, %6 ], [ %40, %16 ]
%47 = phi i64 [ 1, %6 ], [ %42, %16 ]
%48 = icmp eq i64 %12, 0
br i1 %48, label %61, label %49
49: ; preds = %45
%50 = add nsw i64 %47, -1
%51 = fpext float %46 to double
%52 = getelementptr inbounds float, float* %0, i64 %50
%53 = load float, float* %52
%54 = getelementptr inbounds float, float* %0, i64 %47
%55 = load float, float* %54
%56 = fadd float %53, %55
%57 = fpext float %56 to double
%58 = tail call double @llvm.fmuladd.f64(double %8, double %57, double %51)
%59 = fptrunc double %58 to float
%60 = getelementptr inbounds float, float* %4, i64 %47
store float %59, float* %60
br label %61
61: ; preds = %49, %45, %3
ret float* %4
}