void solve_stokes_momentum(double *restrict u, double *restrict v,
                           double *restrict un, double *restrict vn,
                           double *restrict p, double *restrict rho,
                           double *restrict nu, double dt, double dx,
                           double dy) {
  int i, j;
  unsigned int n = 0;
  unsigned int s = 0;
  unsigned int m = 0;
  unsigned int e = 0;
  unsigned int w = 0;
  double dtodx2 = dt / (dx * dx);
  double dtody2 = dt / (dy * dy);
#pragma omp parallel for
  for (j = 1; j < (51 - 1); j++) {
#pragma omp simd safelen(3)
    for (i = 1; i < (101 - 1); i++) {
      m = 101 * j + i;
      n = 101 * (j - 1) + i;
      s = 101 * (j + 1) + i;
      e = 101 * j + (i + 1);
      w = 101 * j + (i - 1);
      u[m] = un[m] - (dt / (rho[m] * 2. * dx)) * (p[e] - p[w]) +
             nu[m] * ((dtodx2 * (un[e] - 2 * un[m] + un[w])) +
                      (dtody2 * (un[s] - 2 * un[m] + un[n])));
      v[m] = (vn[m] - (dt / (rho[m] * 2. * dy)) * (p[s] - p[n]) +
              nu[m] * ((dtodx2 * (vn[e] - 2 * vn[m] + vn[w])) +
                       (dtody2 * (vn[s] - 2 * vn[m] + vn[n])))) -
             (0.0003 * rho[m]);
    }
  }
}


x86 assembly (clang, -O0):
<solve_stokes_momentum>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
movsd  %xmm1,-0x40(%rbp)
movsd  %xmm2,-0x48(%rbp)
movl   $0x0,-0x54(%rbp)
movl   $0x0,-0x58(%rbp)
movl   $0x0,-0x5c(%rbp)
movl   $0x0,-0x60(%rbp)
movl   $0x0,-0x64(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x1,-0x50(%rbp)
cmpl   $0x32,-0x50(%rbp)
jge    140d <solve_stokes_momentum+0x30d>
movl   $0x1,-0x4c(%rbp)
cmpl   $0x64,-0x4c(%rbp)
jge    13fa <solve_stokes_momentum+0x2fa>
imul   $0x65,-0x50(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x50(%rbp),%eax
sub    $0x1,%eax
imul   $0x65,%eax,%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
imul   $0x65,%eax,%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x58(%rbp)
imul   $0x65,-0x50(%rbp),%eax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x60(%rbp)
imul   $0x65,-0x50(%rbp),%eax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x64(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x38(%rbp),%xmm1
mov    -0x30(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  0xdfe(%rip),%xmm2
mulsd  (%rax,%rcx,8),%xmm2
mulsd  -0x40(%rbp),%xmm2
divsd  %xmm2,%xmm1
mov    -0x28(%rbp),%rax
mov    -0x60(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x64(%rbp),%ecx
subsd  (%rax,%rcx,8),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x70(%rbp),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x60(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm4
mov    -0x18(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm3
addsd  %xmm3,%xmm3
subsd  %xmm3,%xmm4
mov    -0x18(%rbp),%rax
mov    -0x64(%rbp),%ecx
addsd  (%rax,%rcx,8),%xmm4
movsd  -0x78(%rbp),%xmm3
mov    -0x18(%rbp),%rax
mov    -0x58(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm5
mov    -0x18(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm6
addsd  %xmm6,%xmm6
subsd  %xmm6,%xmm5
mov    -0x18(%rbp),%rax
mov    -0x54(%rbp),%ecx
addsd  (%rax,%rcx,8),%xmm5
mulsd  %xmm5,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  0xd0b(%rip),%xmm2
mulsd  (%rax,%rcx,8),%xmm2
mulsd  -0x48(%rbp),%xmm2
divsd  %xmm2,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x58(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x54(%rbp),%ecx
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
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x70(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm4
mov    -0x20(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm3
addsd  %xmm3,%xmm3
subsd  %xmm3,%xmm4
mov    -0x20(%rbp),%rax
mov    -0x64(%rbp),%ecx
addsd  (%rax,%rcx,8),%xmm4
movsd  -0x78(%rbp),%xmm3
mov    -0x20(%rbp),%rax
mov    -0x58(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm5
mov    -0x20(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm6
addsd  %xmm6,%xmm6
subsd  %xmm6,%xmm5
mov    -0x20(%rbp),%rax
mov    -0x54(%rbp),%ecx
addsd  (%rax,%rcx,8),%xmm5
mulsd  %xmm5,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x30(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xc28(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    119a <solve_stokes_momentum+0x9a>
jmp    13ff <solve_stokes_momentum+0x2ff>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1189 <solve_stokes_momentum+0x89>
pop    %rbp
ret

x86 assembly (clang, -O2):
<solve_stokes_momentum>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
movapd %xmm2,%xmm14
mov    %rdi,%rbp
mov    0x38(%rsp),%r10
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm0,%xmm3
divsd  %xmm2,%xmm3
movapd %xmm14,%xmm2
mulsd  %xmm14,%xmm2
movapd %xmm0,%xmm7
divsd  %xmm2,%xmm7
xorpd  0xebf(%rip),%xmm0
movapd %xmm1,-0x28(%rsp)
movapd %xmm1,%xmm10
unpcklpd %xmm1,%xmm10
movapd %xmm0,-0x58(%rsp)
pshufd $0x44,%xmm0,%xmm8
movapd %xmm7,-0x48(%rsp)
unpcklpd %xmm7,%xmm7
movapd %xmm3,-0x38(%rsp)
unpcklpd %xmm3,%xmm3
movapd %xmm14,-0x18(%rsp)
unpcklpd %xmm14,%xmm14
lea    0x8(%rcx),%r11
lea    0x8(%rdx),%rax
lea    0x8(%r8),%rbx
mov    %rsi,-0x68(%rsp)
lea    0x330(%rsi),%r14
lea    0x330(%rdi),%r15
lea    0x330(%r10),%r12
mov    %r9,%rsi
lea    0x330(%r9),%r13
mov    $0x1,%r9d
movapd 0xe54(%rip),%xmm15
movapd 0xe5b(%rip),%xmm11
movapd -0x38(%rsp),%xmm9
nopl   0x0(%rax)
imul   $0x65,%r9,%r10
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
movupd 0x328(%rax,%rbp,1),%xmm2
movupd 0x0(%r13,%rbp,1),%xmm13
movapd %xmm13,%xmm4
addpd  %xmm13,%xmm4
movapd %xmm4,%xmm5
mulpd  %xmm10,%xmm5
movupd 0x320(%rbx,%rbp,1),%xmm1
movupd 0x330(%rbx,%rbp,1),%xmm6
subpd  %xmm1,%xmm6
movdqa %xmm8,%xmm1
divpd  %xmm5,%xmm1
mulpd  %xmm6,%xmm1
addpd  %xmm2,%xmm1
movupd (%r12,%rbp,1),%xmm5
movupd 0x330(%rax,%rbp,1),%xmm6
mulpd  %xmm15,%xmm2
addpd  %xmm2,%xmm6
movupd 0x320(%rax,%rbp,1),%xmm0
addpd  %xmm6,%xmm0
movupd 0x650(%rax,%rbp,1),%xmm6
addpd  %xmm2,%xmm6
movupd (%rax,%rbp,1),%xmm2
addpd  %xmm6,%xmm2
mulpd  %xmm7,%xmm2
mulpd  %xmm3,%xmm0
addpd  %xmm2,%xmm0
mulpd  %xmm5,%xmm0
addpd  %xmm1,%xmm0
movupd %xmm0,(%r15,%rbp,1)
movupd 0x328(%r11,%rbp,1),%xmm0
mulpd  %xmm14,%xmm4
movupd 0x650(%rbx,%rbp,1),%xmm1
movupd (%rbx,%rbp,1),%xmm2
subpd  %xmm2,%xmm1
movdqa %xmm8,%xmm2
divpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm2
addpd  %xmm0,%xmm2
movupd 0x330(%r11,%rbp,1),%xmm1
mulpd  %xmm15,%xmm0
addpd  %xmm0,%xmm1
movupd 0x320(%r11,%rbp,1),%xmm4
addpd  %xmm1,%xmm4
movupd 0x650(%r11,%rbp,1),%xmm1
addpd  %xmm0,%xmm1
movupd (%r11,%rbp,1),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm7,%xmm0
mulpd  %xmm3,%xmm4
addpd  %xmm0,%xmm4
mulpd  %xmm5,%xmm4
addpd  %xmm2,%xmm4
mulpd  %xmm11,%xmm13
addpd  %xmm4,%xmm13
movupd %xmm13,(%r14,%rbp,1)
add    $0x10,%rbp
cmp    $0x310,%rbp
jne    11e0 <solve_stokes_momentum+0xe0>
movsd  0x318(%rdx,%r10,8),%xmm0
movsd  0x318(%rsi,%r10,8),%xmm5
movapd %xmm5,%xmm4
addsd  %xmm5,%xmm4
movapd %xmm4,%xmm1
mulsd  -0x28(%rsp),%xmm1
movsd  0x320(%r8,%r10,8),%xmm2
subsd  0x310(%r8,%r10,8),%xmm2
movapd -0x58(%rsp),%xmm6
divsd  %xmm1,%xmm6
mulsd  %xmm2,%xmm6
addsd  %xmm0,%xmm6
mov    0x38(%rsp),%rbp
movsd  0x318(%rbp,%r10,8),%xmm2
movsd  %xmm2,-0x60(%rsp)
movsd  0xca8(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0x320(%rdx,%r10,8),%xmm12
addsd  %xmm0,%xmm12
addsd  0x310(%rdx,%r10,8),%xmm12
addsd  0x640(%rdx,%r10,8),%xmm0
addsd  -0x10(%rdx,%r10,8),%xmm0
movapd -0x48(%rsp),%xmm13
mulsd  %xmm13,%xmm0
mulsd  %xmm9,%xmm12
addsd  %xmm0,%xmm12
mulsd  %xmm2,%xmm12
addsd  %xmm6,%xmm12
movsd  0x318(%rcx,%r10,8),%xmm0
mulsd  -0x18(%rsp),%xmm4
movapd -0x58(%rsp),%xmm1
divsd  %xmm4,%xmm1
movsd  0x640(%r8,%r10,8),%xmm2
subsd  -0x10(%r8,%r10,8),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mulsd  0xc1f(%rip),%xmm0
movsd  0x320(%rcx,%r10,8),%xmm2
addsd  %xmm0,%xmm2
addsd  0x310(%rcx,%r10,8),%xmm2
addsd  0x640(%rcx,%r10,8),%xmm0
movsd  %xmm12,0x318(%rdi,%r10,8)
addsd  -0x10(%rcx,%r10,8),%xmm0
mulsd  %xmm13,%xmm0
mulsd  %xmm9,%xmm2
addsd  %xmm0,%xmm2
mulsd  -0x60(%rsp),%xmm2
addsd  %xmm1,%xmm2
mulsd  0xbd4(%rip),%xmm5
addsd  %xmm2,%xmm5
mov    -0x68(%rsp),%rbp
movsd  %xmm5,0x318(%rbp,%r10,8)
add    $0x1,%r9
add    $0x328,%r11
add    $0x328,%rax
add    $0x328,%rbx
add    $0x328,%r14
add    $0x328,%r15
add    $0x328,%r12
add    $0x328,%r13
cmp    $0x32,%r9
jne    11d0 <solve_stokes_momentum+0xd0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @solve_stokes_momentum(double* noalias nocapture writeonly %0, double* noalias nocapture writeonly %1, double* noalias nocapture %2, double* noalias nocapture %3, double* noalias nocapture %4, double* noalias nocapture %5, double* noalias nocapture %6, double %7, double %8, double %9) {
%11 = fmul double %8, %8
%12 = fdiv double %7, %11
%13 = fmul double %9, %9
%14 = fdiv double %7, %13
%15 = fneg double %7
%16 = insertelement <2 x double> poison, double %8, i64 0
%17 = shufflevector <2 x double> %16, <2 x double> poison, <2 x i32> zeroinitializer
%18 = insertelement <2 x double> poison, double %15, i64 0
%19 = shufflevector <2 x double> %18, <2 x double> poison, <2 x i32> zeroinitializer
%20 = insertelement <2 x double> poison, double %14, i64 0
%21 = shufflevector <2 x double> %20, <2 x double> poison, <2 x i32> zeroinitializer
%22 = insertelement <2 x double> poison, double %12, i64 0
%23 = shufflevector <2 x double> %22, <2 x double> poison, <2 x i32> zeroinitializer
%24 = insertelement <2 x double> poison, double %9, i64 0
%25 = shufflevector <2 x double> %24, <2 x double> poison, <2 x i32> zeroinitializer
br label %26
26: ; preds = %10, %116
%27 = phi i64 [ 1, %10 ], [ %180, %116 ]
%28 = mul nuw nsw i64 %27, 101
br label %29
29: ; preds = %29, %26
%30 = phi i64 [ 0, %26 ], [ %114, %29 ]
%31 = or i64 %30, 1
%32 = add nuw nsw i64 %31, %28
%33 = add nsw i64 %32, -101
%34 = add nuw nsw i64 %32, 101
%35 = add i64 %30, 2
%36 = add nuw nsw i64 %35, %28
%37 = add nsw i64 %32, -1
%38 = getelementptr inbounds double, double* %2, i64 %32
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = getelementptr inbounds double, double* %5, i64 %32
%42 = bitcast double* %41 to <2 x double>*
%43 = load <2 x double>, <2 x double>* %42
%44 = fmul <2 x double> %43, <double 2.000000e+00, double 2.000000e+00>
%45 = fmul <2 x double> %44, %17
%46 = getelementptr inbounds double, double* %4, i64 %36
%47 = bitcast double* %46 to <2 x double>*
%48 = load <2 x double>, <2 x double>* %47
%49 = getelementptr inbounds double, double* %4, i64 %37
%50 = bitcast double* %49 to <2 x double>*
%51 = load <2 x double>, <2 x double>* %50
%52 = fsub <2 x double> %48, %51
%53 = fdiv <2 x double> %19, %45
%54 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %53, <2 x double> %52, <2 x double> %40)
%55 = getelementptr inbounds double, double* %6, i64 %32
%56 = bitcast double* %55 to <2 x double>*
%57 = load <2 x double>, <2 x double>* %56
%58 = getelementptr inbounds double, double* %2, i64 %36
%59 = bitcast double* %58 to <2 x double>*
%60 = load <2 x double>, <2 x double>* %59
%61 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %40, <2 x double> <double -2.000000e+00, double -2.000000e+00>, <2 x double> %60)
%62 = getelementptr inbounds double, double* %2, i64 %37
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = fadd <2 x double> %61, %64
%66 = getelementptr inbounds double, double* %2, i64 %34
%67 = bitcast double* %66 to <2 x double>*
%68 = load <2 x double>, <2 x double>* %67
%69 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %40, <2 x double> <double -2.000000e+00, double -2.000000e+00>, <2 x double> %68)
%70 = getelementptr inbounds double, double* %2, i64 %33
%71 = bitcast double* %70 to <2 x double>*
%72 = load <2 x double>, <2 x double>* %71
%73 = fadd <2 x double> %69, %72
%74 = fmul <2 x double> %21, %73
%75 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %23, <2 x double> %65, <2 x double> %74)
%76 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %57, <2 x double> %75, <2 x double> %54)
%77 = getelementptr inbounds double, double* %0, i64 %32
%78 = bitcast double* %77 to <2 x double>*
store <2 x double> %76, <2 x double>* %78
%79 = getelementptr inbounds double, double* %3, i64 %32
%80 = bitcast double* %79 to <2 x double>*
%81 = load <2 x double>, <2 x double>* %80
%82 = fmul <2 x double> %44, %25
%83 = getelementptr inbounds double, double* %4, i64 %34
%84 = bitcast double* %83 to <2 x double>*
%85 = load <2 x double>, <2 x double>* %84
%86 = getelementptr inbounds double, double* %4, i64 %33
%87 = bitcast double* %86 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = fsub <2 x double> %85, %88
%90 = fdiv <2 x double> %19, %82
%91 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %90, <2 x double> %89, <2 x double> %81)
%92 = getelementptr inbounds double, double* %3, i64 %36
%93 = bitcast double* %92 to <2 x double>*
%94 = load <2 x double>, <2 x double>* %93
%95 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %81, <2 x double> <double -2.000000e+00, double -2.000000e+00>, <2 x double> %94)
%96 = getelementptr inbounds double, double* %3, i64 %37
%97 = bitcast double* %96 to <2 x double>*
%98 = load <2 x double>, <2 x double>* %97
%99 = fadd <2 x double> %95, %98
%100 = getelementptr inbounds double, double* %3, i64 %34
%101 = bitcast double* %100 to <2 x double>*
%102 = load <2 x double>, <2 x double>* %101
%103 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %81, <2 x double> <double -2.000000e+00, double -2.000000e+00>, <2 x double> %102)
%104 = getelementptr inbounds double, double* %3, i64 %33
%105 = bitcast double* %104 to <2 x double>*
%106 = load <2 x double>, <2 x double>* %105
%107 = fadd <2 x double> %103, %106
%108 = fmul <2 x double> %21, %107
%109 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %23, <2 x double> %99, <2 x double> %108)
%110 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %57, <2 x double> %109, <2 x double> %91)
%111 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %43, <2 x double> <double -3.000000e-04, double -3.000000e-04>, <2 x double> %110)
%112 = getelementptr inbounds double, double* %1, i64 %32
%113 = bitcast double* %112 to <2 x double>*
store <2 x double> %111, <2 x double>* %113
%114 = add nuw i64 %30, 2
%115 = icmp eq i64 %114, 98
br i1 %115, label %116, label %29
116: ; preds = %29
%117 = add nuw nsw i64 %28, 99
%118 = add nsw i64 %28, -2
%119 = add nuw nsw i64 %28, 200
%120 = add nuw nsw i64 %28, 100
%121 = add nuw nsw i64 %28, 98
%122 = getelementptr inbounds double, double* %2, i64 %117
%123 = load double, double* %122
%124 = getelementptr inbounds double, double* %5, i64 %117
%125 = load double, double* %124
%126 = fmul double %125, 2.000000e+00
%127 = fmul double %126, %8
%128 = getelementptr inbounds double, double* %4, i64 %120
%129 = load double, double* %128
%130 = getelementptr inbounds double, double* %4, i64 %121
%131 = load double, double* %130
%132 = fsub double %129, %131
%133 = fdiv double %15, %127
%134 = tail call double @llvm.fmuladd.f64(double %133, double %132, double %123)
%135 = getelementptr inbounds double, double* %6, i64 %117
%136 = load double, double* %135
%137 = getelementptr inbounds double, double* %2, i64 %120
%138 = load double, double* %137
%139 = tail call double @llvm.fmuladd.f64(double %123, double -2.000000e+00, double %138)
%140 = getelementptr inbounds double, double* %2, i64 %121
%141 = load double, double* %140
%142 = fadd double %139, %141
%143 = getelementptr inbounds double, double* %2, i64 %119
%144 = load double, double* %143
%145 = tail call double @llvm.fmuladd.f64(double %123, double -2.000000e+00, double %144)
%146 = getelementptr inbounds double, double* %2, i64 %118
%147 = load double, double* %146
%148 = fadd double %145, %147
%149 = fmul double %14, %148
%150 = tail call double @llvm.fmuladd.f64(double %12, double %142, double %149)
%151 = tail call double @llvm.fmuladd.f64(double %136, double %150, double %134)
%152 = getelementptr inbounds double, double* %0, i64 %117
store double %151, double* %152
%153 = getelementptr inbounds double, double* %3, i64 %117
%154 = load double, double* %153
%155 = fmul double %126, %9
%156 = getelementptr inbounds double, double* %4, i64 %119
%157 = load double, double* %156
%158 = getelementptr inbounds double, double* %4, i64 %118
%159 = load double, double* %158
%160 = fsub double %157, %159
%161 = fdiv double %15, %155
%162 = tail call double @llvm.fmuladd.f64(double %161, double %160, double %154)
%163 = getelementptr inbounds double, double* %3, i64 %120
%164 = load double, double* %163
%165 = tail call double @llvm.fmuladd.f64(double %154, double -2.000000e+00, double %164)
%166 = getelementptr inbounds double, double* %3, i64 %121
%167 = load double, double* %166
%168 = fadd double %165, %167
%169 = getelementptr inbounds double, double* %3, i64 %119
%170 = load double, double* %169
%171 = tail call double @llvm.fmuladd.f64(double %154, double -2.000000e+00, double %170)
%172 = getelementptr inbounds double, double* %3, i64 %118
%173 = load double, double* %172
%174 = fadd double %171, %173
%175 = fmul double %14, %174
%176 = tail call double @llvm.fmuladd.f64(double %12, double %168, double %175)
%177 = tail call double @llvm.fmuladd.f64(double %136, double %176, double %162)
%178 = tail call double @llvm.fmuladd.f64(double %125, double -3.000000e-04, double %177)
%179 = getelementptr inbounds double, double* %1, i64 %117
store double %178, double* %179
%180 = add nuw nsw i64 %27, 1
%181 = icmp eq i64 %180, 50
br i1 %181, label %182, label %26
182: ; preds = %116
ret void
}