// File: /workspace/github_c_makefile_star2_O1/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 %r8d, %r12d
movl %ecx, %r15d
movl %edx, %r14d
movq %rsi, %rbx
movq %rdi, %r13
testl %ecx, %ecx
je 0x9ff99
movl $0x0, 0x10(%rsp)
jmp 0x9ffb2
cmpl $0x2, 0x40(%rbx)
setge %al
cmpq $0x0, 0x1080(%rsp)
setne %dl
andb %al, %dl
movl %edx, 0x10(%rsp)
cmpl $0x1, %r15d
jne 0x9ffc1
cmpl $0x0, (%rbx)
je 0xa0383
cmpl $0x2, %r15d
jne 0x9ffd2
cmpw $0x0, 0x3a(%rbx)
je 0xa03de
cmpl $0x3, %r15d
jne 0x9ffe3
cmpw $0x0, 0x38(%rbx)
je 0xa049a
cmpl $0x3, %r15d
movq %rbx, 0x8(%rsp)
je 0xa0062
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0x9fffa
jnp 0xa0062
cmpl $0x1, 0x90(%rbx)
jne 0xa0439
leaq 0x7e4dc(%rip), %rsi     # 0x11e4ea
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movb $0x1, %al
callq 0x30c2c
leaq 0x5105e(%rip), %rsi     # 0xf1082
leaq 0x7e573(%rip), %rdx     # 0x11e59e
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x50fe6(%rip), %rsi     # 0xf1034
leaq 0x7e350(%rip), %rdx     # 0x11e3a5
movq %rbx, %rdi
movl $0x8e, %ecx
jmp 0xa04f0
movslq 0x6c(%rbx), %rax
testq %rax, %rax
jle 0xa01ac
movslq %r12d, %rcx
xorl %edx, %edx
movl %r15d, %esi
leaq 0x7e282(%rip), %rdi     # 0x11e300
cmpl $0x3, %r15d
ja 0xa01a0
movslq (%rdi,%rsi,4), %r8
addq %rdi, %r8
jmpq *%r8
movq 0x8(%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, 0x10(%rsp)
je 0xa01a0
movq 0x8(%rsp), %r8
movslq 0x44(%r8), %r8
testq %r8, %r8
jle 0xa01a0
movq 0x8(%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 0xa0117
jmp 0xa01a0
movq 0x8(%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 0xa0195
leaq (,%rdx,8), %r8
movq 0x8(%rsp), %r9
addq 0x118(%r9), %r8
jmp 0xa0195
movq 0x8(%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 0xa007e
xorl %eax, %eax
testl %r14d, %r14d
jne 0xa04ff
movq 0x8(%rsp), %rcx
cmpl $0x0, 0x6c(%rcx)
jle 0xa04ff
xorl %r15d, %r15d
movq %rbp, 0x30(%rsp)
movsd (%rbp,%r15,8), %xmm0
callq 0x8670
movq 0x8(%rsp), %rdi
movsd %xmm0, (%rbp,%r15,8)
cmpb $0x0, 0x10(%rsp)
movq 0x1080(%rsp), %r8
je 0xa0367
movl 0x40(%rdi), %ebx
testl %ebx, %ebx
jle 0xa024e
xorl %eax, %eax
incl %ebx
imull %eax, %ebx
leal 0x1(%rax), %ebp
imull %ebp, %eax
shrl %eax
subl %eax, %ebx
movl 0x44(%rdi), %eax
imull %r15d, %eax
addl %ebx, %eax
cltq
movsd (%r8,%rax,8), %xmm0
callq 0x8670
movq 0x1080(%rsp), %r8
movq 0x8(%rsp), %rdi
movl 0x44(%rdi), %eax
imull %r15d, %eax
addl %ebx, %eax
cltq
movsd %xmm0, (%r8,%rax,8)
movl 0x40(%rdi), %ebx
movl %ebp, %eax
cmpl %ebx, %ebp
jl 0xa0202
cmpl $0x0, 0x40(%rdi)
jle 0xa0367
movl $0x2, %eax
movl $0x0, 0x24(%rsp)
xorl %ecx, %ecx
movl %eax, 0x2c(%rsp)
movl %ecx, %ebx
leal 0x1(%rbx), %ecx
movl 0x40(%rdi), %eax
movl %ecx, 0x28(%rsp)
cmpl %eax, %ecx
jge 0xa034f
movl 0x24(%rsp), %ecx
movl %ecx, %ebp
shrl %ebp
negl %ebp
movl 0x28(%rsp), %ecx
movl %ecx, %r13d
imull %ebx, %r13d
shrl %r13d
leal 0x1(%rcx), %r12d
movl %r12d, %edx
imull %ecx, %edx
shrl %edx
cmpl %ecx, %ebx
jbe 0xa02b6
movl %eax, %r14d
imull %ecx, %r14d
addl %ebx, %r14d
subl %edx, %r14d
jmp 0xa02c3
movl %ebx, %r14d
imull %eax, %r14d
addl %ebp, %r14d
addl %ecx, %r14d
incl %eax
movl %eax, %esi
imull %ebx, %esi
subl %r13d, %esi
imull %ecx, %eax
subl %edx, %eax
movl 0x44(%rdi), %ecx
imull %r15d, %ecx
leal (%rcx,%r14), %edx
movslq %edx, %rdx
movsd (%r8,%rdx,8), %xmm1
addl %ecx, %esi
movslq %esi, %rdx
movsd (%r8,%rdx,8), %xmm0
addl %eax, %ecx
movslq %ecx, %rax
mulsd (%r8,%rax,8), %xmm0
ucomisd 0x50f54(%rip), %xmm0    # 0xf1258
jb 0xa030c
sqrtsd %xmm0, %xmm0
jmp 0xa032a
movsd %xmm1, 0x18(%rsp)
callq 0x82b0
movsd 0x18(%rsp), %xmm1
movq 0x1080(%rsp), %r8
movq 0x8(%rsp), %rdi
mulsd %xmm0, %xmm1
movl 0x44(%rdi), %eax
imull %r15d, %eax
addl %r14d, %eax
cltq
movsd %xmm1, (%r8,%rax,8)
movl 0x40(%rdi), %eax
movl %r12d, %ecx
cmpl %eax, %r12d
jl 0xa0297
movl 0x2c(%rsp), %eax
addl %eax, 0x24(%rsp)
addl $0x2, %eax
movl 0x28(%rsp), %ecx
cmpl 0x40(%rdi), %ecx
jl 0xa0267
incq %r15
movslq 0x6c(%rdi), %rax
cmpq %rax, %r15
movq 0x30(%rsp), %rbp
jl 0xa01ce
xorl %eax, %eax
jmp 0xa04ff
leaq 0x7df9e(%rip), %rsi     # 0x11e328
leaq 0x40(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x30c2c
leaq 0x50ce2(%rip), %rsi     # 0xf1082
leaq 0x7dfca(%rip), %rdx     # 0x11e371
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x50c6a(%rip), %rsi     # 0xf1034
leaq 0x7dfd4(%rip), %rdx     # 0x11e3a5
movq %rbx, %rdi
movl $0x60, %ecx
jmp 0xa04f0
leaq 0x7dfd0(%rip), %rsi     # 0x11e3b5
leaq 0x40(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x30c2c
leaq 0x50c87(%rip), %rsi     # 0xf1082
leaq 0x7dffa(%rip), %rdx     # 0x11e3fc
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x50c0f(%rip), %rsi     # 0xf1034
leaq 0x7df79(%rip), %rdx     # 0x11e3a5
movq %rbx, %rdi
movl $0x62, %ecx
jmp 0xa04f0
leaq 0x38(%rsp), %rsi
movq %r13, %rdi
callq 0x5efb7
cmpl $0x1, %eax
jne 0xa0511
addq $0xd1c, %r13            # imm = 0xD1C
leaq 0x50f93(%rip), %rsi     # 0xf13f0
leaq 0x7498b(%rip), %rdx     # 0x114def
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movq %r13, %rcx
xorl %eax, %eax
callq 0x30c2c
movq 0x8(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x50bab(%rip), %rsi     # 0xf1034
leaq 0x7df15(%rip), %rdx     # 0x11e3a5
movl $0x95, %ecx
movq %r14, %r8
jmp 0xa04f3
leaq 0x7dfab(%rip), %rsi     # 0x11e44c
leaq 0x40(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x30c2c
leaq 0x50bcb(%rip), %rsi     # 0xf1082
leaq 0x7dfde(%rip), %rdx     # 0x11e49c
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x50b53(%rip), %rsi     # 0xf1034
leaq 0x7debd(%rip), %rdx     # 0x11e3a5
movq %rbx, %rdi
movl $0x64, %ecx
movq %r15, %r8
xorl %eax, %eax
callq 0x30c2c
movl $0x1, %eax
addq $0x1048, %rsp           # imm = 0x1048
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movsd 0x38(%rsp), %xmm0
callq 0x82c0
movq 0x8(%rsp), %r10
movslq 0x90(%r10), %rsi
leal -0x1(%rsi), %eax
movl %eax, 0x14(%rsp)
movq 0x88(%r10), %rdi
movsd (%rdi), %xmm1
ucomisd %xmm0, %xmm1
jae 0xa0603
movsd -0x8(%rdi,%rsi,8), %xmm1
ucomisd %xmm1, %xmm0
jae 0xa0677
testl %r15d, %r15d
je 0xa08f5
cmpl $0x1, %r15d
jne 0xa097a
movl 0x94(%r10), %eax
movl %esi, %ecx
subl %eax, %ecx
movslq %ecx, %rcx
movsd (%rdi,%rcx,8), %xmm1
ucomisd %xmm0, %xmm1
movslq %r12d, %r9
jbe 0xa0a98
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 0x1c(%rsp), %r9
pushq %r15
pushq %r8
pushq %rbp
callq 0x17ee6
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa01ac
leaq 0x50e21(%rip), %rsi     # 0xf13f0
leaq 0x7e133(%rip), %rdx     # 0x11e709
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x50a45(%rip), %rsi     # 0xf1034
leaq 0x7ddaf(%rip), %rdx     # 0x11e3a5
movq %r15, %rdi
movl $0x107, %ecx            # imm = 0x107
jmp 0xa0b23
addsd 0x7dcc5(%rip), %xmm1    # 0x11e2d0
ucomisd %xmm0, %xmm1
jbe 0xa06ed
movsd (%rdi), %xmm1
leaq 0x7df94(%rip), %rsi     # 0x11e5b4
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movb $0x2, %al
movq %r10, %rbx
callq 0x30c2c
leaq 0x50a49(%rip), %rsi     # 0xf1082
leaq 0x7e01a(%rip), %rdx     # 0x11e65a
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x509d1(%rip), %rsi     # 0xf1034
leaq 0x7dd3b(%rip), %rdx     # 0x11e3a5
movq %rbx, %rdi
movl $0xa1, %ecx
jmp 0xa04f0
addsd 0x555c9(%rip), %xmm1    # 0xf5c48
ucomisd %xmm1, %xmm0
jbe 0xa07c8
movsd -0x8(%rdi,%rsi,8), %xmm1
leaq 0x7dfe9(%rip), %rsi     # 0x11e67f
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movb $0x2, %al
movq %r10, %rbx
callq 0x30c2c
leaq 0x509d3(%rip), %rsi     # 0xf1082
leaq 0x7e022(%rip), %rdx     # 0x11e6d8
leaq 0x840(%rsp), %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
addq $0x192, %rbx            # imm = 0x192
leaq 0x5095b(%rip), %rsi     # 0xf1034
leaq 0x7dcc5(%rip), %rdx     # 0x11e3a5
movq %rbx, %rdi
movl $0xbd, %ecx
jmp 0xa04f0
movslq 0x6c(%r10), %rax
testq %rax, %rax
jle 0xa01ac
movq 0x8(%rsp), %rcx
addq $0x130, %rcx            # imm = 0x130
movslq %r12d, %rdx
xorl %esi, %esi
leaq (,%rdx,8), %rdi
testl %r15d, %r15d
je 0xa0741
movq %rcx, %r8
cmpl $0x1, %r15d
jne 0xa0730
movq 0x8(%rsp), %r8
movq 0xb8(%r8), %r8
addq %rdi, %r8
movq (%r8), %r8
movsd (%r8,%rsi,8), %xmm0
movsd %xmm0, (%rbp,%rsi,8)
jmp 0xa07b7
movq 0x8(%rsp), %r8
movq 0xa8(%r8), %r8
movq (%r8,%rdx,8), %r8
movsd (%r8,%rsi,8), %xmm0
movsd %xmm0, (%rbp,%rsi,8)
cmpb $0x0, 0x10(%rsp)
je 0xa07b7
movq 0x8(%rsp), %r8
movslq 0x44(%r8), %r8
testq %r8, %r8
jle 0xa07b7
movq 0x8(%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 0xa07a3
incq %rsi
cmpq %rax, %rsi
jne 0xa0713
jmp 0xa01ac
movslq 0x6c(%r10), %rcx
testq %rcx, %rcx
jle 0xa01ac
movq 0x8(%rsp), %rdx
addq $0x130, %rdx            # imm = 0x130
movq %rdx, 0x18(%rsp)
movslq %r12d, %rsi
movl %ecx, %edi
imull %eax, %edi
movslq %edi, %rdi
leaq (,%rdi,8), %r8
xorl %r9d, %r9d
leaq (,%rsi,8), %r10
xorl %r11d, %r11d
movl %r14d, %edx
testl %r15d, %r15d
je 0xa0832
movq %rbp, %r13
cmpl $0x1, %r15d
jne 0xa08a8
movq 0x8(%rsp), %rbx
movq 0xb8(%rbx), %rbx
addq %r10, %rbx
movl %ecx, %ebp
jmp 0xa08b8
movq 0x8(%rsp), %rbx
movq 0xa8(%rbx), %rbx
movq (%rbx,%rsi,8), %rbx
leaq (%r11,%rdi), %r14
movsd (%rbx,%r14,8), %xmm0
movsd %xmm0, (%rbp,%r11,8)
cmpb $0x0, 0x10(%rsp)
je 0xa08d9
movq 0x8(%rsp), %rbx
movslq 0x44(%rbx), %rbx
testq %rbx, %rbx
jle 0xa08d9
movq 0x8(%rsp), %r14
movq 0x98(%r14), %r13
movq %r9, %r14
imulq %rbx, %r14
addq 0x1080(%rsp), %r14
movq %r8, %r12
imulq %rbx, %r12
addq (%r13,%rsi,8), %r12
xorl %r13d, %r13d
movsd (%r12,%r13,8), %xmm0
movsd %xmm0, (%r14,%r13,8)
incq %r13
cmpq %r13, %rbx
jne 0xa0892
jmp 0xa08d9
movq 0x8(%rsp), %rbx
movl 0x104(%rbx), %ebp
movq 0x18(%rsp), %rbx
imull %eax, %ebp
movslq %ebp, %r14
leaq (,%r11,8), %r12
addq (%rbx), %r12
movsd (%r12,%r14,8), %xmm0
movq %r13, %rbp
movsd %xmm0, (%r13,%r11,8)
incq %r11
addq $0x8, %r9
addq $0x8, %r8
cmpq %rcx, %r11
movl %edx, %r14d
jne 0xa080a
jmp 0xa01ac
movq 0xa8(%r10), %rax
movslq %r12d, %rbx
movq (%rax,%rbx,8), %rdx
movq 0xb0(%r10), %rax
movq (%rax,%rbx,8), %rcx
movl 0x6c(%r10), %r8d
leaq 0x192(%r10), %r15
subq $0x8, %rsp
leaq 0x1c(%rsp), %r9
movsd %xmm0, 0x20(%rsp)
pushq %r15
pushq %r8
pushq %rbp
callq 0x17ee6
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa09f4
leaq 0x50aaa(%rip), %rsi     # 0xf13f0
leaq 0x7ddbc(%rip), %rdx     # 0x11e709
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x506ce(%rip), %rsi     # 0xf1034
leaq 0x7da38(%rip), %rdx     # 0x11e3a5
movq %r15, %rdi
movl $0xe5, %ecx
jmp 0xa0b23
movq 0x130(%r10), %rdx
movq 0x138(%r10), %rcx
movl 0x6c(%r10), %eax
movl 0x104(%r10), %r8d
leaq 0x192(%r10), %r15
subq $0x8, %rsp
leaq 0x1c(%rsp), %r9
pushq %r15
pushq %rax
pushq %rbp
callq 0x17ee6
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa01ac
leaq 0x50a30(%rip), %rsi     # 0xf13f0
leaq 0x7e01c(%rip), %rdx     # 0x11e9e3
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x50654(%rip), %rsi     # 0xf1034
leaq 0x7d9be(%rip), %rdx     # 0x11e3a5
movq %r15, %rdi
movl $0x128, %ecx            # imm = 0x128
jmp 0xa0b23
cmpb $0x0, 0x10(%rsp)
je 0xa01ac
movq 0x8(%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 0x1c(%rsp), %r9
movsd 0x20(%rsp), %xmm0
pushq %r15
pushq %r8
pushq 0x1098(%rsp)
callq 0x17ee6
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa01ac
leaq 0x5098c(%rip), %rsi     # 0xf13f0
leaq 0x7dd53(%rip), %rdx     # 0x11e7be
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x505b0(%rip), %rsi     # 0xf1034
leaq 0x7d91a(%rip), %rdx     # 0x11e3a5
movq %r15, %rdi
movl $0xf4, %ecx
jmp 0xa0b23
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 0x1c(%rsp), %r9
movl %eax, %esi
pushq %r15
pushq %r8
pushq %rbp
callq 0x17ee6
addq $0x20, %rsp
cmpl $0x1, %eax
jne 0xa01ac
leaq 0x508fc(%rip), %rsi     # 0xf13f0
leaq 0x7dda1(%rip), %rdx     # 0x11e89c
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r15, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x50520(%rip), %rsi     # 0xf1034
leaq 0x7d88a(%rip), %rdx     # 0x11e3a5
movq %r15, %rdi
movl $0x117, %ecx            # imm = 0x117
movq %rbx, %r8
jmp 0xa04f3