// File: /workspace/github_c_makefile_star2_O1/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 0x5efb7
cmpl $0x1, %eax
jne 0xa2c66
addq $0xd1c, %r12            # imm = 0xD1C
leaq 0x4e7d4(%rip), %rsi     # 0xf13f0
leaq 0x721cc(%rip), %rdx     # 0x114def
leaq 0x10(%rsp), %r14
movq %r14, %rdi
movq %r12, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x4e3f1(%rip), %rsi     # 0xf1034
leaq 0x7c83a(%rip), %rdx     # 0x11f484
movq %rbx, %rdi
movl $0x4ad, %ecx            # imm = 0x4AD
movq %r14, %r8
xorl %eax, %eax
callq 0x30c2c
movl $0x1, %ebp
jmp 0xa2d5f
cmpw $0x1, 0x50(%rbx)
jne 0xa2d2a
movl 0x140(%rbx), %r9d
cmpl $0x1, %r9d
jne 0xa2c9b
movq 0x158(%rbx), %rax
movslq 0x54(%rbx), %rcx
movq (%rax,%rcx,8), %rax
movsd (%rax), %xmm0
movsd %xmm0, (%r15)
jmp 0xa2d2a
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 0x19f67
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa2d2a
leaq 0x4e701(%rip), %rsi     # 0xf13f0
leaq 0x7c7a0(%rip), %rdx     # 0x11f496
leaq 0x10(%rsp), %rbx
movq %rbx, %rdi
movq %r12, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x4e325(%rip), %rsi     # 0xf1034
leaq 0x7c76e(%rip), %rdx     # 0x11f484
movq %r12, %rdi
movl $0x4c2, %ecx            # imm = 0x4C2
movq %rbx, %r8
xorl %eax, %eax
callq 0x30c2c
jmp 0xa2d5f
cmpw $0x1, 0x52(%rbx)
jne 0xa2d53
movl 0x140(%rbx), %r9d
cmpl $0x1, %r9d
jne 0xa2d71
movq 0x158(%rbx), %rax
movslq 0x58(%rbx), %rcx
movq (%rax,%rcx,8), %rax
movsd (%rax), %xmm0
jmp 0xa2d58
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 0x19f67
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa2d5f
leaq 0x4e62f(%rip), %rsi     # 0xf13f0
leaq 0x7c744(%rip), %rdx     # 0x11f50c
leaq 0x10(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x4e253(%rip), %rsi     # 0xf1034
leaq 0x7c69c(%rip), %rdx     # 0x11f484
movq %rbx, %rdi
movl $0x4d9, %ecx            # imm = 0x4D9
jmp 0xa2c52