// File: /workspace/github_c_makefile_star2_O3/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 $0x848, %rsp            # imm = 0x848
movq %rcx, %r13
movapd %xmm0, 0x30(%rsp)
movq %rdx, %r14
movq %rsi, %r12
movq 0x80(%rcx), %rbx
movl $0x30, %edi
callq 0x8460
testq %rax, %rax
je 0x6b461
movq %rax, %r15
movl $0x0, 0xc(%rax)
movl 0x74(%rbx), %eax
movl 0x90(%rbx), %ebp
movl $0x1, %ecx
movl %eax, 0x1c(%rsp)
cmpl %eax, %ebp
je 0x6b4c8
cmpl 0x78(%rbx), %ebp
je 0x6b4c8
cmpl 0x7c(%rbx), %ebp
je 0x6b4c8
cmpl 0x80(%rbx), %ebp
je 0x6b4c8
cmpl 0x84(%rbx), %ebp
jne 0x6b4a5
movl $0x1, 0x8(%r15)
movl $0x2, %ecx
jmp 0x6b4c8
leaq 0x82b9c(%rip), %rsi     # 0xee004
leaq 0xa45c1(%rip), %rdx     # 0x10fa30
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movl $0x30, %ecx
xorl %eax, %eax
callq 0x30af0
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82ba3(%rip), %rsi     # 0xee034
leaq 0xa459c(%rip), %rdx     # 0x10fa34
movq %r14, %rdi
movl $0x794, %ecx            # imm = 0x794
jmp 0x6b685
cmpl 0x88(%rbx), %ebp
je 0x6b4b5
cmpl 0x8c(%rbx), %ebp
jne 0x6b4c8
movabsq $0x100000002, %rax      # imm = 0x100000002
movq %rax, 0x4(%r15)
movl $0x3, %ecx
movq %r12, 0x10(%rsp)
movzwl 0x4b6(%r12), %edx
movl %ecx, %eax
movw %dx, 0x6(%rsp)
cmpw $0x1, %dx
jne 0x6b4f9
cmpw $0x0, 0x60(%r13)
je 0x6b4f2
movl %ecx, %eax
cmpl 0x70(%rbx), %ebp
je 0x6b4f9
movl %ecx, 0x10(%r15)
leal 0x1(%rcx), %eax
movq %rcx, 0x20(%rsp)
movq %rbx, 0x28(%rsp)
movl %eax, (%r15)
movl %eax, %ebx
movl $0x8, %esi
movq %rbx, %rdi
callq 0x8250
movq %rax, 0x18(%r15)
testq %rax, %rax
je 0x6b646
movq %rax, %r12
movq %r14, 0x8(%rsp)
leal (,%rbx,8), %edi
callq 0x8460
movq %rax, 0x20(%r15)
testq %rax, %rax
je 0x6b68d
movq %rax, %r14
leal (,%rbx,4), %edi
callq 0x8460
movq %rax, 0x28(%r15)
testq %rax, %rax
je 0x6b6d0
decq %rbx
movq %rbx, %xmm0
pshufd $0x44, %xmm0, %xmm1     # xmm1 = xmm0[0,1,0,1]
xorl %ecx, %ecx
testb %cl, %cl
jne 0x6b578
movl $0x1, (%rax)
pxor 0x83710(%rip), %xmm1    # 0xeec90
movdqa 0xa2228(%rip), %xmm0    # 0x10d7b0
pshufd $0xf5, %xmm0, %xmm4     # xmm4 = xmm0[1,1,3,3]
pcmpgtd %xmm1, %xmm0
pshufd $0xf5, %xmm1, %xmm2     # xmm2 = xmm1[1,1,3,3]
pcmpeqd %xmm2, %xmm4
pand %xmm0, %xmm4
pcmpeqd %xmm3, %xmm3
pxor %xmm3, %xmm4
pextrw $0x4, %xmm4, %ecx
testb $0x1, %cl
movq 0x28(%rsp), %rbx
movapd 0x30(%rsp), %xmm6
je 0x6b5c2
movl $0x1, 0x4(%rax)
movapd 0x82c66(%rip), %xmm0    # 0xee230
xorpd %xmm6, %xmm0
movdqa 0xa21ea(%rip), %xmm5    # 0x10d7c0
pshufd $0xf5, %xmm5, %xmm4     # xmm4 = xmm5[1,1,3,3]
pcmpgtd %xmm1, %xmm5
pcmpeqd %xmm2, %xmm4
pand %xmm5, %xmm4
pxor %xmm4, %xmm3
pextrw $0x0, %xmm3, %ecx
testb $0x1, %cl
je 0x6b5fc
movl $0x1, 0x8(%rax)
pcmpeqd %xmm1, %xmm1
pxor %xmm1, %xmm4
pextrw $0x4, %xmm4, %ecx
testb $0x1, %cl
je 0x6b615
movl $0x1, 0xc(%rax)
cmpw $0x1, 0x60(%r13)
jne 0x6b734
cmpl 0x70(%rbx), %ebp
jne 0x6b7d0
movsd 0x82c2e(%rip), %xmm0    # 0xee260
subsd %xmm6, %xmm0
mulsd 0x38(%r13), %xmm0
movsd %xmm0, 0x28(%rbx)
jmp 0x6bac9
leaq 0x829b7(%rip), %rsi     # 0xee004
leaq 0xa43fb(%rip), %rdx     # 0x10fa4f
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movl $0x8, %ecx
xorl %eax, %eax
callq 0x30af0
addq $0x1064, %r14           # imm = 0x1064
leaq 0x829be(%rip), %rsi     # 0xee034
leaq 0xa43b7(%rip), %rdx     # 0x10fa34
movq %r14, %rdi
movl $0x7c1, %ecx            # imm = 0x7C1
movq %rbx, %r8
jmp 0x6b714
shll $0x3, %ebx
leaq 0x8296d(%rip), %rsi     # 0xee004
leaq 0xa43b8(%rip), %rdx     # 0x10fa56
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movl %ebx, %ecx
xorl %eax, %eax
callq 0x30af0
movq 0x8(%rsp), %rdi
addq $0x1064, %rdi           # imm = 0x1064
leaq 0x82972(%rip), %rsi     # 0xee034
leaq 0xa436b(%rip), %rdx     # 0x10fa34
movl $0x7c2, %ecx            # imm = 0x7C2
jmp 0x6b711
shll $0x2, %ebx
leaq 0x8292a(%rip), %rsi     # 0xee004
leaq 0xa437d(%rip), %rdx     # 0x10fa5e
leaq 0x40(%rsp), %r14
movq %r14, %rdi
movl %ebx, %ecx
xorl %eax, %eax
callq 0x30af0
movq 0x8(%rsp), %rdi
addq $0x1064, %rdi           # imm = 0x1064
leaq 0x8292f(%rip), %rsi     # 0xee034
leaq 0xa4328(%rip), %rdx     # 0x10fa34
movl $0x7c3, %ecx            # imm = 0x7C3
movq %r14, %r8
xorl %eax, %eax
callq 0x30af0
movl $0x1, %ebp
movl %ebp, %eax
addq $0x848, %rsp            # imm = 0x848
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
cmpl 0x1c(%rsp), %ebp
jne 0x6b85e
xorl %ebp, %ebp
movsd 0x82b18(%rip), %xmm1    # 0xee260
subsd %xmm6, %xmm1
mulsd 0x38(%r13), %xmm1
movsd %xmm1, 0x28(%rbx)
movq %r15, 0xb8(%rbx)
movq $0x0, 0xb0(%rbx)
cmpw $0x1, 0x6(%rsp)
jne 0x6b720
movq 0x10(%rsp), %rdi
movq 0x8(%rsp), %r14
movq %r14, %rsi
movq %rbx, %rdx
callq 0x6faac
testl %eax, %eax
je 0x6ba6c
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82c54(%rip), %rsi     # 0xee3f0
leaq 0xa42cf(%rip), %rdx     # 0x10fa72
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x82878(%rip), %rsi     # 0xee034
leaq 0xa4271(%rip), %rdx     # 0x10fa34
movq %r14, %rdi
movl $0x7ff, %ecx            # imm = 0x7FF
jmp 0x6b685
cmpl 0x1c(%rsp), %ebp
jne 0x6b85e
movq 0xb8(%rbx), %rdi
movq 0x18(%rdi), %rax
movslq 0xc(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%r12)
movq 0x20(%rdi), %rax
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%r14)
movsd 0x82a56(%rip), %xmm0    # 0xee260
subsd %xmm6, %xmm0
movq 0x10(%rsp), %rax
movsd 0x18(%rax), %xmm1
mulsd %xmm0, %xmm1
movl 0x20(%rsp), %eax
movsd %xmm1, (%r12,%rax,8)
movq 0x8(%rsp), %rcx
cmpw $0x1, 0x38(%rcx)
jne 0x6bac4
movsd 0x1028(%rcx), %xmm1
pxor %xmm2, %xmm2
ucomisd %xmm2, %xmm1
jbe 0x6bac4
movq 0x10(%rsp), %rcx
mulsd 0x100(%rcx), %xmm0
jmp 0x6babe
cmpl 0x84(%rbx), %ebp
jne 0x6b8fc
movq 0xb8(%rbx), %rax
movq 0x18(%rax), %rdi
movslq 0xc(%rax), %rax
movsd 0x829df(%rip), %xmm1    # 0xee260
subsd %xmm6, %xmm1
mulsd 0x38(%r13), %xmm1
addsd (%rdi,%rax,8), %xmm1
movsd %xmm1, 0x28(%rbx)
decl %ebp
movq 0x10(%rsp), %r12
movq %r12, %rsi
movq 0x8(%rsp), %r14
movq %r14, %rdx
movq %r13, %rcx
movl %ebp, %r8d
callq 0x6fe5b
testl %eax, %eax
je 0x6ba04
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82b28(%rip), %rsi     # 0xee3f0
leaq 0xa41dd(%rip), %rdx     # 0x10faac
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x8274c(%rip), %rsi     # 0xee034
leaq 0xa4145(%rip), %rdx     # 0x10fa34
movq %r14, %rdi
movl $0x80b, %ecx            # imm = 0x80B
jmp 0x6b685
cmpl 0x88(%rbx), %ebp
movq 0x8(%rsp), %r14
jne 0x6b99a
movq 0xb8(%rbx), %rax
movq 0x18(%rax), %rdi
movslq 0xc(%rax), %rax
movsd 0x8293c(%rip), %xmm1    # 0xee260
subsd %xmm6, %xmm1
mulsd 0x38(%r13), %xmm1
addsd (%rdi,%rax,8), %xmm1
movsd %xmm1, 0x28(%rbx)
decl %ebp
movq 0x10(%rsp), %r12
movq %r12, %rsi
movq %r14, %rdx
movq %r13, %rcx
movl %ebp, %r8d
callq 0x6fe5b
testl %eax, %eax
je 0x6badd
addq $0x1064, %r14           # imm = 0x1064
leaq 0x82a8a(%rip), %rsi     # 0xee3f0
leaq 0xa413f(%rip), %rdx     # 0x10faac
leaq 0x40(%rsp), %rbx
movq %rbx, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x30af0
leaq 0x826ae(%rip), %rsi     # 0xee034
leaq 0xa40a7(%rip), %rdx     # 0x10fa34
movq %r14, %rdi
movl $0x827, %ecx            # imm = 0x827
jmp 0x6b685
movq 0xb8(%rbx), %rdi
movq 0x18(%rdi), %rax
movslq 0xc(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
cmpl 0x8c(%rbx), %ebp
jne 0x6ba83
movsd %xmm0, (%r12)
movslq 0x4(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, 0x10(%r12)
movslq 0x8(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, 0x8(%r12)
cmpw $0x1, 0x6(%rsp)
jne 0x6bb3b
movslq 0x10(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movl 0x20(%rsp), %eax
movsd %xmm0, (%r12,%rax,8)
jmp 0x6bb3b
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 0x6ba51
movslq 0x10(%rdi), %rdx
movsd (%rax,%rdx,8), %xmm0
movslq 0x10(%r15), %rax
movsd %xmm0, (%rcx,%rax,8)
callq 0x6dafb
movq %r15, 0xb8(%rbx)
movabsq $0x100000000, %rax      # imm = 0x100000000
jmp 0x6bb51
movsd 0x38(%rbx), %xmm0
movq 0x18(%r15), %rax
movslq 0x10(%r15), %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0x6b720
movsd 0x827d5(%rip), %xmm1    # 0xee260
subsd %xmm6, %xmm1
mulsd 0x38(%r13), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x28(%rbx)
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, (%r12)
cmpw $0x1, 0x6(%rsp)
jne 0x6bac4
movslq 0x10(%rdi), %rcx
movsd (%rax,%rcx,8), %xmm0
movl 0x20(%rsp), %eax
movsd %xmm0, (%r12,%rax,8)
callq 0x6dafb
movq %r15, 0xb8(%rbx)
movq $0x0, 0xb0(%rbx)
jmp 0x6bb58
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 0x6bb3b
movslq 0x10(%rdi), %rdx
movsd (%rax,%rdx,8), %xmm0
movslq 0x10(%r15), %rax
movsd %xmm0, (%rcx,%rax,8)
callq 0x6dafb
movq %r15, 0xb8(%rbx)
movabsq $0x100000001, %rax      # imm = 0x100000001
movq %rax, 0xb0(%rbx)
xorl %ebp, %ebp
jmp 0x6b720