// File: /workspace/github_c_makefile_star2_O0/PoulinV[P]AxiCLASS/build/../source/thermodynamics.c
// Function: thermodynamics_output_data
int thermodynamics_output_data(
                               struct background * pba,
                               struct thermodynamics *pth,
                               int number_of_titles,
                               double *data
                               ) {

  int index_z, storeidx;
  double *dataptr, *pvecthermo;
  double z,tau;

  // pth->number_of_thermodynamics_titles = get_number_of_titles(pth->thermodynamics_titles);
  // pth->size_thermodynamics_data = pth->number_of_thermodynamics_titles*pth->tt_size;

  /* Store quantities: */
  for (index_z=0; index_z<pth->tt_size; index_z++) {
    dataptr = data + index_z*number_of_titles;
    pvecthermo = pth->thermodynamics_table+index_z*pth->th_size;
    z = pth->z_table[index_z];
    storeidx=0;

    class_call(background_tau_of_z(pba, z, &tau),
               pba->error_message,
               pth->error_message);

    class_store_double(dataptr,1./(1.+z),_TRUE_,storeidx);
    class_store_double(dataptr,z,_TRUE_,storeidx);
    class_store_double(dataptr,tau,_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_xe],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_dkappa],_TRUE_,storeidx);
    //class_store_double(dataptr,pvecthermo[pth->index_th_ddkappa],_TRUE_,storeidx);
    //class_store_double(dataptr,pvecthermo[pth->index_th_dddkappa],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_exp_m_kappa],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_g],_TRUE_,storeidx);
    //class_store_double(dataptr,pvecthermo[pth->index_th_dg],_TRUE_,storeidx);
    //class_store_double(dataptr,pvecthermo[pth->index_th_ddg],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_Tb],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_dTb],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_wb],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_cb2],_TRUE_,storeidx);
    if (pba->has_idm == _TRUE_) {
      class_store_double(dataptr,pvecthermo[pth->index_th_T_idm],_TRUE_,storeidx);
      class_store_double(dataptr,pvecthermo[pth->index_th_c2_idm],_TRUE_,storeidx);
      if (pth->has_idm_g == _TRUE_) {
        class_store_double(dataptr,pvecthermo[pth->index_th_dmu_idm_g],_TRUE_,storeidx);
        class_store_double(dataptr,pvecthermo[pth->index_th_ddmu_idm_g],_TRUE_,storeidx);
      }
      if (pth->has_idm_b == _TRUE_){
        class_store_double(dataptr,pvecthermo[pth->index_th_R_idm_b],_TRUE_,storeidx);
      }
      if (pth->has_idm_dr == _TRUE_){
        class_store_double(dataptr,pvecthermo[pth->index_th_dmu_idm_dr],_TRUE_,storeidx);
        class_store_double(dataptr,pvecthermo[pth->index_th_tau_idm_dr],_TRUE_,storeidx);
        class_store_double(dataptr,pvecthermo[pth->index_th_tau_idr],_TRUE_,storeidx);
        class_store_double(dataptr,pvecthermo[pth->index_th_g_idm_dr],_TRUE_,storeidx);
      }
    }
    if (pba->has_idr == _TRUE_) {
      class_store_double(dataptr, pvecthermo[pth->index_th_T_idr],_TRUE_,storeidx);
      class_store_double(dataptr,pvecthermo[pth->index_th_dmu_idr],_TRUE_,storeidx);
    }
    class_store_double(dataptr,pvecthermo[pth->index_th_tau_d],_TRUE_,storeidx);
    //class_store_double(dataptr,pvecthermo[pth->index_th_rate],_TRUE_,storeidx);
    class_store_double(dataptr,pvecthermo[pth->index_th_r_d],pth->compute_damping_scale,storeidx);

  }

  return _SUCCESS_;
}

// Assembly:
subq $0x858, %rsp            # imm = 0x858
movq %rdi, 0x848(%rsp)
movq %rsi, 0x840(%rsp)
movl %edx, 0x83c(%rsp)
movq %rcx, 0x830(%rsp)
movl $0x0, 0x82c(%rsp)
movl 0x82c(%rsp), %eax
movq 0x840(%rsp), %rcx
cmpl 0xf08(%rcx), %eax
jge 0xc7703
movq 0x830(%rsp), %rax
movl 0x82c(%rsp), %ecx
imull 0x83c(%rsp), %ecx
movslq %ecx, %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, 0x820(%rsp)
movq 0x840(%rsp), %rax
movq 0xf20(%rax), %rax
movl 0x82c(%rsp), %ecx
movq 0x840(%rsp), %rdx
imull 0xf04(%rdx), %ecx
movslq %ecx, %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, 0x818(%rsp)
movq 0x840(%rsp), %rax
movq 0xf10(%rax), %rax
movslq 0x82c(%rsp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, 0x810(%rsp)
movl $0x0, 0x828(%rsp)
movq 0x848(%rsp), %rdi
movsd 0x810(%rsp), %xmm0
leaq 0x808(%rsp), %rsi
callq 0x9ee00
cmpl $0x1, %eax
jne 0xc70b2
movq %rsp, %rdi
movq 0x848(%rsp), %rcx
addq $0xd1c, %rcx            # imm = 0xD1C
leaq 0x10b370(%rip), %rsi    # 0x1d23d8
leaq 0x12ecf6(%rip), %rdx    # 0x1f5d65
movb $0x0, %al
callq 0x5b9c0
movq 0x840(%rsp), %rdi
addq $0x1064, %rdi           # imm = 0x1064
movq %rsp, %r8
leaq 0x10afa5(%rip), %rsi    # 0x1d2034
leaq 0x12ecf1(%rip), %rdx    # 0x1f5d87
movl $0x115c, %ecx           # imm = 0x115C
movb $0x0, %al
callq 0x5b9c0
movl $0x1, 0x854(%rsp)
jmp 0xc770e
movsd 0x10b18e(%rip), %xmm1   # 0x1d2248
addsd 0x810(%rsp), %xmm1
movsd 0x10b17d(%rip), %xmm0   # 0x1d2248
divsd %xmm1, %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movsd 0x810(%rsp), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movsd 0x808(%rsp), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xe7c(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xe80(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xe90(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xe94(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xee4(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xee8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xeec(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xef0(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x848(%rsp), %rax
movswl 0x4b6(%rax), %eax
cmpl $0x1, %eax
jne 0xc75ca
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xea0(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xea4(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x840(%rsp), %rax
movswl 0x36(%rax), %eax
cmpl $0x1, %eax
jne 0xc7467
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xec8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xecc(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x840(%rsp), %rax
movswl 0x34(%rax), %eax
cmpl $0x1, %eax
jne 0xc74b7
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xed8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x840(%rsp), %rax
movswl 0x38(%rax), %eax
cmpl $0x1, %eax
jne 0xc75c8
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xeac(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xebc(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xec0(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xec4(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0xc75ca
movq 0x848(%rsp), %rax
movswl 0x4c6(%rax), %eax
cmpl $0x1, %eax
jne 0xc765c
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xea8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xeb8(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xe84(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
movq 0x840(%rsp), %rax
movswl 0x32(%rax), %eax
cmpl $0x1, %eax
jne 0xc76eb
movq 0x818(%rsp), %rax
movq 0x840(%rsp), %rcx
movslq 0xf00(%rcx), %rcx
movsd (%rax,%rcx,8), %xmm0
movq 0x820(%rsp), %rax
movl 0x828(%rsp), %ecx
movl %ecx, %edx
addl $0x1, %edx
movl %edx, 0x828(%rsp)
movslq %ecx, %rcx
movsd %xmm0, (%rax,%rcx,8)
jmp 0xc76ed
movl 0x82c(%rsp), %eax
addl $0x1, %eax
movl %eax, 0x82c(%rsp)
jmp 0xc6f81
movl $0x0, 0x854(%rsp)
movl 0x854(%rsp), %eax
addq $0x858, %rsp            # imm = 0x858
retq
nopl (%rax)