// File: /workspace/github_cpp_makefile_star2_O2/j-faria[P]kima/examples/trends/../../celerite/cpp/include/celerite/solver/cholesky.h
// Function: celerite::solver::CholeskySolver<double, -1>::dot_L(Eigen::Matrix<double, -1, -1, 0, -1, -1> const&) const
matrix_t dot_L (const Eigen::MatrixXd& z) const {
  if (z.rows() != this->N_) throw dimension_mismatch();
  if (!(this->computed_)) throw compute_exception();

  T tmp;
  int N = z.rows(), nrhs = z.cols();
  Eigen::Array<T, Eigen::Dynamic, 1> D = sqrt(D_.array());
  vector_j_t f(J_);
  matrix_t y(N, nrhs);

  for (int k = 0; k < nrhs; ++k) {
    f.setZero();
    tmp = z(0, k) * D(0);
    y(0, k) = tmp;
    for (int n = 1; n < N; ++n) {
      f = phi_.col(n-1).asDiagonal() * (f + W_.col(n-1) * tmp);
      tmp = D(n) * z(n, k);
      y(n, k) = tmp + u_.col(n-1).transpose().dot(f);
    }
  }

  return y;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x278, %rsp            # imm = 0x278
movq %rdi, 0x10(%rsp)
movq %rdx, 0x20(%rsp)
movq 0x8(%rdx), %rax
movslq 0xc(%rsi), %rcx
cmpq %rcx, %rax
jne 0x2c894
movq %rsi, %rbx
cmpb $0x0, 0x8(%rsi)
je 0x2c8af
movl %eax, 0x1c(%rsp)
movq 0x20(%rsp), %rax
movl 0x10(%rax), %eax
movl %eax, 0x18(%rsp)
leaq 0x68(%rbx), %rax
leaq 0x180(%rsp), %rsi
movq %rax, (%rsi)
leaq 0x30(%rsp), %rdi
callq 0x32e64
leaq 0x18(%rbx), %rsi
leaq 0x60(%rsp), %rdi
callq 0x11e0a
leaq 0x1c(%rsp), %rsi
leaq 0x18(%rsp), %rdx
movq 0x10(%rsp), %rdi
callq 0x15068
leaq 0x38(%rbx), %rax
movq %rax, 0x50(%rsp)
leaq 0x50(%rbx), %rax
movq %rax, 0x48(%rsp)
addq $0x20, %rbx
movq %rbx, 0x40(%rsp)
xorl %ebp, %ebp
leaq 0x60(%rsp), %rbx
leaq 0x220(%rsp), %r14
leaq 0x118(%rsp), %r12
leaq 0x180(%rsp), %r13
movslq 0x18(%rsp), %rax
cmpq %rax, %rbp
jge 0x2c869
movq %rbx, %rdi
callq 0x1e92a
movq 0x20(%rsp), %rdi
xorl %esi, %esi
movq %rbp, %rdx
callq 0x1e570
movsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
leaq 0x30(%rsp), %rdi
xorl %esi, %esi
callq 0x32e2e
movsd 0x8(%rsp), %xmm0
mulsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd %xmm0, 0x28(%rsp)
movq 0x10(%rsp), %rdi
xorl %esi, %esi
movq %rbp, %rdx
callq 0x11dc0
movsd 0x8(%rsp), %xmm0
movsd %xmm0, (%rax)
pushq $0x1
popq %r15
movslq 0x1c(%rsp), %rax
cmpq %rax, %r15
jge 0x2c861
decq %r15
leaq 0xa8(%rsp), %rdi
movq 0x50(%rsp), %rsi
movq %r15, %rdx
callq 0x203fe
movups 0xa8(%rsp), %xmm0
movaps %xmm0, 0x70(%rsp)
leaq 0xc0(%rsp), %rax
movupd (%rax), %xmm0
movups 0x10(%rax), %xmm1
leaq 0x88(%rsp), %rax
movups %xmm1, 0x10(%rax)
movupd %xmm0, (%rax)
leaq 0xe0(%rsp), %rdi
movq 0x48(%rsp), %rsi
movq %r15, %rdx
callq 0x203fe
movq %r14, %rdi
leaq 0xe0(%rsp), %rsi
leaq 0x28(%rsp), %rdx
callq 0x1e956
movq %r12, %rdi
movq %rbx, %rsi
movq %r14, %rdx
callq 0x1e93a
movq %r13, %rdi
leaq 0x70(%rsp), %rsi
movq %r12, %rdx
callq 0x204d8
movq %rbx, %rdi
movq %r13, %rsi
callq 0x205e6
incq %r15
leaq 0x30(%rsp), %rdi
movq %r15, %rsi
callq 0x32e2e
movsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movq 0x20(%rsp), %rdi
movq %r15, %rsi
movq %rbp, %rdx
callq 0x1e570
movsd 0x8(%rsp), %xmm0
mulsd (%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd %xmm0, 0x28(%rsp)
decq %r15
movq %r12, %rdi
movq 0x40(%rsp), %rsi
movq %r15, %rdx
callq 0x203fe
movups 0x118(%rsp), %xmm0
movaps %xmm0, 0x180(%rsp)
leaq 0x130(%rsp), %rax
movupd (%rax), %xmm0
movups 0x10(%rax), %xmm1
leaq 0x198(%rsp), %rax
movups %xmm1, 0x10(%rax)
movupd %xmm0, (%rax)
movq %r13, %rdi
movq %rbx, %rsi
callq 0x1e9a6
movsd %xmm0, 0x58(%rsp)
incq %r15
movq 0x10(%rsp), %rdi
movq %r15, %rsi
movq %rbp, %rdx
callq 0x11dc0
movsd 0x8(%rsp), %xmm0
addsd 0x58(%rsp), %xmm0
movsd %xmm0, (%rax)
incq %r15
jmp 0x2c6fc
incq %rbp
jmp 0x2c68e
movq 0x60(%rsp), %rdi
callq 0x64a0
movq 0x30(%rsp), %rdi
callq 0x64a0
movq 0x10(%rsp), %rax
addq $0x278, %rsp            # imm = 0x278
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
pushq $0x8
popq %rdi
callq 0x61d0
leaq 0x36f3d(%rip), %rcx     # 0x637e0
movq %rcx, (%rax)
leaq 0x36ef3(%rip), %rsi     # 0x637a0
jmp 0x2c8c8
pushq $0x8
popq %rdi
callq 0x61d0
leaq 0x36f4a(%rip), %rcx     # 0x63808
movq %rcx, (%rax)
leaq 0x36ef0(%rip), %rsi     # 0x637b8
movq 0x37701(%rip), %rdx     # 0x63fd0
movq %rax, %rdi
callq 0x6690
movq %rax, %rbx
jmp 0x2c8f6
movq %rax, %rbx
jmp 0x2c900
jmp 0x2c8e9
jmp 0x2c8e9
jmp 0x2c8e9
jmp 0x2c8e9
movq %rax, %rbx
movq 0x10(%rsp), %rdi
callq 0x144ac
movq 0x60(%rsp), %rdi
callq 0x64a0
movq 0x30(%rsp), %rdi
callq 0x64a0
movq %rbx, %rdi
callq 0x66e0