// File: /workspace/github_c_makefile_star2_O3/PoulinV[P]AxiCLASS/build/../source/fourier.c
// Function: fourier_pk_at_z
int fourier_pk_at_z(
                    struct background * pba,
                    struct fourier *pfo,
                    enum linear_or_logarithmic mode,
                    enum pk_outputs pk_output,
                    double z,
                    int index_pk,
                    double * out_pk, // array out_pk[index_k]
                    double * out_pk_ic // array out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic2]
                    ) {
  double tau;
  double ln_tau;
  int index_k;
  int index_ic1;
  int index_ic2;
  int index_ic1_ic1;
  int index_ic2_ic2;
  int index_ic1_ic2;
  int last_index;
  short do_ic = _FALSE_;

  /** - check whether we need the decomposition into contributions from each initial condition */

  if ((pk_output == pk_linear) && (pfo->ic_size > 1) && (out_pk_ic != NULL))
    do_ic = _TRUE_;

  class_test(pk_output == pk_nonlinear && pfo->method == nl_none, pfo->error_message, "Cannot get nonlinear power spectrum when no nonlinear method is employed");

  class_test(pk_output == pk_numerical_nowiggle && pfo->has_pk_numerical_nowiggle == _FALSE_, pfo->error_message, "Cannot get nowiggle power spectrum since it was not requested in input");

  class_test(pk_output == pk_analytic_nowiggle && pfo->has_pk_analytic_nowiggle == _FALSE_, pfo->error_message, "Cannot get analytic nowiggle power spectrum since it was not requested in input");

  /** - case z=0 requiring no interpolation in z. The
        pk_analytic_nowiggle can also be computed here because, in the
        current implementation, it is only computed at z=0. */
  if ((z == 0) || (pk_output == pk_analytic_nowiggle)) {

    for (index_k=0; index_k<pfo->k_size; index_k++) {

      switch (pk_output) {

      case pk_linear:
        out_pk[index_k] = pfo->ln_pk_l[index_pk][(pfo->ln_tau_size-1)*pfo->k_size+index_k];
        if (do_ic == _TRUE_) {
          for (index_ic1_ic2 = 0; index_ic1_ic2 < pfo->ic_ic_size; index_ic1_ic2++) {
            out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic2] =
              pfo->ln_pk_ic_l[index_pk][((pfo->ln_tau_size-1)*pfo->k_size+index_k)*pfo->ic_ic_size+index_ic1_ic2];
          }
        }
        break;

      case pk_nonlinear:
        out_pk[index_k] = pfo->ln_pk_nl[index_pk][(pfo->ln_tau_size-1)*pfo->k_size+index_k];
        break;

      case pk_numerical_nowiggle:
        out_pk[index_k] = pfo->ln_pk_l_nw_extra[(pfo->ln_tau_size-1)*pfo->k_size_extra+index_k];
        break;

      case pk_analytic_nowiggle:
        out_pk[index_k] = pfo->ln_pk_l_an_extra[index_k];
        break;

      }
    }
  }

  /** - interpolation in z */
  else {

    class_test(pfo->ln_tau_size == 1,
               pfo->error_message,
               "You are asking for the matter power spectrum at z=%e but the code was asked to store it only at z=0. You probably forgot to pass the input parameter z_max_pk (see explanatory.ini)",z);

    /** --> get value of contormal time tau */
    class_call(background_tau_of_z(pba,
                                   z,
                                   &tau),
               pba->error_message,
               pfo->error_message);

    ln_tau = log(tau);
    last_index = pfo->ln_tau_size-1;

    /** -> check that tau is in pre-computed table */

    if (ln_tau <= pfo->ln_tau[0]) {

      /** --> if ln(tau) much too small, raise an error */
      class_test(ln_tau<pfo->ln_tau[0]-100.*_EPSILON_,
                 pfo->error_message,
                 "requested z was not inside of tau tabulation range (Requested ln(tau_=%.10e, Min %.10e). Solution might be to increase input parameter z_max_pk (see explanatory.ini)",ln_tau,pfo->ln_tau[0]);

      /** --> if ln(tau) too small but within tolerance, round it and get right values without interpolating */
      ln_tau = pfo->ln_tau[0];

      for (index_k = 0 ; index_k < pfo->k_size; index_k++) {
        if (pk_output == pk_linear) {
          out_pk[index_k] = pfo->ln_pk_l[index_pk][index_k];
          if (do_ic == _TRUE_) {
            for (index_ic1_ic2 = 0; index_ic1_ic2 < pfo->ic_ic_size; index_ic1_ic2++) {
              out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic2] = pfo->ln_pk_ic_l[index_pk][index_k * pfo->ic_ic_size + index_ic1_ic2];
            }
          }
        }
        else if (pk_output == pk_nonlinear) {
          out_pk[index_k] = pfo->ln_pk_nl[index_pk][index_k];
        }
        else {
          out_pk[index_k] = pfo->ln_pk_l_nw_extra[index_k];
        }
      }
    }

    else if (ln_tau >= pfo->ln_tau[pfo->ln_tau_size-1]) {

      /** --> if ln(tau) much too large, raise an error */
      class_test(ln_tau>pfo->ln_tau[pfo->ln_tau_size-1]+_EPSILON_,
                 pfo->error_message,
                 "requested z was not inside of tau tabulation range (Requested ln(tau_=%.10e, Max %.10e) ",ln_tau,pfo->ln_tau[pfo->ln_tau_size-1]);

      /** --> if ln(tau) too large but within tolerance, round it and get right values without interpolating */
      ln_tau = pfo->ln_tau[pfo->ln_tau_size-1];

      for (index_k = 0 ; index_k < pfo->k_size; index_k++) {
        if (pk_output == pk_linear) {
          out_pk[index_k] = pfo->ln_pk_l[index_pk][(pfo->ln_tau_size-1) * pfo->k_size + index_k];
          if (do_ic == _TRUE_) {
            for (index_ic1_ic2 = 0; index_ic1_ic2 < pfo->ic_ic_size; index_ic1_ic2++) {
              out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic2] = pfo->ln_pk_ic_l[index_pk][((pfo->ln_tau_size-1) * pfo->k_size + index_k) * pfo->ic_ic_size + index_ic1_ic2];
            }
          }
        }
        else if (pk_output == pk_nonlinear) {
          out_pk[index_k] = pfo->ln_pk_nl[index_pk][(pfo->ln_tau_size-1) * pfo->k_size + index_k];
        }
        else {
          out_pk[index_k] = pfo->ln_pk_l_nw_extra[(pfo->ln_tau_size-1) * pfo->k_size_extra + index_k];
        }
      }
    }

    /** -> tau is in pre-computed table: interpolate */
    else {

      if (pk_output == pk_linear) {

        /** --> interpolate P_l(k) at tau from pre-computed array */
        class_call(array_interpolate_spline(pfo->ln_tau,
                                            pfo->ln_tau_size,
                                            pfo->ln_pk_l[index_pk],
                                            pfo->ddln_pk_l[index_pk],
                                            pfo->k_size,
                                            ln_tau,
                                            &last_index,
                                            out_pk,
                                            pfo->k_size,
                                            pfo->error_message),
                   pfo->error_message,
                   pfo->error_message);

        /** --> interpolate P_ic_l(k) at tau from pre-computed array */
        if (do_ic == _TRUE_) {
          class_call(array_interpolate_spline(pfo->ln_tau,
                                              pfo->ln_tau_size,
                                              pfo->ln_pk_ic_l[index_pk],
                                              pfo->ddln_pk_ic_l[index_pk],
                                              pfo->k_size*pfo->ic_ic_size,
                                              ln_tau,
                                              &last_index,
                                              out_pk_ic,
                                              pfo->k_size*pfo->ic_ic_size,
                                              pfo->error_message),
                     pfo->error_message,
                     pfo->error_message);
        }
      }

      else if (pk_output == pk_nonlinear) {

        if (ln_tau < pfo->ln_tau[pfo->ln_tau_size-pfo->ln_tau_size_nl]) {
          /** --> we requested P_nl(k) at tau where P_l is computed but the NL correction is NOT -> Return P_l */
          class_call(array_interpolate_spline(pfo->ln_tau,
                                              pfo->ln_tau_size,
                                              pfo->ln_pk_l[index_pk],
                                              pfo->ddln_pk_l[index_pk],
                                              pfo->k_size,
                                              ln_tau,
                                              &last_index,
                                              out_pk,
                                              pfo->k_size,
                                              pfo->error_message),
                     pfo->error_message,
                     pfo->error_message);

        }
        else {
          /** --> interpolate P_nl(k) at tau from pre-computed array */
          class_call(array_interpolate_spline(pfo->ln_tau+(pfo->ln_tau_size-pfo->ln_tau_size_nl),
                                              pfo->ln_tau_size_nl,
                                              pfo->ln_pk_nl[index_pk]+(pfo->ln_tau_size-pfo->ln_tau_size_nl)*pfo->k_size,
                                              pfo->ddln_pk_nl[index_pk]+(pfo->ln_tau_size-pfo->ln_tau_size_nl)*pfo->k_size,
                                              pfo->k_size,
                                              ln_tau,
                                              &last_index,
                                              out_pk,
                                              pfo->k_size,
                                              pfo->error_message),
                     pfo->error_message,
                     pfo->error_message);
        }
      }
      else {

        /** --> interpolate P_l_nw(k) at tau from pre-computed array */
        class_call(array_interpolate_spline(pfo->ln_tau,
                                            pfo->ln_tau_size,
                                            pfo->ln_pk_l_nw_extra,
                                            pfo->ddln_pk_l_nw_extra,
                                            pfo->k_size_extra,
                                            ln_tau,
                                            &last_index,
                                            out_pk,
                                            pfo->k_size,
                                            pfo->error_message),
                   pfo->error_message,
                   pfo->error_message);
      }
    }
  }

  /** - so far, all output stored in logarithmic format. Eventually, convert to linear one. */

  if (mode == linear) {

    /** --> loop over k */
    for (index_k=0; index_k<pfo->k_size; index_k++) {

      /** --> convert total spectrum */
      out_pk[index_k] = exp(out_pk[index_k]);

      if (do_ic == _TRUE_) {
        /** --> convert contribution of each ic (diagonal elements) */
        for (index_ic1=0; index_ic1 < pfo->ic_size; index_ic1++) {
          index_ic1_ic1 = index_symmetric_matrix(index_ic1,index_ic1,pfo->ic_size);

          out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic1] = exp(out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic1]);
        }

        /** --> convert contribution of each ic (non-diagonal elements) */
        for (index_ic1=0; index_ic1 < pfo->ic_size; index_ic1++) {
          for (index_ic2=index_ic1+1; index_ic2 < pfo->ic_size; index_ic2++) {
            index_ic1_ic1 = index_symmetric_matrix(index_ic1,index_ic1,pfo->ic_size);
            index_ic2_ic2 = index_symmetric_matrix(index_ic2,index_ic2,pfo->ic_size);
            index_ic1_ic2 = index_symmetric_matrix(index_ic1,index_ic2,pfo->ic_size);

            /* P_ic1xic2 = cos(angle) * sqrt(P_ic1 * P_ic2) */
            out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic2]
              = out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic2]
              *sqrt(out_pk_ic[index_k * pfo->ic_ic_size + index_ic1_ic1]
                    *out_pk_ic[index_k * pfo->ic_ic_size + index_ic2_ic2]);
          }
        }
      }
    }
  }

  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x1048, %rsp           # imm = 0x1048
movq %r9, %rbp
movl %ecx, %r15d
movl %edx, %r14d
movq %rdi, %r13
testl %ecx, %ecx
movq %rsi, (%rsp)
je 0x9d095
cmpl $0x3, %r15d
je 0x9d462
cmpl $0x2, %r15d
je 0x9d0b3
cmpl $0x1, %r15d
jne 0x9d0be
cmpl $0x0, (%rsi)
je 0x9d68e
movb $0x1, %r12b
jmp 0x9d0c1
cmpl $0x2, 0x40(%rsi)
setge %al
cmpq $0x0, 0x1080(%rsp)
setne %cl
andb %al, %cl
movl %ecx, 0x8(%rsp)
xorl %r12d, %r12d
jmp 0x9d0c9
cmpw $0x0, 0x3a(%rsi)
je 0x9d5d2
xorl %r12d, %r12d
movl $0x0, 0x8(%rsp)
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0x9d3fb
jp 0x9d3fb
movslq 0x6c(%rsi), %rax
testq %rax, %rax
jle 0x9d221
movslq %r8d, %rcx
xorl %edx, %edx
movl %r15d, %esi
leaq 0x7e247(%rip), %rdi     # 0x11b340
cmpl $0x3, %r15d
ja 0x9d215
movslq (%rdi,%rsi,4), %r8
addq %rdi, %r8
jmpq *%r8
movq (%rsp), %r9
movq 0xa8(%r9), %r8
movq (%r8,%rcx,8), %r8
movl 0x90(%r9), %r9d
decl %r9d
imull %eax, %r9d
leal (%rdx,%r9), %r10d
movslq %r10d, %r10
movsd (%r8,%r10,8), %xmm0
movsd %xmm0, (%rbp,%rdx,8)
cmpb $0x0, 0x8(%rsp)
je 0x9d215
movq (%rsp), %r8
movslq 0x44(%r8), %r8
testq %r8, %r8
jle 0x9d215
movq (%rsp), %r10
movq 0x98(%r10), %r11
movl %r8d, %r10d
imull %edx, %r10d
movslq %r10d, %r10
movq 0x1080(%rsp), %rbx
leaq (%rbx,%r10,8), %r10
addl %edx, %r9d
imull %r8d, %r9d
movslq %r9d, %r9
shlq $0x3, %r9
addq (%r11,%rcx,8), %r9
xorl %r11d, %r11d
movsd (%r9,%r11,8), %xmm0
movsd %xmm0, (%r10,%r11,8)
incq %r11
cmpq %r11, %r8
jne 0x9d18f
jmp 0x9d215
movq (%rsp), %r9
movl 0x90(%r9), %r8d
decl %r8d
imull 0x104(%r9), %r8d
addl %edx, %r8d
movslq %r8d, %r8
shlq $0x3, %r8
addq 0x130(%r9), %r8
jmp 0x9d20a
leaq (,%rdx,8), %r8
movq (%rsp), %r9
addq 0x118(%r9), %r8
jmp 0x9d20a
movq (%rsp), %r8
movq 0xb8(%r8), %r9
movl 0x90(%r8), %r8d
decl %r8d
imull %eax, %r8d
addl %edx, %r8d
movslq %r8d, %r8
shlq $0x3, %r8
addq (%r9,%rcx,8), %r8
movsd (%r8), %xmm0
movsd %xmm0, (%rbp,%rdx,8)
incq %rdx
cmpq %rax, %rdx
jne 0x9d0f9
xorl %eax, %eax
testl %r14d, %r14d
jne 0x9d7d7
movq (%rsp), %rcx
cmpl $0x0, 0x6c(%rcx)
jle 0x9d7d7
xorl %r15d, %r15d
movq %rbp, 0x30(%rsp)
movsd (%rbp,%r15,8), %xmm0
callq 0x8670
movq 0x1080(%rsp), %rdi
movq (%rsp), %rsi
movsd %xmm0, (%rbp,%r15,8)
cmpb $0x0, 0x8(%rsp)
je 0x9d3df
movl 0x40(%rsi), %r13d
testl %r13d, %r13d
jle 0x9d3df
movl 0x44(%rsi), %eax
xorl %ecx, %ecx
incl %r13d
imull %ecx, %r13d
leal 0x1(%rcx), %ebx
imull %ebx, %ecx
shrl %ecx
subl %ecx, %r13d
imull %r15d, %eax
addl %r13d, %eax
cltq
movsd (%rdi,%rax,8), %xmm0
callq 0x8670
movq 0x1080(%rsp), %rdi
movq (%rsp), %rsi
movl 0x44(%rsi), %eax
movl %eax, %ecx
imull %r15d, %ecx
addl %r13d, %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rdi,%rcx,8)
movl 0x40(%rsi), %r13d
movl %ebx, %ecx
cmpl %r13d, %ebx
jl 0x9d27e
testl %r13d, %r13d
jle 0x9d3df
movl $0x2, %eax
xorl %ecx, %ecx
xorl %edx, %edx
movl %eax, 0x2c(%rsp)
movl %edx, %r14d
leal 0x1(%r14), %eax
cmpl %r13d, %eax
movl %ecx, 0x28(%rsp)
movl %eax, 0x24(%rsp)
jge 0x9d3c5
shrl %ecx
negl %ecx
movq %rcx, 0x18(%rsp)
movl %eax, %r12d
imull %r14d, %r12d
shrl %r12d
leal 0x1(%rax), %ebx
movl %ebx, %ecx
imull %eax, %ecx
shrl %ecx
cmpl %eax, %r14d
jbe 0x9d32b
movl %r13d, %ebp
imull %eax, %ebp
addl %r14d, %ebp
subl %ecx, %ebp
jmp 0x9d33b
movl %r14d, %ebp
imull %r13d, %ebp
movq 0x18(%rsp), %rdx
addl %eax, %edx
addl %edx, %ebp
incl %r13d
movl %r13d, %edx
imull %r14d, %edx
subl %r12d, %edx
imull %eax, %r13d
subl %ecx, %r13d
movl 0x44(%rsi), %eax
imull %r15d, %eax
leal (%rax,%rbp), %ecx
movslq %ecx, %rcx
movsd (%rdi,%rcx,8), %xmm1
addl %eax, %edx
movslq %edx, %rcx
movsd (%rdi,%rcx,8), %xmm0
addl %r13d, %eax
cltq
mulsd (%rdi,%rax,8), %xmm0
ucomisd 0x50edb(%rip), %xmm0    # 0xee258
jb 0x9d385
sqrtsd %xmm0, %xmm0
jmp 0x9d3a2
movsd %xmm1, 0x10(%rsp)
callq 0x82b0
movsd 0x10(%rsp), %xmm1
movq 0x1080(%rsp), %rdi
movq (%rsp), %rsi
mulsd %xmm0, %xmm1
movl 0x44(%rsi), %eax
imull %r15d, %eax
addl %ebp, %eax
cltq
movsd %xmm1, (%rdi,%rax,8)
movl 0x40(%rsi), %r13d
movl %ebx, %eax
cmpl %r13d, %ebx
jl 0x9d30f
movl 0x2c(%rsp), %eax
movl 0x28(%rsp), %ecx
addl %eax, %ecx
addl $0x2, %eax
movl 0x24(%rsp), %edx
cmpl %r13d, %edx
jl 0x9d2e0
incq %r15
movslq 0x6c(%rsi), %rax
cmpq %rax, %r15
movq 0x30(%rsp), %rbp
jl 0x9d242
xorl %eax, %eax
jmp 0x9d7d7
cmpl $0x1, 0x90(%rsi)
jne 0x9d47a
movq %rsi, %rbx
leaq 0x7e11c(%rip), %rsi     # 0x11b52a
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movb $0x1, %al
callq 0x30af0
leaq 0x50c5e(%rip), %rsi     # 0xee082
leaq 0x7e1b3(%rip), %rdx     # 0x11b5de
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
addq $0x192, %rbx            # imm = 0x192
leaq 0x50be6(%rip), %rsi     # 0xee034
leaq 0x7df90(%rip), %rdx     # 0x11b3e5
movq %rbx, %rdi
movl $0x8e, %ecx
jmp 0x9d7c8
cmpw $0x0, 0x38(%rsi)
je 0x9d630
movl $0x0, 0x8(%rsp)
jmp 0x9d0dd
movl %r8d, %ebx
leaq 0x38(%rsp), %rsi
movq %r13, %rdi
callq 0x5c896
cmpl $0x1, %eax
jne 0x9d4dc
addq $0xd1c, %r13            # imm = 0xD1C
leaq 0x50f53(%rip), %rsi     # 0xee3f0
leaq 0x7499b(%rip), %rdx     # 0x111e3f
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movq %r13, %rcx
xorl %eax, %eax
callq 0x30af0
movq (%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x50b6c(%rip), %rsi     # 0xee034
leaq 0x7df16(%rip), %rdx     # 0x11b3e5
movl $0x95, %ecx
movq %r14, %r8
jmp 0x9d7cb
movsd 0x38(%rsp), %xmm0
callq 0x82c0
movq (%rsp), %rax
movslq 0x90(%rax), %rsi
leal -0x1(%rsi), %r10d
movl %r10d, 0xc(%rsp)
movq 0x88(%rax), %rdi
movsd (%rdi), %xmm1
ucomisd %xmm0, %xmm1
jae 0x9d6ec
movsd -0x8(%rdi,%rsi,8), %xmm1
ucomisd %xmm1, %xmm0
jae 0x9d75b
testl %r15d, %r15d
je 0x9d9ef
testb %r12b, %r12b
je 0x9da78
movq (%rsp), %r10
movl 0x94(%r10), %eax
movl %esi, %ecx
subl %eax, %ecx
movslq %ecx, %rcx
movsd (%rdi,%rcx,8), %xmm1
ucomisd %xmm0, %xmm1
movslq %ebx, %r9
jbe 0x9db9f
movq 0xa8(%r10), %rax
movq (%rax,%r9,8), %rdx
movq 0xb0(%r10), %rax
movq (%rax,%r9,8), %rcx
movl 0x6c(%r10), %r8d
leaq 0x192(%r10), %r15
subq $0x8, %rsp
leaq 0x14(%rsp), %r9
pushq %r15
pushq %r8
pushq %rbp
callq 0x18071
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0x9d221
leaq 0x50e52(%rip), %rsi     # 0xee3f0
leaq 0x7e1a4(%rip), %rdx     # 0x11b749
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x50a76(%rip), %rsi     # 0xee034
leaq 0x7de20(%rip), %rdx     # 0x11b3e5
movq %r15, %rdi
movl $0x107, %ecx            # imm = 0x107
jmp 0x9dc2a
movq %rsi, %rbx
leaq 0x7de19(%rip), %rsi     # 0x11b3f5
leaq 0x40(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x30af0
leaq 0x50a90(%rip), %rsi     # 0xee082
leaq 0x7de43(%rip), %rdx     # 0x11b43c
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
addq $0x192, %rbx            # imm = 0x192
leaq 0x50a18(%rip), %rsi     # 0xee034
leaq 0x7ddc2(%rip), %rdx     # 0x11b3e5
movq %rbx, %rdi
movl $0x62, %ecx
jmp 0x9d7c8
movq %rsi, %rbx
leaq 0x7de52(%rip), %rsi     # 0x11b48c
leaq 0x40(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x30af0
leaq 0x50a32(%rip), %rsi     # 0xee082
leaq 0x7de85(%rip), %rdx     # 0x11b4dc
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
addq $0x192, %rbx            # imm = 0x192
leaq 0x509ba(%rip), %rsi     # 0xee034
leaq 0x7dd64(%rip), %rdx     # 0x11b3e5
movq %rbx, %rdi
movl $0x64, %ecx
jmp 0x9d7c8
movq %rsi, %rbx
leaq 0x7dcd0(%rip), %rsi     # 0x11b368
leaq 0x40(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x30af0
leaq 0x509d4(%rip), %rsi     # 0xee082
leaq 0x7dcfc(%rip), %rdx     # 0x11b3b1
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
addq $0x192, %rbx            # imm = 0x192
leaq 0x5095c(%rip), %rsi     # 0xee034
leaq 0x7dd06(%rip), %rdx     # 0x11b3e5
movq %rbx, %rdi
movl $0x60, %ecx
jmp 0x9d7c8
movsd 0x7dc1c(%rip), %xmm2    # 0x11b310
addsd %xmm1, %xmm2
ucomisd %xmm0, %xmm2
jbe 0x9d7e9
leaq 0x7deeb(%rip), %rsi     # 0x11b5f4
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movb $0x2, %al
callq 0x30af0
leaq 0x50963(%rip), %rsi     # 0xee082
leaq 0x7df74(%rip), %rdx     # 0x11b69a
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
movq (%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x508e7(%rip), %rsi     # 0xee034
leaq 0x7dc91(%rip), %rdx     # 0x11b3e5
movl $0xa1, %ecx
jmp 0x9d7c8
movsd 0x55515(%rip), %xmm2    # 0xf2c78
addsd %xmm1, %xmm2
ucomisd %xmm2, %xmm0
jbe 0x9d8be
leaq 0x7df47(%rip), %rsi     # 0x11b6bf
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movb $0x2, %al
callq 0x30af0
leaq 0x508f4(%rip), %rsi     # 0xee082
leaq 0x7df83(%rip), %rdx     # 0x11b718
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
movq (%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x50878(%rip), %rsi     # 0xee034
leaq 0x7dc22(%rip), %rdx     # 0x11b3e5
movl $0xbd, %ecx
movq %r15, %r8
xorl %eax, %eax
callq 0x30af0
movl $0x1, %eax
addq $0x1048, %rsp           # imm = 0x1048
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq (%rsp), %rcx
movslq 0x6c(%rcx), %rax
testq %rax, %rax
jle 0x9d221
addq $0x130, %rcx            # imm = 0x130
movslq %ebx, %rdx
xorl %esi, %esi
leaq (,%rdx,8), %rdi
testl %r15d, %r15d
je 0x9d83a
movq %rcx, %r8
testb %r12b, %r12b
je 0x9d829
movq (%rsp), %r8
movq 0xb8(%r8), %r8
addq %rdi, %r8
movq (%r8), %r8
movsd (%r8,%rsi,8), %xmm0
movsd %xmm0, (%rbp,%rsi,8)
jmp 0x9d8ad
movq (%rsp), %r8
movq 0xa8(%r8), %r8
movq (%r8,%rdx,8), %r8
movsd (%r8,%rsi,8), %xmm0
movsd %xmm0, (%rbp,%rsi,8)
cmpb $0x0, 0x8(%rsp)
je 0x9d8ad
movq (%rsp), %r8
movslq 0x44(%r8), %r8
testq %r8, %r8
jle 0x9d8ad
movq (%rsp), %r9
movq 0x98(%r9), %r10
movl %r8d, %r9d
imull %esi, %r9d
movslq %r9d, %r11
leaq (,%r11,8), %r9
addq (%r10,%rdx,8), %r9
movq 0x1080(%rsp), %r10
leaq (%r10,%r11,8), %r10
xorl %r11d, %r11d
movsd (%r9,%r11,8), %xmm0
movsd %xmm0, (%r10,%r11,8)
incq %r11
cmpq %r11, %r8
jne 0x9d899
incq %rsi
cmpq %rax, %rsi
jne 0x9d80e
jmp 0x9d221
movq (%rsp), %rcx
movslq 0x6c(%rcx), %rax
testq %rax, %rax
jle 0x9d221
addq $0x130, %rcx            # imm = 0x130
movq %rcx, 0x10(%rsp)
movslq %ebx, %rdx
movl %eax, %esi
imull %r10d, %esi
movslq %esi, %rsi
leaq (,%rsi,8), %rdi
xorl %r8d, %r8d
leaq (,%rdx,8), %rcx
movq %rcx, 0x18(%rsp)
movslq %r10d, %r10
xorl %r11d, %r11d
movl %r14d, %ecx
movl %r12d, %r9d
testl %r15d, %r15d
je 0x9d931
testb %r12b, %r12b
je 0x9d9a4
movq (%rsp), %rbx
movq 0xb8(%rbx), %rbx
addq 0x18(%rsp), %rbx
movq %rsi, %r14
jmp 0x9d9b8
movq (%rsp), %rbx
movq 0xa8(%rbx), %rbx
movq (%rbx,%rdx,8), %rbx
leaq (%r11,%rsi), %r14
movsd (%rbx,%r14,8), %xmm0
movsd %xmm0, (%rbp,%r11,8)
cmpb $0x0, 0x8(%rsp)
je 0x9d9d0
movq (%rsp), %rbx
movslq 0x44(%rbx), %rbx
testq %rbx, %rbx
jle 0x9d9d0
movq (%rsp), %r14
movq 0x98(%r14), %r13
movq %r8, %r14
imulq %rbx, %r14
addq 0x1080(%rsp), %r14
movq %rdi, %r12
imulq %rbx, %r12
addq (%r13,%rdx,8), %r12
xorl %r13d, %r13d
movsd (%r12,%r13,8), %xmm0
movsd %xmm0, (%r14,%r13,8)
incq %r13
cmpq %r13, %rbx
jne 0x9d98e
jmp 0x9d9d0
movq (%rsp), %rbx
movslq 0x104(%rbx), %r14
imulq %r10, %r14
movq 0x10(%rsp), %rbx
leaq (,%r11,8), %r12
addq (%rbx), %r12
movsd (%r12,%r14,8), %xmm0
movsd %xmm0, (%rbp,%r11,8)
incq %r11
addq $0x8, %r8
addq $0x8, %rdi
cmpq %rax, %r11
movl %r9d, %r12d
jne 0x9d90b
movl %ecx, %r14d
jmp 0x9d221
movq (%rsp), %r9
movq 0xa8(%r9), %rax
movslq %ebx, %rbx
movq (%rax,%rbx,8), %rdx
movq 0xb0(%r9), %rax
movq (%rax,%rbx,8), %rcx
movl 0x6c(%r9), %r8d
leaq 0x192(%r9), %r15
subq $0x8, %rsp
leaq 0x14(%rsp), %r9
movsd %xmm0, 0x18(%rsp)
pushq %r15
pushq %r8
pushq %rbp
callq 0x18071
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0x9daf6
leaq 0x509ac(%rip), %rsi     # 0xee3f0
leaq 0x7dcfe(%rip), %rdx     # 0x11b749
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x505d0(%rip), %rsi     # 0xee034
leaq 0x7d97a(%rip), %rdx     # 0x11b3e5
movq %r15, %rdi
movl $0xe5, %ecx
jmp 0x9dc2a
movq (%rsp), %r9
movq 0x130(%r9), %rdx
movq 0x138(%r9), %rcx
movl 0x6c(%r9), %eax
movl 0x104(%r9), %r8d
leaq 0x192(%r9), %r15
subq $0x8, %rsp
leaq 0x14(%rsp), %r9
pushq %r15
pushq %rax
pushq %rbp
callq 0x18071
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0x9d221
leaq 0x5092e(%rip), %rsi     # 0xee3f0
leaq 0x7df5a(%rip), %rdx     # 0x11ba23
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x50552(%rip), %rsi     # 0xee034
leaq 0x7d8fc(%rip), %rdx     # 0x11b3e5
movq %r15, %rdi
movl $0x128, %ecx            # imm = 0x128
jmp 0x9dc2a
cmpb $0x0, 0x8(%rsp)
je 0x9dc32
movq (%rsp), %r9
movq 0x88(%r9), %rdi
movl 0x90(%r9), %esi
movq 0x98(%r9), %rax
movq (%rax,%rbx,8), %rdx
movq 0xa0(%r9), %rax
movq (%rax,%rbx,8), %rcx
movl 0x44(%r9), %r8d
imull 0x6c(%r9), %r8d
subq $0x8, %rsp
leaq 0x14(%rsp), %r9
movsd 0x18(%rsp), %xmm0
pushq %r15
pushq %r8
pushq 0x1098(%rsp)
callq 0x18071
addq $0x20, %rsp
movb $0x1, %cl
movl %ecx, 0x8(%rsp)
cmpl $0x1, %eax
jne 0x9d221
leaq 0x50885(%rip), %rsi     # 0xee3f0
leaq 0x7dc8c(%rip), %rdx     # 0x11b7fe
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x504a9(%rip), %rsi     # 0xee034
leaq 0x7d853(%rip), %rdx     # 0x11b3e5
movq %r15, %rdi
movl $0xf4, %ecx
jmp 0x9dc2a
leaq (%rdi,%rcx,8), %rdi
movq 0xb8(%r10), %rdx
movl 0x6c(%r10), %r8d
imull %r8d, %ecx
movslq %ecx, %rcx
shlq $0x3, %rcx
movq (%rdx,%r9,8), %rdx
addq %rcx, %rdx
movq 0xc0(%r10), %rsi
addq (%rsi,%r9,8), %rcx
leaq 0x192(%r10), %r15
subq $0x8, %rsp
leaq 0x14(%rsp), %r9
movl %eax, %esi
pushq %r15
pushq %r8
pushq %rbp
callq 0x18071
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0x9d221
leaq 0x507f5(%rip), %rsi     # 0xee3f0
leaq 0x7dcda(%rip), %rdx     # 0x11b8dc
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x50419(%rip), %rsi     # 0xee034
leaq 0x7d7c3(%rip), %rdx     # 0x11b3e5
movq %r15, %rdi
movl $0x117, %ecx            # imm = 0x117
movq %rbx, %r8
jmp 0x9d7cb
movl $0x0, 0x8(%rsp)
jmp 0x9d221