#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 */
int /*<<< orphan*/ BasisFuns(int, double, int, double *, double *, double *,
                             double *);
int FindSpan(int, int, double, double *);
void computeGNmatPnt2Pln(int degU, int degV, double *cpw, double *cp, int ncp,
                         int kcp, double *knotU, double *knotV, double *us,
                         int nus, double *dirOffset, double *cpDir,
                         double *sgnDistDiff, double *HessApprx,
                         double *gradIent, int *matInds, int *derivInds,
                         double *derivVals, double *leftU, double *rightU,
                         double *NU, double *leftV, double *rightV,
                         double *NV) {
  int i, j, ii, twoii, threeii, spanU, spanV, ind, ind2, dind, fourncp, myint,
      numCP, numRows, orderU, orderV, nnz, rowInd, colInd, cpInd;
  double bspnt[4];
  fourncp = 4 * ncp;
  orderU = degU + 1;
  orderV = degV + 1;
  nnz = orderU * orderV;
  numCP = ncp * kcp;
  numRows = 2 * (numCP - 4);
  for (ii = 0; ii < nus; ii++) {
    twoii = 2 * ii;
    threeii = 3 * ii;
    spanU = FindSpan(ncp, degU, us[twoii], knotU);
    BasisFuns(spanU, us[twoii], degU, knotU, NU, leftU, rightU);
    ind = spanU - degU;
    spanV = FindSpan(kcp, degV, us[twoii + 1], knotV);
    BasisFuns(spanV, us[twoii + 1], degV, knotV, NV, leftV, rightV);
    ind2 = spanV - degV;
    myint = ind2 * fourncp + ind * 4;
    bspnt[0] = 0.0;
    bspnt[1] = 0.0;
    bspnt[2] = 0.0;
    bspnt[3] = 0.0;
    for (i = 0; i <= degV; i++) {
      for (j = 0; j <= degU; j++) {
        bspnt[0] += NV[i] * NU[j] * cpw[i * fourncp + j * 4 + myint];
        bspnt[1] += NV[i] * NU[j] * cpw[i * fourncp + j * 4 + myint + 1];
        bspnt[2] += NV[i] * NU[j] * cpw[i * fourncp + j * 4 + myint + 2];
        bspnt[3] += NV[i] * NU[j] * cpw[i * fourncp + j * 4 + myint + 3];
      }
    }
    dind = 0;
    for (i = 0; i <= degV; i++) {
      for (j = 0; j <= degU; j++) {
        cpInd = (ind2 + i) * ncp + ind + j;
        derivInds[dind] = cpInd;
        derivVals[2 * dind] = NV[i] * NU[j] * cpw[4 * cpInd + 3] *
                              (dirOffset[threeii] * cpDir[3 * cpInd] +
                               dirOffset[threeii + 1] * cpDir[3 * cpInd + 1] +
                               dirOffset[threeii + 2] * cpDir[3 * cpInd + 2]) /
                              bspnt[3];
        derivVals[2 * dind + 1] =
            NV[i] * NU[j] *
            (dirOffset[threeii] * (bspnt[3] * cp[3 * cpInd] - bspnt[0]) +
             dirOffset[threeii + 1] *
                 (bspnt[3] * cp[3 * cpInd + 1] - bspnt[1]) +
             dirOffset[threeii + 2] *
                 (bspnt[3] * cp[3 * cpInd + 2] - bspnt[2])) /
            (bspnt[3] * bspnt[3]);
        dind++;
      }
    }
    for (i = 0; i < nnz; i++) {
      colInd = matInds[derivInds[i]];
      if (colInd < numCP) {
        gradIent[2 * colInd] += sgnDistDiff[ii] * derivVals[2 * i];
        gradIent[2 * colInd + 1] += sgnDistDiff[ii] * derivVals[2 * i + 1];
        for (j = 0; j < nnz; j++) {
          rowInd = matInds[derivInds[j]];
          if (rowInd < numCP) {
            if (rowInd > colInd) {
              HessApprx[(2 * colInd) * numRows + 2 * rowInd] +=
                  derivVals[2 * i] * derivVals[2 * j];
              HessApprx[(2 * colInd) * numRows + 2 * rowInd + 1] +=
                  derivVals[2 * i] * derivVals[2 * j + 1];
              HessApprx[(2 * colInd + 1) * numRows + 2 * rowInd] +=
                  derivVals[2 * i + 1] * derivVals[2 * j];
              HessApprx[(2 * colInd + 1) * numRows + 2 * rowInd + 1] +=
                  derivVals[2 * i + 1] * derivVals[2 * j + 1];
            } else if (rowInd == colInd) {
              HessApprx[(2 * colInd) * numRows + 2 * rowInd] +=
                  derivVals[2 * i] * derivVals[2 * j];
              HessApprx[(2 * colInd) * numRows + 2 * rowInd + 1] +=
                  derivVals[2 * i] * derivVals[2 * j + 1];
              HessApprx[(2 * colInd + 1) * numRows + 2 * rowInd + 1] +=
                  derivVals[2 * i + 1] * derivVals[2 * j + 1];
            }
          }
        }
      }
    }
  }
  for (i = 1; i < numRows; i++) {
    for (j = 0; j < i; j++) {
      HessApprx[i * numRows + j] = HessApprx[j * numRows + i];
    }
  }
}


x86 assembly (clang, -O0):
<computeGNmatPnt2Pln>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
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),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
mov    -0x5c(%rbp),%eax
imul   -0x60(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x1c(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
sub    $0x4,%eax
shl    $0x1,%eax
mov    %eax,-0x58(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    1aa2 <computeGNmatPnt2Pln+0x982>
mov    -0x2c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x30(%rbp)
imul   $0x3,-0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x4(%rbp),%esi
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rdx
call   1040 <FindSpan@plt>
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%edi
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x4(%rbp),%esi
mov    0x10(%rbp),%rdx
mov    0x80(%rbp),%rcx
mov    0x70(%rbp),%r8
mov    0x78(%rbp),%r9
call   1030 <BasisFuns@plt>
mov    -0x38(%rbp),%eax
sub    -0x4(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x20(%rbp),%edi
mov    -0x8(%rbp),%esi
mov    0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rdx
call   1040 <FindSpan@plt>
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%edi
mov    0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%esi
mov    0x18(%rbp),%rdx
mov    0x98(%rbp),%rcx
mov    0x88(%rbp),%r8
mov    0x90(%rbp),%r9
call   1030 <BasisFuns@plt>
mov    -0x3c(%rbp),%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
imul   -0x4c(%rbp),%eax
mov    -0x40(%rbp),%ecx
shl    $0x2,%ecx
add    %ecx,%eax
mov    %eax,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     147c <computeGNmatPnt2Pln+0x35c>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1469 <computeGNmatPnt2Pln+0x349>
mov    0x98(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    0x98(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
add    -0x50(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    0x98(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
add    -0x50(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x80(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    0x98(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
add    -0x50(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12fa <computeGNmatPnt2Pln+0x1da>
jmp    146e <computeGNmatPnt2Pln+0x34e>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12e7 <computeGNmatPnt2Pln+0x1c7>
movl   $0x0,-0x48(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     16fe <computeGNmatPnt2Pln+0x5de>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     16eb <computeGNmatPnt2Pln+0x5cb>
mov    -0x44(%rbp),%eax
add    -0x24(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x40(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x70(%rbp)
mov    -0x70(%rbp),%edx
mov    0x60(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    0x98(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x70(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x38(%rbp),%rax
imul   $0x3,-0x70(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm3
mov    0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x38(%rbp),%rax
imul   $0x3,-0x70(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
mov    0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x38(%rbp),%rax
imul   $0x3,-0x70(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
divsd  -0x78(%rbp),%xmm0
mov    0x68(%rbp),%rax
mov    -0x48(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x98(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x78(%rbp),%xmm3
mov    -0x18(%rbp),%rax
imul   $0x3,-0x70(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm4
movsd  -0x90(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm4,%xmm3
addsd  %xmm1,%xmm3
mov    0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x78(%rbp),%xmm4
mov    -0x18(%rbp),%rax
imul   $0x3,-0x70(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm6
movsd  -0x88(%rbp),%xmm5
ff
movq   %xmm5,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm5
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm1
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
mov    0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x78(%rbp),%xmm3
mov    -0x18(%rbp),%rax
imul   $0x3,-0x70(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm5
movsd  -0x80(%rbp),%xmm4
movq   %xmm4,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm4
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x78(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    0x68(%rbp),%rax
mov    -0x48(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    149d <computeGNmatPnt2Pln+0x37d>
jmp    16f0 <computeGNmatPnt2Pln+0x5d0>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    148a <computeGNmatPnt2Pln+0x36a>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x64(%rbp),%eax
jge    1a8f <computeGNmatPnt2Pln+0x96f>
mov    0x58(%rbp),%rax
mov    0x60(%rbp),%rcx
movslq -0x24(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x54(%rbp),%eax
jge    1a7c <computeGNmatPnt2Pln+0x95c>
mov    0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x50(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x50(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x64(%rbp),%eax
jge    1a77 <computeGNmatPnt2Pln+0x957>
mov    0x58(%rbp),%rax
mov    0x60(%rbp),%rcx
movslq -0x28(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x54(%rbp),%eax
jge    1a64 <computeGNmatPnt2Pln+0x944>
mov    -0x68(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jle    194f <computeGNmatPnt2Pln+0x82f>
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1a5f <computeGNmatPnt2Pln+0x93f>
mov    -0x68(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jne    1a5a <computeGNmatPnt2Pln+0x93a>
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x68(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x58(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1a5f <computeGNmatPnt2Pln+0x93f>
jmp    1a64 <computeGNmatPnt2Pln+0x944>
jmp    1a69 <computeGNmatPnt2Pln+0x949>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    17bc <computeGNmatPnt2Pln+0x69c>
jmp    1a7c <computeGNmatPnt2Pln+0x95c>
jmp    1a81 <computeGNmatPnt2Pln+0x961>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1705 <computeGNmatPnt2Pln+0x5e5>
jmp    1a94 <computeGNmatPnt2Pln+0x974>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11d6 <computeGNmatPnt2Pln+0xb6>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x58(%rbp),%eax
jge    1b15 <computeGNmatPnt2Pln+0x9f5>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1b02 <computeGNmatPnt2Pln+0x9e2>
mov    0x48(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x58(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x58(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1abc <computeGNmatPnt2Pln+0x99c>
jmp    1b07 <computeGNmatPnt2Pln+0x9e7>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1aa9 <computeGNmatPnt2Pln+0x989>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<computeGNmatPnt2Pln>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1f8,%rsp
mov    %rdx,%r13
mov    %rsi,0x48(%rsp)
mov    0x248(%rsp),%r15d
mov    %r9d,0xac(%rsp)
mov    %r9d,%edx
mov    %r8,%rax
mov    %r8,0x28(%rsp)
imul   %r8d,%edx
mov    %rdx,0x38(%rsp)
lea    (%rdx,%rdx,1),%eax
add    $0xfffffff8,%eax
mov    %eax,0x40(%rsp)
test   %r15d,%r15d
jle    2052 <computeGNmatPnt2Pln+0xf32>
mov    %edi,%esi
mov    0x288(%rsp),%r8
mov    0x258(%rsp),%rdx
lea    0x1(%rsi),%eax
mov    0x48(%rsp),%r9
lea    0x1(%r9),%r11d
mov    %r11d,%r14d
imul   %eax,%r11d
mov    %eax,0xb4(%rsp)
mov    %eax,%ebx
lea    -0x1(%rbx),%rdi
mov    %rbx,%rax
shl    $0x4,%rax
add    %r8,%rax
mov    %rax,0x130(%rsp)
lea    0x0(,%rbx,8),%rax
lea    (%rax,%rax,2),%rax
add    %rax,%rdx
mov    %rdx,0x120(%rsp)
add    %rcx,%rax
mov    %rax,0x128(%rsp)
mov    %rbx,%rbp
shl    $0x5,%rbp
mov    %edi,%eax
add    %eax,%eax
mov    %rax,0x188(%rsp)
setb   0x11(%rsp)
mov    $0x4,%edx
mov    %edi,%eax
mul    %edx
mov    %eax,%r10d
seto   0x13(%rsp)
mov    $0x10,%edx
mov    %rdi,0x168(%rsp)
mov    %rdi,%rax
mul    %rdx
mov    %rax,0x178(%rsp)
mov    %r15d,%eax
mov    %rax,0x108(%rsp)
lea    -0x18(%rbp,%r13,1),%rax
mov    %rax,0x118(%rsp)
seto   0x12(%rsp)
mov    %ebx,%eax
and    $0xfffffffe,%eax
mov    %rax,0x88(%rsp)
lea    0x0(,%rsi,4),%eax
neg    %eax
mov    %rax,0x100(%rsp)
mov    %r9d,%eax
neg    %eax
mov    %eax,0xa4(%rsp)
mov    $0x7,%edx
mov    %r10,%rax
mov    %r10,0x180(%rsp)
sub    %r10d,%edx
mov    %edx,0xa0(%rsp)
mov    0x2a0(%rsp),%rdx
mov    0x250(%rsp),%r9
mov    0x28(%rsp),%rax
lea    0x0(,%rax,4),%eax
mov    %r11d,0xa8(%rsp)
mov    %r11d,%edi
lea    0x8(%r8),%rbp
mov    %rbp,0x170(%rsp)
lea    (%rdx,%rbx,8),%rdx
mov    %rdx,0x110(%rsp)
lea    0x8(%r9),%rdx
mov    %rdx,0xf8(%rsp)
mov    %rbx,0xb8(%rsp)
lea    0x0(,%rbx,4),%rdx
mov    %rdx,0xf0(%rsp)
mov    %eax,0x24(%rsp)
mov    %eax,%eax
mov    %rax,0xe8(%rsp)
mov    %rdi,0x190(%rsp)
lea    0x0(,%rdi,4),%rax
mov    %rax,0x198(%rsp)
xor    %ebx,%ebx
mov    %rcx,0x90(%rsp)
mov    %r14,0xd0(%rsp)
mov    %rsi,0x50(%rsp)
mov    %r13,0xe0(%rsp)
jmp    1357 <computeGNmatPnt2Pln+0x237>
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x80(%rsp),%rbx
add    $0x1,%rbx
cmp    0x108(%rsp),%rbx
mov    0x50(%rsp),%rsi
mov    0xd0(%rsp),%r14
je     2052 <computeGNmatPnt2Pln+0xf32>
mov    %rbx,0x80(%rsp)
shl    $0x4,%rbx
mov    0x240(%rsp),%r12
movsd  (%r12,%rbx,1),%xmm0
mov    0x28(%rsp),%rdi
mov    0x230(%rsp),%rbp
mov    %rbp,%rdx
call   1040 <FindSpan@plt>
movsd  (%r12,%rbx,1),%xmm0
mov    %rax,0x8(%rsp)
mov    %eax,%edi
mov    0x50(%rsp),%rsi
mov    %rbp,%rdx
mov    0x2a0(%rsp),%rcx
mov    0x290(%rsp),%r8
mov    0x298(%rsp),%r9
call   1030 <BasisFuns@plt>
movsd  0x8(%r12,%rbx,1),%xmm0
mov    0xac(%rsp),%edi
mov    0x48(%rsp),%rbp
mov    %ebp,%esi
mov    0x238(%rsp),%rdx
call   1040 <FindSpan@plt>
mov    %eax,%r15d
movsd  0x8(%r12,%rbx,1),%xmm0
mov    0x258(%rsp),%r12
mov    %eax,%edi
mov    %ebp,%esi
mov    0x238(%rsp),%rdx
mov    0x2b8(%rsp),%rbx
mov    %rbx,%rcx
mov    0x2a8(%rsp),%r8
mov    0x2b0(%rsp),%r9
call   1030 <BasisFuns@plt>
mov    %r15,%rdx
mov    0x2a0(%rsp),%r15
mov    0x50(%rsp),%r11
test   %ebp,%ebp
mov    0x38(%rsp),%rsi
mov    0x280(%rsp),%rax
mov    0x250(%rsp),%r8
js     1e30 <computeGNmatPnt2Pln+0xd10>
mov    0x80(%rsp),%r9
lea    (%r9,%r9,2),%rax
lea    0x2(%rax),%ecx
lea    0x1(%rax),%edi
mov    %edx,%ebp
sub    0x48(%rsp),%ebp
mov    %rbp,0x70(%rsp)
lea    (%r8,%rcx,8),%rbp
mov    %rbp,0x160(%rsp)
mov    0xf8(%rsp),%rbp
lea    0x0(%rbp,%rcx,8),%rcx
mov    %rcx,0x158(%rsp)
lea    (%r8,%rdi,8),%rcx
mov    %rcx,0x150(%rsp)
lea    0x0(%rbp,%rdi,8),%rcx
mov    %rcx,0x148(%rsp)
mov    %eax,%eax
lea    (%r8,%rax,8),%rcx
mov    %rcx,0x140(%rsp)
lea    0x0(%rbp,%rax,8),%rax
mov    %rax,0x138(%rsp)
lea    (%r9,%r9,2),%eax
mov    %rax,0x68(%rsp)
add    0xa4(%rsp),%edx
imul   0x24(%rsp),%edx
mov    0x100(%rsp),%rax
mov    %rdx,0x58(%rsp)
lea    (%rax,%rdx,1),%ecx
mov    0x8(%rsp),%r9
lea    0x0(,%r9,4),%eax
mov    %eax,0x18(%rsp)
lea    (%rcx,%r9,4),%ecx
xorpd  %xmm0,%xmm0
xor    %r8d,%r8d
xorpd  %xmm1,%xmm1
mov    0xf0(%rsp),%rdx
mov    0xe8(%rsp),%r10
jmp    1520 <computeGNmatPnt2Pln+0x400>
xchg   %ax,%ax
add    $0x1,%r8
add    %r10,%rcx
cmp    %r14,%r8
je     15c0 <computeGNmatPnt2Pln+0x4a0>
test   %r11d,%r11d
js     1510 <computeGNmatPnt2Pln+0x3f0>
movsd  (%rbx,%r8,8),%xmm2
xor    %esi,%esi
nopl   (%rax)
movapd %xmm1,%xmm3
movapd %xmm0,%xmm4
movsd  (%r15,%rsi,2),%xmm0
mulsd  %xmm2,%xmm0
lea    (%rcx,%rsi,1),%rdi
movslq %edi,%rdi
lea    0x1(%rdi),%ebp
movslq %ebp,%rbp
lea    0x2(%rdi),%ebx
movslq %ebx,%rbx
lea    0x3(%rdi),%eax
cltq
unpcklpd %xmm0,%xmm0
movsd  0x0(%r13,%rax,8),%xmm5
movhpd 0x0(%r13,%rbx,8),%xmm5
movsd  0x0(%r13,%rbp,8),%xmm6
movhpd 0x0(%r13,%rdi,8),%xmm6
mulpd  %xmm0,%xmm6
movapd %xmm6,%xmm1
addpd  %xmm3,%xmm1
mulpd  %xmm0,%xmm5
movapd %xmm5,%xmm0
addpd  %xmm4,%xmm0
add    $0x4,%rsi
cmp    %rsi,%rdx
jne    1530 <computeGNmatPnt2Pln+0x410>
addpd  %xmm4,%xmm5
addpd  %xmm3,%xmm6
movapd %xmm5,%xmm0
movapd %xmm6,%xmm1
mov    0x38(%rsp),%rsi
mov    0x2b8(%rsp),%rbx
jmp    1510 <computeGNmatPnt2Pln+0x3f0>
nopl   (%rax)
cmpl   $0x0,0x48(%rsp)
mov    %r9,%rdi
mov    0x280(%rsp),%rax
mov    0x58(%rsp),%r9
mov    0x68(%rsp),%rcx
js     1e30 <computeGNmatPnt2Pln+0xd10>
mov    %ecx,%edx
mov    %rdx,0xc8(%rsp)
lea    0x1(%rcx),%edx
mov    %rdx,0xc0(%rsp)
add    $0x2,%ecx
mov    %rcx,0x68(%rsp)
movapd %xmm1,%xmm9
unpckhpd %xmm1,%xmm9
movapd 0x19f1(%rip),%xmm2
xorpd  %xmm2,%xmm9
xorpd  %xmm2,%xmm1
movapd %xmm0,%xmm7
unpckhpd %xmm0,%xmm7
xorpd  %xmm2,%xmm7
movapd %xmm0,%xmm6
mulsd  %xmm0,%xmm6
sub    %r11d,%edi
mov    0x70(%rsp),%rcx
imul   0x28(%rsp),%ecx
add    %edi,%ecx
lea    0x3(,%rcx,4),%edx
movapd %xmm0,%xmm8
unpcklpd %xmm0,%xmm8
pshufd $0x44,%xmm9,%xmm2
movdqa %xmm2,0x1d0(%rsp)
pshufd $0x44,%xmm1,%xmm15
pshufd $0x44,%xmm7,%xmm10
movapd %xmm6,%xmm13
unpcklpd %xmm6,%xmm13
add    0xa0(%rsp),%r9d
add    0x18(%rsp),%r9d
xor    %edi,%edi
mov    %edx,0xb0(%rsp)
mov    %edx,%r10d
mov    %rcx,0x70(%rsp)
mov    %ecx,%ebx
xor    %ebp,%ebp
movaps %xmm9,0x1c0(%rsp)
movaps %xmm7,0x1b0(%rsp)
movapd %xmm6,0x1a0(%rsp)
movdqa %xmm10,0x1e0(%rsp)
jmp    1704 <computeGNmatPnt2Pln+0x5e4>
nopl   0x0(%rax)
mov    0x50(%rsp),%r11
mov    0x38(%rsp),%rsi
mov    0xd0(%rsp),%r14
mov    0x58(%rsp),%r9
mov    0x60(%rsp),%rdi
mov    0x44(%rsp),%r10d
mov    0xd8(%rsp),%rbx
add    $0x1,%rdi
add    0x28(%rsp),%ebx
mov    0x24(%rsp),%ecx
add    %ecx,%r9d
add    %ecx,%r10d
cmp    %r14,%rdi
je     1e30 <computeGNmatPnt2Pln+0xd10>
test   %r11d,%r11d
js     16e9 <computeGNmatPnt2Pln+0x5c9>
movslq %ebx,%r14
movslq %ebp,%rcx
cmpl   $0x2,0xb4(%rsp)
02
mov    %rdi,0x60(%rsp)
mov    %r9,0x58(%rsp)
mov    %r10d,0x44(%rsp)
mov    %rbx,0xd8(%rsp)
jae    1740 <computeGNmatPnt2Pln+0x620>
mov    %rcx,%rbp
xor    %ebx,%ebx
jmp    1c60 <computeGNmatPnt2Pln+0xb40>
nopw   0x0(%rax,%rax,1)
mov    %rcx,0x8(%rsp)
mov    0x24(%rsp),%eax
mov    %eax,%r12d
imul   %edi,%r12d
add    0xb0(%rsp),%r12d
mov    %ebp,%edx
add    %ebp,%edx
mov    0x188(%rsp),%rsi
lea    (%rsi,%rbp,2),%ecx
cmp    %edx,%ecx
setl   %al
or     0x11(%rsp),%al
lea    0x1(%rbp,%rbp,1),%ecx
lea    (%rsi,%rbp,2),%esi
add    $0x1,%esi
cmp    %ecx,%esi
setl   %r8b
mov    0x180(%rsp),%rsi
add    %r12d,%esi
cmp    %r12d,%esi
setl   %r11b
movslq %edx,%rsi
mov    0x288(%rsp),%rdx
lea    (%rdx,%rsi,8),%rdx
mov    0x178(%rsp),%rbp
lea    (%rdx,%rbp,1),%rdi
cmp    %rdx,%rdi
setb   %r9b
mov    0x170(%rsp),%rdi
lea    (%rdi,%rsi,8),%rdi
add    %rdi,%rbp
cmp    %rdi,%rbp
setb   %dil
test   %al,%al
jne    1c12 <computeGNmatPnt2Pln+0xaf2>
mov    0x168(%rsp),%rbp
shr    $0x20,%rbp
jne    1c12 <computeGNmatPnt2Pln+0xaf2>
or     0x11(%rsp),%r8b
jne    1c12 <computeGNmatPnt2Pln+0xaf2>
test   %rbp,%rbp
jne    1c12 <computeGNmatPnt2Pln+0xaf2>
or     0x13(%rsp),%r11b
mov    0x280(%rsp),%rax
jne    1c34 <computeGNmatPnt2Pln+0xb14>
test   %rbp,%rbp
jne    1c34 <computeGNmatPnt2Pln+0xb14>
mov    0x12(%rsp),%al
or     %al,%r9b
mov    0x90(%rsp),%rbp
jne    1c12 <computeGNmatPnt2Pln+0xaf2>
or     %al,%dil
jne    1c12 <computeGNmatPnt2Pln+0xaf2>
mov    %r14,0x18(%rsp)
mov    0x60(%rsp),%rbx
mov    %ebx,%eax
imul   0x28(%rsp),%eax
add    0x70(%rsp),%eax
cltq
shl    $0x3,%rax
lea    (%rax,%rax,2),%rdi
mov    0x258(%rsp),%rax
lea    (%rax,%rdi,1),%r10
mov    0x120(%rsp),%rax
add    %rdi,%rax
mov    %rax,0x78(%rsp)
lea    (%rdi,%rbp,1),%r8
add    0x128(%rsp),%rdi
mov    0x2b8(%rsp),%rax
lea    (%rax,%rbx,8),%rbp
lea    (%rax,%rbx,8),%rbx
add    $0x8,%rbx
mov    0x130(%rsp),%rax
lea    (%rax,%rsi,8),%r11
cmp    %r11,%rbp
setb   0x98(%rsp)
cmp    %rbx,%rdx
setb   0x30(%rsp)
cmp    %r15,%r11
seta   %r14b
cmp    0x110(%rsp),%rdx
setb   0x17(%rsp)
cmp    %r11,0x160(%rsp)
setb   %sil
cmp    0x158(%rsp),%rdx
setb   0x16(%rsp)
cmp    %r11,0x150(%rsp)
setb   %r9b
cmp    0x148(%rsp),%rdx
setb   0x15(%rsp)
cmp    %r11,0x140(%rsp)
setb   %bpl
cmp    0x138(%rsp),%rdx
setb   0x14(%rsp)
cmp    %r11,%r10
setb   %r10b
cmp    0x78(%rsp),%rdx
movslq %r12d,%r12
setb   %al
cmp    %r11,%r8
setb   %r8b
cmp    %rdi,%rdx
mov    0x118(%rsp),%rdi
lea    (%rdi,%r12,8),%rdi
setb   0x78(%rsp)
cmp    %rdi,%rdx
lea    0x0(%r13,%r12,8),%rdi
setb   %bl
cmp    %r11,%rdi
setb   %r11b
mov    0x30(%rsp),%dl
test   %dl,0x98(%rsp)
jne    1e0a <computeGNmatPnt2Pln+0xcea>
and    0x17(%rsp),%r14b
mov    0x258(%rsp),%r12
mov    0x18(%rsp),%r14
jne    1c45 <computeGNmatPnt2Pln+0xb25>
and    0x16(%rsp),%sil
jne    1c45 <computeGNmatPnt2Pln+0xb25>
and    0x15(%rsp),%r9b
jne    1c45 <computeGNmatPnt2Pln+0xb25>
and    0x14(%rsp),%bpl
jne    1c45 <computeGNmatPnt2Pln+0xb25>
and    %al,%r10b
jne    1c45 <computeGNmatPnt2Pln+0xb25>
and    0x78(%rsp),%r8b
mov    0x90(%rsp),%rbp
mov    0x280(%rsp),%rax
jne    1c2b <computeGNmatPnt2Pln+0xb0b>
mov    %ebx,%edx
and    %r11b,%dl
jne    1c2b <computeGNmatPnt2Pln+0xb0b>
mov    %rax,%rsi
lea    0x0(,%r14,8),%rax
lea    (%rax,%rax,2),%rdx
mov    0x8(%rsp),%rdi
mov    0x88(%rsp),%rax
add    %rdi,%rax
mov    %rax,0x30(%rsp)
lea    (%rsi,%rdi,4),%rax
mov    %rax,0x8(%rsp)
mov    0x44(%rsp),%eax
mov    %eax,%r11d
mov    0x58(%rsp),%r9
mov    %r12,%rax
xor    %r12d,%r12d
mov    0x250(%rsp),%rbx
mov    0xc8(%rsp),%rdi
mov    0xc0(%rsp),%r10
mov    0x60(%rsp),%r14
mov    0xd8(%rsp),%r8
nopw   0x0(%rax,%rax,1)
lea    (%r8,%r12,1),%esi
add    $0x1,%esi
movd   %esi,%xmm2
lea    (%r8,%r12,1),%esi
movd   %esi,%xmm3
punpckldq %xmm2,%xmm3
mov    0x8(%rsp),%rsi
movq   %xmm3,(%rsi,%r12,4)
mov    0x2b8(%rsp),%rsi
movsd  (%rsi,%r14,8),%xmm2
movslq %r11d,%r11
movslq %r9d,%r9
movupd (%r15,%r12,8),%xmm6
movsd  0x0(%r13,%r11,8),%xmm5
movhpd 0x0(%r13,%r9,8),%xmm5
unpcklpd %xmm2,%xmm2
movsd  (%rbx,%rdi,8),%xmm3
movapd %xmm8,%xmm14
movupd (%rax,%rdx,1),%xmm8
movupd 0x10(%rax,%rdx,1),%xmm9
unpcklpd %xmm3,%xmm3
movupd 0x20(%rax,%rdx,1),%xmm7
movapd %xmm9,%xmm4
movsd  %xmm8,%xmm4
shufpd $0x1,%xmm7,%xmm8
movsd  %xmm9,%xmm7
movupd 0x0(%rbp,%rdx,1),%xmm10
movupd 0x10(%rbp,%rdx,1),%xmm12
movupd 0x20(%rbp,%rdx,1),%xmm9
movapd %xmm12,%xmm11
movsd  %xmm10,%xmm11
shufpd $0x1,%xmm9,%xmm10
movsd  %xmm12,%xmm9
movsd  (%rbx,%r10,8),%xmm12
mulpd  %xmm2,%xmm6
unpcklpd %xmm12,%xmm12
mulpd  %xmm12,%xmm8
mulpd  %xmm3,%xmm4
addpd  %xmm8,%xmm4
movapd %xmm14,%xmm8
movapd 0x1d0(%rsp),%xmm14
mov    0x68(%rsp),%rsi
movsd  (%rbx,%rsi,8),%xmm2
unpcklpd %xmm2,%xmm2
mulpd  %xmm2,%xmm7
mulpd  %xmm6,%xmm5
mulpd  %xmm8,%xmm11
addpd  %xmm14,%xmm11
mulpd  %xmm8,%xmm10
addpd  %xmm4,%xmm7
addpd  %xmm15,%xmm10
mulpd  %xmm12,%xmm10
mulpd  %xmm3,%xmm11
mulpd  %xmm5,%xmm7
addpd  %xmm10,%xmm11
movapd 0x1e0(%rsp),%xmm10
mulpd  %xmm8,%xmm9
addpd  %xmm10,%xmm9
divpd  %xmm8,%xmm7
mulpd  %xmm2,%xmm9
addpd  %xmm11,%xmm9
mulpd  %xmm6,%xmm9
divpd  %xmm13,%xmm9
movslq %ecx,%rcx
movapd %xmm7,%xmm2
unpcklpd %xmm9,%xmm2
movhlps %xmm7,%xmm9
mov    0x288(%rsp),%rsi
movups %xmm9,0x8(%rsi,%rcx,8)
movupd %xmm2,-0x8(%rsi,%rcx,8)
add    $0x2,%r12
add    $0x30,%rbp
add    $0x4,%ecx
add    $0x30,%rax
add    $0x8,%r9d
add    $0x8,%r11d
cmp    %r12,0x88(%rsp)
jne    1a30 <computeGNmatPnt2Pln+0x910>
mov    0x88(%rsp),%rax
mov    %rax,%rbx
cmp    0xb8(%rsp),%rax
mov    0x280(%rsp),%rax
mov    0x258(%rsp),%r12
movapd 0x1c0(%rsp),%xmm9
movapd 0x1b0(%rsp),%xmm7
movapd 0x1a0(%rsp),%xmm6
mov    0x18(%rsp),%r14
mov    0x30(%rsp),%rbp
je     16c0 <computeGNmatPnt2Pln+0x5a0>
jmp    1c60 <computeGNmatPnt2Pln+0xb40>
mov    0x8(%rsp),%rbp
xor    %ebx,%ebx
mov    0x280(%rsp),%rax
mov    0x258(%rsp),%r12
jmp    1c60 <computeGNmatPnt2Pln+0xb40>
mov    0x8(%rsp),%rbp
xor    %ebx,%ebx
jmp    1c60 <computeGNmatPnt2Pln+0xb40>
mov    0x8(%rsp),%rbp
xor    %ebx,%ebx
mov    0x258(%rsp),%r12
jmp    1c60 <computeGNmatPnt2Pln+0xb40>
mov    0x8(%rsp),%rbp
xor    %ebx,%ebx
mov    0x280(%rsp),%rax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
lea    (%rax,%rbp,4),%rax
mov    %rax,0x8(%rsp)
mov    %rbp,0x30(%rsp)
lea    0x1(%rbp,%rbp,1),%ecx
lea    (%r14,%rbx,1),%eax
lea    0x3(,%rax,4),%edx
movslq %ebx,%rax
add    %r14,%rax
mov    %rax,0x18(%rsp)
add    %rbx,%r14
lea    (%r14,%r14,2),%rsi
add    $0x2,%rsi
lea    (%r15,%rbx,8),%r11
mov    0xb8(%rsp),%rax
sub    %rbx,%rax
mov    %rax,0x98(%rsp)
xor    %ebx,%ebx
mov    %r12,%r14
mov    0x250(%rsp),%r12
mov    0x68(%rsp),%r10
mov    0xc8(%rsp),%r9
mov    0xc0(%rsp),%r13
mov    0x60(%rsp),%r15
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x18(%rsp),%rbp
add    %ebx,%ebp
mov    0x8(%rsp),%r8
mov    %ebp,(%r8,%rbx,4)
mov    0x2b8(%rsp),%rdi
movsd  (%rdi,%r15,8),%xmm2
mulsd  (%r11,%rbx,8),%xmm2
movslq %edx,%rdx
mov    0xe0(%rsp),%rax
mulsd  (%rax,%rdx,8),%xmm2
movsd  (%r12,%r9,8),%xmm3
movsd  (%r12,%r13,8),%xmm4
mulsd  -0x8(%r14,%rsi,8),%xmm4
mulsd  -0x10(%r14,%rsi,8),%xmm3
movsd  (%r12,%r10,8),%xmm5
mulsd  (%r14,%rsi,8),%xmm5
addsd  %xmm4,%xmm3
addsd  %xmm3,%xmm5
mulsd  %xmm2,%xmm5
divsd  %xmm0,%xmm5
lea    -0x1(%rcx),%ebp
movslq %ebp,%rbp
mov    0x288(%rsp),%rax
movsd  %xmm5,(%rax,%rbp,8)
mov    0x90(%rsp),%rbp
movsd  (%rdi,%r15,8),%xmm2
mulsd  (%r11,%rbx,8),%xmm2
movsd  -0x10(%rbp,%rsi,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm9,%xmm3
movsd  -0x8(%rbp,%rsi,8),%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm1,%xmm4
mulsd  (%r12,%r13,8),%xmm4
mulsd  (%r12,%r9,8),%xmm3
addsd  %xmm4,%xmm3
movsd  0x0(%rbp,%rsi,8),%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm7,%xmm4
mulsd  (%r12,%r10,8),%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm2,%xmm4
divsd  %xmm6,%xmm4
movslq %ecx,%rcx
movsd  %xmm4,(%rax,%rcx,8)
add    $0x1,%rbx
add    $0x2,%ecx
add    $0x4,%edx
add    $0x3,%rsi
cmp    %rbx,0x98(%rsp)
jne    1ce0 <computeGNmatPnt2Pln+0xbc0>
mov    0x30(%rsp),%rbp
add    %rbx,%rbp
mov    %r14,%r12
mov    0x2a0(%rsp),%r15
mov    0x280(%rsp),%rax
mov    0xe0(%rsp),%r13
jmp    16c0 <computeGNmatPnt2Pln+0x5a0>
mov    0x8(%rsp),%rbp
xor    %ebx,%ebx
mov    0x280(%rsp),%rax
mov    0x258(%rsp),%r12
mov    0x18(%rsp),%r14
jmp    1c60 <computeGNmatPnt2Pln+0xb40>
nopl   0x0(%rax,%rax,1)
cmpl   $0x0,0xa8(%rsp)
jle    1330 <computeGNmatPnt2Pln+0x210>
xor    %r15d,%r15d
jmp    1e62 <computeGNmatPnt2Pln+0xd42>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r15
cmp    0x190(%rsp),%r15
je     1330 <computeGNmatPnt2Pln+0x210>
movslq (%rax,%r15,4),%rcx
mov    0x278(%rsp),%rdx
movslq (%rdx,%rcx,4),%rcx
cmp    %esi,%ecx
jge    1e50 <computeGNmatPnt2Pln+0xd30>
mov    0x260(%rsp),%r8
mov    0x80(%rsp),%rbx
movsd  (%r8,%rbx,8),%xmm0
lea    (%r15,%r15,1),%rdx
mov    %r15,%rsi
shl    $0x4,%rsi
lea    (%rcx,%rcx,1),%rbp
mov    %rcx,%rdi
shl    $0x4,%rdi
mov    0x288(%rsp),%rax
mulsd  (%rax,%rsi,1),%xmm0
mov    0x270(%rsp),%r9
addsd  (%r9,%rdi,1),%xmm0
movsd  %xmm0,(%r9,%rdi,1)
movsd  (%r8,%rbx,8),%xmm0
lea    0x1(%rbp),%esi
movslq %esi,%r14
mov    0x38(%rsp),%rsi
mulsd  0x8(%rax,%rdx,8),%xmm0
mov    0x280(%rsp),%rax
addsd  (%r9,%r14,8),%xmm0
lea    0x1(%r15,%r15,1),%r12
movsd  %xmm0,(%r9,%r14,8)
mov    %ebp,%r11d
mov    0x40(%rsp),%ebx
imul   %ebx,%r11d
imul   %ebx,%r14d
lea    (%r11,%rbp,1),%ebx
movslq %ebx,%rdi
mov    %rdi,0x8(%rsp)
lea    (%r11,%rbp,1),%ebx
add    $0x1,%ebx
movslq %ebx,%r9
add    %r14d,%ebp
or     $0x1,%ebp
movslq %ebp,%r10
xor    %ebx,%ebx
jmp    1fcd <computeGNmatPnt2Pln+0xead>
nopl   0x0(%rax,%rax,1)
mov    0x288(%rsp),%rax
movsd  (%rax,%rdx,8),%xmm0
lea    (%r11,%rbp,2),%esi
mulsd  (%rax,%rbx,4),%xmm0
movslq %esi,%rsi
mov    0x268(%rsp),%rdi
addsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdi,%rsi,8)
movsd  (%rax,%rdx,8),%xmm0
or     $0x1,%esi
mulsd  0x8(%rax,%rbx,4),%xmm0
movslq %esi,%rsi
addsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdi,%rsi,8)
movsd  (%rax,%r12,8),%xmm0
lea    (%r14,%rbp,2),%esi
mulsd  (%rax,%rbx,4),%xmm0
movslq %esi,%rsi
addsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdi,%rsi,8)
movsd  (%rax,%r12,8),%xmm0
or     $0x1,%esi
mulsd  0x8(%rax,%rbx,4),%xmm0
mov    0x280(%rsp),%rax
movslq %esi,%rsi
addsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdi,%rsi,8)
mov    0x38(%rsp),%rsi
add    $0x4,%rbx
cmp    %rbx,0x198(%rsp)
je     1e50 <computeGNmatPnt2Pln+0xd30>
movslq (%rax,%rbx,1),%rbp
mov    0x278(%rsp),%r8
mov    (%r8,%rbp,4),%ebp
cmp    %esi,%ebp
jge    1fbb <computeGNmatPnt2Pln+0xe9b>
cmp    %ecx,%ebp
jg     1f30 <computeGNmatPnt2Pln+0xe10>
cmp    %ecx,%ebp
jne    1fbb <computeGNmatPnt2Pln+0xe9b>
mov    0x288(%rsp),%rax
movsd  (%rax,%rdx,8),%xmm0
mulsd  (%rax,%rbx,4),%xmm0
mov    0x268(%rsp),%rdi
mov    0x8(%rsp),%rbp
addsd  (%rdi,%rbp,8),%xmm0
movsd  %xmm0,(%rdi,%rbp,8)
movsd  (%rax,%rdx,8),%xmm0
mulsd  0x8(%rax,%rbx,4),%xmm0
addsd  (%rdi,%r9,8),%xmm0
movsd  %xmm0,(%rdi,%r9,8)
movsd  (%rax,%r12,8),%xmm0
mulsd  0x8(%rax,%rbx,4),%xmm0
mov    0x280(%rsp),%rax
addsd  (%rdi,%r10,8),%xmm0
movsd  %xmm0,(%rdi,%r10,8)
jmp    1fbb <computeGNmatPnt2Pln+0xe9b>
mov    0x40(%rsp),%eax
cmp    $0x2,%eax
jl     2135 <computeGNmatPnt2Pln+0x1015>
mov    %eax,%r9d
mov    0x268(%rsp),%rax
lea    (%rax,%r9,8),%rcx
add    $0x8,%rcx
mov    %r9,%r14
shl    $0x4,%r14
lea    0x0(,%r9,8),%r8
lea    0x8(%rax),%rdi
mov    $0x1,%r11d
xor    %r10d,%r10d
mov    %rcx,%rsi
jmp    20b8 <computeGNmatPnt2Pln+0xf98>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r11
add    $0x1,%r10
add    $0x8,%rsi
add    %r8,%rcx
add    $0x8,%rdi
cmp    %r9,%r11
je     2135 <computeGNmatPnt2Pln+0x1015>
test   %r10,%r10
je     2100 <computeGNmatPnt2Pln+0xfe0>
mov    %r11,%rbp
and    $0xfffffffffffffffe,%rbp
xor    %ebx,%ebx
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,1),%xmm0
movsd  %xmm0,-0x8(%rcx,%rdx,8)
movsd  (%rsi,%rbx,1),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x2,%rdx
add    %r14,%rbx
cmp    %rdx,%rbp
jne    20d0 <computeGNmatPnt2Pln+0xfb0>
test   $0x1,%r11b
je     20a0 <computeGNmatPnt2Pln+0xf80>
jmp    2108 <computeGNmatPnt2Pln+0xfe8>
nopl   0x0(%rax)
xor    %edx,%edx
test   $0x1,%r11b
je     20a0 <computeGNmatPnt2Pln+0xf80>
mov    %r11,%rbx
imul   %r9,%rbx
mov    %rdx,%rax
imul   %r9,%rax
add    %r11,%rax
mov    0x268(%rsp),%rbp
movsd  0x0(%rbp,%rax,8),%xmm0
add    %rdx,%rbx
movsd  %xmm0,0x0(%rbp,%rbx,8)
jmp    20a0 <computeGNmatPnt2Pln+0xf80>
add    $0x1f8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @computeGNmatPnt2Pln(i32 %0, i32 %1, double* nocapture %2, double* nocapture %3, i32 %4, i32 %5, double* %6, double* %7, double* nocapture %8, i32 %9, double* nocapture %10, double* nocapture %11, double* nocapture %12, double* nocapture %13, double* nocapture %14, i32* nocapture %15, i32* nocapture %16, double* nocapture %17, double* %18, double* %19, double* %20, double* %21, double* %22, double* %23) {
%25 = shl nsw i32 %4, 2
%26 = add i32 %0, 1
%27 = add i32 %1, 1
%28 = mul i32 %27, %26
%29 = mul nsw i32 %5, %4
%30 = shl i32 %29, 1
%31 = add i32 %30, -8
%32 = icmp sgt i32 %9, 0
br i1 %32, label %33, label %84
33: ; preds = %24
%34 = icmp slt i32 %1, 0
%35 = icmp slt i32 %0, 0
%36 = icmp sgt i32 %28, 0
%37 = zext i32 %9 to i64
%38 = zext i32 %27 to i64
%39 = zext i32 %26 to i64
%40 = zext i32 %27 to i64
%41 = zext i32 %26 to i64
%42 = zext i32 %28 to i64
%43 = add nsw i64 %39, -1
%44 = shl i32 %4, 2
%45 = getelementptr double, double* %17, i64 1
%46 = shl nuw nsw i64 %39, 1
%47 = getelementptr double, double* %17, i64 %46
%48 = getelementptr double, double* %20, i64 %39
%49 = getelementptr double, double* %10, i64 1
%50 = getelementptr double, double* %10, i64 1
%51 = getelementptr double, double* %10, i64 1
%52 = mul nuw nsw i64 %39, 3
%53 = getelementptr double, double* %11, i64 %52
%54 = getelementptr double, double* %3, i64 %52
%55 = shl i32 %4, 2
%56 = shl nuw nsw i64 %39, 2
%57 = add nsw i64 %56, -3
%58 = getelementptr double, double* %2, i64 %57
%59 = icmp ult i32 %26, 2
%60 = trunc i64 %43 to i32
%61 = call { i32, i1 } @llvm.umul.with.overflow.i32(i32 %60, i32 2)
%62 = extractvalue { i32, i1 } %61, 0
%63 = extractvalue { i32, i1 } %61, 1
%64 = icmp ugt i64 %43, 4294967295
%65 = trunc i64 %43 to i32
%66 = call { i32, i1 } @llvm.umul.with.overflow.i32(i32 %65, i32 2)
%67 = extractvalue { i32, i1 } %66, 0
%68 = extractvalue { i32, i1 } %66, 1
%69 = icmp ugt i64 %43, 4294967295
%70 = trunc i64 %43 to i32
%71 = call { i32, i1 } @llvm.umul.with.overflow.i32(i32 %70, i32 4)
%72 = extractvalue { i32, i1 } %71, 0
%73 = extractvalue { i32, i1 } %71, 1
%74 = icmp ugt i64 %43, 4294967295
%75 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %43, i64 16)
%76 = extractvalue { i64, i1 } %75, 0
%77 = extractvalue { i64, i1 } %75, 1
%78 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %43, i64 16)
%79 = extractvalue { i64, i1 } %78, 0
%80 = extractvalue { i64, i1 } %78, 1
%81 = and i64 %39, 4294967294
%82 = getelementptr inbounds double, double* %17, i64 -1
%83 = icmp eq i64 %81, %39
br label %89
84: ; preds = %586, %24
%85 = icmp sgt i32 %31, 1
br i1 %85, label %86, label %630
86: ; preds = %84
%87 = zext i32 %31 to i64
%88 = zext i32 %31 to i64
br label %589
89: ; preds = %33, %586
%90 = phi i64 [ 0, %33 ], [ %587, %586 ]
%91 = mul nuw nsw i64 %90, 3
%92 = add nuw i64 %91, 2
%93 = and i64 %92, 4294967295
%94 = getelementptr double, double* %10, i64 %93
%95 = getelementptr double, double* %49, i64 %93
%96 = add nuw i64 %91, 1
%97 = and i64 %96, 4294967295
%98 = getelementptr double, double* %10, i64 %97
%99 = getelementptr double, double* %50, i64 %97
%100 = and i64 %91, 4294967295
%101 = getelementptr double, double* %10, i64 %100
%102 = getelementptr double, double* %51, i64 %100
%103 = shl nuw nsw i64 %90, 1
%104 = trunc i64 %90 to i32
%105 = mul nsw i32 %104, 3
%106 = getelementptr inbounds double, double* %8, i64 %103
%107 = load double, double* %106
%108 = tail call i32 @FindSpan(i32 %4, i32 %0, double %107, double* %6)
%109 = load double, double* %106
%110 = tail call i32 @BasisFuns(i32 %108, double %109, i32 %0, double* %6, double* %20, double* %18, double* %19)
%111 = sub nsw i32 %108, %0
%112 = or i64 %103, 1
%113 = getelementptr inbounds double, double* %8, i64 %112
%114 = load double, double* %113
%115 = tail call i32 @FindSpan(i32 %5, i32 %1, double %114, double* %7)
%116 = load double, double* %113
%117 = tail call i32 @BasisFuns(i32 %115, double %116, i32 %1, double* %7, double* %23, double* %21, double* %22)
%118 = sub i32 %115, %1
br i1 %34, label %204, label %155
119: ; preds = %200
br i1 %34, label %204, label %120
120: ; preds = %119
%121 = zext i32 %105 to i64
%122 = getelementptr inbounds double, double* %10, i64 %121
%123 = add nuw nsw i32 %105, 1
%124 = zext i32 %123 to i64
%125 = getelementptr inbounds double, double* %10, i64 %124
%126 = add nuw nsw i32 %105, 2
%127 = zext i32 %126 to i64
%128 = getelementptr inbounds double, double* %10, i64 %127
%129 = extractelement <4 x double> %201, i64 3
%130 = fneg double %129
%131 = extractelement <4 x double> %201, i64 2
%132 = fneg double %131
%133 = extractelement <4 x double> %201, i64 1
%134 = fneg double %133
%135 = extractelement <4 x double> %201, i64 0
%136 = fmul double %135, %135
%137 = sub i32 %108, %0
%138 = mul i32 %118, %4
%139 = add i32 %137, %138
%140 = shl i32 %139, 2
%141 = or i32 %140, 3
%142 = sub i32 %108, %0
%143 = add i32 %142, %138
%144 = shl i32 %143, 2
%145 = or i32 %144, 3
%146 = shufflevector <4 x double> %201, <4 x double> undef, <2 x i32> zeroinitializer
%147 = insertelement <2 x double> poison, double %130, i64 0
%148 = shufflevector <2 x double> %147, <2 x double> poison, <2 x i32> zeroinitializer
%149 = insertelement <2 x double> poison, double %132, i64 0
%150 = shufflevector <2 x double> %149, <2 x double> poison, <2 x i32> zeroinitializer
%151 = insertelement <2 x double> poison, double %134, i64 0
%152 = shufflevector <2 x double> %151, <2 x double> poison, <2 x i32> zeroinitializer
%153 = insertelement <2 x double> poison, double %136, i64 0
%154 = shufflevector <2 x double> %153, <2 x double> poison, <2 x i32> zeroinitializer
br label %207
155: ; preds = %89, %200
%156 = phi i64 [ %202, %200 ], [ 0, %89 ]
%157 = phi <4 x double> [ %201, %200 ], [ zeroinitializer, %89 ]
br i1 %35, label %200, label %158
158: ; preds = %155
%159 = getelementptr inbounds double, double* %23, i64 %156
%160 = load double, double* %159
%161 = trunc i64 %156 to i32
%162 = add i32 %118, %161
%163 = mul i32 %162, %25
br label %164
164: ; preds = %158, %164
%165 = phi i64 [ 0, %158 ], [ %196, %164 ]
%166 = phi <4 x double> [ %157, %158 ], [ %195, %164 ]
%167 = getelementptr inbounds double, double* %20, i64 %165
%168 = load double, double* %167
%169 = fmul double %160, %168
%170 = trunc i64 %165 to i32
%171 = add i32 %111, %170
%172 = shl i32 %171, 2
%173 = add i32 %172, %163
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %2, i64 %174
%176 = load double, double* %175
%177 = or i32 %173, 1
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds double, double* %2, i64 %178
%180 = load double, double* %179
%181 = or i32 %173, 2
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %2, i64 %182
%184 = load double, double* %183
%185 = or i32 %173, 3
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %2, i64 %186
%188 = load double, double* %187
%189 = insertelement <4 x double> poison, double %169, i64 0
%190 = shufflevector <4 x double> %189, <4 x double> poison, <4 x i32> zeroinitializer
%191 = insertelement <4 x double> poison, double %188, i64 0
%192 = insertelement <4 x double> %191, double %184, i64 1
%193 = insertelement <4 x double> %192, double %180, i64 2
%194 = insertelement <4 x double> %193, double %176, i64 3
%195 = call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %190, <4 x double> %194, <4 x double> %166)
%196 = add nuw nsw i64 %165, 1
%197 = icmp eq i64 %196, %39
br i1 %197, label %198, label %164
198: ; preds = %164
%199 = call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %190, <4 x double> %194, <4 x double> %166)
br label %200
200: ; preds = %198, %155
%201 = phi <4 x double> [ %157, %155 ], [ %199, %198 ]
%202 = add nuw nsw i64 %156, 1
%203 = icmp eq i64 %202, %38
br i1 %203, label %119, label %155
204: ; preds = %474, %89, %119
br i1 %36, label %205, label %586
205: ; preds = %204
%206 = getelementptr inbounds double, double* %12, i64 %90
br label %478
207: ; preds = %120, %474
%208 = phi i64 [ 0, %120 ], [ %476, %474 ]
%209 = phi i32 [ 0, %120 ], [ %475, %474 ]
%210 = getelementptr double, double* %23, i64 %208
%211 = add nuw i64 %208, 1
%212 = getelementptr double, double* %23, i64 %211
%213 = trunc i64 %208 to i32
%214 = mul i32 %213, %4
%215 = add i32 %214, %143
%216 = sext i32 %215 to i64
%217 = mul nsw i64 %216, 3
%218 = getelementptr double, double* %11, i64 %217
%219 = getelementptr double, double* %53, i64 %217
%220 = getelementptr double, double* %3, i64 %217
%221 = getelementptr double, double* %54, i64 %217
%222 = trunc i64 %208 to i32
%223 = mul i32 %55, %222
%224 = add i32 %223, %145
%225 = sext i32 %224 to i64
%226 = getelementptr double, double* %2, i64 %225
%227 = getelementptr double, double* %58, i64 %225
%228 = trunc i64 %208 to i32
%229 = mul i32 %44, %228
%230 = add i32 %229, %141
br i1 %35, label %474, label %231
231: ; preds = %207
%232 = trunc i64 %208 to i32
%233 = add nsw i32 %118, %232
%234 = mul nsw i32 %233, %4
%235 = add nsw i32 %234, %111
%236 = getelementptr inbounds double, double* %23, i64 %208
%237 = sext i32 %235 to i64
%238 = sext i32 %209 to i64
br i1 %59, label %403, label %239
239: ; preds = %231
%240 = shl i32 %209, 1
%241 = add i32 %240, %62
%242 = icmp slt i32 %241, %240
%243 = or i1 %242, %63
%244 = or i1 %243, %64
%245 = or i32 %240, 1
%246 = add i32 %245, %67
%247 = icmp slt i32 %246, %245
%248 = or i1 %247, %68
%249 = or i1 %248, %69
%250 = add i32 %230, %72
%251 = icmp slt i32 %250, %230
%252 = or i1 %251, %73
%253 = or i1 %252, %74
%254 = sext i32 %240 to i64
%255 = getelementptr double, double* %17, i64 %254
%256 = bitcast double* %255 to i8*
%257 = getelementptr i8, i8* %256, i64 %76
%258 = icmp ult i8* %257, %256
%259 = or i1 %258, %77
%260 = getelementptr double, double* %45, i64 %254
%261 = bitcast double* %260 to i8*
%262 = getelementptr i8, i8* %261, i64 %79
%263 = icmp ult i8* %262, %261
%264 = or i1 %263, %80
%265 = or i1 %244, %249
%266 = or i1 %265, %253
%267 = or i1 %266, %259
%268 = or i1 %267, %264
br i1 %268, label %403, label %269
269: ; preds = %239
%270 = shl i32 %209, 1
%271 = sext i32 %270 to i64
%272 = getelementptr double, double* %17, i64 %271
%273 = getelementptr double, double* %47, i64 %271
%274 = icmp ult double* %210, %273
%275 = icmp ult double* %272, %212
%276 = and i1 %274, %275
%277 = icmp ugt double* %273, %20
%278 = icmp ult double* %272, %48
%279 = and i1 %277, %278
%280 = or i1 %276, %279
%281 = icmp ult double* %94, %273
%282 = icmp ult double* %272, %95
%283 = and i1 %281, %282
%284 = or i1 %280, %283
%285 = icmp ult double* %98, %273
%286 = icmp ult double* %272, %99
%287 = and i1 %285, %286
%288 = or i1 %284, %287
%289 = icmp ult double* %101, %273
%290 = icmp ult double* %272, %102
%291 = and i1 %289, %290
%292 = or i1 %288, %291
%293 = icmp ult double* %218, %273
%294 = icmp ult double* %272, %219
%295 = and i1 %293, %294
%296 = or i1 %292, %295
%297 = icmp ult double* %220, %273
%298 = icmp ult double* %272, %221
%299 = and i1 %297, %298
%300 = or i1 %296, %299
%301 = icmp ult double* %272, %227
%302 = icmp ult double* %226, %273
%303 = and i1 %301, %302
%304 = or i1 %300, %303
br i1 %304, label %403, label %305
305: ; preds = %269
%306 = add nsw i64 %81, %238
br label %307
307: ; preds = %307, %305
%308 = phi i64 [ 0, %305 ], [ %400, %307 ]
%309 = add i64 %308, %238
%310 = or i64 %308, 1
%311 = trunc i64 %308 to i32
%312 = or i32 %311, 1
%313 = trunc i64 %308 to i32
%314 = add i32 %209, %313
%315 = add nsw i64 %308, %237
%316 = add nsw i64 %310, %237
%317 = insertelement <2 x i64> poison, i64 %315, i64 0
%318 = insertelement <2 x i64> %317, i64 %316, i64 1
%319 = add nsw i32 %235, %311
%320 = add nsw i32 %235, %312
%321 = getelementptr inbounds i32, i32* %16, i64 %309
%322 = trunc <2 x i64> %318 to <2 x i32>
%323 = bitcast i32* %321 to <2 x i32>*
store <2 x i32> %322, <2 x i32>* %323
%324 = load double, double* %236
%325 = insertelement <2 x double> poison, double %324, i64 0
%326 = shufflevector <2 x double> %325, <2 x double> poison, <2 x i32> zeroinitializer
%327 = getelementptr inbounds double, double* %20, i64 %308
%328 = bitcast double* %327 to <2 x double>*
%329 = load <2 x double>, <2 x double>* %328
%330 = fmul <2 x double> %326, %329
%331 = shl nsw i32 %319, 2
%332 = shl nsw i32 %320, 2
%333 = or i32 %331, 3
%334 = or i32 %332, 3
%335 = sext i32 %333 to i64
%336 = sext i32 %334 to i64
%337 = getelementptr inbounds double, double* %2, i64 %335
%338 = getelementptr inbounds double, double* %2, i64 %336
%339 = load double, double* %337
%340 = load double, double* %338
%341 = insertelement <2 x double> poison, double %339, i64 0
%342 = insertelement <2 x double> %341, double %340, i64 1
%343 = fmul <2 x double> %330, %342
%344 = load double, double* %122
%345 = insertelement <2 x double> poison, double %344, i64 0
%346 = shufflevector <2 x double> %345, <2 x double> poison, <2 x i32> zeroinitializer
%347 = mul nsw i64 %315, 3
%348 = getelementptr inbounds double, double* %11, i64 %347
%349 = bitcast double* %348 to <6 x double>*
%350 = load <6 x double>, <6 x double>* %349
%351 = shufflevector <6 x double> %350, <6 x double> poison, <2 x i32> <i32 0, i32 3>
%352 = shufflevector <6 x double> %350, <6 x double> poison, <2 x i32> <i32 1, i32 4>
%353 = shufflevector <6 x double> %350, <6 x double> poison, <2 x i32> <i32 2, i32 5>
%354 = load double, double* %125
%355 = insertelement <2 x double> poison, double %354, i64 0
%356 = shufflevector <2 x double> %355, <2 x double> poison, <2 x i32> zeroinitializer
%357 = fmul <2 x double> %356, %352
%358 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %346, <2 x double> %351, <2 x double> %357)
%359 = load double, double* %128
%360 = insertelement <2 x double> poison, double %359, i64 0
%361 = shufflevector <2 x double> %360, <2 x double> poison, <2 x i32> zeroinitializer
%362 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %361, <2 x double> %353, <2 x double> %358)
%363 = fmul <2 x double> %343, %362
%364 = fdiv <2 x double> %363, %146
%365 = shl nsw i32 %314, 1
%366 = load double, double* %236
%367 = insertelement <2 x double> poison, double %366, i64 0
%368 = shufflevector <2 x double> %367, <2 x double> poison, <2 x i32> zeroinitializer
%369 = bitcast double* %327 to <2 x double>*
%370 = load <2 x double>, <2 x double>* %369
%371 = fmul <2 x double> %368, %370
%372 = load double, double* %122
%373 = insertelement <2 x double> poison, double %372, i64 0
%374 = shufflevector <2 x double> %373, <2 x double> poison, <2 x i32> zeroinitializer
%375 = getelementptr inbounds double, double* %3, i64 %347
%376 = bitcast double* %375 to <6 x double>*
%377 = load <6 x double>, <6 x double>* %376
%378 = shufflevector <6 x double> %377, <6 x double> poison, <2 x i32> <i32 0, i32 3>
%379 = shufflevector <6 x double> %377, <6 x double> poison, <2 x i32> <i32 1, i32 4>
%380 = shufflevector <6 x double> %377, <6 x double> poison, <2 x i32> <i32 2, i32 5>
%381 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %146, <2 x double> %378, <2 x double> %148)
%382 = load double, double* %125
%383 = insertelement <2 x double> poison, double %382, i64 0
%384 = shufflevector <2 x double> %383, <2 x double> poison, <2 x i32> zeroinitializer
%385 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %146, <2 x double> %379, <2 x double> %150)
%386 = fmul <2 x double> %384, %385
%387 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %374, <2 x double> %381, <2 x double> %386)
%388 = load double, double* %128
%389 = insertelement <2 x double> poison, double %388, i64 0
%390 = shufflevector <2 x double> %389, <2 x double> poison, <2 x i32> zeroinitializer
%391 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %146, <2 x double> %380, <2 x double> %152)
%392 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %390, <2 x double> %391, <2 x double> %387)
%393 = fmul <2 x double> %371, %392
%394 = fdiv <2 x double> %393, %154
%395 = or i32 %365, 1
%396 = sext i32 %395 to i64
%397 = getelementptr inbounds double, double* %82, i64 %396
%398 = bitcast double* %397 to <4 x double>*
%399 = shufflevector <2 x double> %364, <2 x double> %394, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %399, <4 x double>* %398
%400 = add nuw i64 %308, 2
%401 = icmp eq i64 %400, %81
br i1 %401, label %402, label %307
402: ; preds = %307
br i1 %83, label %471, label %403
403: ; preds = %269, %239, %231, %402
%404 = phi i64 [ %238, %269 ], [ %238, %239 ], [ %238, %231 ], [ %306, %402 ]
%405 = phi i64 [ 0, %269 ], [ 0, %239 ], [ 0, %231 ], [ %81, %402 ]
br label %406
406: ; preds = %403, %406
%407 = phi i64 [ %468, %406 ], [ %404, %403 ]
%408 = phi i64 [ %469, %406 ], [ %405, %403 ]
%409 = trunc i64 %408 to i32
%410 = add nsw i64 %408, %237
%411 = add nsw i32 %235, %409
%412 = getelementptr inbounds i32, i32* %16, i64 %407
%413 = trunc i64 %410 to i32
store i32 %413, i32* %412
%414 = load double, double* %236
%415 = getelementptr inbounds double, double* %20, i64 %408
%416 = load double, double* %415
%417 = fmul double %414, %416
%418 = shl nsw i32 %411, 2
%419 = or i32 %418, 3
%420 = sext i32 %419 to i64
%421 = getelementptr inbounds double, double* %2, i64 %420
%422 = load double, double* %421
%423 = fmul double %417, %422
%424 = load double, double* %122
%425 = mul nsw i64 %410, 3
%426 = getelementptr inbounds double, double* %11, i64 %425
%427 = load double, double* %426
%428 = load double, double* %125
%429 = add nsw i64 %425, 1
%430 = getelementptr inbounds double, double* %11, i64 %429
%431 = load double, double* %430
%432 = fmul double %428, %431
%433 = tail call double @llvm.fmuladd.f64(double %424, double %427, double %432)
%434 = load double, double* %128
%435 = add nsw i64 %425, 2
%436 = getelementptr inbounds double, double* %11, i64 %435
%437 = load double, double* %436
%438 = tail call double @llvm.fmuladd.f64(double %434, double %437, double %433)
%439 = fmul double %423, %438
%440 = fdiv double %439, %135
%441 = trunc i64 %407 to i32
%442 = shl nsw i32 %441, 1
%443 = sext i32 %442 to i64
%444 = getelementptr inbounds double, double* %17, i64 %443
store double %440, double* %444
%445 = load double, double* %236
%446 = load double, double* %415
%447 = fmul double %445, %446
%448 = load double, double* %122
%449 = getelementptr inbounds double, double* %3, i64 %425
%450 = load double, double* %449
%451 = tail call double @llvm.fmuladd.f64(double %135, double %450, double %130)
%452 = load double, double* %125
%453 = getelementptr inbounds double, double* %3, i64 %429
%454 = load double, double* %453
%455 = tail call double @llvm.fmuladd.f64(double %135, double %454, double %132)
%456 = fmul double %452, %455
%457 = tail call double @llvm.fmuladd.f64(double %448, double %451, double %456)
%458 = load double, double* %128
%459 = getelementptr inbounds double, double* %3, i64 %435
%460 = load double, double* %459
%461 = tail call double @llvm.fmuladd.f64(double %135, double %460, double %134)
%462 = tail call double @llvm.fmuladd.f64(double %458, double %461, double %457)
%463 = fmul double %447, %462
%464 = fdiv double %463, %136
%465 = or i32 %442, 1
%466 = sext i32 %465 to i64
%467 = getelementptr inbounds double, double* %17, i64 %466
store double %464, double* %467
%468 = add nsw i64 %407, 1
%469 = add nuw nsw i64 %408, 1
%470 = icmp eq i64 %469, %41
br i1 %470, label %471, label %406
471: ; preds = %406, %402
%472 = phi i64 [ %306, %402 ], [ %468, %406 ]
%473 = trunc i64 %472 to i32
br label %474
474: ; preds = %471, %207
%475 = phi i32 [ %209, %207 ], [ %473, %471 ]
%476 = add nuw nsw i64 %208, 1
%477 = icmp eq i64 %476, %40
br i1 %477, label %204, label %207
478: ; preds = %205, %583
%479 = phi i64 [ 0, %205 ], [ %584, %583 ]
%480 = getelementptr inbounds i32, i32* %16, i64 %479
%481 = load i32, i32* %480
%482 = sext i32 %481 to i64
%483 = getelementptr inbounds i32, i32* %15, i64 %482
%484 = load i32, i32* %483
%485 = icmp slt i32 %484, %29
br i1 %485, label %486, label %583
486: ; preds = %478
%487 = load double, double* %206
%488 = shl nuw nsw i64 %479, 1
%489 = getelementptr inbounds double, double* %17, i64 %488
%490 = load double, double* %489
%491 = shl nsw i32 %484, 1
%492 = sext i32 %491 to i64
%493 = getelementptr inbounds double, double* %14, i64 %492
%494 = load double, double* %493
%495 = tail call double @llvm.fmuladd.f64(double %487, double %490, double %494)
store double %495, double* %493
%496 = load double, double* %206
%497 = or i64 %488, 1
%498 = getelementptr inbounds double, double* %17, i64 %497
%499 = load double, double* %498
%500 = or i32 %491, 1
%501 = sext i32 %500 to i64
%502 = getelementptr inbounds double, double* %14, i64 %501
%503 = load double, double* %502
%504 = tail call double @llvm.fmuladd.f64(double %496, double %499, double %503)
store double %504, double* %502
%505 = mul nsw i32 %491, %31
%506 = mul nsw i32 %500, %31
%507 = add nsw i32 %491, %505
%508 = sext i32 %507 to i64
%509 = getelementptr inbounds double, double* %13, i64 %508
%510 = or i32 %507, 1
%511 = sext i32 %510 to i64
%512 = getelementptr inbounds double, double* %13, i64 %511
%513 = add nsw i32 %491, %506
%514 = or i32 %513, 1
%515 = sext i32 %514 to i64
%516 = getelementptr inbounds double, double* %13, i64 %515
br label %517
517: ; preds = %486, %580
%518 = phi i64 [ 0, %486 ], [ %581, %580 ]
%519 = getelementptr inbounds i32, i32* %16, i64 %518
%520 = load i32, i32* %519
%521 = sext i32 %520 to i64
%522 = getelementptr inbounds i32, i32* %15, i64 %521
%523 = load i32, i32* %522
%524 = icmp slt i32 %523, %29
br i1 %524, label %525, label %580
525: ; preds = %517
%526 = icmp sgt i32 %523, %484
br i1 %526, label %527, label %561
527: ; preds = %525
%528 = load double, double* %489
%529 = shl nuw nsw i64 %518, 1
%530 = getelementptr inbounds double, double* %17, i64 %529
%531 = load double, double* %530
%532 = shl nsw i32 %523, 1
%533 = add nsw i32 %532, %505
%534 = sext i32 %533 to i64
%535 = getelementptr inbounds double, double* %13, i64 %534
%536 = load double, double* %535
%537 = tail call double @llvm.fmuladd.f64(double %528, double %531, double %536)
store double %537, double* %535
%538 = load double, double* %489
%539 = or i64 %529, 1
%540 = getelementptr inbounds double, double* %17, i64 %539
%541 = load double, double* %540
%542 = or i32 %533, 1
%543 = sext i32 %542 to i64
%544 = getelementptr inbounds double, double* %13, i64 %543
%545 = load double, double* %544
%546 = tail call double @llvm.fmuladd.f64(double %538, double %541, double %545)
store double %546, double* %544
%547 = load double, double* %498
%548 = load double, double* %530
%549 = add nsw i32 %532, %506
%550 = sext i32 %549 to i64
%551 = getelementptr inbounds double, double* %13, i64 %550
%552 = load double, double* %551
%553 = tail call double @llvm.fmuladd.f64(double %547, double %548, double %552)
store double %553, double* %551
%554 = load double, double* %498
%555 = load double, double* %540
%556 = or i32 %549, 1
%557 = sext i32 %556 to i64
%558 = getelementptr inbounds double, double* %13, i64 %557
%559 = load double, double* %558
%560 = tail call double @llvm.fmuladd.f64(double %554, double %555, double %559)
store double %560, double* %558
br label %580
561: ; preds = %525
%562 = icmp eq i32 %523, %484
br i1 %562, label %563, label %580
563: ; preds = %561
%564 = load double, double* %489
%565 = shl nuw nsw i64 %518, 1
%566 = getelementptr inbounds double, double* %17, i64 %565
%567 = load double, double* %566
%568 = load double, double* %509
%569 = tail call double @llvm.fmuladd.f64(double %564, double %567, double %568)
store double %569, double* %509
%570 = load double, double* %489
%571 = or i64 %565, 1
%572 = getelementptr inbounds double, double* %17, i64 %571
%573 = load double, double* %572
%574 = load double, double* %512
%575 = tail call double @llvm.fmuladd.f64(double %570, double %573, double %574)
store double %575, double* %512
%576 = load double, double* %498
%577 = load double, double* %572
%578 = load double, double* %516
%579 = tail call double @llvm.fmuladd.f64(double %576, double %577, double %578)
store double %579, double* %516
br label %580
580: ; preds = %517, %561, %563, %527
%581 = add nuw nsw i64 %518, 1
%582 = icmp eq i64 %581, %42
br i1 %582, label %583, label %517
583: ; preds = %580, %478
%584 = add nuw nsw i64 %479, 1
%585 = icmp eq i64 %584, %42
br i1 %585, label %586, label %478
586: ; preds = %583, %204
%587 = add nuw nsw i64 %90, 1
%588 = icmp eq i64 %587, %37
br i1 %588, label %84, label %89
589: ; preds = %86, %626
%590 = phi i64 [ 0, %86 ], [ %629, %626 ]
%591 = phi i64 [ 1, %86 ], [ %627, %626 ]
%592 = mul nsw i64 %591, %87
%593 = and i64 %591, 1
%594 = icmp eq i64 %590, 0
br i1 %594, label %616, label %595
595: ; preds = %589
%596 = and i64 %591, 9223372036854775806
br label %597
597: ; preds = %597, %595
%598 = phi i64 [ 0, %595 ], [ %613, %597 ]
%599 = phi i64 [ 0, %595 ], [ %614, %597 ]
%600 = mul nsw i64 %598, %87
%601 = add nuw nsw i64 %600, %591
%602 = getelementptr inbounds double, double* %13, i64 %601
%603 = load double, double* %602
%604 = add nuw nsw i64 %598, %592
%605 = getelementptr inbounds double, double* %13, i64 %604
store double %603, double* %605
%606 = or i64 %598, 1
%607 = mul nsw i64 %606, %87
%608 = add nuw nsw i64 %607, %591
%609 = getelementptr inbounds double, double* %13, i64 %608
%610 = load double, double* %609
%611 = add nuw nsw i64 %606, %592
%612 = getelementptr inbounds double, double* %13, i64 %611
store double %610, double* %612
%613 = add nuw nsw i64 %598, 2
%614 = add i64 %599, 2
%615 = icmp eq i64 %614, %596
br i1 %615, label %616, label %597
616: ; preds = %597, %589
%617 = phi i64 [ 0, %589 ], [ %613, %597 ]
%618 = icmp eq i64 %593, 0
br i1 %618, label %626, label %619
619: ; preds = %616
%620 = mul nsw i64 %617, %87
%621 = add nuw nsw i64 %620, %591
%622 = getelementptr inbounds double, double* %13, i64 %621
%623 = load double, double* %622
%624 = add nuw nsw i64 %617, %592
%625 = getelementptr inbounds double, double* %13, i64 %624
store double %623, double* %625
br label %626
626: ; preds = %616, %619
%627 = add nuw nsw i64 %591, 1
%628 = icmp eq i64 %627, %88
%629 = add i64 %590, 1
br i1 %628, label %630, label %589
630: ; preds = %626, %84
ret void
}