float diff_t(float floatpic[17][512][512], float kernel[5], int x, int y,
             int n) {
  int i;
  float sum;
  sum = 0.0;
  for (i = (-n); i <= n; i++) {
    sum += kernel[i + 2] * floatpic[2 + i][x][y];
  }
  return (sum);
}


x86 assembly (clang, -O0):
<diff_t>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
xor    %eax,%eax
sub    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     118c <diff_t+0x8c>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
shl    $0x14,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
shl    $0xb,%rcx
add    %rcx,%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1126 <diff_t+0x26>
movss  -0x24(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<diff_t>:
test   %r8d,%r8d
js     11a2 <diff_t+0xa2>
mov    %r8d,%eax
neg    %eax
movslq %edx,%r10
movslq %ecx,%r9
cltq
movss  0x8(%rsi,%rax,4),%xmm0
add    $0x2,%rax
shl    $0x14,%rax
add    %rdi,%rax
shl    $0xb,%r10
add    %r10,%rax
mulss  (%rax,%r9,4),%xmm0
xorps  %xmm1,%xmm1
addss  %xmm0,%xmm1
test   $0x7fffffff,%r8d
je     11a6 <diff_t+0xa6>
lea    (%r8,%r8,1),%eax
movslq %r8d,%rcx
mov    $0x4,%edx
sub    %rcx,%rdx
lea    (%rsi,%rdx,4),%rcx
shl    $0x14,%rdx
add    %rdx,%rdi
add    %r10,%rdi
lea    (%rdi,%r9,4),%rdx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movss  -0x4(%rcx,%rsi,4),%xmm2
mulss  -0x100000(%rdx),%xmm2
ff
movss  (%rcx,%rsi,4),%xmm0
addss  %xmm1,%xmm2
mulss  (%rdx),%xmm0
addss  %xmm2,%xmm0
add    $0x2,%rsi
add    $0x200000,%rdx
movaps %xmm0,%xmm1
cmp    %esi,%eax
jne    1170 <diff_t+0x70>
ret
xorps  %xmm0,%xmm0
ret
movaps %xmm1,%xmm0
ret

LLVM IR (-O0):
define float @diff_t([512 x [512 x float]]* %0, float* %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca [512 x [512 x float]]*
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca float
store [512 x [512 x float]]* %0, [512 x [512 x float]]** %6
store float* %1, float** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i32 %4, i32* %10
store float 0.000000e+00, float* %12
%13 = load i32, i32* %10
%14 = sub nsw i32 0, %13
store i32 %14, i32* %11
br label %15
15: ; preds = %40, %5
%16 = load i32, i32* %11
%17 = load i32, i32* %10
%18 = icmp sle i32 %16, %17
br i1 %18, label %19, label %43
19: ; preds = %15
%20 = load float*, float** %7
%21 = load i32, i32* %11
%22 = add nsw i32 %21, 2
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float, float* %20, i64 %23
%25 = load float, float* %24
%26 = load [512 x [512 x float]]*, [512 x [512 x float]]** %6
%27 = load i32, i32* %11
%28 = add nsw i32 2, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds [512 x [512 x float]], [512 x [512 x float]]* %26, i64 %29
%31 = load i32, i32* %8
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds [512 x [512 x float]], [512 x [512 x float]]* %30, i64 0, i64 %32
%34 = load i32, i32* %9
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds [512 x float], [512 x float]* %33, i64 0, i64 %35
%37 = load float, float* %36
%38 = load float, float* %12
%39 = call float @llvm.fmuladd.f32(float %25, float %37, float %38)
store float %39, float* %12
br label %40
40: ; preds = %19
%41 = load i32, i32* %11
%42 = add nsw i32 %41, 1
store i32 %42, i32* %11
br label %15
43: ; preds = %15
%44 = load float, float* %12
ret float %44
}

LLVM IR (-O2):
define float @diff_t([512 x [512 x float]]* nocapture %0, float* nocapture %1, i32 %2, i32 %3, i32 %4) {
%6 = icmp slt i32 %4, 0
br i1 %6, label %40, label %7
7: ; preds = %5
%8 = sub nsw i32 0, %4
%9 = sext i32 %2 to i64
%10 = sext i32 %3 to i64
%11 = sext i32 %8 to i64
%12 = add nuw i32 %4, 1
%13 = add nsw i64 %11, 2
%14 = getelementptr inbounds float, float* %1, i64 %13
%15 = load float, float* %14
%16 = getelementptr inbounds [512 x [512 x float]], [512 x [512 x float]]* %0, i64 %13, i64 %9, i64 %10
%17 = load float, float* %16
%18 = tail call float @llvm.fmuladd.f32(float %15, float %17, float 0.000000e+00)
%19 = add nsw i64 %11, 1
%20 = and i32 %4, 2147483647
%21 = icmp eq i32 %20, 0
br i1 %21, label %40, label %22
22: ; preds = %7, %22
%23 = phi i64 [ %37, %22 ], [ %19, %7 ]
%24 = phi float [ %36, %22 ], [ %18, %7 ]
%25 = add nsw i64 %23, 2
%26 = getelementptr inbounds float, float* %1, i64 %25
%27 = load float, float* %26
%28 = getelementptr inbounds [512 x [512 x float]], [512 x [512 x float]]* %0, i64 %25, i64 %9, i64 %10
%29 = load float, float* %28
%30 = tail call float @llvm.fmuladd.f32(float %27, float %29, float %24)
%31 = add nsw i64 %23, 3
%32 = getelementptr inbounds float, float* %1, i64 %31
%33 = load float, float* %32
%34 = getelementptr inbounds [512 x [512 x float]], [512 x [512 x float]]* %0, i64 %31, i64 %9, i64 %10
%35 = load float, float* %34
%36 = tail call float @llvm.fmuladd.f32(float %33, float %35, float %30)
%37 = add nsw i64 %23, 2
%38 = trunc i64 %37 to i32
%39 = icmp eq i32 %12, %38
br i1 %39, label %40, label %22
40: ; preds = %7, %22, %5
%41 = phi float [ 0.000000e+00, %5 ], [ %18, %7 ], [ %36, %22 ]
ret float %41
}