void cumsum(double y[], double x[], int NumCel) {
  int i;
  y[0] = x[0];
  for (i = 1; i < NumCel; i++) {
    y[i] = y[i - 1] + x[i];
  }
}


x86 assembly (clang, -O0):
<cumsum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    116c <cumsum+0x6c>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1126 <cumsum+0x26>
pop    %rbp
ret

x86 assembly (clang, -O2):
<cumsum>:
movsd  (%rsi),%xmm0
movsd  %xmm0,(%rdi)
cmp    $0x2,%edx
jl     11a3 <cumsum+0xa3>
mov    %edx,%r8d
lea    -0x1(%r8),%rdx
add    $0xfffffffffffffffe,%r8
mov    %edx,%eax
and    $0x3,%eax
mov    $0x1,%ecx
cmp    $0x3,%r8
jb     117d <cumsum+0x7d>
and    $0xfffffffffffffffc,%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  0x8(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rdi,%rcx,8)
addsd  0x10(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rdi,%rcx,8)
addsd  0x18(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rdi,%rcx,8)
addsd  0x20(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x20(%rdi,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1140 <cumsum+0x40>
add    $0x1,%rcx
test   %rax,%rax
je     11a3 <cumsum+0xa3>
lea    (%rdi,%rcx,8),%rdx
lea    (%rsi,%rcx,8),%rcx
xor    %esi,%esi
nopl   0x0(%rax)
addsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1190 <cumsum+0x90>
ret

LLVM IR (-O0):
define void @cumsum(double* %0, double* %1, i32 %2) {
%4 = alloca double*
%5 = alloca double*
%6 = alloca i32
%7 = alloca i32
store double* %0, double** %4
store double* %1, double** %5
store i32 %2, i32* %6
%8 = load double*, double** %5
%9 = getelementptr inbounds double, double* %8, i64 0
%10 = load double, double* %9
%11 = load double*, double** %4
%12 = getelementptr inbounds double, double* %11, i64 0
store double %10, double* %12
store i32 1, i32* %7
br label %13
13: ; preds = %34, %3
%14 = load i32, i32* %7
%15 = load i32, i32* %6
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %37
17: ; preds = %13
%18 = load double*, double** %4
%19 = load i32, i32* %7
%20 = sub nsw i32 %19, 1
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %18, i64 %21
%23 = load double, double* %22
%24 = load double*, double** %5
%25 = load i32, i32* %7
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = load double, double* %27
%29 = fadd double %23, %28
%30 = load double*, double** %4
%31 = load i32, i32* %7
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
store double %29, double* %33
br label %34
34: ; preds = %17
%35 = load i32, i32* %7
%36 = add nsw i32 %35, 1
store i32 %36, i32* %7
br label %13
37: ; preds = %13
ret void
}

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