// File: /workspace/github_c_makefile_star2_O3/PoulinV[P]AxiCLASS/build/../source/fourier.c
// Function: fourier_k_nl_at_z
int fourier_k_nl_at_z(
                      struct background *pba,
                      struct fourier * pfo,
                      double z,
                      double * k_nl,
                      double * k_nl_cb
                      ) {

  double tau;

  /** - convert input redshift into a conformal time */

  class_call(background_tau_of_z(pba,
                                 z,
                                 &tau),
             pba->error_message,
             pfo->error_message);

  /** - interpolate the precomputed k_nl array at the needed valuetime */

  if (pfo->has_pk_m == _TRUE_) {

    if (pfo->tau_size == 1) {
      *k_nl = pfo->k_nl[pfo->index_pk_m][0];
    }
    else {
      class_call(array_interpolate_two(pfo->tau,
                                       1,
                                       0,
                                       pfo->k_nl[pfo->index_pk_m],
                                       1,
                                       pfo->tau_size,
                                       tau,
                                       k_nl,
                                       1,
                                       pfo->error_message),
                 pfo->error_message,
                 pfo->error_message);
    }
  }

  /** - if needed, do the same for the baryon part only */

  if (pfo->has_pk_cb == _TRUE_) {

    if (pfo->tau_size == 1) {
      *k_nl_cb = pfo->k_nl[pfo->index_pk_cb][0];
    }
    else {
      class_call(array_interpolate_two(pfo->tau,
                                       1,
                                       0,
                                       pfo->k_nl[pfo->index_pk_cb],
                                       1,
                                       pfo->tau_size,
                                       tau,
                                       k_nl_cb,
                                       1,
                                       pfo->error_message),
                 pfo->error_message,
                 pfo->error_message);
    }
  }

  /* otherwise, return the same for k_nl_cb as for k_nl */

  else{
    *k_nl_cb = *k_nl;
  }

  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x810, %rsp            # imm = 0x810
movq %rcx, %r14
movq %rdx, %r15
movq %rsi, %rbx
movq %rdi, %r12
leaq 0x8(%rsp), %rsi
callq 0x5c896
cmpl $0x1, %eax
jne 0x9fd9d
addq $0xd1c, %r12            # imm = 0xD1C
leaq 0x4e69d(%rip), %rsi     # 0xee3f0
leaq 0x720e5(%rip), %rdx     # 0x111e3f
leaq 0x10(%rsp), %r14
movq %r14, %rdi
movq %r12, %rcx
xorl %eax, %eax
callq 0x30af0
addq $0x192, %rbx            # imm = 0x192
leaq 0x4e2ba(%rip), %rsi     # 0xee034
leaq 0x7c743(%rip), %rdx     # 0x11c4c4
movq %rbx, %rdi
movl $0x4ad, %ecx            # imm = 0x4AD
movq %r14, %r8
xorl %eax, %eax
callq 0x30af0
movl $0x1, %ebp
jmp 0x9fe96
cmpw $0x1, 0x50(%rbx)
jne 0x9fe61
movl 0x140(%rbx), %r9d
cmpl $0x1, %r9d
jne 0x9fdd2
movq 0x158(%rbx), %rax
movslq 0x54(%rbx), %rcx
movq (%rax,%rcx,8), %rax
movsd (%rax), %xmm0
movsd %xmm0, (%r15)
jmp 0x9fe61
movq 0x148(%rbx), %rdi
movq 0x158(%rbx), %rax
movslq 0x54(%rbx), %rcx
movq (%rax,%rcx,8), %rcx
movsd 0x8(%rsp), %xmm0
leaq 0x192(%rbx), %r12
subq $0x8, %rsp
movl $0x1, %ebp
movl $0x1, %esi
xorl %edx, %edx
movl $0x1, %r8d
pushq %r12
pushq $0x1
pushq %r15
callq 0x1a08a
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0x9fe61
leaq 0x4e5ca(%rip), %rsi     # 0xee3f0
leaq 0x7c6a9(%rip), %rdx     # 0x11c4d6
leaq 0x10(%rsp), %rbx
movq %rbx, %rdi
movq %r12, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x4e1ee(%rip), %rsi     # 0xee034
leaq 0x7c677(%rip), %rdx     # 0x11c4c4
movq %r12, %rdi
movl $0x4c2, %ecx            # imm = 0x4C2
movq %rbx, %r8
xorl %eax, %eax
callq 0x30af0
jmp 0x9fe96
cmpw $0x1, 0x52(%rbx)
jne 0x9fe8a
movl 0x140(%rbx), %r9d
cmpl $0x1, %r9d
jne 0x9fea8
movq 0x158(%rbx), %rax
movslq 0x58(%rbx), %rcx
movq (%rax,%rcx,8), %rax
movsd (%rax), %xmm0
jmp 0x9fe8f
movsd (%r15), %xmm0
movsd %xmm0, (%r14)
xorl %ebp, %ebp
movl %ebp, %eax
addq $0x810, %rsp            # imm = 0x810
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
movq 0x148(%rbx), %rdi
movq 0x158(%rbx), %rax
movslq 0x58(%rbx), %rcx
movq (%rax,%rcx,8), %rcx
movsd 0x8(%rsp), %xmm0
addq $0x192, %rbx            # imm = 0x192
subq $0x8, %rsp
xorl %ebp, %ebp
movl $0x1, %esi
xorl %edx, %edx
movl $0x1, %r8d
pushq %rbx
pushq $0x1
pushq %r14
callq 0x1a08a
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0x9fe96
leaq 0x4e4f8(%rip), %rsi     # 0xee3f0
leaq 0x7c64d(%rip), %rdx     # 0x11c54c
leaq 0x10(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x4e11c(%rip), %rsi     # 0xee034
leaq 0x7c5a5(%rip), %rdx     # 0x11c4c4
movq %rbx, %rdi
movl $0x4d9, %ecx            # imm = 0x4D9
jmp 0x9fd89