void shallowWater_2D_Evaluate(
    const int nPoints, const double h_eps, const double g,
    const double bedFrictionCoefficient, const double bedFrictionPower,
    const double eddyViscosity, const double *x, const double *grad_b,
    const double *h, const double *hu, const double *hv, double *H,
    double *mass_acc, double *dmass_acc_dh, double *mom_u_acc,
    double *dmom_u_acc_dhu, double *mom_v_acc, double *dmom_v_acc_dhv,
    double *mass_adv, double *dmass_adv_dhu, double *dmass_adv_dhv,
    double *mom_u_adv, double *dmom_u_adv_dh, double *dmom_u_adv_dhu,
    double *dmom_u_adv_dhv, double *mom_v_adv, double *dmom_v_adv_dh,
    double *dmom_v_adv_dhu, double *dmom_v_adv_dhv, double *mom_u_diff,
    double *mom_v_diff, double *mom_u_source, double *dmom_u_source_dh,
    double *dmom_u_source_dhu, double *dmom_u_source_dhv, double *mom_v_source,
    double *dmom_v_source_dh, double *dmom_v_source_dhu,
    double *dmom_v_source_dhv) {
  int k;
  double u, du_dh, du_dhu, v, dv_dh, dv_dhv, heval, unorm, hpow;
  for (k = 0; k < nPoints; k++) {
    if (h[k] < h_eps) {
      u = 0.0;
      du_dh = 0.0;
      du_dhu = 0.0;
      v = 0.0;
      dv_dh = 0.0;
      dv_dhv = 0.0;
    } else {
      u = hu[k] / h[k];
      du_dh = -hu[k] / (h[k] * h[k]);
      du_dhu = 1.0 / h[k];
      v = hv[k] / h[k];
      dv_dh = -hv[k] / (h[k] * h[k]);
      dv_dhv = 1.0 / h[k];
    }
    H[k] = h[k] + x[k * 3 + 2];
    mass_acc[k] = h[k];
    dmass_acc_dh[k] = 1.0;
    mass_adv[2 * k + 0] = hu[k];
    dmass_adv_dhu[2 * k + 0] = 1.0;
    mass_adv[2 * k + 1] = hv[k];
    dmass_adv_dhv[2 * k + 1] = 1.0;
    mom_u_acc[k] = hu[k];
    dmom_u_acc_dhu[k] = 1.0;
    mom_u_adv[2 * k + 0] = hu[k] * u + 0.5 * g * h[k] * h[k];
    dmom_u_adv_dh[2 * k + 0] = hu[k] * du_dh + g * h[k];
    dmom_u_adv_dhu[2 * k + 0] = u + hu[k] * du_dhu;
    mom_u_adv[2 * k + 1] = hu[k] * v;
    dmom_u_adv_dh[2 * k + 1] = hu[k] * dv_dh;
    dmom_u_adv_dhu[2 * k + 1] = v;
    dmom_u_adv_dhv[2 * k + 1] = hu[k] * dv_dhv;
    mom_u_diff[2 * k + 0] = eddyViscosity;
    mom_u_diff[2 * k + 1] = eddyViscosity;
    mom_v_acc[k] = hv[k];
    dmom_v_acc_dhv[k] = 1.0;
    mom_v_adv[2 * k + 0] = hv[k] * u;
    dmom_v_adv_dh[2 * k + 0] = hv[k] * du_dh;
    dmom_v_adv_dhu[2 * k + 0] = hv[k] * du_dhu;
    dmom_v_adv_dhv[2 * k + 0] = u;
    mom_v_adv[2 * k + 1] = hv[k] * v + 0.5 * g * h[k] * h[k];
    dmom_v_adv_dh[2 * k + 1] = hv[k] * dv_dh + g * h[k];
    dmom_v_adv_dhv[2 * k + 1] = v + hv[k] * dv_dhv;
    mom_v_diff[2 * k + 0] = eddyViscosity;
    mom_v_diff[2 * k + 1] = eddyViscosity;
    mom_u_source[k] = g * h[k] * grad_b[k * 2 + 0];
    dmom_u_source_dh[k] = g * grad_b[k * 2 + 0];
    mom_v_source[k] = g * h[k] * grad_b[k * 2 + 1];
    dmom_v_source_dh[k] = g * grad_b[k * 2 + 1];
    heval = fmax(1.0e-3, h[k]);
    hpow = pow(heval, -bedFrictionPower - 2.0);
    unorm = sqrt(hu[k] * hu[k] + hv[k] * hv[k]);
    mom_u_source[k] += g * hu[k] * bedFrictionCoefficient * unorm * hpow;
    dmom_u_source_dh[k] += -(bedFrictionPower + 2.0) * g * hu[k] *
                           bedFrictionCoefficient * unorm *
                           pow(heval, -bedFrictionPower - 3.0);
    dmom_u_source_dhu[k] =
        g * bedFrictionCoefficient * unorm * hpow +
        g * bedFrictionCoefficient * hu[k] * hu[k] * hpow / (unorm + h_eps);
    dmom_u_source_dhv[k] =
        g * bedFrictionCoefficient * unorm * hpow +
        g * bedFrictionCoefficient * hu[k] * hv[k] * hpow / (unorm + h_eps);
    mom_v_source[k] += g * hv[k] * bedFrictionCoefficient * unorm * hpow;
    dmom_v_source_dh[k] += -(bedFrictionPower + 2.0) * g * hv[k] *
                           bedFrictionCoefficient * unorm *
                           pow(heval, -bedFrictionPower - 3.0);
    dmom_v_source_dhu[k] =
        g * bedFrictionCoefficient * unorm * hpow +
        g * bedFrictionCoefficient * hu[k] * hv[k] * hpow / (unorm + h_eps);
    dmom_v_source_dhv[k] =
        g * bedFrictionCoefficient * unorm * hpow +
        g * bedFrictionCoefficient * hv[k] * hv[k] * hpow / (unorm + h_eps);
  }
}


x86 assembly (clang, -O0):
<shallowWater_2D_Evaluate>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0xe8(%rbp),%rax
mov    0xe0(%rbp),%rax
mov    0xd8(%rbp),%rax
mov    0xd0(%rbp),%rax
mov    0xc8(%rbp),%rax
mov    0xc0(%rbp),%rax
mov    0xb8(%rbp),%rax
mov    0xb0(%rbp),%rax
mov    0xa8(%rbp),%rax
mov    0xa0(%rbp),%rax
mov    0x98(%rbp),%rax
mov    0x90(%rbp),%rax
mov    0x88(%rbp),%rax
mov    0x80(%rbp),%rax
mov    0x78(%rbp),%rax
mov    0x70(%rbp),%rax
mov    0x68(%rbp),%rax
mov    0x60(%rbp),%rax
mov    0x58(%rbp),%rax
mov    0x50(%rbp),%rax
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)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
movsd  %xmm4,-0x30(%rbp)
mov    %rsi,-0x38(%rbp)
mov    %rdx,-0x40(%rbp)
mov    %rcx,-0x48(%rbp)
mov    %r8,-0x50(%rbp)
mov    %r9,-0x58(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1cf1 <shallowWater_2D_Evaluate+0xbd1>
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    125f <shallowWater_2D_Evaluate+0x13f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
jmp    1365 <shallowWater_2D_Evaluate+0x245>
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  0xd2b(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  0xca8(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
imul   $0x3,-0x5c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x20(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  0xc44(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x50(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  0xc05(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x58(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  0xbc6(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x30(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  0xb97(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  0xb98(%rip),%xmm1
mulsd  -0x18(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x60(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x68(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x78(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x70(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x80(%rbp),%xmm0
mov    0x60(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
mov    0x68(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    0x70(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
mov    0x78(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    0xa0(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    0xa0(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x38(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x40(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  0x9ce(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
mov    0x80(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x70(%rbp),%xmm0
mov    0x88(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x78(%rbp),%xmm0
mov    0x90(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
mov    0x98(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  0x934(%rip),%xmm1
mulsd  -0x18(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x80(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x18(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x88(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm1
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x90(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x98(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    0xa8(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    0xa8(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0xb0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0xb8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0xd0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x5c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0xd8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0x75b(%rip),%xmm1
maxsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0x71b(%rip),%xmm2
subsd  %xmm2,%xmm1
call   1030 <pow@plt>
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
mov    0xb0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x675(%rip),%xmm0
addsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x18(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0x5fe(%rip),%xmm2
subsd  %xmm2,%xmm1
call   1030 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0xb8(%rbp),%xmm0
ff
mov    0xb8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x18(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
movsd  -0xa0(%rbp),%xmm3
ff
addsd  -0x10(%rbp),%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0xc0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x18(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
movsd  -0xa0(%rbp),%xmm3
ff
addsd  -0x10(%rbp),%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0xc8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
mov    0xd0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x4ab(%rip),%xmm0
addsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  -0x18(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0x434(%rip),%xmm2
subsd  %xmm2,%xmm1
call   1030 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0xb0(%rbp),%xmm0
ff
mov    0xd8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x18(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
mov    -0x50(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
movsd  -0xa0(%rbp),%xmm3
ff
addsd  -0x10(%rbp),%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0xe0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x18(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x58(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
movsd  -0xa0(%rbp),%xmm3
ff
addsd  -0x10(%rbp),%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0xe8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    11fc <shallowWater_2D_Evaluate+0xdc>
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<shallowWater_2D_Evaluate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %rdx,0x60(%rsp)
mov    %rsi,0x70(%rsp)
movsd  %xmm4,0x68(%rsp)
movsd  %xmm2,0x18(%rsp)
movsd  %xmm1,(%rsp)
test   %edi,%edi
jle    1798 <shallowWater_2D_Evaluate+0x678>
mov    %r9,%r14
mov    %r8,%rbx
mov    %rcx,%r15
movapd %xmm0,%xmm2
mov    0x198(%rsp),%rdx
movsd  0xe8f(%rip),%xmm1
movsd  (%rsp),%xmm0
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x58(%rsp)
movsd  0xe80(%rip),%xmm1
movsd  0xe80(%rip),%xmm4
addsd  %xmm3,%xmm4
xorpd  0xe94(%rip),%xmm4
subsd  %xmm3,%xmm1
movsd  %xmm1,0x50(%rsp)
mulsd  %xmm0,%xmm4
movapd %xmm4,0x80(%rsp)
movsd  0xe5d(%rip),%xmm1
subsd  %xmm3,%xmm1
movsd  %xmm1,0x30(%rsp)
mulsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x28(%rsp)
mov    %edi,%eax
shl    $0x3,%rax
mov    %rax,0x48(%rsp)
mov    $0x2,%ebp
xor    %r13d,%r13d
movsd  %xmm2,0x38(%rsp)
jmp    1407 <shallowWater_2D_Evaluate+0x2e7>
nop
call   1040 <sqrt@plt>
movsd  0x20(%rsp),%xmm2
movsd  %xmm0,0x8(%rsp)
movsd  (%rbx,%r13,1),%xmm1
mulsd  (%rsp),%xmm1
movsd  0x18(%rsp),%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm1
mov    0x170(%rsp),%rax
addsd  (%rax,%r13,1),%xmm1
movsd  %xmm1,(%rax,%r13,1)
movsd  (%rbx,%r13,1),%xmm1
mulsd  0x80(%rsp),%xmm1
mulsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  0x40(%rsp),%xmm0
movsd  0x30(%rsp),%xmm1
call   1030 <pow@plt>
mulsd  0x10(%rsp),%xmm0
mov    0x178(%rsp),%rax
addsd  (%rax,%r13,1),%xmm0
movsd  %xmm0,(%rax,%r13,1)
movsd  0x28(%rsp),%xmm3
movapd %xmm3,%xmm2
movsd  0x8(%rsp),%xmm5
mulsd  %xmm5,%xmm2
movsd  (%rbx,%r13,1),%xmm0
movapd %xmm3,%xmm1
movapd %xmm3,%xmm4
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x20(%rsp),%xmm0
mulsd  %xmm0,%xmm1
movapd %xmm5,%xmm3
addsd  0x38(%rsp),%xmm3
divsd  %xmm3,%xmm1
movapd %xmm3,%xmm6
movsd  %xmm3,0x78(%rsp)
mulsd  %xmm0,%xmm2
movapd %xmm0,%xmm3
addsd  %xmm2,%xmm1
movsd  %xmm2,0x10(%rsp)
mov    0x180(%rsp),%rax
movsd  %xmm1,(%rax,%r13,1)
movsd  (%rbx,%r13,1),%xmm0
mulsd  %xmm4,%xmm0
mulsd  (%r14,%r13,1),%xmm0
mulsd  %xmm3,%xmm0
divsd  %xmm6,%xmm0
addsd  %xmm2,%xmm0
mov    0x188(%rsp),%rax
movsd  %xmm0,(%rax,%r13,1)
movsd  (%r14,%r13,1),%xmm0
mulsd  (%rsp),%xmm0
movsd  0x18(%rsp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm5,%xmm0
mulsd  %xmm3,%xmm0
mov    0x190(%rsp),%rax
addsd  (%rax,%r13,1),%xmm0
movsd  %xmm0,(%rax,%r13,1)
movsd  (%r14,%r13,1),%xmm0
mulsd  0x80(%rsp),%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm5,%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x40(%rsp),%xmm0
movsd  0x30(%rsp),%xmm1
call   1030 <pow@plt>
movsd  0x38(%rsp),%xmm2
mulsd  0x8(%rsp),%xmm0
mov    %r12,%rdx
addsd  (%r12,%r13,1),%xmm0
movsd  %xmm0,(%r12,%r13,1)
movsd  (%rbx,%r13,1),%xmm0
movsd  0x28(%rsp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  (%r14,%r13,1),%xmm0
movsd  0x20(%rsp),%xmm5
mulsd  %xmm5,%xmm0
movsd  0x78(%rsp),%xmm3
divsd  %xmm3,%xmm0
movsd  0x10(%rsp),%xmm4
addsd  %xmm4,%xmm0
mov    0x1a0(%rsp),%rax
movsd  %xmm0,(%rax,%r13,1)
movsd  (%r14,%r13,1),%xmm0
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm5,%xmm1
divsd  %xmm3,%xmm1
addsd  %xmm4,%xmm1
mov    0x1a8(%rsp),%rax
movsd  %xmm1,(%rax,%r13,1)
add    $0x3,%rbp
add    $0x8,%r13
cmp    %r13,0x48(%rsp)
je     1798 <shallowWater_2D_Evaluate+0x678>
movsd  (%r15,%r13,1),%xmm3
xorpd  %xmm1,%xmm1
ucomisd %xmm3,%xmm2
xorpd  %xmm2,%xmm2
xorpd  %xmm0,%xmm0
ja     146f <shallowWater_2D_Evaluate+0x34f>
movsd  (%rbx,%r13,1),%xmm1
movapd %xmm1,%xmm2
divsd  %xmm3,%xmm2
movapd 0xbfb(%rip),%xmm6
xorpd  %xmm6,%xmm1
movapd 0xbff(%rip),%xmm0
unpcklpd %xmm1,%xmm0
movapd %xmm3,%xmm4
mulsd  %xmm3,%xmm4
movsd  (%r14,%r13,1),%xmm1
movapd %xmm1,%xmm5
xorpd  %xmm6,%xmm5
unpcklpd %xmm5,%xmm1
movapd %xmm3,%xmm5
unpcklpd %xmm4,%xmm5
divpd  %xmm5,%xmm1
divpd  %xmm5,%xmm0
mov    %ebp,%eax
mov    0x70(%rsp),%rcx
addsd  (%rcx,%rax,8),%xmm3
mov    0xd0(%rsp),%rax
movsd  %xmm3,(%rax,%r13,1)
movsd  (%r15,%r13,1),%xmm3
mov    0xd8(%rsp),%rax
movsd  %xmm3,(%rax,%r13,1)
mov    0xe0(%rsp),%rax
movabs $0x3ff0000000000000,%rdi
mov    %rdi,(%rax,%r13,1)
movsd  (%rbx,%r13,1),%xmm3
mov    0x108(%rsp),%rcx
movsd  %xmm3,(%rcx,%r13,2)
mov    0x110(%rsp),%rax
mov    %rdi,(%rax,%r13,2)
movsd  (%r14,%r13,1),%xmm3
movsd  %xmm3,0x8(%rcx,%r13,2)
mov    0x118(%rsp),%rax
mov    %rdi,0x8(%rax,%r13,2)
movsd  (%rbx,%r13,1),%xmm3
mov    0xe8(%rsp),%rax
movsd  %xmm3,(%rax,%r13,1)
mov    0xf0(%rsp),%rax
mov    %rdi,(%rax,%r13,1)
movsd  (%r15,%r13,1),%xmm3
movsd  0x58(%rsp),%xmm8
movapd %xmm8,%xmm4
mulsd  %xmm3,%xmm4
mulsd  %xmm3,%xmm4
movsd  (%rbx,%r13,1),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm4,%xmm3
mov    0x120(%rsp),%rsi
movsd  %xmm3,(%rsi,%r13,2)
movsd  (%r15,%r13,1),%xmm3
movsd  (%rsp),%xmm7
mulsd  %xmm7,%xmm3
movapd %xmm0,%xmm4
unpckhpd %xmm0,%xmm4
movsd  (%rbx,%r13,1),%xmm5
mulsd  %xmm4,%xmm5
addsd  %xmm3,%xmm5
mov    0x128(%rsp),%rcx
movsd  %xmm5,(%rcx,%r13,2)
movsd  (%rbx,%r13,1),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
mov    0x130(%rsp),%rax
movsd  %xmm3,(%rax,%r13,2)
movsd  (%rbx,%r13,1),%xmm3
mulsd  %xmm1,%xmm3
movsd  %xmm3,0x8(%rsi,%r13,2)
movapd %xmm1,%xmm3
unpckhpd %xmm1,%xmm3
movsd  (%rbx,%r13,1),%xmm5
mulsd  %xmm3,%xmm5
movsd  %xmm5,0x8(%rcx,%r13,2)
movlpd %xmm1,0x8(%rax,%r13,2)
movsd  (%rbx,%r13,1),%xmm5
mulsd  %xmm0,%xmm5
mov    0x138(%rsp),%rax
movsd  %xmm5,0x8(%rax,%r13,2)
movsd  0x68(%rsp),%xmm6
mov    0x160(%rsp),%rax
movsd  %xmm6,(%rax,%r13,2)
movsd  %xmm6,0x8(%rax,%r13,2)
movsd  (%r14,%r13,1),%xmm5
mov    0xf8(%rsp),%rax
movsd  %xmm5,(%rax,%r13,1)
mov    0x100(%rsp),%rax
mov    %rdi,(%rax,%r13,1)
movsd  (%r14,%r13,1),%xmm5
mulsd  %xmm2,%xmm5
mov    0x140(%rsp),%rsi
movsd  %xmm5,(%rsi,%r13,2)
mulsd  (%r14,%r13,1),%xmm4
mov    0x148(%rsp),%rcx
movsd  %xmm4,(%rcx,%r13,2)
movsd  (%r14,%r13,1),%xmm4
mulsd  %xmm0,%xmm4
mov    0x150(%rsp),%rax
movsd  %xmm4,(%rax,%r13,2)
mov    0x158(%rsp),%rax
movsd  %xmm2,(%rax,%r13,2)
movsd  (%r15,%r13,1),%xmm2
movapd %xmm8,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm2,%xmm4
movsd  (%r14,%r13,1),%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,0x8(%rsi,%r13,2)
movsd  (%r15,%r13,1),%xmm2
mulsd  %xmm7,%xmm2
mulsd  (%r14,%r13,1),%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,0x8(%rcx,%r13,2)
mulsd  (%r14,%r13,1),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rax,%r13,2)
mov    0x168(%rsp),%rax
movsd  %xmm6,(%rax,%r13,2)
movsd  %xmm6,0x8(%rax,%r13,2)
movsd  (%r15,%r13,1),%xmm0
mulsd  %xmm7,%xmm0
mov    0x60(%rsp),%rax
mulsd  (%rax,%r13,2),%xmm0
mov    0x170(%rsp),%rcx
movsd  %xmm0,(%rcx,%r13,1)
movsd  (%rax,%r13,2),%xmm0
mulsd  %xmm7,%xmm0
mov    0x178(%rsp),%rcx
movsd  %xmm0,(%rcx,%r13,1)
movsd  (%r15,%r13,1),%xmm0
mulsd  %xmm7,%xmm0
mulsd  0x8(%rax,%r13,2),%xmm0
mov    0x190(%rsp),%rcx
movsd  %xmm0,(%rcx,%r13,1)
movsd  0x8(%rax,%r13,2),%xmm0
mulsd  %xmm7,%xmm0
mov    %rdx,%r12
movsd  %xmm0,(%rdx,%r13,1)
movsd  (%r15,%r13,1),%xmm0
maxsd  0x8d2(%rip),%xmm0
movsd  %xmm0,0x40(%rsp)
movsd  0x50(%rsp),%xmm1
call   1030 <pow@plt>
movapd %xmm0,%xmm2
movsd  (%rbx,%r13,1),%xmm0
movsd  (%r14,%r13,1),%xmm1
mulsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0x8a5(%rip),%xmm0
movsd  %xmm2,0x20(%rsp)
jb     11f0 <shallowWater_2D_Evaluate+0xd0>
sqrtsd %xmm0,%xmm0
jmp    11fb <shallowWater_2D_Evaluate+0xdb>
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @shallowWater_2D_Evaluate(i32 %0, double %1, double %2, double %3, double %4, double %5, double* %6, double* %7, double* %8, double* %9, double* %10, double* %11, double* %12, double* %13, double* %14, double* %15, double* %16, double* %17, double* %18, double* %19, double* %20, double* %21, double* %22, double* %23, double* %24, double* %25, double* %26, double* %27, double* %28, double* %29, double* %30, double* %31, double* %32, double* %33, double* %34, double* %35, double* %36, double* %37, double* %38) {
%40 = alloca i32
%41 = alloca double
%42 = alloca double
%43 = alloca double
%44 = alloca double
%45 = alloca double
%46 = alloca double*
%47 = alloca double*
%48 = alloca double*
%49 = alloca double*
%50 = alloca double*
%51 = alloca double*
%52 = alloca double*
%53 = alloca double*
%54 = alloca double*
%55 = alloca double*
%56 = alloca double*
%57 = alloca double*
%58 = alloca double*
%59 = alloca double*
%60 = alloca double*
%61 = alloca double*
%62 = alloca double*
%63 = alloca double*
%64 = alloca double*
%65 = alloca double*
%66 = alloca double*
%67 = alloca double*
%68 = alloca double*
%69 = alloca double*
%70 = alloca double*
%71 = alloca double*
%72 = alloca double*
%73 = alloca double*
%74 = alloca double*
%75 = alloca double*
%76 = alloca double*
%77 = alloca double*
%78 = alloca double*
%79 = alloca i32
%80 = alloca double
%81 = alloca double
%82 = alloca double
%83 = alloca double
%84 = alloca double
%85 = alloca double
%86 = alloca double
%87 = alloca double
%88 = alloca double
store i32 %0, i32* %40
store double %1, double* %41
store double %2, double* %42
store double %3, double* %43
store double %4, double* %44
store double %5, double* %45
store double* %6, double** %46
store double* %7, double** %47
store double* %8, double** %48
store double* %9, double** %49
store double* %10, double** %50
store double* %11, double** %51
store double* %12, double** %52
store double* %13, double** %53
store double* %14, double** %54
store double* %15, double** %55
store double* %16, double** %56
store double* %17, double** %57
store double* %18, double** %58
store double* %19, double** %59
store double* %20, double** %60
store double* %21, double** %61
store double* %22, double** %62
store double* %23, double** %63
store double* %24, double** %64
store double* %25, double** %65
store double* %26, double** %66
store double* %27, double** %67
store double* %28, double** %68
store double* %29, double** %69
store double* %30, double** %70
store double* %31, double** %71
store double* %32, double** %72
store double* %33, double** %73
store double* %34, double** %74
store double* %35, double** %75
store double* %36, double** %76
store double* %37, double** %77
store double* %38, double** %78
store i32 0, i32* %79
br label %89
89: ; preds = %820, %39
%90 = load i32, i32* %79
%91 = load i32, i32* %40
%92 = icmp slt i32 %90, %91
br i1 %92, label %93, label %823
93: ; preds = %89
%94 = load double*, double** %48
%95 = load i32, i32* %79
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %94, i64 %96
%98 = load double, double* %97
%99 = load double, double* %41
%100 = fcmp olt double %98, %99
br i1 %100, label %101, label %102
101: ; preds = %93
store double 0.000000e+00, double* %80
store double 0.000000e+00, double* %81
store double 0.000000e+00, double* %82
store double 0.000000e+00, double* %83
store double 0.000000e+00, double* %84
store double 0.000000e+00, double* %85
br label %173
102: ; preds = %93
%103 = load double*, double** %49
%104 = load i32, i32* %79
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds double, double* %103, i64 %105
%107 = load double, double* %106
%108 = load double*, double** %48
%109 = load i32, i32* %79
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %108, i64 %110
%112 = load double, double* %111
%113 = fdiv double %107, %112
store double %113, double* %80
%114 = load double*, double** %49
%115 = load i32, i32* %79
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %114, i64 %116
%118 = load double, double* %117
%119 = fneg double %118
%120 = load double*, double** %48
%121 = load i32, i32* %79
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %120, i64 %122
%124 = load double, double* %123
%125 = load double*, double** %48
%126 = load i32, i32* %79
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %125, i64 %127
%129 = load double, double* %128
%130 = fmul double %124, %129
%131 = fdiv double %119, %130
store double %131, double* %81
%132 = load double*, double** %48
%133 = load i32, i32* %79
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = fdiv double 1.000000e+00, %136
store double %137, double* %82
%138 = load double*, double** %50
%139 = load i32, i32* %79
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %138, i64 %140
%142 = load double, double* %141
%143 = load double*, double** %48
%144 = load i32, i32* %79
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds double, double* %143, i64 %145
%147 = load double, double* %146
%148 = fdiv double %142, %147
store double %148, double* %83
%149 = load double*, double** %50
%150 = load i32, i32* %79
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds double, double* %149, i64 %151
%153 = load double, double* %152
%154 = fneg double %153
%155 = load double*, double** %48
%156 = load i32, i32* %79
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds double, double* %155, i64 %157
%159 = load double, double* %158
%160 = load double*, double** %48
%161 = load i32, i32* %79
%162 = sext i32 %161 to i64
%163 = getelementptr inbounds double, double* %160, i64 %162
%164 = load double, double* %163
%165 = fmul double %159, %164
%166 = fdiv double %154, %165
store double %166, double* %84
%167 = load double*, double** %48
%168 = load i32, i32* %79
%169 = sext i32 %168 to i64
%170 = getelementptr inbounds double, double* %167, i64 %169
%171 = load double, double* %170
%172 = fdiv double 1.000000e+00, %171
store double %172, double* %85
br label %173
173: ; preds = %102, %101
%174 = load double*, double** %48
%175 = load i32, i32* %79
%176 = sext i32 %175 to i64
%177 = getelementptr inbounds double, double* %174, i64 %176
%178 = load double, double* %177
%179 = load double*, double** %46
%180 = load i32, i32* %79
%181 = mul nsw i32 %180, 3
%182 = add nsw i32 %181, 2
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds double, double* %179, i64 %183
%185 = load double, double* %184
%186 = fadd double %178, %185
%187 = load double*, double** %51
%188 = load i32, i32* %79
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds double, double* %187, i64 %189
store double %186, double* %190
%191 = load double*, double** %48
%192 = load i32, i32* %79
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds double, double* %191, i64 %193
%195 = load double, double* %194
%196 = load double*, double** %52
%197 = load i32, i32* %79
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds double, double* %196, i64 %198
store double %195, double* %199
%200 = load double*, double** %53
%201 = load i32, i32* %79
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds double, double* %200, i64 %202
store double 1.000000e+00, double* %203
%204 = load double*, double** %49
%205 = load i32, i32* %79
%206 = sext i32 %205 to i64
%207 = getelementptr inbounds double, double* %204, i64 %206
%208 = load double, double* %207
%209 = load double*, double** %58
%210 = load i32, i32* %79
%211 = mul nsw i32 2, %210
%212 = add nsw i32 %211, 0
%213 = sext i32 %212 to i64
%214 = getelementptr inbounds double, double* %209, i64 %213
store double %208, double* %214
%215 = load double*, double** %59
%216 = load i32, i32* %79
%217 = mul nsw i32 2, %216
%218 = add nsw i32 %217, 0
%219 = sext i32 %218 to i64
%220 = getelementptr inbounds double, double* %215, i64 %219
store double 1.000000e+00, double* %220
%221 = load double*, double** %50
%222 = load i32, i32* %79
%223 = sext i32 %222 to i64
%224 = getelementptr inbounds double, double* %221, i64 %223
%225 = load double, double* %224
%226 = load double*, double** %58
%227 = load i32, i32* %79
%228 = mul nsw i32 2, %227
%229 = add nsw i32 %228, 1
%230 = sext i32 %229 to i64
%231 = getelementptr inbounds double, double* %226, i64 %230
store double %225, double* %231
%232 = load double*, double** %60
%233 = load i32, i32* %79
%234 = mul nsw i32 2, %233
%235 = add nsw i32 %234, 1
%236 = sext i32 %235 to i64
%237 = getelementptr inbounds double, double* %232, i64 %236
store double 1.000000e+00, double* %237
%238 = load double*, double** %49
%239 = load i32, i32* %79
%240 = sext i32 %239 to i64
%241 = getelementptr inbounds double, double* %238, i64 %240
%242 = load double, double* %241
%243 = load double*, double** %54
%244 = load i32, i32* %79
%245 = sext i32 %244 to i64
%246 = getelementptr inbounds double, double* %243, i64 %245
store double %242, double* %246
%247 = load double*, double** %55
%248 = load i32, i32* %79
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds double, double* %247, i64 %249
store double 1.000000e+00, double* %250
%251 = load double*, double** %49
%252 = load i32, i32* %79
%253 = sext i32 %252 to i64
%254 = getelementptr inbounds double, double* %251, i64 %253
%255 = load double, double* %254
%256 = load double, double* %80
%257 = load double, double* %42
%258 = fmul double 5.000000e-01, %257
%259 = load double*, double** %48
%260 = load i32, i32* %79
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double, double* %259, i64 %261
%263 = load double, double* %262
%264 = fmul double %258, %263
%265 = load double*, double** %48
%266 = load i32, i32* %79
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds double, double* %265, i64 %267
%269 = load double, double* %268
%270 = fmul double %264, %269
%271 = call double @llvm.fmuladd.f64(double %255, double %256, double %270)
%272 = load double*, double** %61
%273 = load i32, i32* %79
%274 = mul nsw i32 2, %273
%275 = add nsw i32 %274, 0
%276 = sext i32 %275 to i64
%277 = getelementptr inbounds double, double* %272, i64 %276
store double %271, double* %277
%278 = load double*, double** %49
%279 = load i32, i32* %79
%280 = sext i32 %279 to i64
%281 = getelementptr inbounds double, double* %278, i64 %280
%282 = load double, double* %281
%283 = load double, double* %81
%284 = load double, double* %42
%285 = load double*, double** %48
%286 = load i32, i32* %79
%287 = sext i32 %286 to i64
%288 = getelementptr inbounds double, double* %285, i64 %287
%289 = load double, double* %288
%290 = fmul double %284, %289
%291 = call double @llvm.fmuladd.f64(double %282, double %283, double %290)
%292 = load double*, double** %62
%293 = load i32, i32* %79
%294 = mul nsw i32 2, %293
%295 = add nsw i32 %294, 0
%296 = sext i32 %295 to i64
%297 = getelementptr inbounds double, double* %292, i64 %296
store double %291, double* %297
%298 = load double, double* %80
%299 = load double*, double** %49
%300 = load i32, i32* %79
%301 = sext i32 %300 to i64
%302 = getelementptr inbounds double, double* %299, i64 %301
%303 = load double, double* %302
%304 = load double, double* %82
%305 = call double @llvm.fmuladd.f64(double %303, double %304, double %298)
%306 = load double*, double** %63
%307 = load i32, i32* %79
%308 = mul nsw i32 2, %307
%309 = add nsw i32 %308, 0
%310 = sext i32 %309 to i64
%311 = getelementptr inbounds double, double* %306, i64 %310
store double %305, double* %311
%312 = load double*, double** %49
%313 = load i32, i32* %79
%314 = sext i32 %313 to i64
%315 = getelementptr inbounds double, double* %312, i64 %314
%316 = load double, double* %315
%317 = load double, double* %83
%318 = fmul double %316, %317
%319 = load double*, double** %61
%320 = load i32, i32* %79
%321 = mul nsw i32 2, %320
%322 = add nsw i32 %321, 1
%323 = sext i32 %322 to i64
%324 = getelementptr inbounds double, double* %319, i64 %323
store double %318, double* %324
%325 = load double*, double** %49
%326 = load i32, i32* %79
%327 = sext i32 %326 to i64
%328 = getelementptr inbounds double, double* %325, i64 %327
%329 = load double, double* %328
%330 = load double, double* %84
%331 = fmul double %329, %330
%332 = load double*, double** %62
%333 = load i32, i32* %79
%334 = mul nsw i32 2, %333
%335 = add nsw i32 %334, 1
%336 = sext i32 %335 to i64
%337 = getelementptr inbounds double, double* %332, i64 %336
store double %331, double* %337
%338 = load double, double* %83
%339 = load double*, double** %63
%340 = load i32, i32* %79
%341 = mul nsw i32 2, %340
%342 = add nsw i32 %341, 1
%343 = sext i32 %342 to i64
%344 = getelementptr inbounds double, double* %339, i64 %343
store double %338, double* %344
%345 = load double*, double** %49
%346 = load i32, i32* %79
%347 = sext i32 %346 to i64
%348 = getelementptr inbounds double, double* %345, i64 %347
%349 = load double, double* %348
%350 = load double, double* %85
%351 = fmul double %349, %350
%352 = load double*, double** %64
%353 = load i32, i32* %79
%354 = mul nsw i32 2, %353
%355 = add nsw i32 %354, 1
%356 = sext i32 %355 to i64
%357 = getelementptr inbounds double, double* %352, i64 %356
store double %351, double* %357
%358 = load double, double* %45
%359 = load double*, double** %69
%360 = load i32, i32* %79
%361 = mul nsw i32 2, %360
%362 = add nsw i32 %361, 0
%363 = sext i32 %362 to i64
%364 = getelementptr inbounds double, double* %359, i64 %363
store double %358, double* %364
%365 = load double, double* %45
%366 = load double*, double** %69
%367 = load i32, i32* %79
%368 = mul nsw i32 2, %367
%369 = add nsw i32 %368, 1
%370 = sext i32 %369 to i64
%371 = getelementptr inbounds double, double* %366, i64 %370
store double %365, double* %371
%372 = load double*, double** %50
%373 = load i32, i32* %79
%374 = sext i32 %373 to i64
%375 = getelementptr inbounds double, double* %372, i64 %374
%376 = load double, double* %375
%377 = load double*, double** %56
%378 = load i32, i32* %79
%379 = sext i32 %378 to i64
%380 = getelementptr inbounds double, double* %377, i64 %379
store double %376, double* %380
%381 = load double*, double** %57
%382 = load i32, i32* %79
%383 = sext i32 %382 to i64
%384 = getelementptr inbounds double, double* %381, i64 %383
store double 1.000000e+00, double* %384
%385 = load double*, double** %50
%386 = load i32, i32* %79
%387 = sext i32 %386 to i64
%388 = getelementptr inbounds double, double* %385, i64 %387
%389 = load double, double* %388
%390 = load double, double* %80
%391 = fmul double %389, %390
%392 = load double*, double** %65
%393 = load i32, i32* %79
%394 = mul nsw i32 2, %393
%395 = add nsw i32 %394, 0
%396 = sext i32 %395 to i64
%397 = getelementptr inbounds double, double* %392, i64 %396
store double %391, double* %397
%398 = load double*, double** %50
%399 = load i32, i32* %79
%400 = sext i32 %399 to i64
%401 = getelementptr inbounds double, double* %398, i64 %400
%402 = load double, double* %401
%403 = load double, double* %81
%404 = fmul double %402, %403
%405 = load double*, double** %66
%406 = load i32, i32* %79
%407 = mul nsw i32 2, %406
%408 = add nsw i32 %407, 0
%409 = sext i32 %408 to i64
%410 = getelementptr inbounds double, double* %405, i64 %409
store double %404, double* %410
%411 = load double*, double** %50
%412 = load i32, i32* %79
%413 = sext i32 %412 to i64
%414 = getelementptr inbounds double, double* %411, i64 %413
%415 = load double, double* %414
%416 = load double, double* %82
%417 = fmul double %415, %416
%418 = load double*, double** %67
%419 = load i32, i32* %79
%420 = mul nsw i32 2, %419
%421 = add nsw i32 %420, 0
%422 = sext i32 %421 to i64
%423 = getelementptr inbounds double, double* %418, i64 %422
store double %417, double* %423
%424 = load double, double* %80
%425 = load double*, double** %68
%426 = load i32, i32* %79
%427 = mul nsw i32 2, %426
%428 = add nsw i32 %427, 0
%429 = sext i32 %428 to i64
%430 = getelementptr inbounds double, double* %425, i64 %429
store double %424, double* %430
%431 = load double*, double** %50
%432 = load i32, i32* %79
%433 = sext i32 %432 to i64
%434 = getelementptr inbounds double, double* %431, i64 %433
%435 = load double, double* %434
%436 = load double, double* %83
%437 = load double, double* %42
%438 = fmul double 5.000000e-01, %437
%439 = load double*, double** %48
%440 = load i32, i32* %79
%441 = sext i32 %440 to i64
%442 = getelementptr inbounds double, double* %439, i64 %441
%443 = load double, double* %442
%444 = fmul double %438, %443
%445 = load double*, double** %48
%446 = load i32, i32* %79
%447 = sext i32 %446 to i64
%448 = getelementptr inbounds double, double* %445, i64 %447
%449 = load double, double* %448
%450 = fmul double %444, %449
%451 = call double @llvm.fmuladd.f64(double %435, double %436, double %450)
%452 = load double*, double** %65
%453 = load i32, i32* %79
%454 = mul nsw i32 2, %453
%455 = add nsw i32 %454, 1
%456 = sext i32 %455 to i64
%457 = getelementptr inbounds double, double* %452, i64 %456
store double %451, double* %457
%458 = load double*, double** %50
%459 = load i32, i32* %79
%460 = sext i32 %459 to i64
%461 = getelementptr inbounds double, double* %458, i64 %460
%462 = load double, double* %461
%463 = load double, double* %84
%464 = load double, double* %42
%465 = load double*, double** %48
%466 = load i32, i32* %79
%467 = sext i32 %466 to i64
%468 = getelementptr inbounds double, double* %465, i64 %467
%469 = load double, double* %468
%470 = fmul double %464, %469
%471 = call double @llvm.fmuladd.f64(double %462, double %463, double %470)
%472 = load double*, double** %66
%473 = load i32, i32* %79
%474 = mul nsw i32 2, %473
%475 = add nsw i32 %474, 1
%476 = sext i32 %475 to i64
%477 = getelementptr inbounds double, double* %472, i64 %476
store double %471, double* %477
%478 = load double, double* %83
%479 = load double*, double** %50
%480 = load i32, i32* %79
%481 = sext i32 %480 to i64
%482 = getelementptr inbounds double, double* %479, i64 %481
%483 = load double, double* %482
%484 = load double, double* %85
%485 = call double @llvm.fmuladd.f64(double %483, double %484, double %478)
%486 = load double*, double** %68
%487 = load i32, i32* %79
%488 = mul nsw i32 2, %487
%489 = add nsw i32 %488, 1
%490 = sext i32 %489 to i64
%491 = getelementptr inbounds double, double* %486, i64 %490
store double %485, double* %491
%492 = load double, double* %45
%493 = load double*, double** %70
%494 = load i32, i32* %79
%495 = mul nsw i32 2, %494
%496 = add nsw i32 %495, 0
%497 = sext i32 %496 to i64
%498 = getelementptr inbounds double, double* %493, i64 %497
store double %492, double* %498
%499 = load double, double* %45
%500 = load double*, double** %70
%501 = load i32, i32* %79
%502 = mul nsw i32 2, %501
%503 = add nsw i32 %502, 1
%504 = sext i32 %503 to i64
%505 = getelementptr inbounds double, double* %500, i64 %504
store double %499, double* %505
%506 = load double, double* %42
%507 = load double*, double** %48
%508 = load i32, i32* %79
%509 = sext i32 %508 to i64
%510 = getelementptr inbounds double, double* %507, i64 %509
%511 = load double, double* %510
%512 = fmul double %506, %511
%513 = load double*, double** %47
%514 = load i32, i32* %79
%515 = mul nsw i32 %514, 2
%516 = add nsw i32 %515, 0
%517 = sext i32 %516 to i64
%518 = getelementptr inbounds double, double* %513, i64 %517
%519 = load double, double* %518
%520 = fmul double %512, %519
%521 = load double*, double** %71
%522 = load i32, i32* %79
%523 = sext i32 %522 to i64
%524 = getelementptr inbounds double, double* %521, i64 %523
store double %520, double* %524
%525 = load double, double* %42
%526 = load double*, double** %47
%527 = load i32, i32* %79
%528 = mul nsw i32 %527, 2
%529 = add nsw i32 %528, 0
%530 = sext i32 %529 to i64
%531 = getelementptr inbounds double, double* %526, i64 %530
%532 = load double, double* %531
%533 = fmul double %525, %532
%534 = load double*, double** %72
%535 = load i32, i32* %79
%536 = sext i32 %535 to i64
%537 = getelementptr inbounds double, double* %534, i64 %536
store double %533, double* %537
%538 = load double, double* %42
%539 = load double*, double** %48
%540 = load i32, i32* %79
%541 = sext i32 %540 to i64
%542 = getelementptr inbounds double, double* %539, i64 %541
%543 = load double, double* %542
%544 = fmul double %538, %543
%545 = load double*, double** %47
%546 = load i32, i32* %79
%547 = mul nsw i32 %546, 2
%548 = add nsw i32 %547, 1
%549 = sext i32 %548 to i64
%550 = getelementptr inbounds double, double* %545, i64 %549
%551 = load double, double* %550
%552 = fmul double %544, %551
%553 = load double*, double** %75
%554 = load i32, i32* %79
%555 = sext i32 %554 to i64
%556 = getelementptr inbounds double, double* %553, i64 %555
store double %552, double* %556
%557 = load double, double* %42
%558 = load double*, double** %47
%559 = load i32, i32* %79
%560 = mul nsw i32 %559, 2
%561 = add nsw i32 %560, 1
%562 = sext i32 %561 to i64
%563 = getelementptr inbounds double, double* %558, i64 %562
%564 = load double, double* %563
%565 = fmul double %557, %564
%566 = load double*, double** %76
%567 = load i32, i32* %79
%568 = sext i32 %567 to i64
%569 = getelementptr inbounds double, double* %566, i64 %568
store double %565, double* %569
%570 = load double*, double** %48
%571 = load i32, i32* %79
%572 = sext i32 %571 to i64
%573 = getelementptr inbounds double, double* %570, i64 %572
%574 = load double, double* %573
%575 = call double @llvm.maxnum.f64(double 1.000000e-03, double %574)
store double %575, double* %86
%576 = load double, double* %86
%577 = load double, double* %44
%578 = fneg double %577
%579 = fsub double %578, 2.000000e+00
%580 = call double @pow(double %576, double %579)
store double %580, double* %88
%581 = load double*, double** %49
%582 = load i32, i32* %79
%583 = sext i32 %582 to i64
%584 = getelementptr inbounds double, double* %581, i64 %583
%585 = load double, double* %584
%586 = load double*, double** %49
%587 = load i32, i32* %79
%588 = sext i32 %587 to i64
%589 = getelementptr inbounds double, double* %586, i64 %588
%590 = load double, double* %589
%591 = load double*, double** %50
%592 = load i32, i32* %79
%593 = sext i32 %592 to i64
%594 = getelementptr inbounds double, double* %591, i64 %593
%595 = load double, double* %594
%596 = load double*, double** %50
%597 = load i32, i32* %79
%598 = sext i32 %597 to i64
%599 = getelementptr inbounds double, double* %596, i64 %598
%600 = load double, double* %599
%601 = fmul double %595, %600
%602 = call double @llvm.fmuladd.f64(double %585, double %590, double %601)
%603 = call double @sqrt(double %602)
store double %603, double* %87
%604 = load double, double* %42
%605 = load double*, double** %49
%606 = load i32, i32* %79
%607 = sext i32 %606 to i64
%608 = getelementptr inbounds double, double* %605, i64 %607
%609 = load double, double* %608
%610 = fmul double %604, %609
%611 = load double, double* %43
%612 = fmul double %610, %611
%613 = load double, double* %87
%614 = fmul double %612, %613
%615 = load double, double* %88
%616 = load double*, double** %71
%617 = load i32, i32* %79
%618 = sext i32 %617 to i64
%619 = getelementptr inbounds double, double* %616, i64 %618
%620 = load double, double* %619
%621 = call double @llvm.fmuladd.f64(double %614, double %615, double %620)
store double %621, double* %619
%622 = load double, double* %44
%623 = fadd double %622, 2.000000e+00
%624 = fneg double %623
%625 = load double, double* %42
%626 = fmul double %624, %625
%627 = load double*, double** %49
%628 = load i32, i32* %79
%629 = sext i32 %628 to i64
%630 = getelementptr inbounds double, double* %627, i64 %629
%631 = load double, double* %630
%632 = fmul double %626, %631
%633 = load double, double* %43
%634 = fmul double %632, %633
%635 = load double, double* %87
%636 = fmul double %634, %635
%637 = load double, double* %86
%638 = load double, double* %44
%639 = fneg double %638
%640 = fsub double %639, 3.000000e+00
%641 = call double @pow(double %637, double %640)
%642 = load double*, double** %72
%643 = load i32, i32* %79
%644 = sext i32 %643 to i64
%645 = getelementptr inbounds double, double* %642, i64 %644
%646 = load double, double* %645
%647 = call double @llvm.fmuladd.f64(double %636, double %641, double %646)
store double %647, double* %645
%648 = load double, double* %42
%649 = load double, double* %43
%650 = fmul double %648, %649
%651 = load double, double* %87
%652 = fmul double %650, %651
%653 = load double, double* %88
%654 = load double, double* %42
%655 = load double, double* %43
%656 = fmul double %654, %655
%657 = load double*, double** %49
%658 = load i32, i32* %79
%659 = sext i32 %658 to i64
%660 = getelementptr inbounds double, double* %657, i64 %659
%661 = load double, double* %660
%662 = fmul double %656, %661
%663 = load double*, double** %49
%664 = load i32, i32* %79
%665 = sext i32 %664 to i64
%666 = getelementptr inbounds double, double* %663, i64 %665
%667 = load double, double* %666
%668 = fmul double %662, %667
%669 = load double, double* %88
%670 = fmul double %668, %669
%671 = load double, double* %87
%672 = load double, double* %41
%673 = fadd double %671, %672
%674 = fdiv double %670, %673
%675 = call double @llvm.fmuladd.f64(double %652, double %653, double %674)
%676 = load double*, double** %73
%677 = load i32, i32* %79
%678 = sext i32 %677 to i64
%679 = getelementptr inbounds double, double* %676, i64 %678
store double %675, double* %679
%680 = load double, double* %42
%681 = load double, double* %43
%682 = fmul double %680, %681
%683 = load double, double* %87
%684 = fmul double %682, %683
%685 = load double, double* %88
%686 = load double, double* %42
%687 = load double, double* %43
%688 = fmul double %686, %687
%689 = load double*, double** %49
%690 = load i32, i32* %79
%691 = sext i32 %690 to i64
%692 = getelementptr inbounds double, double* %689, i64 %691
%693 = load double, double* %692
%694 = fmul double %688, %693
%695 = load double*, double** %50
%696 = load i32, i32* %79
%697 = sext i32 %696 to i64
%698 = getelementptr inbounds double, double* %695, i64 %697
%699 = load double, double* %698
%700 = fmul double %694, %699
%701 = load double, double* %88
%702 = fmul double %700, %701
%703 = load double, double* %87
%704 = load double, double* %41
%705 = fadd double %703, %704
%706 = fdiv double %702, %705
%707 = call double @llvm.fmuladd.f64(double %684, double %685, double %706)
%708 = load double*, double** %74
%709 = load i32, i32* %79
%710 = sext i32 %709 to i64
%711 = getelementptr inbounds double, double* %708, i64 %710
store double %707, double* %711
%712 = load double, double* %42
%713 = load double*, double** %50
%714 = load i32, i32* %79
%715 = sext i32 %714 to i64
%716 = getelementptr inbounds double, double* %713, i64 %715
%717 = load double, double* %716
%718 = fmul double %712, %717
%719 = load double, double* %43
%720 = fmul double %718, %719
%721 = load double, double* %87
%722 = fmul double %720, %721
%723 = load double, double* %88
%724 = load double*, double** %75
%725 = load i32, i32* %79
%726 = sext i32 %725 to i64
%727 = getelementptr inbounds double, double* %724, i64 %726
%728 = load double, double* %727
%729 = call double @llvm.fmuladd.f64(double %722, double %723, double %728)
store double %729, double* %727
%730 = load double, double* %44
%731 = fadd double %730, 2.000000e+00
%732 = fneg double %731
%733 = load double, double* %42
%734 = fmul double %732, %733
%735 = load double*, double** %50
%736 = load i32, i32* %79
%737 = sext i32 %736 to i64
%738 = getelementptr inbounds double, double* %735, i64 %737
%739 = load double, double* %738
%740 = fmul double %734, %739
%741 = load double, double* %43
%742 = fmul double %740, %741
%743 = load double, double* %87
%744 = fmul double %742, %743
%745 = load double, double* %86
%746 = load double, double* %44
%747 = fneg double %746
%748 = fsub double %747, 3.000000e+00
%749 = call double @pow(double %745, double %748)
%750 = load double*, double** %76
%751 = load i32, i32* %79
%752 = sext i32 %751 to i64
%753 = getelementptr inbounds double, double* %750, i64 %752
%754 = load double, double* %753
%755 = call double @llvm.fmuladd.f64(double %744, double %749, double %754)
store double %755, double* %753
%756 = load double, double* %42
%757 = load double, double* %43
%758 = fmul double %756, %757
%759 = load double, double* %87
%760 = fmul double %758, %759
%761 = load double, double* %88
%762 = load double, double* %42
%763 = load double, double* %43
%764 = fmul double %762, %763
%765 = load double*, double** %49
%766 = load i32, i32* %79
%767 = sext i32 %766 to i64
%768 = getelementptr inbounds double, double* %765, i64 %767
%769 = load double, double* %768
%770 = fmul double %764, %769
%771 = load double*, double** %50
%772 = load i32, i32* %79
%773 = sext i32 %772 to i64
%774 = getelementptr inbounds double, double* %771, i64 %773
%775 = load double, double* %774
%776 = fmul double %770, %775
%777 = load double, double* %88
%778 = fmul double %776, %777
%779 = load double, double* %87
%780 = load double, double* %41
%781 = fadd double %779, %780
%782 = fdiv double %778, %781
%783 = call double @llvm.fmuladd.f64(double %760, double %761, double %782)
%784 = load double*, double** %77
%785 = load i32, i32* %79
%786 = sext i32 %785 to i64
%787 = getelementptr inbounds double, double* %784, i64 %786
store double %783, double* %787
%788 = load double, double* %42
%789 = load double, double* %43
%790 = fmul double %788, %789
%791 = load double, double* %87
%792 = fmul double %790, %791
%793 = load double, double* %88
%794 = load double, double* %42
%795 = load double, double* %43
%796 = fmul double %794, %795
%797 = load double*, double** %50
%798 = load i32, i32* %79
%799 = sext i32 %798 to i64
%800 = getelementptr inbounds double, double* %797, i64 %799
%801 = load double, double* %800
%802 = fmul double %796, %801
%803 = load double*, double** %50
%804 = load i32, i32* %79
%805 = sext i32 %804 to i64
%806 = getelementptr inbounds double, double* %803, i64 %805
%807 = load double, double* %806
%808 = fmul double %802, %807
%809 = load double, double* %88
%810 = fmul double %808, %809
%811 = load double, double* %87
%812 = load double, double* %41
%813 = fadd double %811, %812
%814 = fdiv double %810, %813
%815 = call double @llvm.fmuladd.f64(double %792, double %793, double %814)
%816 = load double*, double** %78
%817 = load i32, i32* %79
%818 = sext i32 %817 to i64
%819 = getelementptr inbounds double, double* %816, i64 %818
store double %815, double* %819
br label %820
820: ; preds = %173
%821 = load i32, i32* %79
%822 = add nsw i32 %821, 1
store i32 %822, i32* %79
br label %89
823: ; preds = %89
ret void
}

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