void calculate_rs(double dt, double dx, double dy, double dz, int imax,
                  int jmax, int kmax, double ***F, double ***G, double ***H,
                  double ***RS) {
  int i, j, k;
  for (i = 1; i < imax + 1; i++) {
    for (j = 1; j < jmax + 1; j++) {
      for (k = 1; k < kmax + 1; k++) {
        RS[i][j][k] = ((F[i][j][k] - F[i - 1][j][k]) / dx +
                       (G[i][j][k] - G[i][j - 1][k]) / dy +
                       (H[i][j][k] - H[i][j][k - 1]) / dz) /
                      dt;
      }
    }
  }
}


x86 assembly (clang, -O0):
<calculate_rs>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
mov    %edi,-0x24(%rbp)
mov    %esi,-0x28(%rbp)
mov    %edx,-0x2c(%rbp)
mov    %rcx,-0x38(%rbp)
mov    %r8,-0x40(%rbp)
mov    %r9,-0x48(%rbp)
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    12a3 <calculate_rs+0x1a3>
movl   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1290 <calculate_rs+0x190>
movl   $0x1,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    127d <calculate_rs+0x17d>
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x10(%rbp),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x40(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  -0x20(%rbp),%xmm1
addsd  %xmm1,%xmm0
divsd  -0x8(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1168 <calculate_rs+0x68>
jmp    1282 <calculate_rs+0x182>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1150 <calculate_rs+0x50>
jmp    1295 <calculate_rs+0x195>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1138 <calculate_rs+0x38>
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculate_rs>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8,-0x10(%rsp)
mov    %rcx,-0x18(%rsp)
mov    %rdx,-0x28(%rsp)
mov    %rsi,-0x40(%rsp)
mov    %rdi,-0x60(%rsp)
test   %edi,%edi
jle    141c <calculate_rs+0x31c>
mov    -0x28(%rsp),%rax
lea    0x1(%rax),%ecx
mov    -0x40(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x8(%rsp)
mov    -0x60(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x60(%rsp)
mov    %rcx,-0x30(%rsp)
lea    -0x1(%rcx),%rax
mov    %rax,%r8
and    $0xfffffffffffffffe,%r8
mov    %rax,-0x38(%rsp)
or     $0x1,%rax
mov    %rax,-0x20(%rsp)
movapd %xmm1,%xmm8
unpcklpd %xmm1,%xmm8
movapd %xmm2,%xmm9
unpcklpd %xmm2,%xmm9
movapd %xmm3,%xmm10
unpcklpd %xmm3,%xmm10
movapd %xmm0,%xmm7
unpcklpd %xmm0,%xmm7
mov    $0x1,%eax
mov    %rax,-0x48(%rsp)
mov    %r9,-0x70(%rsp)
jmp    11c9 <calculate_rs+0xc9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x48(%rsp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rsp)
cmp    -0x60(%rsp),%rax
je     141c <calculate_rs+0x31c>
cmpl   $0x0,-0x40(%rsp)
jle    11b0 <calculate_rs+0xb0>
mov    $0x1,%r12d
jmp    11eb <calculate_rs+0xeb>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r12
cmp    -0x8(%rsp),%r12
je     11b0 <calculate_rs+0xb0>
cmpl   $0x0,-0x28(%rsp)
jle    11e0 <calculate_rs+0xe0>
mov    -0x18(%rsp),%rcx
mov    -0x48(%rsp),%rdx
mov    -0x8(%rcx,%rdx,8),%rax
mov    (%rcx,%rdx,8),%rcx
mov    (%rcx,%r12,8),%r11
mov    (%rax,%r12,8),%r15
mov    -0x10(%rsp),%rax
mov    (%rax,%rdx,8),%rax
mov    -0x8(%rax,%r12,8),%rdi
mov    (%rax,%r12,8),%rsi
mov    (%r9,%rdx,8),%rax
mov    (%rax,%r12,8),%r10
mov    0x38(%rsp),%rax
mov    (%rax,%rdx,8),%rax
mov    (%rax,%r12,8),%r14
mov    $0x1,%r13d
cmpq   $0x2,-0x38(%rsp)
jb     13a0 <calculate_rs+0x2a0>
lea    0x8(%r14),%rbx
mov    -0x30(%rsp),%rax
lea    (%r14,%rax,8),%rbp
lea    0x8(%r11),%rcx
mov    %rcx,-0x50(%rsp)
lea    (%r11,%rax,8),%rcx
lea    0x8(%r15),%rdx
mov    %rdx,-0x68(%rsp)
lea    (%r15,%rax,8),%rdx
lea    0x8(%rsi),%r9
mov    %r9,-0x58(%rsp)
cmp    %rcx,%rbx
setb   -0x71(%rsp)
cmp    %rbp,-0x50(%rsp)
lea    (%rsi,%rax,8),%rcx
setb   -0x50(%rsp)
cmp    %rdx,%rbx
setb   %r9b
cmp    %rbp,-0x68(%rsp)
lea    0x8(%rdi),%rdx
mov    %rdx,-0x68(%rsp)
setb   -0x72(%rsp)
cmp    %rcx,%rbx
setb   %cl
cmp    %rbp,-0x58(%rsp)
lea    (%rdi,%rax,8),%rdx
setb   -0x58(%rsp)
cmp    %rdx,%rbx
setb   -0x73(%rsp)
cmp    %rbp,-0x68(%rsp)
lea    (%r10,%rax,8),%rdx
setb   -0x68(%rsp)
cmp    %rdx,%rbx
setb   %bl
cmp    %rbp,%r10
setb   %dl
mov    -0x50(%rsp),%al
test   %al,-0x71(%rsp)
jne    1390 <calculate_rs+0x290>
and    -0x72(%rsp),%r9b
jne    1390 <calculate_rs+0x290>
and    -0x58(%rsp),%cl
jne    1390 <calculate_rs+0x290>
mov    -0x68(%rsp),%al
and    %al,-0x73(%rsp)
jne    1390 <calculate_rs+0x290>
and    %dl,%bl
mov    -0x70(%rsp),%r9
jne    13a0 <calculate_rs+0x2a0>
xor    %eax,%eax
nopl   0x0(%rax)
movupd 0x8(%r11,%rax,8),%xmm4
movupd 0x8(%r15,%rax,8),%xmm5
subpd  %xmm5,%xmm4
divpd  %xmm8,%xmm4
movupd 0x8(%rsi,%rax,8),%xmm5
movupd 0x8(%rdi,%rax,8),%xmm6
subpd  %xmm6,%xmm5
divpd  %xmm9,%xmm5
addpd  %xmm4,%xmm5
movupd (%r10,%rax,8),%xmm4
movupd 0x8(%r10,%rax,8),%xmm6
subpd  %xmm4,%xmm6
divpd  %xmm10,%xmm6
addpd  %xmm5,%xmm6
divpd  %xmm7,%xmm6
movupd %xmm6,0x8(%r14,%rax,8)
add    $0x2,%rax
cmp    %rax,%r8
jne    1320 <calculate_rs+0x220>
mov    -0x20(%rsp),%r13
cmp    %r8,-0x38(%rsp)
je     11e0 <calculate_rs+0xe0>
jmp    13a0 <calculate_rs+0x2a0>
mov    -0x70(%rsp),%r9
cs nopw 0x0(%rax,%rax,1)
nop
lea    (%r14,%r13,8),%r14
lea    (%r10,%r13,8),%rbp
lea    (%rdi,%r13,8),%rdi
lea    (%rsi,%r13,8),%rsi
lea    (%r15,%r13,8),%rdx
lea    (%r11,%r13,8),%rax
mov    -0x30(%rsp),%rcx
sub    %r13,%rcx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rax,%rbx,8),%xmm4
subsd  (%rdx,%rbx,8),%xmm4
movsd  (%rsi,%rbx,8),%xmm5
subsd  (%rdi,%rbx,8),%xmm5
divsd  %xmm1,%xmm4
divsd  %xmm2,%xmm5
addsd  %xmm4,%xmm5
movsd  0x0(%rbp,%rbx,8),%xmm4
subsd  -0x8(%rbp,%rbx,8),%xmm4
divsd  %xmm3,%xmm4
addsd  %xmm5,%xmm4
divsd  %xmm0,%xmm4
movsd  %xmm4,(%r14,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%rcx
jne    13d0 <calculate_rs+0x2d0>
jmp    11e0 <calculate_rs+0xe0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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