// File: /workspace/github_c_makefile_star2_O1/PoulinV[P]AxiCLASS/build/../source/thermodynamics.c
// Function: thermodynamics_vector_init
int thermodynamics_vector_init(
                               struct precision * ppr,
                               struct background * pba,
                               struct thermodynamics * pth,
                               double mz,
                               struct thermo_workspace * ptw
                               ) {

  /** Summary: */

  /** Define local variables */
  int index_ti;
  /* ptdw->ptv unallocated if ap_current == index_ap_brec, allocated and filled otherwise */
  struct thermo_vector * ptv;
  struct thermo_diffeq_workspace * ptdw = ptw->ptdw;
  double z;

  class_alloc(ptv,sizeof(struct thermo_vector),pth->error_message);

  /* mz = Minus z is inverted*/
  z = -mz;

  /* Start from no component */
  index_ti = 0;

  /* Add common indices (Have to be added before) */
  class_define_index(ptv->index_ti_D_Tmat,_TRUE_,index_ti,1);

  /* Add all components that should be evolved */
  if (ptdw->ap_current == ptdw->index_ap_brec) {
    /* Nothing else to add */
  }
  else if (ptdw->ap_current == ptdw->index_ap_He1) {
    /* Nothing else to add */
  }
  else if (ptdw->ap_current == ptdw->index_ap_He1f) {
    /* Nothing else to add */
  }
  else if (ptdw->ap_current == ptdw->index_ap_He2) {
    /* Nothing else to add */
  }
  else if (ptdw->ap_current == ptdw->index_ap_H) {
    class_define_index(ptv->index_ti_x_He,_TRUE_,index_ti,1);
  }
  else if (ptdw->ap_current == ptdw->index_ap_frec) {
    class_define_index(ptv->index_ti_x_He,_TRUE_,index_ti,1);
    class_define_index(ptv->index_ti_x_H,_TRUE_,index_ti,1);
  }
  else if (ptdw->ap_current == ptdw->index_ap_reio) {
    class_define_index(ptv->index_ti_x_He,_TRUE_,index_ti,1);
    class_define_index(ptv->index_ti_x_H,_TRUE_,index_ti,1);
  }

  /* Only in the case where we in idm case but not idm TCA define T_idm */
  if (pba->has_idm == _TRUE_){
    class_define_index(ptv->index_ti_T_idm,(ptw->has_ap_idmtca == _FALSE_ || ptdw->ap_current != ptdw->index_ap_idmtca),index_ti,1);
  }

  /* We have now obtained the full size */
  ptv->ti_size = index_ti;

  /* Allocate all arrays used during the evolution */
  class_calloc(ptv->y,ptv->ti_size,sizeof(double),pth->error_message);
  class_alloc(ptv->dy,ptv->ti_size*sizeof(double),pth->error_message);
  class_alloc(ptv->used_in_output,ptv->ti_size*sizeof(int),pth->error_message);

  for (index_ti=0; index_ti<ptv->ti_size; index_ti++) {
    ptv->used_in_output[index_ti] = _TRUE_;
  }

  if (ptw->has_ap_idmtca == _TRUE_ && ptdw->ap_current == ptdw->index_ap_idmtca) {
    /* Store Tmat in workspace for later use */
    ptdw->Tmat = ptw->Tcmb*(1.+z);

    /* Set the new vector and its indices */
    ptdw->ptv = ptv;

    ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat] = 0.;

    /* set other initial conditions */
    ptdw->require_H = _FALSE_;
    ptdw->require_He = _FALSE_;
  }
  else if (ptw->has_ap_idmtca == _TRUE_ && ptdw->ap_current == ptdw->index_ap_brec) {
    /* Copy old values */
    ptv->y[ptv->index_ti_D_Tmat] = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat];
    ptv->dy[ptv->index_ti_D_Tmat] = ptdw->ptv->dy[ptdw->ptv->index_ti_D_Tmat];

    /* Set initial conditions */
    ptv->y[ptv->index_ti_T_idm] = pba->T_cmb * (1.+z);
    /* If we instead have idr tight coupling, we choose instead */
    if ((pth->has_idm_dr == _TRUE_) && pth->n_index_idm_dr > 0)
      ptv->y[ptv->index_ti_T_idm] = pba->T_idr * (1.+z);

    /* Free the old vector and its indices */
    class_call(thermodynamics_vector_free(ptdw->ptv),
               pth->error_message,
               pth->error_message);

    /* Copy the new vector into the position of the old one */
    ptdw->ptv = ptv;
    ptdw->require_H = _FALSE_;
    ptdw->require_He = _FALSE_;
  }
  else if (ptdw->ap_current == ptdw->index_ap_brec) {

    /* Store Tmat in workspace for later use */
    ptdw->Tmat = ptw->Tcmb*(1.+z);

    /* Set the new vector and its indices */
    ptdw->ptv = ptv;

    ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat] = 0.;

    ptdw->require_H = _FALSE_;
    ptdw->require_He = _FALSE_;

    /* Set idm initial temperature if not coupled */
    if (pba->has_idm == _TRUE_) {
      class_call(thermodynamics_idm_initial_temperature(pba,
                                                        pth,
                                                        z,
                                                        ptdw),
                 pth->error_message,
                 pth->error_message);
      ptv->y[ptv->index_ti_T_idm] = ptdw->T_idm;
    }
  }
  /* - in this scheme we start to evolve Helium and thus need to set its initial condition via the analytic function */
  else if (ptdw->ap_current == ptdw->index_ap_H) {
    /* Store Tmat in workspace for later use */
    ptdw->Tmat = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat] + ptw->Tcmb*(1.+z);

    /* Obtain initial contents of new vector analytically, especially x_He */
    class_call(thermodynamics_ionization_fractions(z,ptdw->ptv->y,pba,pth,ptw,ptdw->ap_current-1),
               pth->error_message,
               pth->error_message);

    /* Set the new vector and its indices */
    ptv->y[ptv->index_ti_D_Tmat] = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat];
    ptv->y[ptv->index_ti_x_He] = ptdw->x_He;
    if (pba->has_idm == _TRUE_)
      ptv->y[ptv->index_ti_T_idm] = ptdw->ptv->y[ptdw->ptv->index_ti_T_idm];

    /* Free the old vector and its indices */
    class_call(thermodynamics_vector_free(ptdw->ptv),
               pth->error_message,
               pth->error_message);

    /* Copy the new vector into the position of the old one*/
    ptdw->ptv = ptv;

    ptdw->require_H = _FALSE_;
    ptdw->require_He = _TRUE_;
  }
  /* - in the scheme of full recombination (=frec) we evolve all quantities and thus need to set their initial conditions.
     Tmat and x_He are solely taken from the previous scheme, x_H is set via the analytic function */
  else if (ptdw->ap_current == ptdw->index_ap_frec) {
    /* Store Tmat in workspace for later use */
    ptdw->Tmat = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat] + ptw->Tcmb*(1.+z);

    /* Obtain initial contents of new vector analytically, especially x_H */
    class_call(thermodynamics_ionization_fractions(z,ptdw->ptv->y,pba,pth,ptw,ptdw->ap_current-1),
               pth->error_message,
               pth->error_message);

    /* Set the new vector and its indices */
    ptv->y[ptv->index_ti_D_Tmat] = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat];
    ptv->y[ptv->index_ti_x_H] = ptdw->x_H;
    ptv->y[ptv->index_ti_x_He] = ptdw->ptv->y[ptdw->ptv->index_ti_x_He];
    if (pba->has_idm == _TRUE_)
      ptv->y[ptv->index_ti_T_idm] = ptdw->ptv->y[ptdw->ptv->index_ti_T_idm];

    /* Free the old vector and its indices */
    class_call(thermodynamics_vector_free(ptdw->ptv),
               pth->error_message,
               pth->error_message);

    /* Copy the new vector into the position of the old one*/

    ptdw->ptv = ptv;

    ptdw->require_H = _TRUE_;
    ptdw->require_He = _TRUE_;
  }
  /* - during reionization we continue to evolve all quantities. Now all three intial conditions are just taken from the previous scheme */
  else if (ptdw->ap_current == ptdw->index_ap_reio) {

    /* Set the new vector and its indices */
    ptv->y[ptv->index_ti_D_Tmat] = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat];
    ptv->y[ptv->index_ti_x_H] = ptdw->ptv->y[ptdw->ptv->index_ti_x_H];
    ptv->y[ptv->index_ti_x_He] = ptdw->ptv->y[ptdw->ptv->index_ti_x_He];
    if (pba->has_idm == _TRUE_)
      ptv->y[ptv->index_ti_T_idm] = ptdw->ptv->y[ptdw->ptv->index_ti_T_idm];

    /* Free the old vector and its indices */
    class_call(thermodynamics_vector_free(ptdw->ptv),
               pth->error_message,
               pth->error_message);

    /* Copy the new vector into the position of the old one*/

    ptdw->ptv = ptv;

    ptdw->require_H = _TRUE_;
    ptdw->require_He = _TRUE_;
  }
  /* - in all other approximations we only evolve Tmat and set its initial conditions from the previous scheme */
  else{
    /* Store Tmat in workspace for later use */
    ptdw->Tmat = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat] + ptw->Tcmb*(1.+z);

    /* Set the new vector and its indices */
    ptv->y[ptv->index_ti_D_Tmat] = ptdw->ptv->y[ptdw->ptv->index_ti_D_Tmat];
    if (pba->has_idm == _TRUE_)
      ptv->y[ptv->index_ti_T_idm] = ptdw->ptv->y[ptdw->ptv->index_ti_T_idm];

    /* Free the old vector and its indices */
    class_call(thermodynamics_vector_free(ptdw->ptv),
               pth->error_message,
               pth->error_message);

    /* Copy the new vector into the position of the old one*/
    ptdw->ptv = ptv;

    ptdw->require_H = _FALSE_;
    ptdw->require_He = _FALSE_;
  }

  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x838, %rsp            # imm = 0x838
movq %rcx, %r13
movdqa %xmm0, 0x20(%rsp)
movq %rdx, %r14
movq %rsi, %r12
movq 0x80(%rcx), %rbx
movl $0x30, %edi
callq 0x8460
testq %rax, %rax
je 0x6ddc3
movq %rax, %r15
movl $0x0, 0xc(%rax)
movl 0x90(%rbx), %ecx
movl $0x1, %eax
cmpl 0x74(%rbx), %ecx
je 0x6de2a
cmpl 0x78(%rbx), %ecx
je 0x6de2a
cmpl 0x7c(%rbx), %ecx
je 0x6de2a
cmpl 0x80(%rbx), %ecx
je 0x6de2a
cmpl 0x84(%rbx), %ecx
jne 0x6de07
movl $0x1, 0x8(%r15)
movl $0x2, %eax
jmp 0x6de2a
leaq 0x8323a(%rip), %rsi     # 0xf1004
leaq 0xa4c0f(%rip), %rdx     # 0x1129e0
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movl $0x30, %ecx
xorl %eax, %eax
callq 0x30c2c
addq $0x1064, %r14           # imm = 0x1064
leaq 0x83241(%rip), %rsi     # 0xf1034
leaq 0xa4bea(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x794, %ecx            # imm = 0x794
jmp 0x6e0d7
cmpl 0x88(%rbx), %ecx
je 0x6de17
cmpl 0x8c(%rbx), %ecx
jne 0x6de2a
movabsq $0x100000002, %rax      # imm = 0x100000002
movq %rax, 0x4(%r15)
movl $0x3, %eax
movzwl 0x4b6(%r12), %edx
cmpw $0x1, %dx
movq %r12, 0x10(%rsp)
movw %dx, 0xe(%rsp)
jne 0x6de56
cmpw $0x0, 0x60(%r13)
je 0x6de50
cmpl 0x70(%rbx), %ecx
je 0x6de56
movl %eax, 0x10(%r15)
incl %eax
movl %eax, (%r15)
movl %eax, %r12d
movl $0x8, %esi
movq %r12, %rdi
callq 0x8250
movq %rax, 0x18(%r15)
testq %rax, %rax
je 0x6e00f
movq %rax, %rbp
leal (,%r12,8), %edi
callq 0x8460
movq %rax, 0x20(%r15)
testq %rax, %rax
je 0x6e053
movq %rax, 0x18(%rsp)
leal (,%r12,4), %edi
callq 0x8460
movq %rax, 0x28(%r15)
testq %rax, %rax
je 0x6e096
movslq (%r15), %rdx
testq %rdx, %rdx
jle 0x6dfba
leal 0x3(%rdx), %ecx
decq %rdx
movq %rdx, %xmm0
pshufd $0x44, %xmm0, %xmm0     # xmm0 = xmm0[0,1,0,1]
shrl $0x2, %ecx
shlq $0x4, %rcx
movdqa 0x87b43(%rip), %xmm1    # 0xf5a20
movdqa 0x83d9b(%rip), %xmm2    # 0xf1c80
xorl %edx, %edx
movdqa 0x83da1(%rip), %xmm3    # 0xf1c90
pxor %xmm3, %xmm0
pcmpeqd %xmm4, %xmm4
movdqa 0x87b31(%rip), %xmm5    # 0xf5a30
movdqa %xmm2, %xmm6
pxor %xmm3, %xmm6
movdqa %xmm6, %xmm7
pcmpgtd %xmm0, %xmm7
pcmpeqd %xmm0, %xmm6
pshufd $0xf5, %xmm6, %xmm8     # xmm8 = xmm6[1,1,3,3]
pand %xmm7, %xmm8
pshufd $0xf5, %xmm7, %xmm6     # xmm6 = xmm7[1,1,3,3]
por %xmm8, %xmm6
movd %xmm6, %esi
notl %esi
testb $0x1, %sil
je 0x6df3b
movl $0x1, (%rax,%rdx)
pxor %xmm4, %xmm6
pextrw $0x4, %xmm6, %esi
testb $0x1, %sil
je 0x6df52
movl $0x1, 0x4(%rax,%rdx)
movdqa %xmm1, %xmm6
pxor %xmm3, %xmm6
movdqa %xmm6, %xmm7
pcmpgtd %xmm0, %xmm7
pcmpeqd %xmm0, %xmm6
pshufd $0xf5, %xmm6, %xmm8     # xmm8 = xmm6[1,1,3,3]
pand %xmm7, %xmm8
pshufd $0xf5, %xmm7, %xmm6     # xmm6 = xmm7[1,1,3,3]
por %xmm8, %xmm6
pxor %xmm4, %xmm6
pextrw $0x0, %xmm6, %esi
testb $0x1, %sil
je 0x6df92
movl $0x1, 0x8(%rax,%rdx)
pextrw $0x4, %xmm6, %esi
testb $0x1, %sil
je 0x6dfa5
movl $0x1, 0xc(%rax,%rdx)
paddq %xmm5, %xmm2
paddq %xmm5, %xmm1
addq $0x10, %rdx
cmpq %rdx, %rcx
jne 0x6deff
movzwl 0x60(%r13), %eax
cmpw $0x1, %ax
movapd 0x20(%rsp), %xmm2
jne 0x6e199
movl 0x90(%rbx), %ecx
cmpl 0x70(%rbx), %ecx
jne 0x6e0fc
movsd 0x8327a(%rip), %xmm0    # 0xf1260
subsd %xmm2, %xmm0
mulsd 0x38(%r13), %xmm0
movsd %xmm0, 0x28(%rbx)
movq %r15, 0xb8(%rbx)
xorl %r12d, %r12d
movq %r12, (%rbp)
movq %r12, 0xb0(%rbx)
jmp 0x6e0e7
leaq 0x82fee(%rip), %rsi     # 0xf1004
leaq 0xa49e2(%rip), %rdx     # 0x1129ff
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movl $0x8, %ecx
xorl %eax, %eax
callq 0x30c2c
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82ff5(%rip), %rsi     # 0xf1034
leaq 0xa499e(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x7c1, %ecx            # imm = 0x7C1
jmp 0x6e0d7
shll $0x3, %r12d
leaq 0x82fa6(%rip), %rsi     # 0xf1004
leaq 0xa49a1(%rip), %rdx     # 0x112a06
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movl %r12d, %ecx
xorl %eax, %eax
callq 0x30c2c
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82faf(%rip), %rsi     # 0xf1034
leaq 0xa4958(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x7c2, %ecx            # imm = 0x7C2
jmp 0x6e0d7
shll $0x2, %r12d
leaq 0x82f63(%rip), %rsi     # 0xf1004
leaq 0xa4966(%rip), %rdx     # 0x112a0e
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movl %r12d, %ecx
xorl %eax, %eax
callq 0x30c2c
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82f6c(%rip), %rsi     # 0xf1034
leaq 0xa4915(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x7c3, %ecx            # imm = 0x7C3
movq %rbx, %r8
xorl %eax, %eax
callq 0x30c2c
movl $0x1, %r12d
movl %r12d, %eax
addq $0x838, %rsp            # imm = 0x838
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
cmpw $0x1, %ax
jne 0x6e199
movl 0x90(%rbx), %eax
cmpl 0x74(%rbx), %eax
jne 0x6e199
movq 0xb8(%rbx), %rdi
movq 0x18(%rdi), %rax
movslq 0xc(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rbp)
movq 0x20(%rdi), %rax
movsd (%rax,%rcx,8), %xmm0
movq 0x18(%rsp), %rax
movsd %xmm0, (%rax)
movsd 0x83118(%rip), %xmm0    # 0xf1260
subsd %xmm2, %xmm0
movq 0x10(%rsp), %rax
movsd 0x18(%rax), %xmm1
mulsd %xmm0, %xmm1
movslq 0x10(%r15), %rax
movsd %xmm1, (%rbp,%rax,8)
cmpw $0x1, 0x38(%r14)
jne 0x6e4a7
movsd 0x1028(%r14), %xmm1
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm1
jbe 0x6e4a7
movq 0x10(%rsp), %rcx
mulsd 0x100(%rcx), %xmm0
jmp 0x6e4a1
movapd 0x8308f(%rip), %xmm0    # 0xf1230
xorpd %xmm2, %xmm0
movl 0x90(%rbx), %r8d
cmpl 0x74(%rbx), %r8d
jne 0x6e248
xorl %r12d, %r12d
movsd 0x8309f(%rip), %xmm1    # 0xf1260
subsd %xmm2, %xmm1
mulsd 0x38(%r13), %xmm1
movsd %xmm1, 0x28(%rbx)
movq %r15, 0xb8(%rbx)
movq %r12, (%rbp)
movq %r12, 0xb0(%rbx)
cmpw $0x1, 0xe(%rsp)
jne 0x6e0e7
movq 0x10(%rsp), %rdi
movq %r14, %rsi
movq %rbx, %rdx
callq 0x7259a
testl %eax, %eax
je 0x6e3e8
addq $0x1064, %r14           # imm = 0x1064
leaq 0x831dc(%rip), %rsi     # 0xf13f0
leaq 0xa4807(%rip), %rdx     # 0x112a22
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x82e00(%rip), %rsi     # 0xf1034
leaq 0xa47a9(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x7ff, %ecx            # imm = 0x7FF
jmp 0x6e0d7
cmpl 0x84(%rbx), %r8d
jne 0x6e2e0
movq 0xb8(%rbx), %rax
movq 0x18(%rax), %rdi
movslq 0xc(%rax), %rax
movsd 0x82ff4(%rip), %xmm1    # 0xf1260
subsd %xmm2, %xmm1
mulsd 0x38(%r13), %xmm1
addsd (%rdi,%rax,8), %xmm1
movsd %xmm1, 0x28(%rbx)
decl %r8d
movq 0x10(%rsp), %r12
movq %r12, %rsi
movq %r14, %rdx
movq %r13, %rcx
callq 0x72949
testl %eax, %eax
je 0x6e3ff
addq $0x1064, %r14           # imm = 0x1064
leaq 0x83144(%rip), %rsi     # 0xf13f0
leaq 0xa47a9(%rip), %rdx     # 0x112a5c
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x82d68(%rip), %rsi     # 0xf1034
leaq 0xa4711(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x80b, %ecx            # imm = 0x80B
jmp 0x6e0d7
cmpl 0x88(%rbx), %r8d
jne 0x6e378
movq 0xb8(%rbx), %rax
movq 0x18(%rax), %rdi
movslq 0xc(%rax), %rax
movsd 0x82f5c(%rip), %xmm1    # 0xf1260
subsd %xmm2, %xmm1
mulsd 0x38(%r13), %xmm1
addsd (%rdi,%rax,8), %xmm1
movsd %xmm1, 0x28(%rbx)
decl %r8d
movq 0x10(%rsp), %r12
movq %r12, %rsi
movq %r14, %rdx
movq %r13, %rcx
callq 0x72949
testl %eax, %eax
je 0x6e4c0
addq $0x1064, %r14           # imm = 0x1064
leaq 0x830ac(%rip), %rsi     # 0xf13f0
leaq 0xa4711(%rip), %rdx     # 0x112a5c
leaq 0x30(%rsp), %rbx
movq %rbx, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30c2c
leaq 0x82cd0(%rip), %rsi     # 0xf1034
leaq 0xa4679(%rip), %rdx     # 0x1129e4
movq %r14, %rdi
movl $0x827, %ecx            # imm = 0x827
jmp 0x6e0d7
movq 0xb8(%rbx), %rdi
movq 0x18(%rdi), %rax
movslq 0xc(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
cmpl 0x8c(%rbx), %r8d
jne 0x6e467
movsd %xmm0, (%rbp)
movslq 0x4(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movslq 0x4(%r15), %rcx
movsd %xmm0, (%rbp,%rcx,8)
movslq 0x8(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movslq 0x8(%r15), %rcx
movsd %xmm0, (%rbp,%rcx,8)
cmpw $0x1, 0xe(%rsp)
jne 0x6e51e
movslq 0x10(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movslq 0x10(%r15), %rax
movsd %xmm0, (%rbp,%rax,8)
jmp 0x6e51e
movsd 0x38(%rbx), %xmm0
movq 0x18(%r15), %rax
movslq 0x10(%r15), %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0x6e0e7
movq 0xb8(%rbx), %rdi
movq 0x18(%rdi), %rax
movslq 0xc(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x18(%r15), %rcx
movslq 0xc(%r15), %rdx
movsd %xmm0, (%rcx,%rdx,8)
movsd 0x8(%rbx), %xmm0
movslq 0x8(%r15), %rdx
movsd %xmm0, (%rcx,%rdx,8)
cmpw $0x1, 0x4b6(%r12)
jne 0x6e44c
movslq 0x10(%rdi), %rdx
movsd (%rax,%rdx,8), %xmm0
movslq 0x10(%r15), %rax
movsd %xmm0, (%rcx,%rax,8)
callq 0x70586
movq %r15, 0xb8(%rbx)
movabsq $0x100000000, %rax      # imm = 0x100000000
jmp 0x6e534
movsd 0x82df1(%rip), %xmm1    # 0xf1260
subsd %xmm2, %xmm1
mulsd 0x38(%r13), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x28(%rbx)
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%rbp)
cmpw $0x1, 0xe(%rsp)
jne 0x6e4a7
movslq 0x10(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movslq 0x10(%r15), %rax
movsd %xmm0, (%rbp,%rax,8)
callq 0x70586
movq %r15, 0xb8(%rbx)
movq $0x0, 0xb0(%rbx)
jmp 0x6e53b
movq 0xb8(%rbx), %rdi
movq 0x18(%rdi), %rax
movslq 0xc(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x18(%r15), %rcx
movslq 0xc(%r15), %rdx
movsd %xmm0, (%rcx,%rdx,8)
movsd (%rbx), %xmm0
movslq 0x4(%r15), %rdx
movsd %xmm0, (%rcx,%rdx,8)
movslq 0x8(%rdi), %rdx
movsd (%rax,%rdx,8), %xmm0
movslq 0x8(%r15), %rdx
movsd %xmm0, (%rcx,%rdx,8)
cmpw $0x1, 0x4b6(%r12)
jne 0x6e51e
movslq 0x10(%rdi), %rdx
movsd (%rax,%rdx,8), %xmm0
movslq 0x10(%r15), %rax
movsd %xmm0, (%rcx,%rax,8)
callq 0x70586
movq %r15, 0xb8(%rbx)
movabsq $0x100000001, %rax      # imm = 0x100000001
movq %rax, 0xb0(%rbx)
xorl %r12d, %r12d
jmp 0x6e0e7