void tricntr(double *pxn, int kig1s, int kig1e, int kig2s, int kig2e, int knd,
             int kjd, int kni) {
  int id1, id2, id3, ioffset;
  double r1, r2, r3;
  int j, mi1, mi2;
  double zxnorm;
  id1 = kig1e - kig1s + 1;
  id2 = id1 * (kig2e - kig2s + 1);
  id3 = id2 * 3;
  ioffset = -(id1 + id2 + id3);
  for (j = 1; j <= 2; ++j) {
    mi1 = j * kni / 3;
    mi2 = (j - 1) * kni + 1;
    pxn[mi1 + id1 * (mi1 + 1) + id2 * 1 + id3 * kjd + ioffset] =
        pxn[mi2 - 1 + id1 * (mi2) + id2 * 1 + id3 * kjd + ioffset] +
        pxn[kni + id1 * 1 + id2 * 1 + id3 * kjd + ioffset] +
        pxn[0 + id1 * (kni + 1) + id2 * 1 + id3 * kjd + ioffset];
    pxn[mi1 + id1 * (mi1 + 1) + id2 * 2 + id3 * kjd + ioffset] =
        pxn[mi2 - 1 + id1 * (mi2) + id2 * 2 + id3 * kjd + ioffset] +
        pxn[kni + id1 * 1 + id2 * 2 + id3 * kjd + ioffset] +
        pxn[0 + id1 * (kni + 1) + id2 * 2 + id3 * kjd + ioffset];
    pxn[mi1 + id1 * (mi1 + 1) + id2 * 3 + id3 * kjd + ioffset] =
        pxn[mi2 - 1 + id1 * (mi2) + id2 * 3 + id3 * kjd + ioffset] +
        pxn[kni + id1 * 1 + id2 * 3 + id3 * kjd + ioffset] +
        pxn[0 + id1 * (kni + 1) + id2 * 3 + id3 * kjd + ioffset];
    r1 = pxn[mi1 + id1 * (mi1 + 1) + id2 * 1 + id3 * kjd + ioffset];
    r2 = pxn[mi1 + id1 * (mi1 + 1) + id2 * 2 + id3 * kjd + ioffset];
    r3 = pxn[mi1 + id1 * (mi1 + 1) + id2 * 3 + id3 * kjd + ioffset];
    zxnorm = 1.0 / sqrt(r1 * r1 + r2 * r2 + r3 * r3);
    pxn[mi1 + id1 * (mi1 + 1) + id2 * 1 + id3 * kjd + ioffset] =
        zxnorm * pxn[mi1 + id1 * (mi1 + 1) + id2 * 1 + id3 * kjd + ioffset];
    pxn[mi1 + id1 * (mi1 + 1) + id2 * 2 + id3 * kjd + ioffset] =
        zxnorm * pxn[mi1 + id1 * (mi1 + 1) + id2 * 2 + id3 * kjd + ioffset];
    pxn[mi1 + id1 * (mi1 + 1) + id2 * 3 + id3 * kjd + ioffset] =
        zxnorm * pxn[mi1 + id1 * (mi1 + 1) + id2 * 3 + id3 * kjd + ioffset];
  }
  return;
}


x86 assembly (clang, -O0):
<tricntr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x24(%rbp)
imul   $0x3,-0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x20(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x2c(%rbp)
movl   $0x1,-0x4c(%rbp)
cmpl   $0x2,-0x4c(%rbp)
jg     15e7 <tricntr+0x4d7>
mov    -0x4c(%rbp),%eax
imul   0x18(%rbp),%eax
mov    $0x3,%ecx
cltd
idiv   %ecx
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
sub    $0x1,%eax
imul   0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x20(%rbp),%edx
imul   -0x54(%rbp),%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x20(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    0x18(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
add    $0x0,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x20(%rbp),%edx
imul   -0x54(%rbp),%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x20(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    0x18(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
add    $0x0,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x20(%rbp),%edx
imul   -0x54(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x20(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    0x18(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
add    $0x0,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xb61(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x20(%rbp),%edx
mov    -0x50(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
imul   $0x3,-0x24(%rbp),%edx
add    %edx,%ecx
mov    -0x28(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    116f <tricntr+0x5f>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tricntr>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rdi,%rbx
movslq 0x98(%rsp),%r11
sub    %esi,%edx
mov    %edx,%esi
not    %esi
lea    0x1(%rdx),%r9d
sub    %ecx,%r8d
add    $0x1,%r8d
imul   %r9d,%r8d
lea    (%r8,%r8,2),%eax
mov    %eax,0x8(%rsp)
lea    0x0(,%r8,4),%ecx
sub    %ecx,%esi
mov    0x90(%rsp),%r14d
imul   %eax,%r14d
add    %esi,%r14d
lea    (%r14,%r8,2),%ecx
mov    %r8d,%r15d
add    %r14d,%r15d
lea    (%rdx,%r11,1),%r8d
add    $0x1,%r8d
lea    (%r15,%r8,1),%esi
movslq %esi,%r13
lea    0x1(%r11),%ebp
imul   %r9d,%ebp
lea    (%r15,%rbp,1),%esi
movslq %esi,%r12
imul   $0x55555556,%r11,%rsi
mov    %rsi,%rdi
shr    $0x3f,%rdi
shr    $0x20,%rsi
lea    (%rdx,%r15,1),%eax
add    $0x1,%eax
cltq
movsd  (%rbx,%rax,8),%xmm0
mov    %r13,0x48(%rsp)
addsd  (%rbx,%r13,8),%xmm0
lea    (%rsi,%rdi,1),%eax
mov    %r12,0x50(%rsp)
addsd  (%rbx,%r12,8),%xmm0
add    %edi,%esi
add    $0x1,%esi
mov    %r9d,0x14(%rsp)
imul   %r9d,%esi
add    %eax,%esi
lea    (%r15,%rsi,1),%eax
movslq %eax,%r10
movsd  %xmm0,(%rbx,%r10,8)
lea    (%rdx,%rcx,1),%eax
add    $0x1,%eax
cltq
movsd  (%rbx,%rax,8),%xmm0
lea    (%rcx,%r8,1),%eax
movslq %eax,%rdi
lea    (%rcx,%rbp,1),%eax
movslq %eax,%r9
mov    %rdi,0x40(%rsp)
addsd  (%rbx,%rdi,8),%xmm0
addsd  (%rbx,%r9,8),%xmm0
add    0x8(%rsp),%r14d
mov    %rcx,0x8(%rsp)
lea    (%rcx,%rsi,1),%eax
movslq %eax,%r13
movsd  %xmm0,(%rbx,%r13,8)
lea    (%rdx,%r14,1),%eax
add    $0x1,%eax
cltq
movsd  (%rbx,%rax,8),%xmm0
add    %r14d,%r8d
movslq %r8d,%r8
addsd  (%rbx,%r8,8),%xmm0
lea    (%r14,%rbp,1),%eax
movslq %eax,%rdi
addsd  (%rbx,%rdi,8),%xmm0
add    %r14d,%esi
movslq %esi,%r12
movsd  %xmm0,(%rbx,%r12,8)
mov    %r10,0x38(%rsp)
movsd  (%rbx,%r10,8),%xmm1
movsd  (%rbx,%r13,8),%xmm2
mulsd  %xmm2,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm0
jb     1284 <tricntr+0x174>
sqrtsd %xmm0,%xmm0
jmp    12b5 <tricntr+0x1a5>
mov    %r11,0x30(%rsp)
mov    %r9,0x28(%rsp)
mov    %r8,0x20(%rsp)
mov    %rdi,0x18(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm4,%xmm4
mov    0x18(%rsp),%rdi
mov    0x20(%rsp),%r8
mov    0x28(%rsp),%r9
mov    0x30(%rsp),%r11
movsd  0xd43(%rip),%xmm3
movapd %xmm3,%xmm1
divsd  %xmm0,%xmm1
mov    0x38(%rsp),%rax
movsd  (%rbx,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%rax,8)
movsd  (%rbx,%r13,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%r13,8)
mulsd  (%rbx,%r12,8),%xmm1
movsd  %xmm1,(%rbx,%r12,8)
add    %r11d,%ebp
lea    (%r15,%rbp,1),%eax
cltq
movsd  (%rbx,%rax,8),%xmm0
lea    (%r11,%r11,1),%eax
cltq
imul   $0x55555556,%rax,%rax
mov    0x48(%rsp),%rcx
addsd  (%rbx,%rcx,8),%xmm0
mov    %rax,%rcx
shr    $0x3f,%rcx
shr    $0x20,%rax
mov    0x50(%rsp),%rdx
addsd  (%rbx,%rdx,8),%xmm0
lea    (%rax,%rcx,1),%edx
add    %ecx,%eax
add    $0x1,%eax
imul   0x14(%rsp),%eax
add    %edx,%eax
add    %eax,%r15d
movslq %r15d,%r15
movsd  %xmm0,(%rbx,%r15,8)
mov    0x8(%rsp),%rdx
lea    (%rdx,%rbp,1),%ecx
movslq %ecx,%rcx
movsd  (%rbx,%rcx,8),%xmm0
mov    0x40(%rsp),%rcx
addsd  (%rbx,%rcx,8),%xmm0
addsd  (%rbx,%r9,8),%xmm0
add    %eax,%edx
movslq %edx,%r12
movsd  %xmm0,(%rbx,%r12,8)
add    %r14d,%ebp
movslq %ebp,%rcx
movsd  (%rbx,%rcx,8),%xmm0
addsd  (%rbx,%r8,8),%xmm0
addsd  (%rbx,%rdi,8),%xmm0
add    %r14d,%eax
movslq %eax,%rbp
movsd  %xmm0,(%rbx,%rbp,8)
movsd  (%rbx,%r15,8),%xmm1
movsd  (%rbx,%r12,8),%xmm2
mulsd  %xmm2,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
ucomisd %xmm4,%xmm0
jb     13c1 <tricntr+0x2b1>
sqrtsd %xmm0,%xmm0
jmp    13ce <tricntr+0x2be>
call   1030 <sqrt@plt>
movsd  0xc32(%rip),%xmm3
divsd  %xmm0,%xmm3
movsd  (%rbx,%r15,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%r15,8)
movsd  (%rbx,%r12,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%r12,8)
mulsd  (%rbx,%rbp,8),%xmm3
movsd  %xmm3,(%rbx,%rbp,8)
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @tricntr(double* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = sub nsw i32 %2, %1
%10 = xor i32 %9, -1
%11 = add nsw i32 %9, 1
%12 = sub i32 1, %3
%13 = add i32 %12, %4
%14 = mul nsw i32 %13, %11
%15 = mul nsw i32 %14, 3
%16 = mul i32 %14, -4
%17 = add i32 %16, %10
%18 = mul nsw i32 %15, %6
%19 = add i32 %17, %18
%20 = add i32 %19, %14
%21 = add nsw i32 %11, %7
%22 = add i32 %20, %21
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %0, i64 %23
%25 = add nsw i32 %7, 1
%26 = mul nsw i32 %25, %11
%27 = add i32 %20, %26
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %0, i64 %28
%30 = shl nsw i32 %14, 1
%31 = add i32 %19, %30
%32 = add i32 %31, %21
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %0, i64 %33
%35 = add i32 %31, %26
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %0, i64 %36
%38 = add i32 %19, %15
%39 = add i32 %38, %21
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %0, i64 %40
%42 = add i32 %38, %26
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %0, i64 %43
%45 = sdiv i32 %7, 3
%46 = add i32 %20, %11
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = load double, double* %24
%51 = fadd double %49, %50
%52 = load double, double* %29
%53 = fadd double %51, %52
%54 = add nsw i32 %45, 1
%55 = mul nsw i32 %54, %11
%56 = add nsw i32 %55, %45
%57 = add i32 %20, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %0, i64 %58
store double %53, double* %59
%60 = add i32 %31, %11
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = load double, double* %62
%64 = load double, double* %34
%65 = fadd double %63, %64
%66 = load double, double* %37
%67 = fadd double %65, %66
%68 = add i32 %31, %56
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %0, i64 %69
store double %67, double* %70
%71 = add i32 %38, %11
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %0, i64 %72
%74 = load double, double* %73
%75 = load double, double* %41
%76 = fadd double %74, %75
%77 = load double, double* %44
%78 = fadd double %76, %77
%79 = add i32 %38, %56
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %0, i64 %80
store double %78, double* %81
%82 = load double, double* %59
%83 = load double, double* %70
%84 = fmul double %83, %83
%85 = tail call double @llvm.fmuladd.f64(double %82, double %82, double %84)
%86 = tail call double @llvm.fmuladd.f64(double %78, double %78, double %85)
%87 = tail call double @sqrt(double %86)
%88 = fdiv double 1.000000e+00, %87
%89 = load double, double* %59
%90 = fmul double %88, %89
store double %90, double* %59
%91 = load double, double* %70
%92 = fmul double %88, %91
store double %92, double* %70
%93 = load double, double* %81
%94 = fmul double %88, %93
store double %94, double* %81
%95 = shl i32 %7, 1
%96 = sdiv i32 %95, 3
%97 = add i32 %26, %7
%98 = add i32 %20, %97
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %0, i64 %99
%101 = load double, double* %100
%102 = load double, double* %24
%103 = fadd double %101, %102
%104 = load double, double* %29
%105 = fadd double %103, %104
%106 = add nsw i32 %96, 1
%107 = mul nsw i32 %106, %11
%108 = add nsw i32 %107, %96
%109 = add i32 %20, %108
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %0, i64 %110
store double %105, double* %111
%112 = add i32 %31, %97
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %0, i64 %113
%115 = load double, double* %114
%116 = load double, double* %34
%117 = fadd double %115, %116
%118 = load double, double* %37
%119 = fadd double %117, %118
%120 = add i32 %31, %108
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %0, i64 %121
store double %119, double* %122
%123 = add i32 %38, %97
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %0, i64 %124
%126 = load double, double* %125
%127 = load double, double* %41
%128 = fadd double %126, %127
%129 = load double, double* %44
%130 = fadd double %128, %129
%131 = add i32 %38, %108
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %0, i64 %132
store double %130, double* %133
%134 = load double, double* %111
%135 = load double, double* %122
%136 = fmul double %135, %135
%137 = tail call double @llvm.fmuladd.f64(double %134, double %134, double %136)
%138 = tail call double @llvm.fmuladd.f64(double %130, double %130, double %137)
%139 = tail call double @sqrt(double %138)
%140 = fdiv double 1.000000e+00, %139
%141 = load double, double* %111
%142 = fmul double %140, %141
store double %142, double* %111
%143 = load double, double* %122
%144 = fmul double %140, %143
store double %144, double* %122
%145 = load double, double* %133
%146 = fmul double %140, %145
store double %146, double* %133
ret void
}