void sms_filterArray(double *pFArray, int size1, int size2,
                     double *pFOutArray) {
  int i, j, iPoint, iFrame, size2_2 = size2 - 2, size2_1 = size2 - 1;
  double *pFCurrentArray = pFArray + (size2_1)*size1;
  double fVal, fWeighting, fTotalWeighting, fTmpVal;
  for (i = 0; i < size1; i++) {
    fVal = pFCurrentArray[i];
    fTotalWeighting = 1;
    for (j = 1; j < (size2_2); j++) {
      fWeighting = (double)size2 / (1 + j);
      iPoint = i - (size2_1) + j;
      if (iPoint >= 0) {
        fVal += pFCurrentArray[iPoint] * fWeighting;
        fTotalWeighting += fWeighting;
      }
      iPoint = i + (size2_1)-j;
      if (iPoint < size1) {
        fVal += pFCurrentArray[iPoint] * fWeighting;
        fTotalWeighting += fWeighting;
      }
      iFrame = j;
      fTmpVal = pFArray[iFrame * size1 + i];
      if (fTmpVal) {
        fVal += fTmpVal * fWeighting;
        fTotalWeighting += fWeighting;
      }
    }
    pFOutArray[i] = fVal / fTotalWeighting;
  }
}


x86 assembly (clang, -O0):
<sms_filterArray>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    -0x10(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12cf <sms_filterArray+0x1cf>
mov    -0x38(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe96(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    12aa <sms_filterArray+0x1aa>
cvtsi2sdl -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x1c(%rbp),%eax
sub    -0x30(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     11e3 <sms_filterArray+0xe3>
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x30(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    122e <sms_filterArray+0x12e>
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    126c <sms_filterArray+0x16c>
jp     126c <sms_filterArray+0x16c>
jmp    1297 <sms_filterArray+0x197>
movsd  -0x58(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    129c <sms_filterArray+0x19c>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1176 <sms_filterArray+0x76>
movsd  -0x40(%rbp),%xmm0
divsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1144 <sms_filterArray+0x44>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sms_filterArray>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x8(%rsp)
mov    %rdx,-0x10(%rsp)
test   %esi,%esi
jle    1261 <sms_filterArray+0x161>
mov    -0x10(%rsp),%rcx
lea    -0x1(%rcx),%eax
imul   %esi,%eax
cltq
lea    (%rdi,%rax,8),%r13
cvtsi2sd %ecx,%xmm0
lea    -0x2(%rcx),%r8d
mov    %esi,%edx
mov    $0x2,%r10d
sub    %ecx,%r10d
lea    (%rdi,%rdx,8),%r11
lea    0x0(,%rdx,8),%r15
mov    $0x1,%r12d
sub    %r8,%r12
xor    %r14d,%r14d
movsd  0xea0(%rip),%xmm1
xorpd  %xmm2,%xmm2
movsd  0xe9c(%rip),%xmm3
jmp    119c <sms_filterArray+0x9c>
xchg   %ax,%ax
movapd %xmm1,%xmm4
divsd  %xmm4,%xmm5
mov    -0x8(%rsp),%rax
movsd  %xmm5,(%rax,%r14,8)
add    $0x1,%r14
add    $0x1,%r8d
add    $0x1,%r10d
add    $0x8,%r11
cmp    %rdx,%r14
je     1261 <sms_filterArray+0x161>
movsd  0x0(%r13,%r14,8),%xmm5
cmpl   $0x4,-0x10(%rsp)
jl     1170 <sms_filterArray+0x70>
mov    $0x2,%edi
mov    %r11,%rbx
mov    %r10d,%ebp
mov    %r8d,%eax
movapd %xmm1,%xmm4
jmp    11df <sms_filterArray+0xdf>
xchg   %ax,%ax
addsd  %xmm6,%xmm4
add    $0xffffffff,%eax
add    $0x1,%ebp
add    %r15,%rbx
lea    (%r12,%rdi,1),%rcx
add    $0x1,%rcx
add    $0x1,%rdi
cmp    $0x2,%rcx
je     1174 <sms_filterArray+0x74>
xorps  %xmm7,%xmm7
cvtsi2sd %edi,%xmm7
movapd %xmm0,%xmm6
divsd  %xmm7,%xmm6
test   %ebp,%ebp
js     1208 <sms_filterArray+0x108>
mov    %ebp,%r9d
movsd  0x0(%r13,%r9,8),%xmm7
mulsd  %xmm6,%xmm7
addsd  %xmm7,%xmm5
addsd  %xmm6,%xmm4
cmp    %esi,%eax
jge    1230 <sms_filterArray+0x130>
movslq %eax,%rcx
movsd  0x0(%r13,%rcx,8),%xmm7
mulsd  %xmm6,%xmm7
addsd  %xmm5,%xmm7
addsd  %xmm6,%xmm4
jmp    1234 <sms_filterArray+0x134>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm5,%xmm7
movsd  (%rbx),%xmm5
ucomisd %xmm5,%xmm2
mulsd  %xmm6,%xmm5
addsd  %xmm7,%xmm5
jne    124c <sms_filterArray+0x14c>
jp     124c <sms_filterArray+0x14c>
movapd %xmm7,%xmm5
jne    11c0 <sms_filterArray+0xc0>
jp     11c0 <sms_filterArray+0xc0>
movapd %xmm3,%xmm6
jmp    11c0 <sms_filterArray+0xc0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @sms_filterArray(double* nocapture %0, i32 %1, i32 %2, double* nocapture writeonly %3) {
%5 = add nsw i32 %2, -1
%6 = mul nsw i32 %5, %1
%7 = sext i32 %6 to i64
%8 = getelementptr inbounds double, double* %0, i64 %7
%9 = icmp sgt i32 %1, 0
br i1 %9, label %10, label %74
10: ; preds = %4
%11 = add i32 %2, -2
%12 = icmp sgt i32 %2, 3
%13 = sitofp i32 %2 to double
%14 = zext i32 %1 to i64
%15 = zext i32 %1 to i64
%16 = zext i32 %11 to i64
br label %17
17: ; preds = %10, %67
%18 = phi i64 [ 0, %10 ], [ %72, %67 ]
%19 = getelementptr inbounds double, double* %8, i64 %18
%20 = load double, double* %19
br i1 %12, label %21, label %67
21: ; preds = %17
%22 = trunc i64 %18 to i32
%23 = sub i32 %22, %5
%24 = trunc i64 %18 to i32
%25 = add i32 %5, %24
br label %26
26: ; preds = %21, %54
%27 = phi i64 [ 1, %21 ], [ %31, %54 ]
%28 = phi double [ 1.000000e+00, %21 ], [ %65, %54 ]
%29 = phi double [ %20, %21 ], [ %63, %54 ]
%30 = trunc i64 %27 to i32
%31 = add nuw nsw i64 %27, 1
%32 = trunc i64 %31 to i32
%33 = sitofp i32 %32 to double
%34 = fdiv double %13, %33
%35 = add nsw i32 %23, %30
%36 = icmp sgt i32 %35, -1
br i1 %36, label %37, label %43
37: ; preds = %26
%38 = zext i32 %35 to i64
%39 = getelementptr inbounds double, double* %8, i64 %38
%40 = load double, double* %39
%41 = tail call double @llvm.fmuladd.f64(double %40, double %34, double %29)
%42 = fadd double %28, %34
br label %43
43: ; preds = %37, %26
%44 = phi double [ %41, %37 ], [ %29, %26 ]
%45 = phi double [ %42, %37 ], [ %28, %26 ]
%46 = sub i32 %25, %30
%47 = icmp slt i32 %46, %1
br i1 %47, label %48, label %54
48: ; preds = %43
%49 = sext i32 %46 to i64
%50 = getelementptr inbounds double, double* %8, i64 %49
%51 = load double, double* %50
%52 = tail call double @llvm.fmuladd.f64(double %51, double %34, double %44)
%53 = fadd double %34, %45
br label %54
54: ; preds = %48, %43
%55 = phi double [ %52, %48 ], [ %44, %43 ]
%56 = phi double [ %53, %48 ], [ %45, %43 ]
%57 = mul nsw i64 %27, %14
%58 = add nuw nsw i64 %57, %18
%59 = getelementptr inbounds double, double* %0, i64 %58
%60 = load double, double* %59
%61 = fcmp une double %60, 0.000000e+00
%62 = tail call double @llvm.fmuladd.f64(double %60, double %34, double %55)
%63 = select i1 %61, double %62, double %55
%64 = select i1 %61, double %34, double -0.000000e+00
%65 = fadd double %56, %64
%66 = icmp eq i64 %31, %16
br i1 %66, label %67, label %26
67: ; preds = %54, %17
%68 = phi double [ %20, %17 ], [ %63, %54 ]
%69 = phi double [ 1.000000e+00, %17 ], [ %65, %54 ]
%70 = fdiv double %68, %69
%71 = getelementptr inbounds double, double* %3, i64 %18
store double %70, double* %71
%72 = add nuw nsw i64 %18, 1
%73 = icmp eq i64 %72, %15
br i1 %73, label %74, label %17
74: ; preds = %67, %4
ret void
}