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

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double const fabs(double);
void palDmat(int n, double *a, double *y, double *d, int *jf, int *iw) {
  const double SFA = 1e-20;
  int k;
  double *aoff;
  *jf = 0;
  *d = 1.0;
  for (k = 0, aoff = a; k < n; k++, aoff += n) {
    int imx;
    double *aoff2 = aoff;
    double amx = fabs(aoff[k]);
    imx = k;
    if (k != n) {
      int i;
      double *apos2;
      for (i = k + 1, apos2 = aoff + n; i < n; i++, apos2 += n) {
        double t = fabs(apos2[k]);
        if (t > amx) {
          amx = t;
          imx = i;
          aoff2 = apos2;
        }
      }
    }
    if (amx < SFA) {
      *jf = -1;
    } else {
      if (imx != k) {
        double t;
        int j;
        for (j = 0; j < n; j++) {
          t = aoff[j];
          aoff[j] = aoff2[j];
          aoff2[j] = t;
        }
        t = y[k];
        y[k] = y[imx];
        y[imx] = t;
        *d = -*d;
      }
      iw[k] = imx;
      *d *= aoff[k];
      if (fabs(*d) < SFA) {
        *jf = -1;
      } else {
        double yk;
        double *apos2;
        int i, j;
        aoff[k] = 1.0 / aoff[k];
        for (j = 0; j < n; j++) {
          if (j != k) {
            aoff[j] *= aoff[k];
          }
        }
        yk = y[k] * aoff[k];
        y[k] = yk;
        for (i = 0, apos2 = a; i < n; i++, apos2 += n) {
          if (i != k) {
            for (j = 0; j < n; j++) {
              if (j != k) {
                apos2[j] -= apos2[k] * aoff[j];
              }
            }
            y[i] -= apos2[k] * yk;
          }
        }
        for (i = 0, apos2 = a; i < n; i++, apos2 += n) {
          if (i != k) {
            apos2[k] *= -aoff[k];
          }
        }
      }
    }
  }
  if (*jf != 0) {
    *d = 0.0;
  } else {
    for (k = n; k-- > 0;) {
      int ki = iw[k];
      if (k != ki) {
        int i;
        double *apos = a;
        for (i = 0; i < n; i++, apos += n) {
          double t = apos[k];
          apos[k] = apos[ki];
          apos[ki] = t;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<palDmat>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  0xece(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%rax
movl   $0x0,(%rax)
mov    -0x20(%rbp),%rax
movsd  0xeab(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1674 <palDmat+0x564>
mov    -0x48(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
movsd  %xmm0,-0x60(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
je     1239 <palDmat+0x129>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
mov    -0x48(%rbp),%rax
movslq -0x4(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1234 <palDmat+0x124>
mov    -0x70(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
ucomisd -0x60(%rbp),%xmm0
jbe    120c <palDmat+0xfc>
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x70(%rbp),%rax
mov    %rax,-0x58(%rbp)
jmp    1211 <palDmat+0x101>
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
mov    -0x4(%rbp),%ecx
mov    -0x70(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
jmp    11c1 <palDmat+0xb1>
jmp    1239 <palDmat+0x129>
movsd  0xdc7(%rip),%xmm0
ucomisd -0x60(%rbp),%xmm0
jbe    125b <palDmat+0x14b>
mov    -0x28(%rbp),%rax
movl   $0xffffffff,(%rax)
jmp    164c <palDmat+0x53c>
mov    -0x4c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
je     1343 <palDmat+0x233>
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12de <palDmat+0x1ce>
mov    -0x48(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x58(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    1271 <palDmat+0x161>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x4c(%rbp),%edx
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
call   1030 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0xc86(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    139b <palDmat+0x28b>
mov    -0x28(%rbp),%rax
movl   $0xffffffff,(%rax)
jmp    1647 <palDmat+0x537>
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xc55(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1420 <palDmat+0x310>
mov    -0xa0(%rbp),%eax
cmp    -0x3c(%rbp),%eax
je     1407 <palDmat+0x2f7>
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0xa0(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    140c <palDmat+0x2fc>
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa0(%rbp)
jmp    13c7 <palDmat+0x2b7>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x9c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x98(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15a2 <palDmat+0x492>
mov    -0x9c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
je     156e <palDmat+0x45e>
movl   $0x0,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1522 <palDmat+0x412>
mov    -0xa0(%rbp),%eax
cmp    -0x3c(%rbp),%eax
je     1509 <palDmat+0x3f9>
mov    -0x98(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0xa0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x98(%rbp),%rax
movslq -0xa0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    150e <palDmat+0x3fe>
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa0(%rbp)
jmp    1494 <palDmat+0x384>
mov    -0x98(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x90(%rbp),%xmm2
ff
mov    -0x18(%rbp),%rax
movslq -0x9c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1573 <palDmat+0x463>
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
mov    -0x4(%rbp),%ecx
mov    -0x98(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x98(%rbp)
jmp    146c <palDmat+0x35c>
movl   $0x0,-0x9c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x98(%rbp)
mov    -0x9c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1642 <palDmat+0x532>
mov    -0x9c(%rbp),%eax
cmp    -0x3c(%rbp),%eax
je     160e <palDmat+0x4fe>
mov    -0x48(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x98(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1613 <palDmat+0x503>
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x9c(%rbp)
mov    -0x4(%rbp),%ecx
mov    -0x98(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x98(%rbp)
jmp    15b7 <palDmat+0x4a7>
jmp    1647 <palDmat+0x537>
jmp    164c <palDmat+0x53c>
jmp    1651 <palDmat+0x541>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4(%rbp),%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
jmp    1168 <palDmat+0x58>
mov    -0x28(%rbp),%rax
cmpl   $0x0,(%rax)
je     1691 <palDmat+0x581>
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
jmp    1783 <palDmat+0x673>
mov    -0x4(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x3c(%rbp)
cmp    $0x0,%eax
jle    177e <palDmat+0x66e>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xa4(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xa4(%rbp),%eax
je     1779 <palDmat+0x669>
mov    -0x10(%rbp),%rax
mov    %rax,-0xb0(%rbp)
movl   $0x0,-0xa8(%rbp)
mov    -0xa8(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1774 <palDmat+0x664>
mov    -0xb0(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0xb0(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xb0(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
mov    -0xb0(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xa8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa8(%rbp)
mov    -0x4(%rbp),%ecx
mov    -0xb0(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0xb0(%rbp)
jmp    16e0 <palDmat+0x5d0>
jmp    1779 <palDmat+0x669>
jmp    1697 <palDmat+0x587>
jmp    1783 <palDmat+0x673>
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<palDmat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x48(%rsp)
mov    %rsi,-0x68(%rsp)
mov    %r8,-0x40(%rsp)
movl   $0x0,(%r8)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rcx)
test   %edi,%edi
jle    1733 <palDmat+0x633>
mov    %edi,%r15d
lea    -0x1(%rdi),%eax
mov    %eax,-0x6c(%rsp)
lea    -0x2(%rdi),%eax
mov    %eax,-0x4c(%rsp)
mov    %r15d,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x20(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x28(%rsp)
mov    %rax,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
lea    -0x1(%r15),%rax
mov    %rax,-0x58(%rsp)
mov    %rsi,%rax
mov    %rsi,-0x30(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x38(%rsp)
mov    %r15d,%r14d
and    $0xfffffffe,%r14d
mov    %edi,%eax
and    $0xfffffffe,%eax
mov    %eax,-0x5c(%rsp)
lea    (%r15,%r15,1),%rax
mov    %rax,-0x18(%rsp)
mov    %r15,%r8
shl    $0x4,%r8
lea    0x0(,%r15,8),%r12
xor    %r13d,%r13d
movapd 0xe4e(%rip),%xmm0
movapd 0xe46(%rip),%xmm1
movsd  0xe5e(%rip),%xmm2
movsd  0xe5d(%rip),%xmm8
movapd 0xe3d(%rip),%xmm4
mov    -0x68(%rsp),%r9
mov    %rcx,-0x10(%rsp)
xor    %ebp,%ebp
jmp    1210 <palDmat+0x110>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    -0x40(%rsp),%rax
movl   $0xffffffff,(%rax)
add    %r12,%r9
add    $0x8,%r13
mov    -0x8(%rsp),%rbp
cmp    %r15,%rbp
je     1722 <palDmat+0x622>
movsd  (%r9,%rbp,8),%xmm5
movapd %xmm5,%xmm7
andpd  %xmm0,%xmm7
lea    0x1(%rbp),%rcx
cmp    %r15,%rcx
mov    %rcx,-0x8(%rsp)
jae    1250 <palDmat+0x150>
mov    -0x6c(%rsp),%eax
sub    %ebp,%eax
test   $0x1,%al
jne    1260 <palDmat+0x160>
mov    %ecx,%eax
mov    %r9,%r11
mov    %ebp,%r10d
mov    %r9,%rsi
cmp    %ebp,-0x4c(%rsp)
jne    1291 <palDmat+0x191>
jmp    12f4 <palDmat+0x1f4>
nopl   0x0(%rax)
movapd %xmm7,%xmm6
mov    %r9,%rsi
mov    %ebp,%r10d
jmp    12f4 <palDmat+0x1f4>
nop
lea    (%r9,%r15,8),%r11
movsd  (%r11,%rbp,8),%xmm6
andpd  %xmm0,%xmm6
ucomisd %xmm7,%xmm6
maxsd  %xmm7,%xmm6
mov    %r9,%rsi
cmova  %r11,%rsi
mov    %ebp,%r10d
cmova  %ecx,%r10d
lea    0x1(%rcx),%eax
movapd %xmm6,%xmm7
cmp    %ebp,-0x4c(%rsp)
je     12f4 <palDmat+0x1f4>
mov    -0x18(%rsp),%rcx
lea    (%r11,%rcx,8),%rcx
lea    (%r11,%r15,8),%r11
movapd %xmm7,%xmm6
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r11,%r13,1),%xmm7
andpd  %xmm1,%xmm7
ucomisd %xmm6,%xmm7
cmova  %r11,%rsi
maxsd  %xmm6,%xmm7
cmova  %eax,%r10d
movsd  (%rcx,%r13,1),%xmm6
andpd  %xmm1,%xmm6
ucomisd %xmm7,%xmm6
lea    0x1(%rax),%ebx
cmova  %rcx,%rsi
maxsd  %xmm7,%xmm6
cmova  %ebx,%r10d
add    %r8,%rcx
add    %r8,%r11
add    $0x2,%eax
cmp    %edi,%eax
jne    12b0 <palDmat+0x1b0>
ucomisd %xmm6,%xmm2
mov    -0x10(%rsp),%r11
ja     11f0 <palDmat+0xf0>
mov    %r10d,%eax
cmp    %rax,%rbp
jne    1320 <palDmat+0x220>
movsd  (%r11),%xmm6
jmp    14bb <palDmat+0x3bb>
cs nopw 0x0(%rax,%rax,1)
nop
cmp    $0x4,%edi
jb     134a <palDmat+0x24a>
lea    (%rsi,%r15,8),%rax
cmp    %rax,%r9
jae    13be <palDmat+0x2be>
mov    %rbp,%rax
imul   %r15,%rax
add    %r15,%rax
mov    -0x68(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
cmp    %rax,%rsi
jae    13be <palDmat+0x2be>
xor    %ecx,%ecx
mov    %rcx,%rax
test   $0x1,%r15b
je     1372 <palDmat+0x272>
movsd  (%r9,%rcx,8),%xmm3
movsd  (%rsi,%rcx,8),%xmm5
movsd  %xmm5,(%r9,%rcx,8)
movsd  %xmm3,(%rsi,%rcx,8)
mov    %rcx,%rax
or     $0x1,%rax
cmp    %rcx,-0x58(%rsp)
je     1490 <palDmat+0x390>
nopl   (%rax)
movsd  (%r9,%rax,8),%xmm3
movsd  (%rsi,%rax,8),%xmm5
movsd  %xmm5,(%r9,%rax,8)
movsd  %xmm3,(%rsi,%rax,8)
movsd  0x8(%r9,%rax,8),%xmm3
movsd  0x8(%rsi,%rax,8),%xmm5
movsd  %xmm5,0x8(%r9,%rax,8)
movsd  %xmm3,0x8(%rsi,%rax,8)
add    $0x2,%rax
cmp    %rax,%r15
jne    1380 <palDmat+0x280>
jmp    1490 <palDmat+0x390>
cmpq   $0x0,-0x28(%rsp)
je     1710 <palDmat+0x610>
mov    -0x38(%rsp),%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups (%r9,%rax,8),%xmm5
movups 0x10(%r9,%rax,8),%xmm6
movups (%rsi,%rax,8),%xmm7
movups 0x10(%rsi,%rax,8),%xmm3
movups %xmm7,(%r9,%rax,8)
movups %xmm3,0x10(%r9,%rax,8)
movups %xmm5,(%rsi,%rax,8)
movups %xmm6,0x10(%rsi,%rax,8)
movupd 0x20(%r9,%rax,8),%xmm3
movupd 0x30(%r9,%rax,8),%xmm5
movupd 0x20(%rsi,%rax,8),%xmm6
movupd 0x30(%rsi,%rax,8),%xmm7
movupd %xmm6,0x20(%r9,%rax,8)
movupd %xmm7,0x30(%r9,%rax,8)
movupd %xmm3,0x20(%rsi,%rax,8)
movupd %xmm5,0x30(%rsi,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    13e0 <palDmat+0x2e0>
testb  $0x1,-0x30(%rsp)
je     147d <palDmat+0x37d>
movupd (%r9,%rax,8),%xmm3
movupd 0x10(%r9,%rax,8),%xmm5
movupd (%rsi,%rax,8),%xmm6
movupd 0x10(%rsi,%rax,8),%xmm7
movupd %xmm6,(%r9,%rax,8)
movupd %xmm7,0x10(%r9,%rax,8)
movupd %xmm3,(%rsi,%rax,8)
movupd %xmm5,0x10(%rsi,%rax,8)
mov    -0x20(%rsp),%rax
mov    %rax,%rcx
cmp    %r15,%rax
jne    134c <palDmat+0x24c>
xchg   %ax,%ax
movsd  (%rdx,%rbp,8),%xmm3
movslq %r10d,%rax
movsd  (%rdx,%rax,8),%xmm5
movsd  %xmm5,(%rdx,%rbp,8)
movsd  %xmm3,(%rdx,%rax,8)
movsd  (%r11),%xmm6
xorpd  %xmm4,%xmm6
movlpd %xmm6,(%r11)
movsd  (%r9,%rbp,8),%xmm5
mov    -0x48(%rsp),%rax
mov    %r10d,(%rax,%rbp,4)
mulsd  %xmm5,%xmm6
movsd  %xmm6,(%r11)
andpd  %xmm0,%xmm6
ucomisd %xmm6,%xmm2
ja     11f0 <palDmat+0xf0>
movapd %xmm8,%xmm3
divsd  (%r9,%rbp,8),%xmm3
movsd  %xmm3,(%r9,%rbp,8)
cmpq   $0x0,-0x58(%rsp)
je     1548 <palDmat+0x448>
xor    %ecx,%ecx
jmp    1509 <palDmat+0x409>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rcx
cmp    %r14,%rcx
je     1542 <palDmat+0x442>
mov    %rcx,%rax
cmp    %rcx,%rbp
je     1523 <palDmat+0x423>
movsd  (%r9,%rbp,8),%xmm3
mulsd  (%r9,%rax,8),%xmm3
movsd  %xmm3,(%r9,%rax,8)
lea    0x1(%rax),%rcx
cmp    %rbp,%rcx
je     1500 <palDmat+0x400>
movsd  (%r9,%rbp,8),%xmm3
mulsd  0x8(%r9,%rax,8),%xmm3
movsd  %xmm3,0x8(%r9,%rax,8)
jmp    1500 <palDmat+0x400>
add    $0x2,%rax
jmp    154a <palDmat+0x44a>
xor    %eax,%eax
mov    -0x5c(%rsp),%r10d
test   $0x1,%r15b
je     156c <palDmat+0x46c>
cmp    %rbp,%rax
je     156c <palDmat+0x46c>
movsd  (%r9,%rbp,8),%xmm3
mulsd  (%r9,%rax,8),%xmm3
movsd  %xmm3,(%r9,%rax,8)
movsd  (%rdx,%rbp,8),%xmm5
mulsd  (%r9,%rbp,8),%xmm5
movsd  %xmm5,(%rdx,%rbp,8)
mov    -0x68(%rsp),%rsi
xor    %eax,%eax
jmp    15b7 <palDmat+0x4b7>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx,%rax,8),%xmm3
movsd  (%rsi,%rbp,8),%xmm6
mulsd  %xmm5,%xmm6
subsd  %xmm6,%xmm3
movsd  %xmm3,(%rdx,%rax,8)
add    $0x1,%rax
add    %r12,%rsi
cmp    %r15,%rax
je     1680 <palDmat+0x580>
cmp    %rbp,%rax
je     15a7 <palDmat+0x4a7>
cmpq   $0x0,-0x58(%rsp)
je     163f <palDmat+0x53f>
xor    %ecx,%ecx
jmp    15d9 <palDmat+0x4d9>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rcx
cmp    %r14,%rcx
je     1630 <palDmat+0x530>
mov    %rcx,%rbx
cmp    %rcx,%rbp
je     15fa <palDmat+0x4fa>
movsd  (%rsi,%rbp,8),%xmm3
movsd  (%rsi,%rbx,8),%xmm6
mulsd  (%r9,%rbx,8),%xmm3
subsd  %xmm3,%xmm6
movsd  %xmm6,(%rsi,%rbx,8)
lea    0x1(%rbx),%rcx
cmp    %rbp,%rcx
je     15d0 <palDmat+0x4d0>
movsd  (%rsi,%rbp,8),%xmm3
movsd  0x8(%rsi,%rbx,8),%xmm6
mulsd  0x8(%r9,%rbx,8),%xmm3
subsd  %xmm3,%xmm6
movsd  %xmm6,0x8(%rsi,%rbx,8)
jmp    15d0 <palDmat+0x4d0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x2,%rbx
test   $0x1,%r15b
jne    164b <palDmat+0x54b>
jmp    1590 <palDmat+0x490>
xor    %ebx,%ebx
test   $0x1,%r15b
je     1590 <palDmat+0x490>
cmp    %rbp,%rbx
je     1590 <palDmat+0x490>
movsd  (%rsi,%rbp,8),%xmm3
movsd  (%rsi,%rbx,8),%xmm6
mulsd  (%r9,%rbx,8),%xmm3
subsd  %xmm3,%xmm6
movsd  %xmm6,(%rsi,%rbx,8)
jmp    1590 <palDmat+0x490>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cmpl   $0x0,-0x6c(%rsp)
je     16db <palDmat+0x5db>
xor    %eax,%eax
mov    -0x68(%rsp),%rcx
jmp    169c <palDmat+0x59c>
add    $0x2,%eax
lea    (%rcx,%r15,8),%rcx
cmp    %r10d,%eax
je     16e2 <palDmat+0x5e2>
mov    %eax,%esi
cmp    %rbp,%rsi
je     16b7 <palDmat+0x5b7>
movsd  (%r9,%rbp,8),%xmm3
xorpd  %xmm4,%xmm3
mulsd  (%rcx,%rbp,8),%xmm3
movsd  %xmm3,(%rcx,%rbp,8)
mov    %eax,%esi
or     $0x1,%esi
lea    (%rcx,%r15,8),%rcx
cmp    %rbp,%rsi
je     1690 <palDmat+0x590>
movsd  (%r9,%rbp,8),%xmm3
xorpd  %xmm4,%xmm3
mulsd  (%rcx,%rbp,8),%xmm3
movsd  %xmm3,(%rcx,%rbp,8)
jmp    1690 <palDmat+0x590>
xor    %eax,%eax
mov    -0x68(%rsp),%rcx
test   $0x1,%dil
je     11fb <palDmat+0xfb>
mov    %eax,%eax
cmp    %rax,%rbp
je     11fb <palDmat+0xfb>
movsd  (%r9,%rbp,8),%xmm3
xorpd  %xmm4,%xmm3
mulsd  (%rcx,%rbp,8),%xmm3
movsd  %xmm3,(%rcx,%rbp,8)
jmp    11fb <palDmat+0xfb>
xor    %eax,%eax
testb  $0x1,-0x30(%rsp)
jne    144d <palDmat+0x34d>
jmp    147d <palDmat+0x37d>
mov    -0x40(%rsp),%rax
cmpl   $0x0,(%rax)
je     173e <palDmat+0x63e>
movq   $0x0,(%r11)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
test   %edi,%edi
jle    1733 <palDmat+0x633>
mov    %r15,%r12
shl    $0x4,%r12
lea    0x0(,%r15,8),%r8
mov    -0x6c(%rsp),%r14d
jmp    176a <palDmat+0x66a>
nopl   0x0(%rax,%rax,1)
add    $0xffffffff,%r14d
cmp    $0x1,%r10
jle    1733 <palDmat+0x633>
mov    %r15,%r10
add    $0xffffffffffffffff,%r15
mov    %r15d,%r9d
mov    -0x48(%rsp),%rax
movslq (%rax,%r9,4),%r11
cmp    %r15d,%r11d
je     1760 <palDmat+0x660>
mov    -0x68(%rsp),%rcx
cmpl   $0x0,-0x6c(%rsp)
je     17e0 <palDmat+0x6e0>
mov    %r14d,%ebp
shl    $0x3,%rbp
lea    0x0(,%r11,8),%rsi
lea    (%r8,%r11,8),%rbx
lea    (%r8,%rbp,1),%rdx
mov    -0x5c(%rsp),%eax
mov    -0x68(%rsp),%rcx
xchg   %ax,%ax
movsd  (%rcx,%rbp,1),%xmm0
movsd  (%rcx,%rsi,1),%xmm1
movsd  %xmm1,(%rcx,%rbp,1)
movsd  %xmm0,(%rcx,%rsi,1)
movsd  (%rcx,%rdx,1),%xmm0
movsd  (%rcx,%rbx,1),%xmm1
movsd  %xmm1,(%rcx,%rdx,1)
movsd  %xmm0,(%rcx,%rbx,1)
add    %r12,%rcx
add    $0xfffffffe,%eax
jne    17b0 <palDmat+0x6b0>
test   $0x1,%dil
je     1760 <palDmat+0x660>
movsd  (%rcx,%r9,8),%xmm0
movsd  (%rcx,%r11,8),%xmm1
movsd  %xmm1,(%rcx,%r9,8)
movsd  %xmm0,(%rcx,%r11,8)
jmp    1760 <palDmat+0x660>

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

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