void laplacian2(image, lap, nx, ny) float *lap, *image;
int nx, ny;
{
  int index, i, j;
  for (i = 1; i < ny - 1; i++) {
    for (j = 1; j < nx - 1; j++) {
      index = (i * nx) + j;
      lap[index] = (((4.0 * image[index]) - image[index + 1] -
                     image[index + nx] - image[index - 1] - image[index - nx]) /
                    16.0);
    }
  }
  for (j = 0; j < nx; j++) {
    lap[(ny - 1) * nx + j] = 0.0;
    lap[j] = 0.0;
  }
  for (i = 0; i < ny; i++) {
    lap[(i * nx) + nx - 1] = 0.0;
    lap[i * nx] = 0.0;
  }
}


x86 assembly (clang, -O0):
<laplacian2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1229 <laplacian2+0x129>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1216 <laplacian2+0x116>
mov    -0x20(%rbp),%eax
imul   -0x14(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xe73(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
subsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
subsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  0xe0d(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1131 <laplacian2+0x31>
jmp    121b <laplacian2+0x11b>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1119 <laplacian2+0x19>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1276 <laplacian2+0x176>
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x14(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1230 <laplacian2+0x130>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    12c9 <laplacian2+0x1c9>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    127d <laplacian2+0x17d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<laplacian2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %rsi,-0x78(%rsp)
mov    %rdi,-0x58(%rsp)
lea    -0x1(%rcx),%eax
mov    %eax,-0x7c(%rsp)
movslq %edx,%rbp
mov    %rcx,-0x38(%rsp)
cmp    $0x3,%ecx
mov    %rbp,-0x70(%rsp)
jl     1476 <laplacian2+0x376>
lea    -0x1(%rbp),%r11d
mov    -0x7c(%rsp),%eax
mov    %rax,0x10(%rsp)
lea    0x1(%rbp),%r8
lea    (%r11,%rbp,1),%rax
mov    %rax,-0x8(%rsp)
lea    0x1(%rbp,%r11,1),%rax
mov    %rax,-0x10(%rsp)
lea    0x0(%rbp,%rbp,1),%rax
mov    %rax,-0x40(%rsp)
lea    0x1(%rbp,%rbp,1),%rax
mov    %rax,-0x18(%rsp)
lea    (%r11,%rbp,2),%rax
mov    %rax,-0x20(%rsp)
lea    -0x1(%r11),%rax
mov    %rax,-0x50(%rsp)
and    $0xfffffffffffffffc,%rax
mov    %rax,-0x28(%rsp)
add    $0x1,%rax
mov    %rax,-0x30(%rsp)
mov    -0x78(%rsp),%rax
lea    0x4(%rax,%rbp,4),%rcx
lea    0x0(,%rbp,4),%rax
mov    %rax,0x8(%rsp)
mov    -0x58(%rsp),%rdi
lea    0x8(%rdi,%rbp,4),%r15
lea    0x4(%rdi,%rbp,8),%r12
movsd  0xe61(%rip),%xmm0
movsd  0xe61(%rip),%xmm1
movapd 0xe31(%rip),%xmm2
movapd 0xe39(%rip),%xmm3
mov    $0x1,%r13d
xor    %eax,%eax
mov    %rax,-0x48(%rsp)
mov    %r8,(%rsp)
mov    %rbp,%rax
mov    %rdx,0x18(%rsp)
mov    %r11,-0x68(%rsp)
jmp    123b <laplacian2+0x13b>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
addq   $0x1,-0x48(%rsp)
mov    0x8(%rsp),%rdx
add    %rdx,%rcx
add    %rdx,%rdi
add    %rdx,%r15
add    %rdx,%r12
mov    -0x70(%rsp),%rbp
add    %rbp,%rax
add    %rbp,%r8
add    %rbp,-0x40(%rsp)
cmp    0x10(%rsp),%r13
mov    0x18(%rsp),%rdx
je     1476 <laplacian2+0x376>
cmp    $0x3,%edx
jl     1200 <laplacian2+0x100>
mov    $0x1,%r10d
cmpq   $0x4,-0x50(%rsp)
jb     13d0 <laplacian2+0x2d0>
mov    -0x48(%rsp),%rbp
mov    -0x70(%rsp),%rdx
imul   %rdx,%rbp
mov    (%rsp),%rsi
lea    (%rsi,%rbp,1),%rbx
mov    -0x8(%rsp),%rsi
add    %rbp,%rsi
lea    0x0(%rbp,%r11,1),%r14
add    %rbp,%rdx
mov    %rdx,-0x60(%rsp)
mov    -0x78(%rsp),%rdx
lea    (%rdx,%rbx,4),%rbx
lea    (%rdx,%rsi,4),%rsi
mov    -0x58(%rsp),%r11
lea    (%r11,%rbp,4),%r9
add    $0x4,%r9
lea    (%r11,%r14,4),%rdx
cmp    %rdx,%rbx
setb   -0x7d(%rsp)
cmp    %rsi,%r9
mov    -0x10(%rsp),%rdx
lea    (%rdx,%rbp,1),%r14
mov    -0x60(%rsp),%rdx
lea    (%r11,%rdx,4),%r9
lea    (%r11,%r14,4),%rdx
setb   -0x60(%rsp)
cmp    %rdx,%rbx
setb   %r14b
cmp    %rsi,%r9
mov    -0x20(%rsp),%rdx
lea    (%rdx,%rbp,1),%rdx
lea    (%r11,%rdx,4),%rdx
setb   %r9b
cmp    %rdx,%rbx
mov    -0x18(%rsp),%rdx
lea    (%rdx,%rbp,1),%rdx
lea    (%r11,%rdx,4),%rdx
setb   %bl
cmp    %rsi,%rdx
setb   %sil
mov    -0x60(%rsp),%dl
test   %dl,-0x7d(%rsp)
jne    13be <laplacian2+0x2be>
and    %r9b,%r14b
jne    13be <laplacian2+0x2be>
and    %sil,%bl
mov    -0x68(%rsp),%r11
jne    13d0 <laplacian2+0x2d0>
xor    %edx,%edx
mov    -0x28(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cvtps2pd 0x4(%r15,%rdx,4),%xmm5
cvtps2pd -0x4(%r15,%rdx,4),%xmm4
cvtps2pd 0x8(%r15,%rdx,4),%xmm6
cvtps2pd (%r15,%rdx,4),%xmm7
mulpd  %xmm2,%xmm4
subpd  %xmm7,%xmm4
mulpd  %xmm2,%xmm5
subpd  %xmm6,%xmm5
cvtps2pd (%r12,%rdx,4),%xmm6
subpd  %xmm6,%xmm4
cvtps2pd 0x8(%r12,%rdx,4),%xmm6
subpd  %xmm6,%xmm5
cvtps2pd -0x8(%r15,%rdx,4),%xmm6
subpd  %xmm6,%xmm4
subpd  %xmm7,%xmm5
cvtps2pd 0x4(%rdi,%rdx,4),%xmm6
cvtps2pd 0xc(%rdi,%rdx,4),%xmm7
subpd  %xmm6,%xmm4
subpd  %xmm7,%xmm5
mulpd  %xmm3,%xmm4
mulpd  %xmm3,%xmm5
cvtpd2ps %xmm5,%xmm5
cvtpd2ps %xmm4,%xmm4
unpcklpd %xmm5,%xmm4
movupd %xmm4,(%rcx,%rdx,4)
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1330 <laplacian2+0x230>
mov    -0x30(%rsp),%r10
cmp    %rsi,-0x50(%rsp)
je     1200 <laplacian2+0x100>
jmp    13d0 <laplacian2+0x2d0>
mov    -0x68(%rsp),%r11
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r11,%rdx
sub    %r10,%rdx
lea    (%rdi,%r10,4),%r9
lea    (%r10,%rax,1),%rsi
mov    -0x78(%rsp),%rbx
lea    (%rbx,%rsi,4),%rsi
lea    (%r10,%r8,1),%rbx
mov    -0x58(%rsp),%rbp
lea    0x0(,%rbx,4),%r14
add    %rbp,%r14
add    -0x40(%rsp),%r10
lea    0x0(%rbp,%r10,4),%rbp
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movss  -0x8(%r14,%rbx,4),%xmm4
movss  -0x4(%r14,%rbx,4),%xmm5
cvtss2sd %xmm5,%xmm5
movss  (%r14,%rbx,4),%xmm6
cvtss2sd %xmm6,%xmm6
mulsd  %xmm0,%xmm5
subsd  %xmm6,%xmm5
movss  0x0(%rbp,%rbx,4),%xmm6
cvtss2sd %xmm6,%xmm6
subsd  %xmm6,%xmm5
cvtss2sd %xmm4,%xmm4
subsd  %xmm4,%xmm5
movss  (%r9,%rbx,4),%xmm4
cvtss2sd %xmm4,%xmm4
subsd  %xmm4,%xmm5
mulsd  %xmm1,%xmm5
xorps  %xmm4,%xmm4
cvtsd2ss %xmm5,%xmm4
movss  %xmm4,(%rsi,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%rdx
jne    1410 <laplacian2+0x310>
jmp    1200 <laplacian2+0x100>
test   %edx,%edx
mov    -0x78(%rsp),%r15
jle    164b <laplacian2+0x54b>
mov    -0x7c(%rsp),%eax
imul   %edx,%eax
movslq %eax,%r8
mov    %edx,%eax
cmp    $0x8,%edx
jb     14b6 <laplacian2+0x3b6>
lea    (%r15,%r8,4),%rcx
lea    (%r15,%rax,4),%rsi
cmp    %rsi,%rcx
jae    154c <laplacian2+0x44c>
lea    (%r8,%rax,1),%rcx
lea    (%r15,%rcx,4),%rcx
cmp    %r15,%rcx
jbe    154c <laplacian2+0x44c>
xor    %edi,%edi
mov    %rdi,%rbx
not    %rbx
add    %rax,%rbx
mov    %rax,%rsi
and    $0x3,%rsi
je     14e9 <laplacian2+0x3e9>
lea    (%r15,%r8,4),%rcx
xchg   %ax,%ax
movl   $0x0,(%rcx,%rdi,4)
movl   $0x0,(%r15,%rdi,4)
add    $0x1,%rdi
add    $0xffffffffffffffff,%rsi
jne    14d0 <laplacian2+0x3d0>
cmp    $0x3,%rbx
jb     164b <laplacian2+0x54b>
lea    0xc(%r15),%rsi
lea    (%r15,%r8,4),%rcx
add    $0xc,%rcx
nop
movl   $0x0,-0xc(%rcx,%rdi,4)
movl   $0x0,-0xc(%rsi,%rdi,4)
movl   $0x0,-0x8(%rcx,%rdi,4)
movl   $0x0,-0x8(%rsi,%rdi,4)
movl   $0x0,-0x4(%rcx,%rdi,4)
movl   $0x0,-0x4(%rsi,%rdi,4)
movl   $0x0,(%rcx,%rdi,4)
movl   $0x0,(%rsi,%rdi,4)
add    $0x4,%rdi
cmp    %rdi,%rax
jne    1500 <laplacian2+0x400>
jmp    164b <laplacian2+0x54b>
mov    %rdx,%r9
mov    %eax,%edi
and    $0xfffffff8,%edi
lea    -0x8(%rdi),%rsi
mov    %rsi,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x18,%rsi
jae    1575 <laplacian2+0x475>
xor    %esi,%esi
jmp    1600 <laplacian2+0x500>
and    $0xfffffffffffffffc,%rdx
lea    (%r15,%r8,4),%rbx
add    $0x70,%rbx
xor    %esi,%esi
xorpd  %xmm0,%xmm0
nopw   0x0(%rax,%rax,1)
movupd %xmm0,-0x70(%rbx,%rsi,4)
movupd %xmm0,-0x60(%rbx,%rsi,4)
movupd %xmm0,(%r15,%rsi,4)
movupd %xmm0,0x10(%r15,%rsi,4)
movupd %xmm0,-0x50(%rbx,%rsi,4)
movupd %xmm0,-0x40(%rbx,%rsi,4)
movupd %xmm0,0x20(%r15,%rsi,4)
movupd %xmm0,0x30(%r15,%rsi,4)
movupd %xmm0,-0x30(%rbx,%rsi,4)
movupd %xmm0,-0x20(%rbx,%rsi,4)
movupd %xmm0,0x40(%r15,%rsi,4)
movupd %xmm0,0x50(%r15,%rsi,4)
movupd %xmm0,-0x10(%rbx,%rsi,4)
movupd %xmm0,(%rbx,%rsi,4)
movupd %xmm0,0x60(%r15,%rsi,4)
movupd %xmm0,0x70(%r15,%rsi,4)
add    $0x20,%rsi
add    $0xfffffffffffffffc,%rdx
jne    1590 <laplacian2+0x490>
test   %rcx,%rcx
je     163f <laplacian2+0x53f>
lea    (%r15,%rsi,4),%rdx
add    $0x10,%rdx
lea    (%rdx,%r8,4),%rsi
shl    $0x5,%rcx
xor    %ebx,%ebx
xorpd  %xmm0,%xmm0
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%rsi,%rbx,1)
movupd %xmm0,(%rsi,%rbx,1)
movupd %xmm0,-0x10(%rdx,%rbx,1)
movupd %xmm0,(%rdx,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%rcx
jne    1620 <laplacian2+0x520>
cmp    %rax,%rdi
mov    %r9,%rdx
jne    14b8 <laplacian2+0x3b8>
mov    -0x38(%rsp),%rcx
test   %ecx,%ecx
jle    1803 <laplacian2+0x703>
lea    -0x1(%rdx),%r8d
mov    %ecx,%eax
xor    %r10d,%r10d
cmp    $0x8,%ecx
jb     1760 <laplacian2+0x660>
cmp    $0x1,%edx
jne    1760 <laplacian2+0x660>
mov    %eax,%r10d
and    $0xfffffff8,%r10d
lea    -0x8(%r10),%rcx
mov    %rcx,%r9
shr    $0x3,%r9
add    $0x1,%r9
test   %rcx,%rcx
je     171b <laplacian2+0x61b>
movabs $0x3ffffffffffffffc,%rdi
ff ff 3f
or     $0x2,%rdi
and    %r9,%rdi
mov    -0x70(%rsp),%r11
mov    %r11,%r14
shl    $0x5,%r14
lea    0x10(%r15),%rbp
shl    $0x6,%r11
lea    (%rdx,%rdx,8),%ecx
add    $0xffffffff,%ecx
mov    %rdx,%r12
shl    $0x4,%edx
xor    %ebx,%ebx
xorpd  %xmm0,%xmm0
mov    %r8d,%esi
nopl   0x0(%rax)
movslq %esi,%rsi
movupd %xmm0,(%r15,%rsi,4)
movupd %xmm0,0x10(%r15,%rsi,4)
movupd %xmm0,-0x10(%rbp)
movupd %xmm0,0x0(%rbp)
movslq %ecx,%rcx
movupd %xmm0,(%r15,%rcx,4)
movupd %xmm0,0x10(%r15,%rcx,4)
movupd %xmm0,-0x10(%rbp,%r14,1)
movupd %xmm0,0x0(%rbp,%r14,1)
add    $0x10,%rbx
add    %r11,%rbp
add    %edx,%ecx
add    %edx,%esi
add    $0xfffffffffffffffe,%rdi
jne    16d0 <laplacian2+0x5d0>
jmp    1720 <laplacian2+0x620>
mov    %rdx,%r12
xor    %ebx,%ebx
test   $0x1,%r9b
mov    -0x70(%rsp),%rbp
je     1754 <laplacian2+0x654>
imul   %rbp,%rbx
lea    (%r8,%rbx,1),%ecx
movslq %ecx,%rcx
xorpd  %xmm0,%xmm0
movupd %xmm0,(%r15,%rcx,4)
movupd %xmm0,0x10(%r15,%rcx,4)
movupd %xmm0,(%r15,%rbx,4)
movupd %xmm0,0x10(%r15,%rbx,4)
cmp    %rax,%r10
mov    %r12,%rdx
je     1803 <laplacian2+0x703>
mov    %r10,%rcx
not    %rcx
test   $0x1,%al
je     178b <laplacian2+0x68b>
mov    %r10,%rdi
imul   %rbp,%rdi
add    %edi,%r8d
movslq %r8d,%rsi
movl   $0x0,(%r15,%rsi,4)
movl   $0x0,(%r15,%rdi,4)
or     $0x1,%r10
add    %rax,%rcx
je     1803 <laplacian2+0x703>
sub    %r10,%rax
lea    0x1(%r10),%edi
imul   %edx,%edi
add    $0xffffffff,%edi
lea    (%rdx,%rdx,1),%ecx
mov    %r10,%r8
imul   %rbp,%r8
mov    %r10,%r9
add    $0x1,%r9
imul   %rbp,%r9
shl    $0x3,%rbp
lea    0x2(%r10),%esi
imul   %edx,%esi
add    $0xffffffff,%esi
mov    %r15,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movslq %edi,%rdi
movl   $0x0,(%r15,%rdi,4)
movl   $0x0,(%rbx,%r8,4)
movslq %esi,%rsi
movl   $0x0,(%r15,%rsi,4)
movl   $0x0,(%rbx,%r9,4)
add    %ecx,%edi
add    %rbp,%rbx
add    %ecx,%esi
add    $0xfffffffffffffffe,%rax
jne    17d0 <laplacian2+0x6d0>
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @laplacian2(float* nocapture %0, float* nocapture writeonly %1, i32 %2, i32 %3) {
%5 = add i32 %3, -1
%6 = icmp sgt i32 %3, 2
br i1 %6, label %7, label %100
7: ; preds = %4
%8 = add i32 %2, -1
%9 = icmp sgt i32 %2, 2
%10 = sext i32 %2 to i64
%11 = zext i32 %5 to i64
%12 = zext i32 %8 to i64
%13 = add nsw i64 %10, 1
%14 = add nsw i64 %10, %12
%15 = add nsw i64 %14, 1
%16 = shl nsw i64 %10, 1
%17 = or i64 %16, 1
%18 = add nsw i64 %16, %12
%19 = add nsw i64 %12, -1
%20 = icmp ult i64 %19, 4
%21 = and i64 %19, -4
%22 = or i64 %21, 1
%23 = icmp eq i64 %19, %21
br label %24
24: ; preds = %7, %239
%25 = phi i64 [ 0, %7 ], [ %242, %239 ]
%26 = phi i64 [ 1, %7 ], [ %240, %239 ]
%27 = mul i64 %25, %10
%28 = add i64 %13, %27
%29 = getelementptr float, float* %1, i64 %28
%30 = add i64 %14, %27
%31 = getelementptr float, float* %1, i64 %30
%32 = add i64 %27, 1
%33 = getelementptr float, float* %0, i64 %32
%34 = add i64 %27, %12
%35 = getelementptr float, float* %0, i64 %34
%36 = add i64 %27, %10
%37 = getelementptr float, float* %0, i64 %36
%38 = add i64 %15, %27
%39 = getelementptr float, float* %0, i64 %38
%40 = add i64 %17, %27
%41 = getelementptr float, float* %0, i64 %40
%42 = add i64 %18, %27
%43 = getelementptr float, float* %0, i64 %42
br i1 %9, label %44, label %239
44: ; preds = %24
%45 = mul nsw i64 %26, %10
br i1 %20, label %98, label %46
46: ; preds = %44
%47 = icmp ult float* %29, %35
%48 = icmp ult float* %33, %31
%49 = and i1 %47, %48
%50 = icmp ult float* %29, %39
%51 = icmp ult float* %37, %31
%52 = and i1 %50, %51
%53 = or i1 %49, %52
%54 = icmp ult float* %29, %43
%55 = icmp ult float* %41, %31
%56 = and i1 %54, %55
%57 = or i1 %53, %56
br i1 %57, label %98, label %58
58: ; preds = %46, %58
%59 = phi i64 [ %95, %58 ], [ 0, %46 ]
%60 = or i64 %59, 1
%61 = add nsw i64 %60, %45
%62 = getelementptr inbounds float, float* %0, i64 %61
%63 = bitcast float* %62 to <4 x float>*
%64 = load <4 x float>, <4 x float>* %63
%65 = fpext <4 x float> %64 to <4 x double>
%66 = add nsw i64 %61, 1
%67 = getelementptr inbounds float, float* %0, i64 %66
%68 = bitcast float* %67 to <4 x float>*
%69 = load <4 x float>, <4 x float>* %68
%70 = fpext <4 x float> %69 to <4 x double>
%71 = fneg <4 x double> %70
%72 = call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %65, <4 x double> <double 4.000000e+00, double 4.000000e+00, double 4.000000e+00, double 4.000000e+00>, <4 x double> %71)
%73 = add nsw i64 %61, %10
%74 = getelementptr inbounds float, float* %0, i64 %73
%75 = bitcast float* %74 to <4 x float>*
%76 = load <4 x float>, <4 x float>* %75
%77 = fpext <4 x float> %76 to <4 x double>
%78 = fsub <4 x double> %72, %77
%79 = add nsw i64 %61, -1
%80 = getelementptr inbounds float, float* %0, i64 %79
%81 = bitcast float* %80 to <4 x float>*
%82 = load <4 x float>, <4 x float>* %81
%83 = fpext <4 x float> %82 to <4 x double>
%84 = fsub <4 x double> %78, %83
%85 = sub nsw i64 %61, %10
%86 = getelementptr inbounds float, float* %0, i64 %85
%87 = bitcast float* %86 to <4 x float>*
%88 = load <4 x float>, <4 x float>* %87
%89 = fpext <4 x float> %88 to <4 x double>
%90 = fsub <4 x double> %84, %89
%91 = fmul <4 x double> %90, <double 6.250000e-02, double 6.250000e-02, double 6.250000e-02, double 6.250000e-02>
%92 = fptrunc <4 x double> %91 to <4 x float>
%93 = getelementptr inbounds float, float* %1, i64 %61
%94 = bitcast float* %93 to <4 x float>*
store <4 x float> %92, <4 x float>* %94
%95 = add nuw i64 %59, 4
%96 = icmp eq i64 %95, %21
br i1 %96, label %97, label %58
97: ; preds = %58
br i1 %23, label %239, label %98
98: ; preds = %46, %44, %97
%99 = phi i64 [ 1, %46 ], [ 1, %44 ], [ %22, %97 ]
br label %207
100: ; preds = %239, %4
%101 = icmp sgt i32 %2, 0
br i1 %101, label %102, label %243
102: ; preds = %100
%103 = mul nsw i32 %5, %2
%104 = sext i32 %103 to i64
%105 = zext i32 %2 to i64
%106 = icmp ult i32 %2, 8
br i1 %106, label %189, label %107
107: ; preds = %102
%108 = getelementptr float, float* %1, i64 %104
%109 = add nsw i64 %104, %105
%110 = getelementptr float, float* %1, i64 %109
%111 = getelementptr float, float* %1, i64 %105
%112 = icmp ult float* %108, %111
%113 = icmp ugt float* %110, %1
%114 = and i1 %112, %113
br i1 %114, label %189, label %115
115: ; preds = %107
%116 = and i64 %105, 4294967288
%117 = add nsw i64 %116, -8
%118 = lshr exact i64 %117, 3
%119 = add nuw nsw i64 %118, 1
%120 = and i64 %119, 3
%121 = icmp ult i64 %117, 24
br i1 %121, label %169, label %122
122: ; preds = %115
%123 = and i64 %119, 4611686018427387900
br label %124
124: ; preds = %124, %122
%125 = phi i64 [ 0, %122 ], [ %166, %124 ]
%126 = phi i64 [ 0, %122 ], [ %167, %124 ]
%127 = add nsw i64 %125, %104
%128 = getelementptr inbounds float, float* %1, i64 %127
%129 = bitcast float* %128 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %129
%130 = getelementptr inbounds float, float* %128, i64 4
%131 = bitcast float* %130 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %131
%132 = getelementptr inbounds float, float* %1, i64 %125
%133 = bitcast float* %132 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %133
%134 = getelementptr inbounds float, float* %132, i64 4
%135 = bitcast float* %134 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %135
%136 = or i64 %125, 8
%137 = add nsw i64 %136, %104
%138 = getelementptr inbounds float, float* %1, i64 %137
%139 = bitcast float* %138 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %139
%140 = getelementptr inbounds float, float* %138, i64 4
%141 = bitcast float* %140 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %141
%142 = getelementptr inbounds float, float* %1, i64 %136
%143 = bitcast float* %142 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %143
%144 = getelementptr inbounds float, float* %142, i64 4
%145 = bitcast float* %144 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %145
%146 = or i64 %125, 16
%147 = add nsw i64 %146, %104
%148 = getelementptr inbounds float, float* %1, i64 %147
%149 = bitcast float* %148 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %149
%150 = getelementptr inbounds float, float* %148, i64 4
%151 = bitcast float* %150 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %151
%152 = getelementptr inbounds float, float* %1, i64 %146
%153 = bitcast float* %152 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %153
%154 = getelementptr inbounds float, float* %152, i64 4
%155 = bitcast float* %154 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %155
%156 = or i64 %125, 24
%157 = add nsw i64 %156, %104
%158 = getelementptr inbounds float, float* %1, i64 %157
%159 = bitcast float* %158 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %159
%160 = getelementptr inbounds float, float* %158, i64 4
%161 = bitcast float* %160 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %161
%162 = getelementptr inbounds float, float* %1, i64 %156
%163 = bitcast float* %162 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %163
%164 = getelementptr inbounds float, float* %162, i64 4
%165 = bitcast float* %164 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %165
%166 = add nuw i64 %125, 32
%167 = add i64 %126, 4
%168 = icmp eq i64 %167, %123
br i1 %168, label %169, label %124
169: ; preds = %124, %115
%170 = phi i64 [ 0, %115 ], [ %166, %124 ]
%171 = icmp eq i64 %120, 0
br i1 %171, label %187, label %172
172: ; preds = %169, %172
%173 = phi i64 [ %184, %172 ], [ %170, %169 ]
%174 = phi i64 [ %185, %172 ], [ 0, %169 ]
%175 = add nsw i64 %173, %104
%176 = getelementptr inbounds float, float* %1, i64 %175
%177 = bitcast float* %176 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %177
%178 = getelementptr inbounds float, float* %176, i64 4
%179 = bitcast float* %178 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %179
%180 = getelementptr inbounds float, float* %1, i64 %173
%181 = bitcast float* %180 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %181
%182 = getelementptr inbounds float, float* %180, i64 4
%183 = bitcast float* %182 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %183
%184 = add nuw i64 %173, 8
%185 = add i64 %174, 1
%186 = icmp eq i64 %185, %120
br i1 %186, label %187, label %172
187: ; preds = %172, %169
%188 = icmp eq i64 %116, %105
br i1 %188, label %243, label %189
189: ; preds = %107, %102, %187
%190 = phi i64 [ 0, %107 ], [ 0, %102 ], [ %116, %187 ]
%191 = xor i64 %190, -1
%192 = add nsw i64 %191, %105
%193 = and i64 %105, 3
%194 = icmp eq i64 %193, 0
br i1 %194, label %204, label %195
195: ; preds = %189, %195
%196 = phi i64 [ %201, %195 ], [ %190, %189 ]
%197 = phi i64 [ %202, %195 ], [ 0, %189 ]
%198 = add nsw i64 %196, %104
%199 = getelementptr inbounds float, float* %1, i64 %198
store float 0.000000e+00, float* %199
%200 = getelementptr inbounds float, float* %1, i64 %196
store float 0.000000e+00, float* %200
%201 = add nuw nsw i64 %196, 1
%202 = add i64 %197, 1
%203 = icmp eq i64 %202, %193
br i1 %203, label %204, label %195
204: ; preds = %195, %189
%205 = phi i64 [ %190, %189 ], [ %201, %195 ]
%206 = icmp ult i64 %192, 3
br i1 %206, label %243, label %327
207: ; preds = %98, %207
%208 = phi i64 [ %237, %207 ], [ %99, %98 ]
%209 = add nsw i64 %208, %45
%210 = getelementptr inbounds float, float* %0, i64 %209
%211 = load float, float* %210
%212 = fpext float %211 to double
%213 = add nsw i64 %209, 1
%214 = getelementptr inbounds float, float* %0, i64 %213
%215 = load float, float* %214
%216 = fpext float %215 to double
%217 = fneg double %216
%218 = tail call double @llvm.fmuladd.f64(double %212, double 4.000000e+00, double %217)
%219 = add nsw i64 %209, %10
%220 = getelementptr inbounds float, float* %0, i64 %219
%221 = load float, float* %220
%222 = fpext float %221 to double
%223 = fsub double %218, %222
%224 = add nsw i64 %209, -1
%225 = getelementptr inbounds float, float* %0, i64 %224
%226 = load float, float* %225
%227 = fpext float %226 to double
%228 = fsub double %223, %227
%229 = sub nsw i64 %209, %10
%230 = getelementptr inbounds float, float* %0, i64 %229
%231 = load float, float* %230
%232 = fpext float %231 to double
%233 = fsub double %228, %232
%234 = fmul double %233, 6.250000e-02
%235 = fptrunc double %234 to float
%236 = getelementptr inbounds float, float* %1, i64 %209
store float %235, float* %236
%237 = add nuw nsw i64 %208, 1
%238 = icmp eq i64 %237, %12
br i1 %238, label %239, label %207
239: ; preds = %207, %97, %24
%240 = add nuw nsw i64 %26, 1
%241 = icmp eq i64 %240, %11
%242 = add i64 %25, 1
br i1 %241, label %100, label %24
243: ; preds = %204, %327, %187, %100
%244 = icmp sgt i32 %3, 0
br i1 %244, label %245, label %363
245: ; preds = %243
%246 = add i32 %2, -1
%247 = sext i32 %2 to i64
%248 = zext i32 %3 to i64
%249 = icmp ugt i32 %3, 7
%250 = icmp eq i32 %2, 1
%251 = and i1 %249, %250
br i1 %251, label %252, label %310
252: ; preds = %245
%253 = and i64 %248, 4294967288
%254 = add nsw i64 %253, -8
%255 = lshr exact i64 %254, 3
%256 = add nuw nsw i64 %255, 1
%257 = and i64 %256, 1
%258 = icmp eq i64 %254, 0
br i1 %258, label %292, label %259
259: ; preds = %252
%260 = and i64 %256, 4611686018427387902
br label %261
261: ; preds = %261, %259
%262 = phi i64 [ 0, %259 ], [ %289, %261 ]
%263 = phi i64 [ 0, %259 ], [ %290, %261 ]
%264 = mul nsw i64 %262, %247
%265 = trunc i64 %264 to i32
%266 = add i32 %246, %265
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds float, float* %1, i64 %267
%269 = bitcast float* %268 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %269
%270 = getelementptr inbounds float, float* %268, i64 4
%271 = bitcast float* %270 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %271
%272 = getelementptr inbounds float, float* %1, i64 %264
%273 = bitcast float* %272 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %273
%274 = getelementptr inbounds float, float* %272, i64 4
%275 = bitcast float* %274 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %275
%276 = or i64 %262, 8
%277 = mul nsw i64 %276, %247
%278 = trunc i64 %277 to i32
%279 = add i32 %246, %278
%280 = sext i32 %279 to i64
%281 = getelementptr inbounds float, float* %1, i64 %280
%282 = bitcast float* %281 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %282
%283 = getelementptr inbounds float, float* %281, i64 4
%284 = bitcast float* %283 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %284
%285 = getelementptr inbounds float, float* %1, i64 %277
%286 = bitcast float* %285 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %286
%287 = getelementptr inbounds float, float* %285, i64 4
%288 = bitcast float* %287 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %288
%289 = add nuw i64 %262, 16
%290 = add i64 %263, 2
%291 = icmp eq i64 %290, %260
br i1 %291, label %292, label %261
292: ; preds = %261, %252
%293 = phi i64 [ 0, %252 ], [ %289, %261 ]
%294 = icmp eq i64 %257, 0
br i1 %294, label %308, label %295
295: ; preds = %292
%296 = mul nsw i64 %293, %247
%297 = trunc i64 %296 to i32
%298 = add i32 %246, %297
%299 = sext i32 %298 to i64
%300 = getelementptr inbounds float, float* %1, i64 %299
%301 = bitcast float* %300 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %301
%302 = getelementptr inbounds float, float* %300, i64 4
%303 = bitcast float* %302 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %303
%304 = getelementptr inbounds float, float* %1, i64 %296
%305 = bitcast float* %304 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %305
%306 = getelementptr inbounds float, float* %304, i64 4
%307 = bitcast float* %306 to <4 x float>*
store <4 x float> zeroinitializer, <4 x float>* %307
br label %308
308: ; preds = %292, %295
%309 = icmp eq i64 %253, %248
br i1 %309, label %363, label %310
310: ; preds = %245, %308
%311 = phi i64 [ 0, %245 ], [ %253, %308 ]
%312 = xor i64 %311, -1
%313 = and i64 %248, 1
%314 = icmp eq i64 %313, 0
br i1 %314, label %323, label %315
315: ; preds = %310
%316 = mul nsw i64 %311, %247
%317 = trunc i64 %316 to i32
%318 = add i32 %246, %317
%319 = sext i32 %318 to i64
%320 = getelementptr inbounds float, float* %1, i64 %319
store float 0.000000e+00, float* %320
%321 = getelementptr inbounds float, float* %1, i64 %316
store float 0.000000e+00, float* %321
%322 = or i64 %311, 1
br label %323
323: ; preds = %315, %310
%324 = phi i64 [ %311, %310 ], [ %322, %315 ]
%325 = sub nsw i64 0, %248
%326 = icmp eq i64 %312, %325
br i1 %326, label %363, label %346
327: ; preds = %204, %327
%328 = phi i64 [ %344, %327 ], [ %205, %204 ]
%329 = add nsw i64 %328, %104
%330 = getelementptr inbounds float, float* %1, i64 %329
store float 0.000000e+00, float* %330
%331 = getelementptr inbounds float, float* %1, i64 %328
store float 0.000000e+00, float* %331
%332 = add nuw nsw i64 %328, 1
%333 = add nsw i64 %332, %104
%334 = getelementptr inbounds float, float* %1, i64 %333
store float 0.000000e+00, float* %334
%335 = getelementptr inbounds float, float* %1, i64 %332
store float 0.000000e+00, float* %335
%336 = add nuw nsw i64 %328, 2
%337 = add nsw i64 %336, %104
%338 = getelementptr inbounds float, float* %1, i64 %337
store float 0.000000e+00, float* %338
%339 = getelementptr inbounds float, float* %1, i64 %336
store float 0.000000e+00, float* %339
%340 = add nuw nsw i64 %328, 3
%341 = add nsw i64 %340, %104
%342 = getelementptr inbounds float, float* %1, i64 %341
store float 0.000000e+00, float* %342
%343 = getelementptr inbounds float, float* %1, i64 %340
store float 0.000000e+00, float* %343
%344 = add nuw nsw i64 %328, 4
%345 = icmp eq i64 %344, %105
br i1 %345, label %243, label %327
346: ; preds = %323, %346
%347 = phi i64 [ %361, %346 ], [ %324, %323 ]
%348 = mul nsw i64 %347, %247
%349 = trunc i64 %348 to i32
%350 = add i32 %246, %349
%351 = sext i32 %350 to i64
%352 = getelementptr inbounds float, float* %1, i64 %351
store float 0.000000e+00, float* %352
%353 = getelementptr inbounds float, float* %1, i64 %348
store float 0.000000e+00, float* %353
%354 = add nuw nsw i64 %347, 1
%355 = mul nsw i64 %354, %247
%356 = trunc i64 %355 to i32
%357 = add i32 %246, %356
%358 = sext i32 %357 to i64
%359 = getelementptr inbounds float, float* %1, i64 %358
store float 0.000000e+00, float* %359
%360 = getelementptr inbounds float, float* %1, i64 %355
store float 0.000000e+00, float* %360
%361 = add nuw nsw i64 %347, 2
%362 = icmp eq i64 %361, %248
br i1 %362, label %363, label %346
363: ; preds = %323, %346, %308, %243
ret void
}