// File: /workspace/github_c_makefile_star2_O1/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: background_output_budget
int background_output_budget(
                             struct background* pba
                             ) {

  double budget_matter, budget_radiation, budget_other,budget_neutrino;
  int index_ncdm;

  budget_matter = 0;
  budget_radiation = 0;
  budget_other = 0;
  budget_neutrino = 0;

  //The name for the class_print_species macro can be at most 30 characters total
  if (pba->background_verbose > 1) {

    printf(" ---------------------------- Budget equation ----------------------- \n");

    printf(" ---> Nonrelativistic Species \n");
    class_print_species("Bayrons",b);
    budget_matter+=pba->Omega0_b;
    if (pba->has_cdm == _TRUE_) {
      class_print_species("Cold Dark Matter",cdm);
      budget_matter+=pba->Omega0_cdm;
    }
    if (pba->has_idm == _TRUE_){
      class_print_species("Interacting DM - idr,b,g",idm);
      budget_matter+=pba->Omega0_idm;
    }
    if (pba->has_dcdm == _TRUE_) {
      class_print_species("Decaying Cold Dark Matter",dcdm);
      budget_matter+=pba->Omega0_dcdm;
    }
   if (pba->has_scf == _TRUE_ && (pba->scf_potential == axion && pba->n_axion == 1) ) {
      class_print_species("Axion DM",scf);
      budget_matter+=pba->Omega0_scf;
    }



    if (pba->N_ncdm > 0) {
      printf(" ---> Non-Cold Dark Matter Species (incl. massive neutrinos)\n");
    }
    if (pba->N_ncdm > 0) {
      for (index_ncdm=0;index_ncdm<pba->N_ncdm;++index_ncdm) {
        printf("-> %-26s%-4d Omega = %-15g , omega = %-15g\n","Non-Cold Species Nr.",index_ncdm+1,pba->Omega0_ncdm[index_ncdm],pba->Omega0_ncdm[index_ncdm]*pba->h*pba->h);
        budget_neutrino+=pba->Omega0_ncdm[index_ncdm];
        budget_matter+=pba->Omega0_ncdm[index_ncdm];
      }
    }

    printf(" ---> Relativistic Species \n");
    class_print_species("Photons",g);
    budget_radiation+=pba->Omega0_g;
    if (pba->has_ur == _TRUE_) {
      class_print_species("Ultra-relativistic relics",ur);
      budget_radiation+=pba->Omega0_ur;
    }
    if (pba->has_dr == _TRUE_) {
      class_print_species("Dark Radiation (from decay)",dr);
      budget_radiation+=pba->Omega0_dr;
    }
    if (pba->has_idr == _TRUE_) {
      class_print_species("Interacting Dark Radiation",idr);
      budget_radiation+=pba->Omega0_idr;
    }

    if ((pba->has_lambda == _TRUE_) || (pba->has_fld == _TRUE_) || (pba->has_scf == _TRUE_) || (pba->has_curvature == _TRUE_)) {
      printf(" ---> Other Content \n");
    }
    if (pba->has_lambda == _TRUE_) {
      class_print_species("Cosmological Constant",lambda);
      budget_other+=pba->Omega0_lambda;
    }
    if (pba->has_fld == _TRUE_) {
      class_print_species("Dark Energy Fluid",fld);
      budget_other+=pba->Omega0_fld;
    }
    if (pba->has_scf == _TRUE_ && !(pba->scf_potential == axion && pba->n_axion == 1) ) {
      class_print_species("Scalar Field",scf);
      budget_other+=pba->Omega0_scf;
      // printf("pba->Omega0_axion %e\n", pba->Omega0_axion);
      if(pba->Omega0_axion!=0 && pba->Omega0_scf == 0.0){
        class_print_species("Scalar Axion",axion);
        budget_other+=pba->Omega0_axion;
      }
    }
    // if(pba->has_scf && (pba->scf_potential == axion || pba->scf_potential == phi_2n)){
    //   _class_print_species_("Axion",axion);
    //   budget_other+=pba->Omega0_axion;
    // }
    if (pba->has_curvature == _TRUE_) {
      class_print_species("Spatial Curvature",k);
      budget_other+=pba->Omega0_k;
    }

    printf(" ---> Total budgets \n");
    printf(" Radiation                        Omega = %-15g , omega = %-15g \n",budget_radiation,budget_radiation*pba->h*pba->h);
    printf(" Non-relativistic                 Omega = %-15g , omega = %-15g \n",budget_matter,budget_matter*pba->h*pba->h);
    if (pba->N_ncdm > 0) {
      printf(" - Non-Free-Streaming Matter      Omega = %-15g , omega = %-15g \n",pba->Omega0_nfsm,pba->Omega0_nfsm*pba->h*pba->h);
      printf(" - Non-Cold Dark Matter           Omega = %-15g , omega = %-15g \n",budget_neutrino,budget_neutrino*pba->h*pba->h);
    }
    if ((pba->has_lambda == _TRUE_) || (pba->has_fld == _TRUE_) || (pba->has_scf == _TRUE_) || (pba->has_curvature == _TRUE_)) {
      printf(" Other Content                    Omega = %-15g , omega = %-15g \n",budget_other,budget_other*pba->h*pba->h);
    }
    printf(" TOTAL                            Omega = %-15g , omega = %-15g \n",budget_radiation+budget_matter+budget_other,(budget_radiation+budget_matter+budget_other)*pba->h*pba->h);
    printf(" -------------------------------------------------------------------- \n");
  }

  return _SUCCESS_;
}

// Assembly:
cmpw $0x2, 0xd1a(%rdi)
jl 0x6474a
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x20, %rsp
movq %rdi, %rbx
leaq 0xac447(%rip), %rdi     # 0x1104f9
callq 0x8510
leaq 0xac482(%rip), %rdi     # 0x110540
callq 0x8510
movsd 0x8(%rbx), %xmm2
movsd 0x20(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xac081(%rip), %rdi     # 0x110161
leaq 0xac0a2(%rip), %rsi     # 0x110189
movb $0x2, %al
callq 0x8030
xorpd %xmm1, %xmm1
addsd 0x20(%rbx), %xmm1
cmpw $0x1, 0x4b4(%rbx)
jne 0x6413b
movsd %xmm1, (%rsp)
movsd 0x8(%rbx), %xmm2
movsd 0x28(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xac03e(%rip), %rdi     # 0x110161
leaq 0xac089(%rip), %rsi     # 0x1101b3
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm1
addsd 0x28(%rbx), %xmm1
cmpw $0x1, 0x4b6(%rbx)
jne 0x64185
movsd %xmm1, (%rsp)
movsd 0x8(%rbx), %xmm2
movsd 0xf0(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabff7(%rip), %rdi     # 0x110161
leaq 0xac020(%rip), %rsi     # 0x110191
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm1
addsd 0xf0(%rbx), %xmm1
cmpw $0x1, 0x4b8(%rbx)
jne 0x641cf
movsd 0x330(%rbx), %xmm0
movsd %xmm1, (%rsp)
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabfad(%rip), %rdi     # 0x110161
leaq 0xabfef(%rip), %rsi     # 0x1101aa
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm1
addsd 0x330(%rbx), %xmm1
cmpw $0x1, 0x4bc(%rbx)
jne 0x64236
cmpl $0x2, 0x160(%rbx)
jne 0x64236
movsd 0x210(%rbx), %xmm0
ucomisd 0x8d06e(%rip), %xmm0    # 0xf1260
jne 0x64236
jp 0x64236
movsd %xmm1, (%rsp)
movsd 0x8(%rbx), %xmm2
movsd 0x128(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabf46(%rip), %rdi     # 0x110161
leaq 0xabfa2(%rip), %rsi     # 0x1101c4
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm1
addsd 0x128(%rbx), %xmm1
cmpl $0x0, 0x248(%rbx)
jle 0x64255
leaq 0xac319(%rip), %rdi     # 0x11055f
movsd %xmm1, (%rsp)
callq 0x8510
movsd (%rsp), %xmm1
cmpl $0x0, 0x248(%rbx)
jle 0x642f2
cmpl $0x0, 0x248(%rbx)
jle 0x642f2
xorpd %xmm2, %xmm2
xorl %r13d, %r13d
leaq 0xabf50(%rip), %r14     # 0x1101cd
leaq 0xabf75(%rip), %r15     # 0x1101f9
movsd %xmm1, (%rsp)
movsd %xmm2, 0x18(%rsp)
leaq 0x1(%r13), %r12
movq 0x280(%rbx), %rax
movsd (%rax,%r13,8), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
movq %r14, %rdi
movq %r15, %rsi
movl %r12d, %edx
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm1
movsd 0x18(%rsp), %xmm2
movq 0x280(%rbx), %rax
movsd (%rax,%r13,8), %xmm0
addsd %xmm0, %xmm2
addsd %xmm0, %xmm1
movslq 0x248(%rbx), %rax
movq %r12, %r13
cmpq %rax, %r12
jl 0x64284
jmp 0x642f6
xorpd %xmm2, %xmm2
movsd %xmm2, 0x18(%rsp)
movsd %xmm1, (%rsp)
leaq 0xac294(%rip), %rdi     # 0x11059c
callq 0x8510
movsd 0x8(%rbx), %xmm2
movsd 0x10(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabe37(%rip), %rdi     # 0x110161
leaq 0xabedd(%rip), %rsi     # 0x11020e
movb $0x2, %al
callq 0x8030
xorpd %xmm0, %xmm0
addsd 0x10(%rbx), %xmm0
cmpw $0x1, 0x4c4(%rbx)
jne 0x6438d
movsd 0x8(%rbx), %xmm2
movsd %xmm0, 0x10(%rsp)
movsd 0xe8(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabdf0(%rip), %rdi     # 0x110161
leaq 0xabe9e(%rip), %rsi     # 0x110216
movb $0x2, %al
callq 0x8030
movsd 0x10(%rsp), %xmm0
addsd 0xe8(%rbx), %xmm0
cmpw $0x1, 0x4ba(%rbx)
jne 0x643d9
movsd %xmm0, 0x10(%rsp)
movsd 0x338(%rbx), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabda4(%rip), %rdi     # 0x110161
leaq 0xabe6c(%rip), %rsi     # 0x110230
movb $0x2, %al
callq 0x8030
movsd 0x10(%rsp), %xmm0
addsd 0x338(%rbx), %xmm0
cmpw $0x1, 0x4c6(%rbx)
jne 0x64425
movsd 0x8(%rbx), %xmm2
movsd %xmm0, 0x10(%rsp)
movsd 0xf8(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabd58(%rip), %rdi     # 0x110161
leaq 0xabe3c(%rip), %rsi     # 0x11024c
movb $0x2, %al
callq 0x8030
movsd 0x10(%rsp), %xmm0
addsd 0xf8(%rbx), %xmm0
movsd %xmm0, 0x10(%rsp)
cmpw $0x1, 0x4c0(%rbx)
je 0x64453
cmpw $0x1, 0x4c2(%rbx)
je 0x64453
cmpw $0x1, 0x4bc(%rbx)
je 0x64453
cmpw $0x1, 0x4c8(%rbx)
jne 0x6445f
leaq 0xac15e(%rip), %rdi     # 0x1105b8
callq 0x8510
xorpd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
cmpw $0x1, 0x4c0(%rbx)
jne 0x644ad
movsd 0x8(%rbx), %xmm2
movsd 0x30(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabcd1(%rip), %rdi     # 0x110161
leaq 0xabdd0(%rip), %rsi     # 0x110267
movb $0x2, %al
callq 0x8030
xorpd %xmm0, %xmm0
addsd 0x30(%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
cmpw $0x1, 0x4c2(%rbx)
jne 0x644f3
movsd 0x8(%rbx), %xmm2
movsd 0x38(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabc8d(%rip), %rdi     # 0x110161
leaq 0xabda2(%rip), %rsi     # 0x11027d
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm0
addsd 0x38(%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
cmpw $0x1, 0x4bc(%rbx)
jne 0x645be
cmpl $0x2, 0x160(%rbx)
jne 0x64522
movsd 0x210(%rbx), %xmm0
ucomisd 0x8cd46(%rip), %xmm0    # 0xf1260
jne 0x64522
jnp 0x645be
movsd 0x8(%rbx), %xmm2
movsd 0x128(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabc1f(%rip), %rdi     # 0x110161
leaq 0xabd46(%rip), %rsi     # 0x11028f
movb $0x2, %al
callq 0x8030
movsd 0x128(%rbx), %xmm1
movsd 0x1e0(%rbx), %xmm0
movsd 0x8(%rsp), %xmm2
addsd %xmm1, %xmm2
movsd %xmm2, 0x8(%rsp)
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm0
jne 0x6457c
jnp 0x645be
ucomisd %xmm2, %xmm1
jne 0x645be
jp 0x645be
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabbc5(%rip), %rdi     # 0x110161
leaq 0xabcf9(%rip), %rsi     # 0x11029c
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm0
addsd 0x1e0(%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
cmpw $0x1, 0x4c8(%rbx)
jne 0x6460a
movsd 0x2d0(%rbx), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabb79(%rip), %rdi     # 0x110161
leaq 0xabcba(%rip), %rsi     # 0x1102a9
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm0
addsd 0x2d0(%rbx), %xmm0
movsd %xmm0, 0x8(%rsp)
leaq 0xabfbc(%rip), %rdi     # 0x1105cd
callq 0x8510
movsd 0x8(%rbx), %xmm2
movsd 0x10(%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabc87(%rip), %rdi     # 0x1102bb
movb $0x2, %al
callq 0x8030
movsd 0x8(%rbx), %xmm2
movsd (%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabca5(%rip), %rdi     # 0x1102fd
movb $0x2, %al
callq 0x8030
cmpl $0x0, 0x248(%rbx)
jle 0x646b4
movsd 0x358(%rbx), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabcb7(%rip), %rdi     # 0x11033f
movb $0x2, %al
callq 0x8030
movsd 0x8(%rbx), %xmm2
movsd 0x18(%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabcd4(%rip), %rdi     # 0x110381
movb $0x2, %al
callq 0x8030
cmpw $0x1, 0x4c0(%rbx)
je 0x646dc
cmpw $0x1, 0x4c2(%rbx)
je 0x646dc
cmpw $0x1, 0x4bc(%rbx)
je 0x646dc
cmpw $0x1, 0x4c8(%rbx)
jne 0x64701
movsd 0x8(%rbx), %xmm2
movsd 0x8(%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabcc9(%rip), %rdi     # 0x1103c3
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
addsd 0x8(%rsp), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xabcdb(%rip), %rdi     # 0x110405
movb $0x2, %al
callq 0x8030
leaq 0xabeaa(%rip), %rdi     # 0x1105e2
callq 0x8510
addq $0x20, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
xorl %eax, %eax
retq