void de_emp(float Sn_de[], float Sn[], float *mem, int Nsam) {
  int i;
  for (i = 0; i < Nsam; i++) {
    Sn_de[i] = Sn[i] + 0.94 * mem[0];
    mem[0] = Sn_de[i];
  }
}


x86 assembly (clang, -O0):
<de_emp>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1187 <de_emp+0x87>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeb5(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    111a <de_emp+0x1a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<de_emp>:
test   %ecx,%ecx
jle    11ab <de_emp+0xab>
mov    %ecx,%r8d
movss  (%rdx),%xmm0
cmp    $0x1,%ecx
jne    1118 <de_emp+0x18>
xor    %eax,%eax
jmp    117f <de_emp+0x7f>
mov    %r8d,%ecx
and    $0xfffffffe,%ecx
xor    %eax,%eax
movsd  0xed8(%rip),%xmm1
nopl   0x0(%rax,%rax,1)
movss  (%rsi,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rdi,%rax,4)
movss  %xmm0,(%rdx)
movss  0x4(%rsi,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rdi,%rax,4)
movss  %xmm0,(%rdx)
add    $0x2,%rax
cmp    %rax,%rcx
jne    1130 <de_emp+0x30>
test   $0x1,%r8b
je     11ab <de_emp+0xab>
movss  (%rsi,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
cvtss2sd %xmm0,%xmm0
mulsd  0xe66(%rip),%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rdi,%rax,4)
movss  %xmm0,(%rdx)
ret

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

LLVM IR (-O2):
define void @de_emp(float* nocapture writeonly %0, float* nocapture %1, float* nocapture %2, i32 %3) {
%5 = icmp sgt i32 %3, 0
br i1 %5, label %6, label %47
6: ; preds = %4
%7 = zext i32 %3 to i64
%8 = load float, float* %2
%9 = and i64 %7, 1
%10 = icmp eq i32 %3, 1
br i1 %10, label %35, label %11
11: ; preds = %6
%12 = and i64 %7, 4294967294
br label %13
13: ; preds = %13, %11
%14 = phi float [ %8, %11 ], [ %30, %13 ]
%15 = phi i64 [ 0, %11 ], [ %32, %13 ]
%16 = phi i64 [ 0, %11 ], [ %33, %13 ]
%17 = getelementptr inbounds float, float* %1, i64 %15
%18 = load float, float* %17
%19 = fpext float %18 to double
%20 = fpext float %14 to double
%21 = tail call double @llvm.fmuladd.f64(double %20, double 0x3FEE147AE147AE14, double %19)
%22 = fptrunc double %21 to float
%23 = getelementptr inbounds float, float* %0, i64 %15
store float %22, float* %23
store float %22, float* %2
%24 = or i64 %15, 1
%25 = getelementptr inbounds float, float* %1, i64 %24
%26 = load float, float* %25
%27 = fpext float %26 to double
%28 = fpext float %22 to double
%29 = tail call double @llvm.fmuladd.f64(double %28, double 0x3FEE147AE147AE14, double %27)
%30 = fptrunc double %29 to float
%31 = getelementptr inbounds float, float* %0, i64 %24
store float %30, float* %31
store float %30, float* %2
%32 = add nuw nsw i64 %15, 2
%33 = add i64 %16, 2
%34 = icmp eq i64 %33, %12
br i1 %34, label %35, label %13
35: ; preds = %13, %6
%36 = phi float [ %8, %6 ], [ %30, %13 ]
%37 = phi i64 [ 0, %6 ], [ %32, %13 ]
%38 = icmp eq i64 %9, 0
br i1 %38, label %47, label %39
39: ; preds = %35
%40 = getelementptr inbounds float, float* %1, i64 %37
%41 = load float, float* %40
%42 = fpext float %41 to double
%43 = fpext float %36 to double
%44 = tail call double @llvm.fmuladd.f64(double %43, double 0x3FEE147AE147AE14, double %42)
%45 = fptrunc double %44 to float
%46 = getelementptr inbounds float, float* %0, i64 %37
store float %45, float* %46
store float %45, float* %2
br label %47
47: ; preds = %39, %35, %4
ret void
}