#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double sqrt(double);
void variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat(
    const int nSimplex, const int nPointsPerSimplex, const int nSpace,
    const double d, const int *materialTypes, const double *theta,
    const double *alpha_L_types, const double *alpha_T_types, const double *v,
    const double *u, double *m, double *dm, double *f, double *df, double *a) {
  int i, j, k, I, J, matID;
  const int nSpace2 = nSpace * nSpace;
  double norm_v;
  for (i = 0; i < nSimplex; i++) {
    matID = materialTypes[i];
    for (j = 0; j < nPointsPerSimplex; j++) {
      k = i * nPointsPerSimplex + j;
      m[k] = theta[k] * u[k];
      dm[k] = theta[k];
      norm_v = 0.0;
      for (I = 0; I < nSpace; I++) {
        f[k * nSpace + I] = v[k * nSpace + I] * u[k];
        df[k * nSpace + I] = v[k * nSpace + I];
        norm_v += v[k * nSpace + I] * v[k * nSpace + I];
      }
      norm_v = sqrt(norm_v);
      if (norm_v > 0.0) {
        for (I = 0; I < nSpace; I++) {
          a[k * nSpace2 + I * nSpace + I] =
              theta[k] * d + alpha_T_types[matID] * norm_v +
              (alpha_L_types[matID] - alpha_T_types[matID]) *
                  v[k * nSpace + I] * v[k * nSpace + I] / norm_v;
          for (J = I + 1; J < nSpace; J++) {
            a[k * nSpace2 + I * nSpace + J] =
                (alpha_L_types[matID] - alpha_T_types[matID]) *
                v[k * nSpace + I] * v[k * nSpace + J] / norm_v;
            a[k * nSpace2 + J * nSpace + I] = a[k * nSpace2 + I * nSpace + J];
          }
        }
      } else
        for (I = 0; I < nSpace; I++)
          a[k * nSpace2 + I * nSpace + I] = theta[k] * d;
    }
  }
}


x86 assembly (clang, -O0):
<variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14d8 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3c8>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x48(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    14c5 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3b5>
mov    -0x34(%rbp),%eax
imul   -0x8(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12aa <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x19a>
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x38(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x40(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11ed <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0xdd>
movsd  -0x58(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    145b <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x34b>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1456 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x346>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x18(%rbp),%xmm2
mov    0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
divsd  -0x58(%rbp),%xmm1
addsd  %xmm1,%xmm0
mov    0x48(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1443 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x333>
mov    -0x30(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x44(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
divsd  -0x58(%rbp),%xmm0
mov    0x48(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x44(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x48(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x44(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x44(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1381 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x271>
jmp    1448 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x338>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12d2 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x1c2>
jmp    14b2 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3a2>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14ad <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x39d>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x18(%rbp),%xmm0
mov    0x48(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x40(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1462 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x352>
jmp    14b2 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3a2>
jmp    14b7 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3a7>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1184 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x74>
jmp    14ca <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3ba>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1163 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x53>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x128,%rsp
mov    %r8,0x20(%rsp)
mov    %rcx,0xc0(%rsp)
mov    %esi,(%rsp)
test   %edi,%edi
jle    1771 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x661>
mov    %r9,%r15
movapd %xmm0,%xmm4
mov    %edx,%r14d
mov    0x198(%rsp),%r11
mov    0x190(%rsp),%r9
mov    0x188(%rsp),%r10
mov    0x170(%rsp),%rsi
mov    %edx,%ebp
imul   %ebp,%ebp
movslq %edx,%rbx
mov    (%rsp),%eax
movslq %eax,%rcx
mov    %edi,%edi
mov    %rdi,0xd8(%rsp)
mov    %ebx,%edi
mov    %edi,%edx
and    $0xfffffffe,%edx
mov    %rdx,0x80(%rsp)
mov    %rbx,%rdx
imul   %rcx,%rdx
shl    $0x3,%rdx
mov    %rdx,0xd0(%rsp)
mov    %ebp,%edx
imul   %eax,%edx
mov    %edx,0x54(%rsp)
mov    %rcx,%rax
imul   %rbp,%rax
shl    $0x3,%rax
mov    %rax,0xc8(%rsp)
xorpd  %xmm5,%xmm5
mov    0x168(%rsp),%rax
mov    0x160(%rsp),%r8
mov    %rcx,0xb8(%rsp)
mov    %ecx,%ecx
mov    %rcx,0xf8(%rsp)
mov    %rdi,0x18(%rsp)
lea    -0x1(%rdi),%rcx
mov    %rcx,0x88(%rsp)
mov    %rax,%rdi
lea    0x8(%r9),%rax
mov    %rax,0x78(%rsp)
lea    0x0(,%rbx,8),%rax
mov    %rax,0xf0(%rsp)
lea    0x8(%r10),%rax
mov    %rax,0x70(%rsp)
lea    0x8(%rdi),%rax
mov    %rax,0x68(%rsp)
lea    0x1(%rbx),%eax
lea    0x2(%rbx,%rbx,1),%ecx
mov    %ecx,0x58(%rsp)
mov    %rbp,0x40(%rsp)
lea    0x0(,%rbp,8),%rcx
mov    %rcx,0xe8(%rsp)
mov    %rbx,0xa8(%rsp)
lea    0x8(,%rbx,8),%rcx
mov    %rcx,0x110(%rsp)
xor    %ecx,%ecx
mov    %rcx,0x38(%rsp)
mov    %rdi,0x60(%rsp)
mov    %r11,0x30(%rsp)
mov    %r14d,0x8(%rsp)
mov    %eax,0x5c(%rsp)
mov    %eax,0x4(%rsp)
xor    %eax,%eax
mov    %rax,0x28(%rsp)
movsd  %xmm0,0xe0(%rsp)
jmp    12f5 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x1e5>
cs nopw 0x0(%rax,%rax,1)
mov    0x28(%rsp),%rbp
add    $0x1,%rbp
mov    0xd0(%rsp),%rax
add    %rax,0x78(%rsp)
add    %rax,0x70(%rsp)
add    %rax,0x68(%rsp)
mov    0x54(%rsp),%ecx
add    %ecx,0x4(%rsp)
mov    0x38(%rsp),%rdx
add    %ecx,%edx
mov    %rdx,0x38(%rsp)
add    %ecx,0x8(%rsp)
mov    0x30(%rsp),%rcx
add    0xc8(%rsp),%rcx
mov    %rcx,0x30(%rsp)
add    %rax,0x60(%rsp)
mov    %rbp,%rax
mov    %rbp,0x28(%rsp)
cmp    0xd8(%rsp),%rbp
je     1771 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x661>
cmpl   $0x0,(%rsp)
jle    1290 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x180>
mov    0xc0(%rsp),%rax
mov    0x28(%rsp),%rcx
movslq (%rax,%rcx,4),%r12
imul   0xb8(%rsp),%rcx
mov    %rcx,0x100(%rsp)
mov    0x60(%rsp),%rax
mov    %rax,0x90(%rsp)
mov    0x30(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    0x8(%rsp),%eax
mov    %eax,0x14(%rsp)
mov    0x38(%rsp),%rax
mov    %eax,0x10(%rsp)
mov    0x4(%rsp),%eax
mov    %eax,0xc(%rsp)
mov    0x68(%rsp),%r13
mov    0x70(%rsp),%rdx
mov    0x78(%rsp),%rbp
xor    %ebx,%ebx
jmp    13ce <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x2be>
mov    0x108(%rsp),%rbx
add    $0x1,%rbx
mov    0xf0(%rsp),%rcx
mov    0x98(%rsp),%rbp
add    %rcx,%rbp
mov    0xa0(%rsp),%rdx
add    %rcx,%rdx
add    %rcx,%r13
mov    0x40(%rsp),%rax
add    %eax,0xc(%rsp)
add    %eax,0x10(%rsp)
add    %eax,0x14(%rsp)
mov    0x48(%rsp),%rax
add    0xe8(%rsp),%rax
mov    %rax,0x48(%rsp)
add    %rcx,0x90(%rsp)
cmp    0xf8(%rsp),%rbx
mov    0x170(%rsp),%rsi
je     1290 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x180>
mov    %rbp,0x98(%rsp)
mov    %rdx,0xa0(%rsp)
mov    0x100(%rsp),%rax
mov    %rbx,0x108(%rsp)
add    %rax,%rbx
mov    0x20(%rsp),%rcx
movsd  (%rcx,%rbx,8),%xmm0
mulsd  (%rsi,%rbx,8),%xmm0
mov    0x178(%rsp),%rax
movsd  %xmm0,(%rax,%rbx,8)
movsd  (%rcx,%rbx,8),%xmm0
mov    0x180(%rsp),%rax
movsd  %xmm0,(%rax,%rbx,8)
test   %r14d,%r14d
jle    1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
cmpq   $0x0,0x88(%rsp)
je     14d0 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3c0>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
mov    0x170(%rsp),%rcx
mov    0x80(%rsp),%rdx
mov    0xa0(%rsp),%rsi
mov    0x98(%rsp),%rbp
nopl   (%rax)
movsd  -0x8(%r13,%rax,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
movsd  %xmm1,-0x8(%rsi,%rax,8)
movsd  -0x8(%r13,%rax,8),%xmm1
movsd  %xmm1,-0x8(%rbp,%rax,8)
movsd  -0x8(%r13,%rax,8),%xmm1
movsd  0x0(%r13,%rax,8),%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
mulsd  (%rcx,%rbx,8),%xmm2
movsd  %xmm2,(%rsi,%rax,8)
movsd  0x0(%r13,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
movsd  0x0(%r13,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%rax
cmp    %rax,%rdx
jne    1460 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x350>
jmp    14de <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x3ce>
cs nopw 0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %eax,%eax
mov    0x170(%rsp),%rcx
mov    %rbx,%rdx
imul   0xa8(%rsp),%rdx
mov    %rdx,0xb0(%rsp)
testb  $0x1,0x18(%rsp)
je     1537 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x427>
add    0xb0(%rsp),%rax
movsd  (%rdi,%rax,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
mov    0x188(%rsp),%rcx
movsd  %xmm1,(%rcx,%rax,8)
movsd  (%rdi,%rax,8),%xmm1
mov    0x190(%rsp),%rcx
movsd  %xmm1,(%rcx,%rax,8)
movsd  (%rdi,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
ucomisd %xmm5,%xmm0
jb     15e0 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x4d0>
sqrtsd %xmm0,%xmm0
ucomisd %xmm5,%xmm0
ja     1612 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x502>
test   %r14d,%r14d
jle    1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
cmpq   $0x0,0x88(%rsp)
je     173a <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x62a>
mov    0x10(%rsp),%eax
mov    %eax,%ecx
mov    0xc(%rsp),%eax
mov    %eax,%edx
xor    %eax,%eax
mov    0x80(%rsp),%r9
mov    0x58(%rsp),%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x20(%rsp),%rbp
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  %xmm4,%xmm0
movslq %ecx,%rcx
movsd  %xmm0,(%r11,%rcx,8)
movsd  0x0(%rbp,%rbx,8),%xmm0
mulsd  %xmm4,%xmm0
movslq %edx,%rdx
movsd  %xmm0,(%r11,%rdx,8)
add    $0x2,%rax
add    %esi,%edx
add    %esi,%ecx
cmp    %rax,%r9
jne    1590 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x480>
testb  $0x1,0x18(%rsp)
je     1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
jmp    1747 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x637>
nopl   0x0(%rax,%rax,1)
mov    %r11,%rbp
call   1030 <sqrt@plt>
mov    0x160(%rsp),%r8
mov    0x168(%rsp),%rdi
xorpd  %xmm5,%xmm5
mov    %rbp,%r11
movsd  0xe0(%rsp),%xmm4
ucomisd %xmm5,%xmm0
jbe    154f <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x43f>
test   %r14d,%r14d
jle    1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
mov    %rbx,%rax
imul   0x40(%rsp),%rax
mov    %rax,0x118(%rsp)
mov    0x18(%rsp),%rcx
mov    0x90(%rsp),%rdx
mov    0x48(%rsp),%rsi
mov    0x14(%rsp),%r9d
xor    %eax,%eax
mov    %rbx,0x120(%rsp)
jmp    167b <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x56b>
nop
add    0x5c(%rsp),%r9d
add    0x110(%rsp),%rsi
add    $0x8,%rdx
add    $0xffffffffffffffff,%rcx
cmp    0x18(%rsp),%rax
mov    %r10,%rdi
mov    0x120(%rsp),%rbx
je     1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
movsd  (%r8,%r12,8),%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
mov    0x20(%rsp),%rbp
movsd  0x0(%rbp,%rbx,8),%xmm3
mulsd  %xmm4,%xmm3
addsd  %xmm2,%xmm3
movsd  (%r15,%r12,8),%xmm2
subsd  %xmm1,%xmm2
mov    0xb0(%rsp),%rbp
lea    (%rax,%rbp,1),%rbx
movsd  (%rdi,%rbx,8),%xmm1
mulsd  %xmm1,%xmm2
mulsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
addsd  %xmm3,%xmm2
mov    %rdi,%r10
mov    %eax,%edi
imul   %r14d,%edi
movslq %edi,%rdi
add    0x118(%rsp),%rdi
add    %rax,%rdi
movsd  %xmm2,(%r11,%rdi,8)
add    $0x1,%rax
cmp    0xa8(%rsp),%rax
jge    1650 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x540>
mov    $0x1,%ebp
mov    %r9d,%edi
xchg   %ax,%ax
movsd  (%r15,%r12,8),%xmm1
subsd  (%r8,%r12,8),%xmm1
mulsd  (%r10,%rbx,8),%xmm1
mulsd  (%rdx,%rbp,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rbp,8)
movslq %edi,%rdi
movsd  %xmm1,(%r11,%rdi,8)
add    %r14d,%edi
add    $0x1,%rbp
cmp    %rbp,%rcx
jne    1700 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x5f0>
jmp    1650 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x540>
xor    %eax,%eax
testb  $0x1,0x18(%rsp)
je     1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
mov    %ebx,%ecx
imul   0x40(%rsp),%ecx
mov    0x20(%rsp),%rdx
movsd  (%rdx,%rbx,8),%xmm0
mulsd  %xmm4,%xmm0
add    %eax,%ecx
imul   %r14d,%eax
add    %ecx,%eax
cltq
movsd  %xmm0,(%r11,%rax,8)
jmp    1360 <variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat+0x250>
add    $0x128,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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