// File: /workspace/github_c_makefile_star2_O0/PoulinV[P]AxiCLASS/build/../source/fourier.c
// Function: fourier_pks_at_kvec_and_zvec
int fourier_pks_at_kvec_and_zvec(
                                 struct background * pba,
                                 struct fourier * pfo,
                                 enum pk_outputs pk_output,
                                 double * kvec, // kvec[index_kvec]
                                 int kvec_size,
                                 double * zvec, // zvec[index_zvec]
                                 int zvec_size,
                                 double * out_pk,   // output_pk[index_zvec*kvec_size+index_kvec],
                                                    // already allocated
                                                    //(or NULL if user knows there is no _m output)
                                 double * out_pk_cb // output_pk[index_zvec*kvec_size+index_kvec],
                                                    //already allocated
                                                    //(or NULL if user knows there is no _cb output)
                                 ) {

  /** Summary: */

  /** - define local variables */

  int index_k, index_kvec, index_zvec;
  double * ln_kvec;
  double * ln_pk_table = NULL;
  double * ddln_pk_table = NULL;
  double * ln_pk_cb_table = NULL;
  double * ddln_pk_cb_table = NULL;
  double h, a, b;

  /** - Allocate arrays */

  class_alloc(ln_kvec, sizeof(double)*kvec_size,
              pfo->error_message);

  if (pfo->has_pk_m == _TRUE_) {
    class_alloc(ln_pk_table, sizeof(double)*pfo->k_size*zvec_size,
                pfo->error_message);
    class_alloc(ddln_pk_table, sizeof(double)*pfo->k_size*zvec_size,
                pfo->error_message);
  }
  if (pfo->has_pk_cb == _TRUE_) {
    class_alloc(ln_pk_cb_table, sizeof(double)*pfo->k_size*zvec_size,
                pfo->error_message);
    class_alloc(ddln_pk_cb_table, sizeof(double)*pfo->k_size*zvec_size,
                pfo->error_message);
  }

  /** - Construct table of log(P(k_n,z_j)) for pre-computed wavenumbers but requested redshifts: */

  for (index_zvec=0; index_zvec<zvec_size; index_zvec++){

    if (pfo->has_pk_m == _TRUE_) {
      class_call(fourier_pk_at_z(pba,
                                 pfo,
                                 logarithmic,
                                 pk_output,
                                 zvec[index_zvec],
                                 pfo->index_pk_m,
                                 &(ln_pk_table[index_zvec * pfo->k_size]),
                                 NULL),
                 pfo->error_message,
                 pfo->error_message);
    }
    if (pfo->has_pk_cb == _TRUE_) {
      class_call(fourier_pk_at_z(pba,
                                 pfo,
                                 logarithmic,
                                 pk_output,
                                 zvec[index_zvec],
                                 pfo->index_pk_cb,
                                 &(ln_pk_cb_table[index_zvec * pfo->k_size]),
                                 NULL),
                 pfo->error_message,
                 pfo->error_message);
    }
  }

  /** - Spline it for interpolation along k */

  if (pfo->has_pk_m == _TRUE_) {

    class_call(array_spline_table_columns2(pfo->ln_k,
                                           pfo->k_size,
                                           ln_pk_table,
                                           zvec_size,
                                           ddln_pk_table,
                                           _SPLINE_NATURAL_,
                                           pfo->error_message),
               pfo->error_message,
               pfo->error_message);
  }
  if (pfo->has_pk_cb == _TRUE_) {

    class_call(array_spline_table_columns2(pfo->ln_k,
                                           pfo->k_size,
                                           ln_pk_cb_table,
                                           zvec_size,
                                           ddln_pk_cb_table,
                                           _SPLINE_NATURAL_,
                                           pfo->error_message),
               pfo->error_message,
               pfo->error_message);
  }

  /** - Construct ln(kvec): */

  for (index_kvec=0; index_kvec<kvec_size; index_kvec++)
    ln_kvec[index_kvec] = log(kvec[index_kvec]);

  /** - Loop over first k values. If k<kmin, fill output with zeros. If not, go to next step. */

  for (index_kvec = 0; index_kvec<kvec_size; index_kvec++){

    /* check whether one should go to next step */
    if (ln_kvec[index_kvec] >= pfo->ln_k[0])
      break;

    /* deal with k<k_min */
    for (index_zvec = 0; index_zvec < zvec_size; index_zvec++) {
      if (pfo->has_pk_m == _TRUE_)     out_pk[index_zvec*kvec_size+index_kvec] = 0.;
      if (pfo->has_pk_cb == _TRUE_) out_pk_cb[index_zvec*kvec_size+index_kvec] = 0.;
      /* (If needed, one could add instead some extrapolation here) */
    }
  }

  /** - Deal with case kmin<=k<=kmax. For better performance, do not
      loop through kvec, but through pre-computed k values. */

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

    /** --> Loop through k_i's that fall in interval [k_n,k_n+1] */

    while ((index_kvec < kvec_size) && (ln_kvec[index_kvec] <= pfo->ln_k[index_k+1])){

      /** --> for each of them, perform spine interpolation */

      h = pfo->ln_k[index_k+1]-pfo->ln_k[index_k];
      b = (ln_kvec[index_kvec] - pfo->ln_k[index_k])/h;
      a = 1.-b;

      for (index_zvec = 0; index_zvec < zvec_size; index_zvec++) {

        if (pfo->has_pk_m == _TRUE_) {

          out_pk[index_zvec*kvec_size+index_kvec] =
            exp(
                array_spline_eval(ln_pk_table,
                                  ddln_pk_table,
                                  (index_zvec * pfo->k_size + index_k),
                                  (index_zvec * pfo->k_size + index_k+1),
                                  h,a,b)
                );
        }
        if (pfo->has_pk_cb == _TRUE_) {

          out_pk_cb[index_zvec*kvec_size+index_kvec] =
            exp(
                array_spline_eval(ln_pk_cb_table,
                                  ddln_pk_cb_table,
                                  (index_zvec * pfo->k_size + index_k),
                                  (index_zvec * pfo->k_size + index_k+1),
                                  h,a,b)
                );
        }
      }
      index_kvec++;
    }
  }

  /** - Loop over possible remaining k values with k > kmax, to fill output with zeros. */

  while (index_kvec < kvec_size) {

    for (index_zvec = 0; index_zvec < zvec_size; index_zvec++) {
      if (pfo->has_pk_m == _TRUE_)     out_pk[index_zvec*kvec_size+index_kvec] = 0.;
      if (pfo->has_pk_cb == _TRUE_) out_pk_cb[index_zvec*kvec_size+index_kvec] = 0.;
      /* (If needed, one could add instead some extrapolation here) */
    }
    index_kvec++;
  }

  free(ln_kvec);
  if (pfo->has_pk_m == _TRUE_) {
    free(ln_pk_table);
    free(ddln_pk_table);
  }
  if (pfo->has_pk_cb == _TRUE_) {
    free(ln_pk_cb_table);
    free(ddln_pk_cb_table);
  }

  return _SUCCESS_;
}

// Assembly:
subq $0x48e8, %rsp           # imm = 0x48E8
movq 0x4900(%rsp), %rax
movq 0x48f8(%rsp), %rax
movl 0x48f0(%rsp), %eax
movq %rdi, 0x48d8(%rsp)
movq %rsi, 0x48d0(%rsp)
movl %edx, 0x48cc(%rsp)
movq %rcx, 0x48c0(%rsp)
movl %r8d, 0x48bc(%rsp)
movq %r9, 0x48b0(%rsp)
movq $0x0, 0x4890(%rsp)
movq $0x0, 0x4888(%rsp)
movq $0x0, 0x4880(%rsp)
movq $0x0, 0x4878(%rsp)
movslq 0x48bc(%rsp), %rdi
shlq $0x3, %rdi
callq 0x8460
movq %rax, 0x4898(%rsp)
cmpq $0x0, 0x4898(%rsp)
jne 0x126839
movslq 0x48bc(%rsp), %rax
shlq $0x3, %rax
movl %eax, 0x485c(%rsp)
leaq 0x4050(%rsp), %rdi
movl 0x485c(%rsp), %ecx
leaq 0xab81a(%rip), %rsi     # 0x1d2004
leaq 0xd9b5c(%rip), %rdx     # 0x20034d
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x4050(%rsp), %r8
leaq 0xab81e(%rip), %rsi     # 0x1d2034
leaq 0xd9b38(%rip), %rdx     # 0x200355
movl $0x353, %ecx            # imm = 0x353
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x1269d0
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rdi
shlq $0x3, %rdi
movslq 0x48f0(%rsp), %rax
imulq %rax, %rdi
callq 0x8460
movq %rax, 0x4890(%rsp)
cmpq $0x0, 0x4890(%rsp)
jne 0x12690e
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rax
shlq $0x3, %rax
movslq 0x48f0(%rsp), %rcx
imulq %rcx, %rax
movl %eax, 0x404c(%rsp)
leaq 0x3840(%rsp), %rdi
movl 0x404c(%rsp), %ecx
leaq 0xab745(%rip), %rsi     # 0x1d2004
leaq 0xd9aae(%rip), %rdx     # 0x200374
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x3840(%rsp), %r8
leaq 0xab749(%rip), %rsi     # 0x1d2034
leaq 0xd9a63(%rip), %rdx     # 0x200355
movl $0x357, %ecx            # imm = 0x357
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rdi
shlq $0x3, %rdi
movslq 0x48f0(%rsp), %rax
imulq %rax, %rdi
callq 0x8460
movq %rax, 0x4888(%rsp)
cmpq $0x0, 0x4888(%rsp)
jne 0x1269ce
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rax
shlq $0x3, %rax
movslq 0x48f0(%rsp), %rcx
imulq %rcx, %rax
movl %eax, 0x383c(%rsp)
leaq 0x3030(%rsp), %rdi
movl 0x383c(%rsp), %ecx
leaq 0xab685(%rip), %rsi     # 0x1d2004
leaq 0xd99ec(%rip), %rdx     # 0x200372
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x3030(%rsp), %r8
leaq 0xab689(%rip), %rsi     # 0x1d2034
leaq 0xd99a3(%rip), %rdx     # 0x200355
movl $0x359, %ecx            # imm = 0x359
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
jmp 0x1269d0
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x126b67
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rdi
shlq $0x3, %rdi
movslq 0x48f0(%rsp), %rax
imulq %rax, %rdi
callq 0x8460
movq %rax, 0x4880(%rsp)
cmpq $0x0, 0x4880(%rsp)
jne 0x126aa5
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rax
shlq $0x3, %rax
movslq 0x48f0(%rsp), %rcx
imulq %rcx, %rax
movl %eax, 0x302c(%rsp)
leaq 0x2820(%rsp), %rdi
movl 0x302c(%rsp), %ecx
leaq 0xab5ae(%rip), %rsi     # 0x1d2004
leaq 0xd9925(%rip), %rdx     # 0x200382
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x2820(%rsp), %r8
leaq 0xab5b2(%rip), %rsi     # 0x1d2034
leaq 0xd98cc(%rip), %rdx     # 0x200355
movl $0x35d, %ecx            # imm = 0x35D
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rdi
shlq $0x3, %rdi
movslq 0x48f0(%rsp), %rax
imulq %rax, %rdi
callq 0x8460
movq %rax, 0x4878(%rsp)
cmpq $0x0, 0x4878(%rsp)
jne 0x126b65
movq 0x48d0(%rsp), %rax
movslq 0x6c(%rax), %rax
shlq $0x3, %rax
movslq 0x48f0(%rsp), %rcx
imulq %rcx, %rax
movl %eax, 0x281c(%rsp)
leaq 0x2010(%rsp), %rdi
movl 0x281c(%rsp), %ecx
leaq 0xab4ee(%rip), %rsi     # 0x1d2004
leaq 0xd9863(%rip), %rdx     # 0x200380
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x2010(%rsp), %r8
leaq 0xab4f2(%rip), %rsi     # 0x1d2034
leaq 0xd980c(%rip), %rdx     # 0x200355
movl $0x35f, %ecx            # imm = 0x35F
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
jmp 0x126b67
movl $0x0, 0x48a4(%rsp)
movl 0x48a4(%rsp), %eax
cmpl 0x48f0(%rsp), %eax
jge 0x126d90
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x126c7f
movq 0x48d8(%rsp), %rdi
movq 0x48d0(%rsp), %rsi
movl 0x48cc(%rsp), %ecx
movq 0x48b0(%rsp), %rax
movslq 0x48a4(%rsp), %rdx
movsd (%rax,%rdx,8), %xmm0
movq 0x48d0(%rsp), %rax
movl 0x54(%rax), %r8d
movq 0x4890(%rsp), %r9
movl 0x48a4(%rsp), %eax
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %eax
cltq
shlq $0x3, %rax
addq %rax, %r9
movl $0x1, %edx
xorl %eax, %eax
movq $0x0, (%rsp)
callq 0x1236f0
cmpl $0x1, %eax
jne 0x126c7d
leaq 0x1810(%rsp), %rdi
movq 0x48d0(%rsp), %rcx
addq $0x192, %rcx            # imm = 0x192
leaq 0xab7aa(%rip), %rsi     # 0x1d23d8
leaq 0xd975c(%rip), %rdx     # 0x200391
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x1810(%rsp), %r8
leaq 0xab3da(%rip), %rsi     # 0x1d2034
leaq 0xd96f4(%rip), %rdx     # 0x200355
movl $0x370, %ecx            # imm = 0x370
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
jmp 0x126c7f
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x126d78
movq 0x48d8(%rsp), %rdi
movq 0x48d0(%rsp), %rsi
movl 0x48cc(%rsp), %ecx
movq 0x48b0(%rsp), %rax
movslq 0x48a4(%rsp), %rdx
movsd (%rax,%rdx,8), %xmm0
movq 0x48d0(%rsp), %rax
movl 0x58(%rax), %r8d
movq 0x4880(%rsp), %r9
movl 0x48a4(%rsp), %eax
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %eax
cltq
shlq $0x3, %rax
addq %rax, %r9
movl $0x1, %edx
xorl %eax, %eax
movq $0x0, (%rsp)
callq 0x1236f0
cmpl $0x1, %eax
jne 0x126d76
leaq 0x1010(%rsp), %rdi
movq 0x48d0(%rsp), %rcx
addq $0x192, %rcx            # imm = 0x192
leaq 0xab6b1(%rip), %rsi     # 0x1d23d8
leaq 0xd96ee(%rip), %rdx     # 0x20041c
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x1010(%rsp), %r8
leaq 0xab2e1(%rip), %rsi     # 0x1d2034
leaq 0xd95fb(%rip), %rdx     # 0x200355
movl $0x37c, %ecx            # imm = 0x37C
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
jmp 0x126d78
jmp 0x126d7a
movl 0x48a4(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a4(%rsp)
jmp 0x126b72
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x126e64
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rdi
movq 0x48d0(%rsp), %rax
movl 0x6c(%rax), %esi
movq 0x4890(%rsp), %rdx
movl 0x48f0(%rsp), %ecx
movq 0x4888(%rsp), %r8
movq 0x48d0(%rsp), %rax
addq $0x192, %rax            # imm = 0x192
xorl %r9d, %r9d
movq %rax, (%rsp)
callq 0x20820
cmpl $0x1, %eax
jne 0x126e62
leaq 0x810(%rsp), %rdi
movq 0x48d0(%rsp), %rcx
addq $0x192, %rcx            # imm = 0x192
leaq 0xab5c5(%rip), %rsi     # 0x1d23d8
leaq 0xd9691(%rip), %rdx     # 0x2004ab
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x810(%rsp), %r8
leaq 0xab1f5(%rip), %rsi     # 0x1d2034
leaq 0xd950f(%rip), %rdx     # 0x200355
movl $0x38c, %ecx            # imm = 0x38C
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
jmp 0x126e64
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x126f32
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rdi
movq 0x48d0(%rsp), %rax
movl 0x6c(%rax), %esi
movq 0x4880(%rsp), %rdx
movl 0x48f0(%rsp), %ecx
movq 0x4878(%rsp), %r8
movq 0x48d0(%rsp), %rax
addq $0x192, %rax            # imm = 0x192
xorl %r9d, %r9d
movq %rax, (%rsp)
callq 0x20820
cmpl $0x1, %eax
jne 0x126f30
leaq 0x10(%rsp), %rdi
movq 0x48d0(%rsp), %rcx
addq $0x192, %rcx            # imm = 0x192
leaq 0xab4f4(%rip), %rsi     # 0x1d23d8
leaq 0xd9632(%rip), %rdx     # 0x20051d
movb $0x0, %al
callq 0x5b9c0
movq 0x48d0(%rsp), %rdi
addq $0x192, %rdi            # imm = 0x192
leaq 0x10(%rsp), %r8
leaq 0xab127(%rip), %rsi     # 0x1d2034
leaq 0xd9441(%rip), %rdx     # 0x200355
movl $0x398, %ecx            # imm = 0x398
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x48e4(%rsp)
jmp 0x1276eb
jmp 0x126f32
movl $0x0, 0x48a8(%rsp)
movl 0x48a8(%rsp), %eax
cmpl 0x48bc(%rsp), %eax
jge 0x126f8f
movq 0x48c0(%rsp), %rax
movslq 0x48a8(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
callq 0x82f0
movq 0x4898(%rsp), %rax
movslq 0x48a8(%rsp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movl 0x48a8(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a8(%rsp)
jmp 0x126f3d
movl $0x0, 0x48a8(%rsp)
movl 0x48a8(%rsp), %eax
cmpl 0x48bc(%rsp), %eax
jge 0x1270a0
movq 0x4898(%rsp), %rax
movslq 0x48a8(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rax
ucomisd (%rax), %xmm0
jb 0x126fdd
jmp 0x1270a0
movl $0x0, 0x48a4(%rsp)
movl 0x48a4(%rsp), %eax
cmpl 0x48f0(%rsp), %eax
jge 0x127088
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x127036
movq 0x48f8(%rsp), %rax
movl 0x48a4(%rsp), %ecx
imull 0x48bc(%rsp), %ecx
addl 0x48a8(%rsp), %ecx
movslq %ecx, %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x127070
movq 0x4900(%rsp), %rax
movl 0x48a4(%rsp), %ecx
imull 0x48bc(%rsp), %ecx
addl 0x48a8(%rsp), %ecx
movslq %ecx, %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
jmp 0x127072
movl 0x48a4(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a4(%rsp)
jmp 0x126fe8
jmp 0x12708a
movl 0x48a8(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a8(%rsp)
jmp 0x126f9a
movl $0x0, 0x48ac(%rsp)
movl 0x48ac(%rsp), %eax
movq 0x48d0(%rsp), %rcx
movl 0x6c(%rcx), %ecx
subl $0x1, %ecx
cmpl %ecx, %eax
jge 0x1275a6
jmp 0x1270ca
movl 0x48a8(%rsp), %ecx
xorl %eax, %eax
cmpl 0x48bc(%rsp), %ecx
movb %al, 0xf(%rsp)
jge 0x127121
movq 0x4898(%rsp), %rax
movslq 0x48a8(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm1
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rax
movl 0x48ac(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
ucomisd %xmm1, %xmm0
setae %al
movb %al, 0xf(%rsp)
movb 0xf(%rsp), %al
testb $0x1, %al
jne 0x12712e
jmp 0x12758e
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rax
movl 0x48ac(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rax
movslq 0x48ac(%rsp), %rcx
subsd (%rax,%rcx,8), %xmm0
movsd %xmm0, 0x4870(%rsp)
movq 0x4898(%rsp), %rax
movslq 0x48a8(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x48d0(%rsp), %rax
movq 0x80(%rax), %rax
movslq 0x48ac(%rsp), %rcx
subsd (%rax,%rcx,8), %xmm0
divsd 0x4870(%rsp), %xmm0
movsd %xmm0, 0x4860(%rsp)
movsd 0xab089(%rip), %xmm0    # 0x1d2248
subsd 0x4860(%rsp), %xmm0
movsd %xmm0, 0x4868(%rsp)
movl $0x0, 0x48a4(%rsp)
movl 0x48a4(%rsp), %eax
cmpl 0x48f0(%rsp), %eax
jge 0x127578
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x1273a8
movsd 0x4868(%rsp), %xmm0
movq 0x4890(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movsd 0x4860(%rsp), %xmm1
movq 0x4890(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
mulsd (%rax,%rcx,8), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x4868(%rsp), %xmm1
mulsd 0x4868(%rsp), %xmm1
movsd 0x4868(%rsp), %xmm3
movsd 0x4868(%rsp), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
movq 0x4888(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm3
movsd 0x4860(%rsp), %xmm2
mulsd 0x4860(%rsp), %xmm2
movsd 0x4860(%rsp), %xmm5
movsd 0x4860(%rsp), %xmm4
movq %xmm4, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm4
mulsd %xmm5, %xmm2
addsd %xmm4, %xmm2
movq 0x4888(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
mulsd (%rax,%rcx,8), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd 0x4870(%rsp), %xmm1
mulsd 0x4870(%rsp), %xmm1
movsd 0xac313(%rip), %xmm2    # 0x1d3688
divsd %xmm2, %xmm1
addsd %xmm1, %xmm0
callq 0x8670
movq 0x48f8(%rsp), %rax
movl 0x48a4(%rsp), %ecx
imull 0x48bc(%rsp), %ecx
addl 0x48a8(%rsp), %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x127560
movsd 0x4868(%rsp), %xmm0
movq 0x4880(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm2
movsd 0x4860(%rsp), %xmm1
movq 0x4880(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
mulsd (%rax,%rcx,8), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd 0x4868(%rsp), %xmm1
mulsd 0x4868(%rsp), %xmm1
movsd 0x4868(%rsp), %xmm3
movsd 0x4868(%rsp), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
movq 0x4878(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
movslq %ecx, %rcx
movsd (%rax,%rcx,8), %xmm3
movsd 0x4860(%rsp), %xmm2
mulsd 0x4860(%rsp), %xmm2
movsd 0x4860(%rsp), %xmm5
movsd 0x4860(%rsp), %xmm4
movq %xmm4, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm4
mulsd %xmm5, %xmm2
addsd %xmm4, %xmm2
movq 0x4878(%rsp), %rax
movl 0x48a4(%rsp), %ecx
movq 0x48d0(%rsp), %rdx
imull 0x6c(%rdx), %ecx
addl 0x48ac(%rsp), %ecx
addl $0x1, %ecx
movslq %ecx, %rcx
mulsd (%rax,%rcx,8), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd 0x4870(%rsp), %xmm1
mulsd 0x4870(%rsp), %xmm1
movsd 0xac15b(%rip), %xmm2    # 0x1d3688
divsd %xmm2, %xmm1
addsd %xmm1, %xmm0
callq 0x8670
movq 0x4900(%rsp), %rax
movl 0x48a4(%rsp), %ecx
imull 0x48bc(%rsp), %ecx
addl 0x48a8(%rsp), %ecx
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0x127562
movl 0x48a4(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a4(%rsp)
jmp 0x1271dc
movl 0x48a8(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a8(%rsp)
jmp 0x1270ca
jmp 0x127590
movl 0x48ac(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48ac(%rsp)
jmp 0x1270ab
jmp 0x1275a8
movl 0x48a8(%rsp), %eax
cmpl 0x48bc(%rsp), %eax
jge 0x12767d
movl $0x0, 0x48a4(%rsp)
movl 0x48a4(%rsp), %eax
cmpl 0x48f0(%rsp), %eax
jge 0x127667
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x127615
movq 0x48f8(%rsp), %rax
movl 0x48a4(%rsp), %ecx
imull 0x48bc(%rsp), %ecx
addl 0x48a8(%rsp), %ecx
movslq %ecx, %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x12764f
movq 0x4900(%rsp), %rax
movl 0x48a4(%rsp), %ecx
imull 0x48bc(%rsp), %ecx
addl 0x48a8(%rsp), %ecx
movslq %ecx, %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
jmp 0x127651
movl 0x48a4(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a4(%rsp)
jmp 0x1275c7
movl 0x48a8(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x48a8(%rsp)
jmp 0x1275a8
movq 0x4898(%rsp), %rdi
callq 0x8420
movq 0x48d0(%rsp), %rax
movswl 0x50(%rax), %eax
cmpl $0x1, %eax
jne 0x1276b5
movq 0x4890(%rsp), %rdi
callq 0x8420
movq 0x4888(%rsp), %rdi
callq 0x8420
movq 0x48d0(%rsp), %rax
movswl 0x52(%rax), %eax
cmpl $0x1, %eax
jne 0x1276e0
movq 0x4880(%rsp), %rdi
callq 0x8420
movq 0x4878(%rsp), %rdi
callq 0x8420
movl $0x0, 0x48e4(%rsp)
movl 0x48e4(%rsp), %eax
addq $0x48e8, %rsp           # imm = 0x48E8
retq
nopw (%rax,%rax)