void smoothed_diff_int(int *data, unsigned int datalen) {
  unsigned int k;
  double memory = 0.0;
  for (k = 0; k < datalen - 1; k++) {
    memory = 0.7 * memory + (1 - 0.7) * (double)(data[k + 1] - data[k]);
    data[k] = (int)memory;
  };
  data[datalen - 1] = data[datalen - 2];
}


x86 assembly (clang, -O0):
<smoothed_diff_int>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jae    118b <smoothed_diff_int+0x8b>
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,%ecx
mov    (%rax,%rcx,4),%eax
mov    -0x8(%rbp),%rcx
mov    -0x10(%rbp),%edx
sub    (%rcx,%rdx,4),%eax
cvtsi2sd %eax,%xmm2
movsd  0xeb3(%rip),%xmm1
mulsd  %xmm2,%xmm1
movsd  0xe9f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
cvttsd2si -0x18(%rbp),%edx
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    111a <smoothed_diff_int+0x1a>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
mov    %ecx,%ecx
mov    (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
mov    %edx,(%rax,%rcx,4)
pop    %rbp
ret

x86 assembly (clang, -O2):
<smoothed_diff_int>:
mov    %esi,%ecx
add    $0xffffffff,%ecx
je     111d <smoothed_diff_int+0x1d>
mov    %ecx,%r8d
mov    (%rdi),%r10d
cmp    $0x1,%ecx
jne    1125 <smoothed_diff_int+0x25>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    11a6 <smoothed_diff_int+0xa6>
xor    %r8d,%r8d
jmp    11d5 <smoothed_diff_int+0xd5>
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
movsd  0xec6(%rip),%xmm1
movsd  0xec6(%rip),%xmm2
mov    %r10d,%edx
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x4(%rdi,%rcx,4),%r11d
mov    0x8(%rdi,%rcx,4),%r10d
mov    %r11d,%eax
sub    %edx,%eax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm1,%xmm3
mulsd  %xmm2,%xmm0
addsd  %xmm3,%xmm0
cvttsd2si %xmm0,%eax
mov    %r10d,%edx
sub    %r11d,%edx
xorps  %xmm3,%xmm3
cvtsi2sd %edx,%xmm3
mulsd  %xmm1,%xmm3
mulsd  %xmm2,%xmm0
addsd  %xmm3,%xmm0
cvttsd2si %xmm0,%edx
mov    %eax,(%rdi,%rcx,4)
mov    %edx,0x4(%rdi,%rcx,4)
add    $0x2,%rcx
mov    %r10d,%edx
cmp    %rcx,%r9
jne    1150 <smoothed_diff_int+0x50>
test   $0x1,%r8b
je     11d5 <smoothed_diff_int+0xd5>
mov    0x4(%rdi,%rcx,4),%eax
sub    %r10d,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  0xe3e(%rip),%xmm1
mulsd  0xe3e(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,(%rdi,%rcx,4)
add    $0xfffffffe,%esi
mov    (%rdi,%rsi,4),%eax
mov    %eax,(%rdi,%r8,4)
ret

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

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