// File: /workspace/github_c_makefile_star2_O2/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 0x5df86
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x20, %rsp
movq %rdi, %rbx
leaq 0xa6aea(%rip), %rdi     # 0x1043d9
callq 0x8540
leaq 0xa6b25(%rip), %rdi     # 0x104420
callq 0x8540
movsd 0x8(%rbx), %xmm2
movsd 0x20(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6724(%rip), %rdi     # 0x104041
leaq 0xa6745(%rip), %rsi     # 0x104069
movb $0x2, %al
callq 0x8030
xorpd %xmm1, %xmm1
addsd 0x20(%rbx), %xmm1
cmpw $0x1, 0x4b4(%rbx)
jne 0x5d97a
movsd %xmm1, 0x8(%rsp)
movsd 0x8(%rbx), %xmm2
movsd 0x28(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa66e0(%rip), %rdi     # 0x104041
leaq 0xa672b(%rip), %rsi     # 0x104093
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm1
addsd 0x28(%rbx), %xmm1
cmpw $0x1, 0x4b6(%rbx)
jne 0x5d9c6
movsd %xmm1, 0x8(%rsp)
movsd 0x8(%rbx), %xmm2
movsd 0xf0(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6697(%rip), %rdi     # 0x104041
leaq 0xa66c0(%rip), %rsi     # 0x104071
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm1
addsd 0xf0(%rbx), %xmm1
cmpw $0x1, 0x4b8(%rbx)
jne 0x5da12
movsd 0x330(%rbx), %xmm0
movsd %xmm1, 0x8(%rsp)
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa664b(%rip), %rdi     # 0x104041
leaq 0xa668d(%rip), %rsi     # 0x10408a
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm1
addsd 0x330(%rbx), %xmm1
cmpw $0x1, 0x4bc(%rbx)
jne 0x5da7b
cmpl $0x2, 0x160(%rbx)
jne 0x5da7b
movsd 0x210(%rbx), %xmm0
ucomisd 0x8782b(%rip), %xmm0    # 0xe5260
jne 0x5da7b
jp 0x5da7b
movsd %xmm1, 0x8(%rsp)
movsd 0x8(%rbx), %xmm2
movsd 0x128(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa65e2(%rip), %rdi     # 0x104041
leaq 0xa663e(%rip), %rsi     # 0x1040a4
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm1
addsd 0x128(%rbx), %xmm1
movsd %xmm1, 0x8(%rsp)
cmpl $0x0, 0x248(%rbx)
jle 0x5db34
leaq 0xa69aa(%rip), %rdi     # 0x10443f
callq 0x8540
movl 0x248(%rbx), %eax
testl %eax, %eax
jle 0x5db34
xorpd %xmm0, %xmm0
movsd %xmm0, 0x18(%rsp)
xorl %r13d, %r13d
leaq 0xa65f1(%rip), %r14     # 0x1040ad
leaq 0xa6616(%rip), %r15     # 0x1040d9
cltq
cmpq %rax, %r13
jge 0x5db3e
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
movq 0x280(%rbx), %rax
movsd (%rax,%r13,8), %xmm0
movsd 0x18(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x18(%rsp)
movsd 0x8(%rsp), %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
movl 0x248(%rbx), %eax
movq %r12, %r13
jmp 0x5dac3
xorpd %xmm0, %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0xa6937(%rip), %rdi     # 0x10447c
callq 0x8540
movsd 0x8(%rbx), %xmm2
movsd 0x10(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa64da(%rip), %rdi     # 0x104041
leaq 0xa6580(%rip), %rsi     # 0x1040ee
movb $0x2, %al
callq 0x8030
xorpd %xmm0, %xmm0
addsd 0x10(%rbx), %xmm0
cmpw $0x1, 0x4c4(%rbx)
jne 0x5dbca
movsd 0x8(%rbx), %xmm2
movsd %xmm0, 0x10(%rsp)
movsd 0xe8(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6493(%rip), %rdi     # 0x104041
leaq 0xa6541(%rip), %rsi     # 0x1040f6
movb $0x2, %al
callq 0x8030
movsd 0x10(%rsp), %xmm0
addsd 0xe8(%rbx), %xmm0
cmpw $0x1, 0x4ba(%rbx)
jne 0x5dc16
movsd %xmm0, 0x10(%rsp)
movsd 0x338(%rbx), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6447(%rip), %rdi     # 0x104041
leaq 0xa650f(%rip), %rsi     # 0x104110
movb $0x2, %al
callq 0x8030
movsd 0x10(%rsp), %xmm0
addsd 0x338(%rbx), %xmm0
cmpw $0x1, 0x4c6(%rbx)
movsd %xmm0, 0x10(%rsp)
jne 0x5dc68
movsd 0x8(%rbx), %xmm2
movsd 0xf8(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa63fb(%rip), %rdi     # 0x104041
leaq 0xa64df(%rip), %rsi     # 0x10412c
movb $0x2, %al
callq 0x8030
movsd 0x10(%rsp), %xmm0
addsd 0xf8(%rbx), %xmm0
movsd %xmm0, 0x10(%rsp)
cmpw $0x1, 0x4c0(%rbx)
je 0x5dc99
cmpw $0x1, 0x4c2(%rbx)
je 0x5dc99
cmpw $0x1, 0x4bc(%rbx)
je 0x5dc99
xorpd %xmm0, %xmm0
movsd %xmm0, (%rsp)
cmpw $0x1, 0x4c8(%rbx)
jne 0x5dcf1
leaq 0xa67f8(%rip), %rdi     # 0x104498
callq 0x8540
xorpd %xmm0, %xmm0
movsd %xmm0, (%rsp)
cmpw $0x1, 0x4c0(%rbx)
jne 0x5dcf1
movsd 0x8(%rbx), %xmm2
movsd 0x30(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa636c(%rip), %rdi     # 0x104041
leaq 0xa646b(%rip), %rsi     # 0x104147
movb $0x2, %al
callq 0x8030
xorpd %xmm0, %xmm0
addsd 0x30(%rbx), %xmm0
movsd %xmm0, (%rsp)
cmpw $0x1, 0x4c2(%rbx)
jne 0x5dd35
movsd 0x8(%rbx), %xmm2
movsd 0x38(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6329(%rip), %rdi     # 0x104041
leaq 0xa643e(%rip), %rsi     # 0x10415d
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm0
addsd 0x38(%rbx), %xmm0
movsd %xmm0, (%rsp)
cmpw $0x1, 0x4bc(%rbx)
jne 0x5ddfc
cmpl $0x2, 0x160(%rbx)
jne 0x5dd64
movsd 0x210(%rbx), %xmm0
ucomisd 0x87504(%rip), %xmm0    # 0xe5260
jne 0x5dd64
jnp 0x5ddfc
movsd 0x8(%rbx), %xmm2
movsd 0x128(%rbx), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa62bd(%rip), %rdi     # 0x104041
leaq 0xa63e4(%rip), %rsi     # 0x10416f
movb $0x2, %al
callq 0x8030
movsd 0x128(%rbx), %xmm1
movsd 0x1e0(%rbx), %xmm0
movsd (%rsp), %xmm2
addsd %xmm1, %xmm2
movsd %xmm2, (%rsp)
xorpd %xmm2, %xmm2
ucomisd %xmm2, %xmm0
jne 0x5ddbc
jnp 0x5ddfc
ucomisd %xmm2, %xmm1
jne 0x5ddfc
jp 0x5ddfc
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6265(%rip), %rdi     # 0x104041
leaq 0xa6399(%rip), %rsi     # 0x10417c
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm0
addsd 0x1e0(%rbx), %xmm0
movsd %xmm0, (%rsp)
cmpw $0x1, 0x4c8(%rbx)
jne 0x5de46
movsd 0x2d0(%rbx), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa621b(%rip), %rdi     # 0x104041
leaq 0xa635c(%rip), %rsi     # 0x104189
movb $0x2, %al
callq 0x8030
movsd (%rsp), %xmm0
addsd 0x2d0(%rbx), %xmm0
movsd %xmm0, (%rsp)
leaq 0xa6660(%rip), %rdi     # 0x1044ad
callq 0x8540
movsd 0x8(%rbx), %xmm2
movsd 0x10(%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa632b(%rip), %rdi     # 0x10419b
movb $0x2, %al
callq 0x8030
movsd 0x8(%rbx), %xmm2
movsd 0x8(%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6348(%rip), %rdi     # 0x1041dd
movb $0x2, %al
callq 0x8030
cmpl $0x0, 0x248(%rbx)
jle 0x5def1
movsd 0x358(%rbx), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa635a(%rip), %rdi     # 0x10421f
movb $0x2, %al
callq 0x8030
movsd 0x8(%rbx), %xmm2
movsd 0x18(%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa6377(%rip), %rdi     # 0x104261
movb $0x2, %al
callq 0x8030
cmpw $0x1, 0x4c0(%rbx)
je 0x5df19
cmpw $0x1, 0x4c2(%rbx)
je 0x5df19
cmpw $0x1, 0x4bc(%rbx)
je 0x5df19
cmpw $0x1, 0x4c8(%rbx)
jne 0x5df3d
movsd 0x8(%rbx), %xmm2
movsd (%rsp), %xmm0
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa636d(%rip), %rdi     # 0x1042a3
movb $0x2, %al
callq 0x8030
movsd 0x8(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
addsd (%rsp), %xmm0
movsd 0x8(%rbx), %xmm2
movapd %xmm0, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm2, %xmm1
leaq 0xa637f(%rip), %rdi     # 0x1042e5
movb $0x2, %al
callq 0x8030
leaq 0xa654e(%rip), %rdi     # 0x1044c2
callq 0x8540
addq $0x20, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
xorl %eax, %eax
retq