void correct_velocities(double **U, double **V, double **F, double **G,
                        double **P, int nx, int ny, double dt, double dx,
                        double dy) {
  int i, j;
  for (i = 1; i < nx; i++)
    for (j = 1; j <= ny; j++)
      U[i][j] = F[i][j] - dt / dx * (P[i + 1][j] - P[i][j]);
  for (i = 1; i <= nx; i++)
    for (j = 1; j < ny; j++)
      V[i][j] = G[i][j] - dt / dy * (P[i][j + 1] - P[i][j]);
}


x86 assembly (clang, -O0):
<correct_velocities>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movsd  %xmm0,-0x38(%rbp)
movsd  %xmm1,-0x40(%rbp)
movsd  %xmm2,-0x48(%rbp)
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    11f7 <correct_velocities+0xf7>
movl   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x10(%rbp),%eax
jg     11e4 <correct_velocities+0xe4>
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1148 <correct_velocities+0x48>
jmp    11e9 <correct_velocities+0xe9>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1135 <correct_velocities+0x35>
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     12c0 <correct_velocities+0x1c0>
movl   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    12ad <correct_velocities+0x1ad>
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1211 <correct_velocities+0x111>
jmp    12b2 <correct_velocities+0x1b2>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11fe <correct_velocities+0xfe>
pop    %rbp
ret

x86 assembly (clang, -O2):
<correct_velocities>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x8(%rsp)
mov    %rsi,-0x10(%rsp)
mov    0x38(%rsp),%r11d
cmp    $0x2,%r9d
mov    %r8,-0x50(%rsp)
jl     1387 <correct_velocities+0x287>
movapd 0xed0(%rip),%xmm3
xorpd  %xmm0,%xmm3
divsd  %xmm1,%xmm3
lea    0x1(%r11),%eax
mov    %r9d,%r14d
lea    -0x1(%rax),%rbx
mov    %rbx,-0x48(%rsp)
and    $0xfffffffffffffffc,%rbx
lea    0x1(%rbx),%rcx
mov    %rcx,-0x18(%rsp)
movapd %xmm3,%xmm8
unpcklpd %xmm3,%xmm8
mov    %rax,%rcx
neg    %rcx
mov    %rcx,-0x28(%rsp)
mov    $0x1,%ecx
mov    %r9,-0x58(%rsp)
mov    %rdx,-0x30(%rsp)
mov    %rdi,-0x38(%rsp)
mov    %r14,-0x40(%rsp)
jmp    11a3 <correct_velocities+0xa3>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rcx
mov    %rcx,%r13
mov    %r13,%rcx
cmp    %r14,%r13
je     1387 <correct_velocities+0x287>
test   %r11d,%r11d
jle    1190 <correct_velocities+0x90>
lea    0x1(%rcx),%r13
mov    (%rdx,%rcx,8),%rbp
mov    (%r8,%rcx,8),%r15
mov    0x8(%r8,%rcx,8),%r10
mov    (%rdi,%rcx,8),%rcx
mov    $0x1,%r12d
cmpq   $0x4,-0x48(%rsp)
jb     1300 <correct_velocities+0x200>
mov    %r13,-0x60(%rsp)
lea    0x8(%rcx),%rsi
lea    (%rcx,%rax,8),%rdi
lea    0x8(%rbp),%r9
lea    0x0(,%rax,8),%r13
add    %rbp,%r13
lea    0x8(%r10),%r14
lea    (%r10,%rax,8),%r11
lea    0x8(%r15),%rdx
mov    %rdx,-0x20(%rsp)
lea    (%r15,%rax,8),%r8
cmp    %r13,%rsi
setb   %r13b
cmp    %rdi,%r9
setb   %r9b
cmp    %r11,%rsi
setb   %dl
cmp    %rdi,%r14
setb   %r11b
cmp    %r8,%rsi
setb   %sil
cmp    %rdi,-0x20(%rsp)
setb   %r14b
test   %r9b,%r13b
jne    12d9 <correct_velocities+0x1d9>
and    %r11b,%dl
jne    12d9 <correct_velocities+0x1d9>
and    %r14b,%sil
mov    -0x58(%rsp),%r9
mov    -0x50(%rsp),%r8
mov    0x38(%rsp),%r11d
mov    -0x30(%rsp),%rdx
mov    -0x38(%rsp),%rdi
mov    -0x40(%rsp),%r14
mov    -0x60(%rsp),%r13
jne    1300 <correct_velocities+0x200>
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movupd 0x8(%rbp,%rsi,8),%xmm4
movupd 0x18(%rbp,%rsi,8),%xmm5
movupd 0x8(%r10,%rsi,8),%xmm6
movupd 0x18(%r10,%rsi,8),%xmm7
movupd 0x8(%r15,%rsi,8),%xmm1
subpd  %xmm1,%xmm6
movupd 0x18(%r15,%rsi,8),%xmm1
subpd  %xmm1,%xmm7
mulpd  %xmm8,%xmm6
addpd  %xmm4,%xmm6
mulpd  %xmm8,%xmm7
addpd  %xmm5,%xmm7
movupd %xmm6,0x8(%rcx,%rsi,8)
movupd %xmm7,0x18(%rcx,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%rbx
jne    1270 <correct_velocities+0x170>
mov    -0x18(%rsp),%r12
cmp    %rbx,-0x48(%rsp)
je     1197 <correct_velocities+0x97>
jmp    1300 <correct_velocities+0x200>
mov    -0x58(%rsp),%r9
mov    -0x50(%rsp),%r8
mov    0x38(%rsp),%r11d
mov    -0x30(%rsp),%rdx
mov    -0x38(%rsp),%rdi
mov    -0x40(%rsp),%r14
mov    -0x60(%rsp),%r13
nopl   0x0(%rax)
mov    %r12,%rsi
test   $0x1,%al
jne    1329 <correct_velocities+0x229>
movsd  (%r10,%r12,8),%xmm1
subsd  (%r15,%r12,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  0x0(%rbp,%r12,8),%xmm1
movsd  %xmm1,(%rcx,%r12,8)
lea    0x1(%r12),%rsi
not    %r12
cmp    -0x28(%rsp),%r12
je     1197 <correct_velocities+0x97>
nopw   0x0(%rax,%rax,1)
movsd  (%r10,%rsi,8),%xmm1
subsd  (%r15,%rsi,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  0x0(%rbp,%rsi,8),%xmm1
movsd  %xmm1,(%rcx,%rsi,8)
movsd  0x8(%r10,%rsi,8),%xmm1
subsd  0x8(%r15,%rsi,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  0x8(%rbp,%rsi,8),%xmm1
movsd  %xmm1,0x8(%rcx,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%rax
jne    1340 <correct_velocities+0x240>
jmp    1197 <correct_velocities+0x97>
test   %r9d,%r9d
jle    1597 <correct_velocities+0x497>
xorpd  0xc68(%rip),%xmm0
divsd  %xmm2,%xmm0
add    $0x1,%r9d
mov    %r11d,%eax
lea    0x1(%rax),%rcx
mov    %rcx,-0x20(%rsp)
lea    -0x1(%rax),%r10
mov    %r10,%r12
and    $0xfffffffffffffffc,%r12
lea    0x1(%r12),%rcx
mov    %rcx,-0x28(%rsp)
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %rax,%rcx
neg    %rcx
mov    %rcx,-0x48(%rsp)
mov    $0x1,%edx
mov    %r9,-0x58(%rsp)
mov    %r10,-0x60(%rsp)
jmp    13fd <correct_velocities+0x2fd>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rdx
cmp    %r9,%rdx
je     1597 <correct_velocities+0x497>
cmp    $0x2,%r11d
jl     13f0 <correct_velocities+0x2f0>
mov    -0x8(%rsp),%rcx
mov    (%rcx,%rdx,8),%rcx
mov    (%r8,%rdx,8),%rbp
mov    -0x10(%rsp),%rsi
mov    (%rsi,%rdx,8),%rsi
mov    $0x1,%edi
cmp    $0x4,%r10
jb     1510 <correct_velocities+0x410>
lea    0x8(%rsi),%r10
lea    (%rsi,%rax,8),%r8
lea    0x8(%rcx),%r9
lea    (%rcx,%rax,8),%r13
lea    0x8(%rbp),%r11
mov    -0x20(%rsp),%rbx
lea    0x0(,%rbx,8),%r14
add    %rbp,%r14
cmp    %r13,%r10
setb   %r13b
cmp    %r8,%r9
setb   %r15b
cmp    %r14,%r10
setb   %bl
cmp    %r8,%r11
setb   %r9b
test   %r15b,%r13b
jne    14f1 <correct_velocities+0x3f1>
and    %r9b,%bl
mov    -0x58(%rsp),%r9
mov    -0x50(%rsp),%r8
mov    0x38(%rsp),%r11d
mov    -0x60(%rsp),%r10
jne    1510 <correct_velocities+0x410>
xor    %edi,%edi
nop
movupd 0x8(%rcx,%rdi,8),%xmm2
movupd 0x18(%rcx,%rdi,8),%xmm3
movupd 0x8(%rbp,%rdi,8),%xmm4
movupd 0x10(%rbp,%rdi,8),%xmm5
subpd  %xmm4,%xmm5
movupd 0x18(%rbp,%rdi,8),%xmm4
movupd 0x20(%rbp,%rdi,8),%xmm6
subpd  %xmm4,%xmm6
mulpd  %xmm1,%xmm5
addpd  %xmm2,%xmm5
mulpd  %xmm1,%xmm6
addpd  %xmm3,%xmm6
movupd %xmm5,0x8(%rsi,%rdi,8)
movupd %xmm6,0x18(%rsi,%rdi,8)
add    $0x4,%rdi
cmp    %rdi,%r12
jne    1490 <correct_velocities+0x390>
mov    -0x28(%rsp),%rdi
cmp    %r12,%r10
je     13f0 <correct_velocities+0x2f0>
jmp    1510 <correct_velocities+0x410>
mov    -0x58(%rsp),%r9
mov    -0x50(%rsp),%r8
mov    0x38(%rsp),%r11d
mov    -0x60(%rsp),%r10
cs nopw 0x0(%rax,%rax,1)
nop
mov    %rdi,%rbx
test   $0x1,%al
jne    1535 <correct_velocities+0x435>
lea    0x1(%rdi),%rbx
movsd  0x8(%rbp,%rdi,8),%xmm2
subsd  0x0(%rbp,%rdi,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rcx,%rdi,8),%xmm2
movsd  %xmm2,(%rsi,%rdi,8)
not    %rdi
cmp    -0x48(%rsp),%rdi
je     13f0 <correct_velocities+0x2f0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x8(%rbp,%rbx,8),%xmm2
subsd  0x0(%rbp,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rcx,%rbx,8),%xmm2
movsd  %xmm2,(%rsi,%rbx,8)
movsd  0x10(%rbp,%rbx,8),%xmm2
subsd  0x8(%rbp,%rbx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  0x8(%rcx,%rbx,8),%xmm2
movsd  %xmm2,0x8(%rsi,%rbx,8)
lea    0x2(%rbx),%rdi
mov    %rdi,%rbx
cmp    %rdi,%rax
jne    1550 <correct_velocities+0x450>
jmp    13f0 <correct_velocities+0x2f0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @correct_velocities(double** nocapture %0, double** nocapture %1, double** nocapture %2, double** nocapture %3, double** nocapture %4, i32 %5, i32 %6, double %7, double %8, double %9) {
%11 = icmp sgt i32 %5, 1
br i1 %11, label %12, label %114
12: ; preds = %10
%13 = icmp slt i32 %6, 1
%14 = fneg double %7
%15 = fdiv double %14, %8
%16 = add i32 %6, 1
%17 = zext i32 %5 to i64
%18 = zext i32 %16 to i64
%19 = add nsw i64 %18, -1
%20 = icmp ult i64 %19, 4
%21 = and i64 %19, -4
%22 = or i64 %21, 1
%23 = insertelement <2 x double> poison, double %15, i64 0
%24 = shufflevector <2 x double> %23, <2 x double> poison, <2 x i32> zeroinitializer
%25 = insertelement <2 x double> poison, double %15, i64 0
%26 = shufflevector <2 x double> %25, <2 x double> poison, <2 x i32> zeroinitializer
%27 = icmp eq i64 %19, %21
%28 = and i64 %18, 1
%29 = icmp eq i64 %28, 0
%30 = sub nsw i64 0, %18
br label %31
31: ; preds = %12, %159
%32 = phi i64 [ 1, %12 ], [ %160, %159 ]
br i1 %13, label %33, label %35
33: ; preds = %31
%34 = add nuw nsw i64 %32, 1
br label %159
35: ; preds = %31
%36 = getelementptr inbounds double*, double** %2, i64 %32
%37 = load double*, double** %36
%38 = add nuw nsw i64 %32, 1
%39 = getelementptr inbounds double*, double** %4, i64 %38
%40 = load double*, double** %39
%41 = getelementptr inbounds double*, double** %4, i64 %32
%42 = load double*, double** %41
%43 = getelementptr inbounds double*, double** %0, i64 %32
%44 = load double*, double** %43
br i1 %20, label %97, label %45
45: ; preds = %35
%46 = getelementptr double, double* %44, i64 1
%47 = getelementptr double, double* %44, i64 %18
%48 = getelementptr double, double* %37, i64 1
%49 = getelementptr double, double* %37, i64 %18
%50 = getelementptr double, double* %40, i64 1
%51 = getelementptr double, double* %40, i64 %18
%52 = getelementptr double, double* %42, i64 1
%53 = getelementptr double, double* %42, i64 %18
%54 = icmp ult double* %46, %49
%55 = icmp ult double* %48, %47
%56 = and i1 %54, %55
%57 = icmp ult double* %46, %51
%58 = icmp ult double* %50, %47
%59 = and i1 %57, %58
%60 = or i1 %56, %59
%61 = icmp ult double* %46, %53
%62 = icmp ult double* %52, %47
%63 = and i1 %61, %62
%64 = or i1 %60, %63
br i1 %64, label %97, label %65
65: ; preds = %45, %65
%66 = phi i64 [ %94, %65 ], [ 0, %45 ]
%67 = or i64 %66, 1
%68 = getelementptr inbounds double, double* %37, i64 %67
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = getelementptr inbounds double, double* %68, i64 2
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = getelementptr inbounds double, double* %40, i64 %67
%75 = bitcast double* %74 to <2 x double>*
%76 = load <2 x double>, <2 x double>* %75
%77 = getelementptr inbounds double, double* %74, i64 2
%78 = bitcast double* %77 to <2 x double>*
%79 = load <2 x double>, <2 x double>* %78
%80 = getelementptr inbounds double, double* %42, i64 %67
%81 = bitcast double* %80 to <2 x double>*
%82 = load <2 x double>, <2 x double>* %81
%83 = getelementptr inbounds double, double* %80, i64 2
%84 = bitcast double* %83 to <2 x double>*
%85 = load <2 x double>, <2 x double>* %84
%86 = fsub <2 x double> %76, %82
%87 = fsub <2 x double> %79, %85
%88 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %24, <2 x double> %86, <2 x double> %70)
%89 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %26, <2 x double> %87, <2 x double> %73)
%90 = getelementptr inbounds double, double* %44, i64 %67
%91 = bitcast double* %90 to <2 x double>*
store <2 x double> %88, <2 x double>* %91
%92 = getelementptr inbounds double, double* %90, i64 2
%93 = bitcast double* %92 to <2 x double>*
store <2 x double> %89, <2 x double>* %93
%94 = add nuw i64 %66, 4
%95 = icmp eq i64 %94, %21
br i1 %95, label %96, label %65
96: ; preds = %65
br i1 %27, label %159, label %97
97: ; preds = %45, %35, %96
%98 = phi i64 [ 1, %45 ], [ 1, %35 ], [ %22, %96 ]
%99 = xor i64 %98, -1
br i1 %29, label %100, label %111
100: ; preds = %97
%101 = getelementptr inbounds double, double* %37, i64 %98
%102 = load double, double* %101
%103 = getelementptr inbounds double, double* %40, i64 %98
%104 = load double, double* %103
%105 = getelementptr inbounds double, double* %42, i64 %98
%106 = load double, double* %105
%107 = fsub double %104, %106
%108 = tail call double @llvm.fmuladd.f64(double %15, double %107, double %102)
%109 = getelementptr inbounds double, double* %44, i64 %98
store double %108, double* %109
%110 = add nuw nsw i64 %98, 1
br label %111
111: ; preds = %100, %97
%112 = phi i64 [ %98, %97 ], [ %110, %100 ]
%113 = icmp eq i64 %99, %30
br i1 %113, label %159, label %136
114: ; preds = %159, %10
%115 = icmp slt i32 %5, 1
br i1 %115, label %261, label %116
116: ; preds = %114
%117 = icmp sgt i32 %6, 1
%118 = fneg double %7
%119 = fdiv double %118, %9
%120 = add nuw i32 %5, 1
%121 = zext i32 %120 to i64
%122 = zext i32 %6 to i64
%123 = add nuw nsw i64 %122, 1
%124 = add nsw i64 %122, -1
%125 = icmp ult i64 %124, 4
%126 = and i64 %124, -4
%127 = or i64 %126, 1
%128 = insertelement <2 x double> poison, double %119, i64 0
%129 = shufflevector <2 x double> %128, <2 x double> poison, <2 x i32> zeroinitializer
%130 = insertelement <2 x double> poison, double %119, i64 0
%131 = shufflevector <2 x double> %130, <2 x double> poison, <2 x i32> zeroinitializer
%132 = icmp eq i64 %124, %126
%133 = and i64 %122, 1
%134 = icmp eq i64 %133, 0
%135 = sub nsw i64 0, %122
br label %162
136: ; preds = %111, %136
%137 = phi i64 [ %157, %136 ], [ %112, %111 ]
%138 = getelementptr inbounds double, double* %37, i64 %137
%139 = load double, double* %138
%140 = getelementptr inbounds double, double* %40, i64 %137
%141 = load double, double* %140
%142 = getelementptr inbounds double, double* %42, i64 %137
%143 = load double, double* %142
%144 = fsub double %141, %143
%145 = tail call double @llvm.fmuladd.f64(double %15, double %144, double %139)
%146 = getelementptr inbounds double, double* %44, i64 %137
store double %145, double* %146
%147 = add nuw nsw i64 %137, 1
%148 = getelementptr inbounds double, double* %37, i64 %147
%149 = load double, double* %148
%150 = getelementptr inbounds double, double* %40, i64 %147
%151 = load double, double* %150
%152 = getelementptr inbounds double, double* %42, i64 %147
%153 = load double, double* %152
%154 = fsub double %151, %153
%155 = tail call double @llvm.fmuladd.f64(double %15, double %154, double %149)
%156 = getelementptr inbounds double, double* %44, i64 %147
store double %155, double* %156
%157 = add nuw nsw i64 %137, 2
%158 = icmp eq i64 %157, %18
br i1 %158, label %159, label %136
159: ; preds = %111, %136, %96, %33
%160 = phi i64 [ %34, %33 ], [ %38, %96 ], [ %38, %136 ], [ %38, %111 ]
%161 = icmp eq i64 %160, %17
br i1 %161, label %114, label %31
162: ; preds = %116, %258
%163 = phi i64 [ 1, %116 ], [ %259, %258 ]
br i1 %117, label %164, label %258
164: ; preds = %162
%165 = getelementptr inbounds double*, double** %3, i64 %163
%166 = load double*, double** %165
%167 = getelementptr inbounds double*, double** %4, i64 %163
%168 = load double*, double** %167
%169 = getelementptr inbounds double*, double** %1, i64 %163
%170 = load double*, double** %169
br i1 %125, label %218, label %171
171: ; preds = %164
%172 = getelementptr double, double* %170, i64 1
%173 = getelementptr double, double* %170, i64 %122
%174 = getelementptr double, double* %166, i64 1
%175 = getelementptr double, double* %166, i64 %122
%176 = getelementptr double, double* %168, i64 1
%177 = getelementptr double, double* %168, i64 %123
%178 = icmp ult double* %172, %175
%179 = icmp ult double* %174, %173
%180 = and i1 %178, %179
%181 = icmp ult double* %172, %177
%182 = icmp ult double* %176, %173
%183 = and i1 %181, %182
%184 = or i1 %180, %183
br i1 %184, label %218, label %185
185: ; preds = %171, %185
%186 = phi i64 [ %215, %185 ], [ 0, %171 ]
%187 = or i64 %186, 1
%188 = getelementptr inbounds double, double* %166, i64 %187
%189 = bitcast double* %188 to <2 x double>*
%190 = load <2 x double>, <2 x double>* %189
%191 = getelementptr inbounds double, double* %188, i64 2
%192 = bitcast double* %191 to <2 x double>*
%193 = load <2 x double>, <2 x double>* %192
%194 = or i64 %186, 2
%195 = getelementptr inbounds double, double* %168, i64 %194
%196 = bitcast double* %195 to <2 x double>*
%197 = load <2 x double>, <2 x double>* %196
%198 = getelementptr inbounds double, double* %195, i64 2
%199 = bitcast double* %198 to <2 x double>*
%200 = load <2 x double>, <2 x double>* %199
%201 = getelementptr inbounds double, double* %168, i64 %187
%202 = bitcast double* %201 to <2 x double>*
%203 = load <2 x double>, <2 x double>* %202
%204 = getelementptr inbounds double, double* %201, i64 2
%205 = bitcast double* %204 to <2 x double>*
%206 = load <2 x double>, <2 x double>* %205
%207 = fsub <2 x double> %197, %203
%208 = fsub <2 x double> %200, %206
%209 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %129, <2 x double> %207, <2 x double> %190)
%210 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %131, <2 x double> %208, <2 x double> %193)
%211 = getelementptr inbounds double, double* %170, i64 %187
%212 = bitcast double* %211 to <2 x double>*
store <2 x double> %209, <2 x double>* %212
%213 = getelementptr inbounds double, double* %211, i64 2
%214 = bitcast double* %213 to <2 x double>*
store <2 x double> %210, <2 x double>* %214
%215 = add nuw i64 %186, 4
%216 = icmp eq i64 %215, %126
br i1 %216, label %217, label %185
217: ; preds = %185
br i1 %132, label %258, label %218
218: ; preds = %171, %164, %217
%219 = phi i64 [ 1, %171 ], [ 1, %164 ], [ %127, %217 ]
%220 = xor i64 %219, -1
br i1 %134, label %221, label %232
221: ; preds = %218
%222 = getelementptr inbounds double, double* %166, i64 %219
%223 = load double, double* %222
%224 = add nuw nsw i64 %219, 1
%225 = getelementptr inbounds double, double* %168, i64 %224
%226 = load double, double* %225
%227 = getelementptr inbounds double, double* %168, i64 %219
%228 = load double, double* %227
%229 = fsub double %226, %228
%230 = tail call double @llvm.fmuladd.f64(double %119, double %229, double %223)
%231 = getelementptr inbounds double, double* %170, i64 %219
store double %230, double* %231
br label %232
232: ; preds = %221, %218
%233 = phi i64 [ %219, %218 ], [ %224, %221 ]
%234 = icmp eq i64 %220, %135
br i1 %234, label %258, label %235
235: ; preds = %232, %235
%236 = phi i64 [ %249, %235 ], [ %233, %232 ]
%237 = getelementptr inbounds double, double* %166, i64 %236
%238 = load double, double* %237
%239 = add nuw nsw i64 %236, 1
%240 = getelementptr inbounds double, double* %168, i64 %239
%241 = load double, double* %240
%242 = getelementptr inbounds double, double* %168, i64 %236
%243 = load double, double* %242
%244 = fsub double %241, %243
%245 = tail call double @llvm.fmuladd.f64(double %119, double %244, double %238)
%246 = getelementptr inbounds double, double* %170, i64 %236
store double %245, double* %246
%247 = getelementptr inbounds double, double* %166, i64 %239
%248 = load double, double* %247
%249 = add nuw nsw i64 %236, 2
%250 = getelementptr inbounds double, double* %168, i64 %249
%251 = load double, double* %250
%252 = getelementptr inbounds double, double* %168, i64 %239
%253 = load double, double* %252
%254 = fsub double %251, %253
%255 = tail call double @llvm.fmuladd.f64(double %119, double %254, double %248)
%256 = getelementptr inbounds double, double* %170, i64 %239
store double %255, double* %256
%257 = icmp eq i64 %249, %122
br i1 %257, label %258, label %235
258: ; preds = %232, %235, %217, %162
%259 = add nuw nsw i64 %163, 1
%260 = icmp eq i64 %259, %121
br i1 %260, label %261, label %162
261: ; preds = %258, %114
ret void
}