LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: ClosedModY
int ClosedModY(int l, int beta, double *y, int * phisign, int * dphisign){

  *phisign = 1;
  *dphisign = 1;


  while (*y > _TWOPI_)
    *y -= _TWOPI_;

  if ((*y) > _PI_){
    *y = 2.0*_PI_-(*y);
    //phisign *= pow(-1,l)
    if (l%2==1) //l is odd
      *phisign = -*phisign;
    else        //l is even
      *dphisign = -*dphisign;
  }
  if ((*y)>0.5*_PI_){
    *y = _PI_-(*y);
    //phisign *= pow(-1,beta-l-1)
    if ((beta-l)%2==0) //beta-l-1 odd
      *phisign = -*phisign;
    else               //beta-l-1 even
      *dphisign = -*dphisign;
  }
  return _SUCCESS_;
}

// Assembly:
pushq $0x1
popq %rax
movl %eax, (%rcx)
movl %eax, (%r8)
movsd (%rdx), %xmm1
movsd 0xc6711(%rip), %xmm0    # 0xe83e8
movsd 0xc6771(%rip), %xmm2    # 0xe8450
ucomisd %xmm0, %xmm1
jbe 0x21cef
addsd %xmm2, %xmm1
movsd %xmm1, (%rdx)
jmp 0x21cdf
ucomisd 0xc62a1(%rip), %xmm1    # 0xe7f98
jbe 0x21d11
subsd %xmm1, %xmm0
movsd %xmm0, (%rdx)
movl %edi, %eax
andl $0x80000001, %eax       # imm = 0x80000001
cmpl $0x1, %eax
jne 0x21d17
negl (%rcx)
jmp 0x21d1b
movapd %xmm1, %xmm0
jmp 0x21d1b
orl $-0x1, (%r8)
ucomisd 0xc6725(%rip), %xmm0    # 0xe8448
jbe 0x21d42
movsd 0xc626b(%rip), %xmm1    # 0xe7f98
subsd %xmm0, %xmm1
subl %edi, %esi
testb $0x1, %sil
cmoveq %rcx, %r8
movsd %xmm1, (%rdx)
negl (%r8)
xorl %eax, %eax
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: hyperspherical_forwards_recurrence
int hyperspherical_forwards_recurrence(int K,
                                       int lmax,
                                       double beta,
                                       double x,
                                       double sinK,
                                       double cotK,
                                       double * __restrict__ sqrtK,
                                       double * __restrict__ one_over_sqrtK,
                                       double * __restrict__ PhiL){
  int l;
  PhiL[0] = 1.0/beta*sin(beta*x)/sinK;
  PhiL[1] = PhiL[0]*(cotK-beta/tan(beta*x))*one_over_sqrtK[1];
  for (l=2; l<=lmax; l++){
    PhiL[l] = ((2*l-1)*cotK*PhiL[l-1]-PhiL[l-2]*sqrtK[l-1])*one_over_sqrtK[l];
  }
  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x28, %rsp
movq %r8, %rbx
movq %rcx, %r14
movq %rdx, %r15
movsd %xmm3, 0x20(%rsp)
movsd %xmm2, (%rsp)
movsd %xmm0, 0x18(%rsp)
movl %esi, %ebp
movsd 0xc34ed(%rip), %xmm2    # 0xe5260
divsd %xmm0, %xmm2
movsd %xmm2, 0x8(%rsp)
mulsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rsp)
movapd %xmm1, %xmm0
callq 0x80a0
mulsd 0x8(%rsp), %xmm0
divsd (%rsp), %xmm0
movsd %xmm0, (%rsp)
movsd %xmm0, (%rbx)
movsd 0x10(%rsp), %xmm0
callq 0x8230
movsd 0x20(%rsp), %xmm2
movsd 0x18(%rsp), %xmm1
divsd %xmm0, %xmm1
movapd %xmm2, %xmm0
subsd %xmm1, %xmm0
mulsd (%rsp), %xmm0
mulsd 0x8(%r14), %xmm0
movsd %xmm0, 0x8(%rbx)
movslq %ebp, %rax
pushq $0x2
popq %rcx
pushq $0x3
popq %rdx
cmpq %rax, %rcx
jg 0x21e17
cvtsi2sd %edx, %xmm0
mulsd %xmm2, %xmm0
movsd -0x8(%r15,%rcx,8), %xmm1
mulsd -0x10(%rbx,%rcx,8), %xmm1
mulsd -0x8(%rbx,%rcx,8), %xmm0
subsd %xmm1, %xmm0
mulsd (%r14,%rcx,8), %xmm0
movsd %xmm0, (%rbx,%rcx,8)
incq %rcx
addl $0x2, %edx
jmp 0x21de0
xorl %eax, %eax
addq $0x28, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: hyperspherical_forwards_recurrence_chunk
int hyperspherical_forwards_recurrence_chunk(int K,
                                             int lmax,
                                             double beta,
                                             double * __restrict__ x,
                                             double * __restrict__ sinK,
                                             double * __restrict__ cotK,
                                             int chunk,
                                             double * __restrict__ sqrtK,
                                             double * __restrict__ one_over_sqrtK,
                                             double * __restrict__ PhiL){
  int l;
  int index_x;
  for (index_x=0; index_x<chunk; index_x++){
    PhiL[index_x] = 1.0/beta*sin(beta*x[index_x])/sinK[index_x];
    PhiL[chunk+index_x] = PhiL[index_x]*
      (cotK[index_x]-beta/tan(beta*x[index_x]))*one_over_sqrtK[1];
  }
  for (l=2; l<=lmax; l++){
    for (index_x=0; index_x<chunk; index_x++)
      PhiL[l*chunk+index_x] =
        ((2*l-1)*cotK[index_x]*PhiL[(l-1)*chunk+index_x]-
         PhiL[(l-2)*chunk+index_x]*sqrtK[l-1])*one_over_sqrtK[l];
  }
  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x58, %rsp
movq %r8, %rbx
movq %rcx, 0x38(%rsp)
movq %rdx, 0x30(%rsp)
movl %esi, 0xc(%rsp)
movq 0xa0(%rsp), %rbp
movq 0x98(%rsp), %r12
movq 0x90(%rsp), %r13
movsd 0xc33fd(%rip), %xmm1    # 0xe5260
movsd %xmm0, 0x10(%rsp)
divsd %xmm0, %xmm1
movsd %xmm1, 0x28(%rsp)
movslq %r9d, %rax
xorl %r14d, %r14d
testl %eax, %eax
movl $0x0, %r15d
cmovgl %eax, %r15d
movq %rax, 0x18(%rsp)
leaq (,%rax,8), %rax
addq %rbp, %rax
movq %rax, 0x20(%rsp)
cmpq %r14, %r15
je 0x21f28
movq 0x30(%rsp), %rax
movsd (%rax,%r14,8), %xmm0
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x48(%rsp)
callq 0x80a0
mulsd 0x28(%rsp), %xmm0
movq 0x38(%rsp), %rax
divsd (%rax,%r14,8), %xmm0
movsd %xmm0, 0x50(%rsp)
movsd %xmm0, (%rbp,%r14,8)
movsd (%rbx,%r14,8), %xmm0
movsd %xmm0, 0x40(%rsp)
movsd 0x48(%rsp), %xmm0
callq 0x8230
movsd 0x10(%rsp), %xmm1
divsd %xmm0, %xmm1
movsd 0x40(%rsp), %xmm0
subsd %xmm1, %xmm0
mulsd 0x50(%rsp), %xmm0
mulsd 0x8(%r12), %xmm0
movq 0x20(%rsp), %rax
movsd %xmm0, (%rax,%r14,8)
incq %r14
jmp 0x21e9a
movslq 0xc(%rsp), %rax
movq 0x18(%rsp), %r8
movq %r8, %rcx
shlq $0x4, %rcx
addq %rbp, %rcx
shlq $0x3, %r8
pushq $0x2
popq %rdx
cmpq %rax, %rdx
jg 0x21f98
leal -0x1(,%rdx,2), %esi
cvtsi2sd %esi, %xmm0
leaq (%r8,%rbp), %rsi
xorl %edi, %edi
cmpq %rdi, %r15
je 0x21f8d
movsd (%rbx,%rdi,8), %xmm1
movsd -0x8(%r13,%rdx,8), %xmm2
mulsd (%rbp,%rdi,8), %xmm2
mulsd %xmm0, %xmm1
mulsd (%rsi,%rdi,8), %xmm1
subsd %xmm2, %xmm1
mulsd (%r12,%rdx,8), %xmm1
movsd %xmm1, (%rcx,%rdi,8)
incq %rdi
jmp 0x21f59
incq %rdx
addq %r8, %rcx
movq %rsi, %rbp
jmp 0x21f43
xorl %eax, %eax
addq $0x58, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: coef1
double coef1(double z){
  const double A[5] = {1.1282427601,-0.6803534e-4,0.16687e-6,-0.128e-8,0.2e-10};
  const double B[5] = {0.7822108673e-1,-0.6895649e-4,0.32857e-6,-0.37e-8,0.7e-10};
  double x,y,t,Ai,zeta,theta,sintheta,costheta,FA,FB;

  x = -z;
  zeta = _TWO_OVER_THREE_*x*sqrt(x);
  theta = zeta+0.25*_PI_;
  sintheta = sin(theta);
  costheta = cos(theta);

  y = pow(7.0/x,3);//y = (7.0/x)*(7.0/x)*(7.0/x);
  FA = cheb(y,5,A);
  FB = cheb(y,5,B)/zeta;

  t = pow(x,-0.25);
  Ai = t*(sintheta*FA-costheta*FB);
  //Bi = t*(costheta*FA+sintheta*FB);
  return Ai;
}

// Assembly:
pushq %rbx
subq $0xa0, %rsp
movapd %xmm0, %xmm2
movaps 0xc5364(%rip), %xmm0    # 0xe8810
movaps %xmm0, 0x70(%rsp)
movaps 0xc5368(%rip), %xmm0    # 0xe8820
movaps %xmm0, 0x80(%rsp)
movabsq $0x3db5fd7fe1796495, %rax # imm = 0x3DB5FD7FE1796495
movq %rax, 0x90(%rsp)
movaps 0xc5357(%rip), %xmm0    # 0xe8830
movaps %xmm0, 0x40(%rsp)
movaps 0xc535b(%rip), %xmm0    # 0xe8840
movaps %xmm0, 0x50(%rsp)
movabsq $0x3dd33dcfe54a3803, %rax # imm = 0x3DD33DCFE54A3803
movq %rax, 0x60(%rsp)
movapd 0xc1d2f(%rip), %xmm0    # 0xe5230
xorpd %xmm2, %xmm0
movsd 0xc4fb3(%rip), %xmm1    # 0xe84c0
mulsd %xmm2, %xmm1
movsd %xmm1, (%rsp)
movsd 0xc4faa(%rip), %xmm1    # 0xe84c8
movapd %xmm2, 0x10(%rsp)
ucomisd %xmm2, %xmm1
movapd %xmm0, 0x30(%rsp)
jb 0x23536
sqrtsd %xmm0, %xmm0
jmp 0x2353b
callq 0x82d0
movsd (%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
movsd 0xc4f7f(%rip), %xmm0    # 0xe84d0
addsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0x80a0
movsd %xmm0, 0x28(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x86d0
movsd %xmm0, 0x8(%rsp)
movsd 0xc4f21(%rip), %xmm0    # 0xe84a0
divsd 0x10(%rsp), %xmm0
movsd 0xc31bb(%rip), %xmm1    # 0xe6748
callq 0x8650
movsd %xmm0, 0x20(%rsp)
pushq $0x5
popq %rbx
leaq 0x70(%rsp), %rsi
movl %ebx, %edi
callq 0x23788
movsd %xmm0, 0x10(%rsp)
leaq 0x40(%rsp), %rsi
movsd 0x20(%rsp), %xmm0
movl %ebx, %edi
callq 0x23788
divsd (%rsp), %xmm0
movsd %xmm0, (%rsp)
movsd 0xc1cb7(%rip), %xmm1    # 0xe5288
movapd 0x30(%rsp), %xmm0
callq 0x8650
movsd (%rsp), %xmm2
mulsd 0x8(%rsp), %xmm2
movsd 0x10(%rsp), %xmm1
mulsd 0x28(%rsp), %xmm1
subsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
addq $0xa0, %rsp
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: coef2
double coef2(double z){
  const double A[17] = {0.11535880704,0.6542816649e-1,0.26091774326,0.21959346500,
              0.12476382168,-0.43476292594,0.28357718605,-0.9751797082e-1,
              0.2182551823e-1,-0.350454097e-2,0.42778312e-3,
              -0.4127564e-4,0.323880e-5,-0.21123e-6,0.1165e-7,
              -0.55e-9,0.2e-10};
  const double B[16] = {0.10888288487,-0.17511655051,0.13887871101,-0.11469998998,
             0.22377807641,-0.18546243714,0.8063565186e-1,
             -0.2208847864e-1,0.422444527e-2,-0.60131028e-3,
             0.6653890e-4,-0.590842e-5,0.43131e-6,-0.2638e-7,
             0.137e-8,-0.6e-10};
  //Ej = 3^(-j/3)/Gamma(j/3);
  double E1 = 0.355028053887817, E2 = 0.258819403792807;
  double x,FA,FB,Ai;
  x = -(z/7.0)*(z/7.0)*(z/7.0);
  FA = E1*cheb(x,17,A);
  FB = E2*z*cheb(x,16,B);
  Ai = FA-FB;
  //Bi = sqrt(3)*(FA+FB);
  return Ai;
}

// Assembly:
subq $0x28, %rsp
movsd %xmm0, (%rsp)
movapd %xmm0, %xmm1
divsd 0xc3127(%rip), %xmm1    # 0xe6740
movapd 0xc1c0f(%rip), %xmm0    # 0xe5230
xorpd %xmm1, %xmm0
mulsd %xmm1, %xmm0
mulsd %xmm1, %xmm0
movapd %xmm0, 0x10(%rsp)
leaq 0xc52d6(%rip), %rsi     # 0xe8910
pushq $0x11
popq %rdi
callq 0x23788
mulsd 0xc4e8e(%rip), %xmm0    # 0xe84d8
movsd %xmm0, 0x8(%rsp)
movsd (%rsp), %xmm0
mulsd 0xc4e83(%rip), %xmm0    # 0xe84e0
movsd %xmm0, (%rsp)
leaq 0xc5337(%rip), %rsi     # 0xe89a0
pushq $0x10
popq %rdi
movapd 0x10(%rsp), %xmm0
callq 0x23788
mulsd (%rsp), %xmm0
addsd 0x8(%rsp), %xmm0
addq $0x28, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: coef4
double coef4(double z){
  const double A[7]={0.56265126169,-0.76136219e-3,0.765252e-5,-0.14228e-6,
            0.380e-8,-0.13e-9,0.1e-10};
  /**  double B[7]={1.1316635302,0.166141673e-02,0.1968882e-04,0.47047e-06,
            0.1769e-7,0.94e-9,0.6e-10};
  */
  double x,Y,t,zeta,EX,EY,Ai;

  Y = z*sqrt(z);
  zeta = 2.0/3.0*Y;
  EX = exp(zeta);
  EY = 1.0/EX;
  x = 7*sqrt(7)/Y;
  t = pow(z,-0.25);

  Ai = t*EY*cheb(x, 7, A);
  //Bi = t*EX*cheb(x, 7, &(B[0]));
  return Ai;
}

// Assembly:
subq $0x18, %rsp
movapd %xmm0, %xmm1
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
movsd %xmm1, (%rsp)
jb 0x236f9
sqrtsd %xmm1, %xmm0
jmp 0x23707
movapd %xmm1, %xmm0
callq 0x82d0
movsd (%rsp), %xmm1
mulsd %xmm1, %xmm0
movapd %xmm0, %xmm1
movsd %xmm0, 0x8(%rsp)
movsd 0xc4d93(%rip), %xmm0    # 0xe84b0
mulsd %xmm1, %xmm0
callq 0x86a0
movsd 0xc1b32(%rip), %xmm1    # 0xe5260
divsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rsp)
movsd 0xc4db0(%rip), %xmm0    # 0xe84f0
divsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0xc1b34(%rip), %xmm1    # 0xe5288
movsd (%rsp), %xmm0
callq 0x8650
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, (%rsp)
leaq 0xc5350(%rip), %rsi     # 0xe8ac0
pushq $0x7
popq %rdi
movsd 0x8(%rsp), %xmm0
callq 0x23788
mulsd (%rsp), %xmm0
addq $0x18, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: cheb
double cheb(double x, int n, const double A[]){
  double b,d,u,y,c,F;
  int j;
  b = 0.0;
  d = A[n-1];
  u = 2*x-1.0;
  y = 2*u;

  for (j=(n-1);j>1;j--){
    c=b;
    b=d;
    d=y*b-c+A[j-1];
  }
  F = u*d-b+0.5*A[0];
  return F;
}

// Assembly:
movslq %edi, %rax
decl %edi
movsd -0x8(%rsi,%rax,8), %xmm4
addsd %xmm0, %xmm0
addsd 0xc2049(%rip), %xmm0    # 0xe57e8
movapd %xmm0, %xmm1
addsd %xmm0, %xmm1
xorpd %xmm2, %xmm2
movapd 0xc1a7d(%rip), %xmm3    # 0xe5230
movapd %xmm2, %xmm5
movapd %xmm4, %xmm2
xorpd %xmm3, %xmm5
cmpl $0x2, %edi
jl 0x237db
movapd %xmm1, %xmm4
mulsd %xmm2, %xmm4
addsd %xmm5, %xmm4
addsd -0x8(%rsi,%rdi,8), %xmm4
decq %rdi
jmp 0x237b3
mulsd %xmm2, %xmm0
addsd %xmm5, %xmm0
movsd (%rsi), %xmm1
mulsd 0xc24c1(%rip), %xmm1    # 0xe5cb0
addsd %xmm0, %xmm1
movapd %xmm1, %xmm0
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/hyperspherical.c
// Function: get_value_at_small_phi
double get_value_at_small_phi(int K,int l,double beta,double Phi){
  double nu, lhs, alpha, xval;

  nu = l+0.5;
  lhs = 1.0/nu*log(2*Phi*nu);
  alpha = -2*lhs/5.0*(1.0+2.0*cosh(1.0/3.0*acosh(1.0+375/(16.0*lhs*lhs))));
  xval = nu/cosh(alpha)/beta;
  //Correct for geometry:
  if (K==1)
    xval *= asin(l/beta)/(l/beta);
  else if(K==-1)
    xval *= asinh(l/beta)/(l/beta);
  return xval;
}

// Assembly:
pushq %rbx
subq $0x20, %rsp
movsd %xmm0, 0x18(%rsp)
cvtsi2sd %esi, %xmm0
movl %edi, %ebx
movsd 0xc249f(%rip), %xmm2    # 0xe5cb0
movsd %xmm0, 0x8(%rsp)
addsd %xmm0, %xmm2
movsd %xmm2, (%rsp)
movsd 0xc1a38(%rip), %xmm0    # 0xe5260
divsd %xmm2, %xmm0
movsd %xmm0, 0x10(%rsp)
addsd %xmm1, %xmm1
mulsd %xmm2, %xmm1
movapd %xmm1, %xmm0
callq 0x82e0
mulsd 0x10(%rsp), %xmm0
movsd 0xc2ee7(%rip), %xmm1    # 0xe6738
mulsd %xmm0, %xmm1
divsd 0xc1a13(%rip), %xmm1    # 0xe5270
movsd %xmm1, 0x10(%rsp)
movsd 0xc1f8d(%rip), %xmm1    # 0xe57f8
mulsd %xmm0, %xmm1
mulsd %xmm0, %xmm1
movsd 0xc4b85(%rip), %xmm0    # 0xe8400
divsd %xmm1, %xmm0
addsd 0xc19d9(%rip), %xmm0    # 0xe5260
callq 0x8660
mulsd 0xc4b74(%rip), %xmm0    # 0xe8408
callq 0x84d0
addsd %xmm0, %xmm0
addsd 0xc19bb(%rip), %xmm0    # 0xe5260
mulsd 0x10(%rsp), %xmm0
callq 0x84d0
movsd 0x18(%rsp), %xmm1
movsd (%rsp), %xmm2
divsd %xmm0, %xmm2
movapd %xmm2, %xmm0
divsd %xmm1, %xmm0
cmpl $-0x1, %ebx
je 0x238f1
cmpl $0x1, %ebx
jne 0x23928
movsd 0x8(%rsp), %xmm2
divsd %xmm1, %xmm2
movsd %xmm2, 0x8(%rsp)
movsd %xmm0, (%rsp)
movapd %xmm2, %xmm0
callq 0x85d0
jmp 0x2390f
movsd 0x8(%rsp), %xmm2
divsd %xmm1, %xmm2
movsd %xmm2, 0x8(%rsp)
movsd %xmm0, (%rsp)
movapd %xmm2, %xmm0
callq 0x84a0
divsd 0x8(%rsp), %xmm0
movsd (%rsp), %xmm1
mulsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
movsd (%rsp), %xmm0
addq $0x20, %rsp
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/common.c
// Function: class_protect_fprintf
void class_protect_fprintf(FILE* stream, char* tpl,...) {
  va_list args;
  char dest[6000];
  va_start(args,tpl);
  vsnprintf(dest, 2048,tpl,args);
  va_end(args);
  fprintf(stream,"%s",dest);
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x1848, %rsp           # imm = 0x1848
movq %rsi, %r10
movq %rdi, %rbx
leaq 0x20(%rsp), %rsi
movq %rdx, 0x10(%rsi)
movq %rcx, 0x18(%rsi)
movq %r8, 0x20(%rsi)
movq %r9, 0x28(%rsi)
testb %al, %al
je 0x2e103
movaps %xmm0, 0x50(%rsp)
movaps %xmm1, 0x60(%rsp)
movaps %xmm2, 0x70(%rsp)
movaps %xmm3, 0x80(%rsp)
movaps %xmm4, 0x90(%rsp)
movaps %xmm5, 0xa0(%rsp)
movaps %xmm6, 0xb0(%rsp)
movaps %xmm7, 0xc0(%rsp)
movabsq $0x3000000010, %rax     # imm = 0x3000000010
movq %rsp, %rcx
movq %rax, (%rcx)
leaq 0x1860(%rsp), %rax
movq %rax, 0x8(%rcx)
movq %rsi, 0x10(%rcx)
leaq 0xd0(%rsp), %r14
movl $0x800, %esi            # imm = 0x800
movq %r14, %rdi
movq %r10, %rdx
callq 0x8260
movq %r14, %rdi
movq %rbx, %rsi
callq 0x8410
addq $0x1848, %rsp           # imm = 0x1848
popq %rbx
popq %r14
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/common.c
// Function: get_number_of_titles
int get_number_of_titles(char * titlestring){
  int i;
  int number_of_titles=0;

  for (i=0; i<strlen(titlestring); i++){
    if (titlestring[i] == '\t')
      number_of_titles++;
  }
  return number_of_titles;
}

// Assembly:
pushq %rbx
movq %rdi, %rbx
callq 0x8190
xorl %edx, %edx
xorl %ecx, %ecx
cmpq %rdx, %rax
je 0x2e179
xorl %esi, %esi
cmpb $0x9, (%rbx,%rdx)
sete %sil
addl %esi, %ecx
incq %rdx
jmp 0x2e163
movl %ecx, %eax
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/common.c
// Function: file_exists
int file_exists(const char *fname){
  FILE *file = fopen(fname, "r");
  if (file != NULL){
    fclose(file);
    return _TRUE_;
  }

  return _FALSE_;

}

// Assembly:
pushq %rax
leaq 0xd5adf(%rip), %rsi     # 0x103c64
callq 0x8420
testq %rax, %rax
je 0x2e19c
movq %rax, %rdi
callq 0x83c0
pushq $0x1
popq %rax
jmp 0x2e19e
xorl %eax, %eax
popq %rcx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../tools/common.c
// Function: string_begins_with
int string_begins_with(char* thestring, char beginchar){

  /** Define temporary variables */
  int int_temp=0;
  int strlength = strlen((thestring));
  int result = _FALSE_;

  /** Check through the beginning of the string to see if the beginchar is met */
  for(int_temp=0;int_temp<strlength;++int_temp){
    /* Skip over whitespaces (very important) */
    if(thestring[int_temp]==' ' || thestring[int_temp]=='\t'){continue;}
    /* If the beginchar is met, everything is good */
    else if(thestring[int_temp]==beginchar){result=_TRUE_;}
    /* If something else is met, cancel */
    else{break;}
  }

  return result;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movl %esi, %ebx
movq %rdi, %r14
callq 0x8190
xorl %edx, %edx
testl %eax, %eax
cmovlel %edx, %eax
pushq $0x1
popq %rsi
xorl %ecx, %ecx
cmpq %rdx, %rax
je 0x2e1fe
movzbl (%r14,%rdx), %r8d
cmpl $0x9, %r8d
je 0x2e1ec
cmpl $0x20, %r8d
jne 0x2e1f0
movl %ecx, %edi
jmp 0x2e1f7
movl %esi, %edi
cmpb %bl, %r8b
jne 0x2e1fe
incq %rdx
movl %edi, %ecx
jmp 0x2e1d6
movl %ecx, %eax
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/input.c
// Function: input_init
int input_init(int argc,
               char **argv,
               struct precision * ppr,
               struct background *pba,
               struct thermodynamics *pth,
               struct perturbations *ppt,
               struct transfer *ptr,
               struct primordial *ppm,
               struct harmonic *phr,
               struct fourier * pfo,
               struct lensing *ple,
               struct distortions *psd,
               struct output *pop,
               ErrorMsg errmsg){

  /** Summary: */

  /** Define local variables */
  struct file_content fc;        // Structure with all parameters

  /** Find and read input file */
  class_call(input_find_file(argc,
                             argv,
                             &fc,
                             errmsg),
             errmsg,
             errmsg);

  /** Initialize all parameters given the input 'file_content' structure.
      If its size is null, all parameters take their default values. */
  class_call(input_read_from_file(&fc,ppr,pba,pth,ppt,ptr,ppm,phr,pfo,ple,psd,pop,
                                  errmsg),
             errmsg,
             errmsg);

  /** Free local struture */
  class_call(parser_free(&fc),
             errmsg,
             errmsg);

  return _SUCCESS_;

}

// Assembly:
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x830, %rsp            # imm = 0x830
movq %r9, %r14
movq %r8, %r15
movq %rcx, %r12
movq %rdx, %r13
movq 0x898(%rsp), %rbx
leaq 0x8(%rsp), %rdx
movq %rbx, %rcx
callq 0x2ea12
testl %eax, %eax
je 0x2e926
leaq 0xb6af8(%rip), %rsi     # 0xe53f0
leaq 0xbb389(%rip), %rdx     # 0xe9c88
leaq 0x30(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x2e014
leaq 0xb671c(%rip), %rsi     # 0xe5034
leaq 0xbb392(%rip), %rdx     # 0xe9cb1
pushq $0x4b
jmp 0x2e9f0
movq 0x860(%rsp), %r9
subq $0x8, %rsp
leaq 0x10(%rsp), %rdi
movq %r13, %rsi
movq %r12, %rdx
movq %r15, %rcx
movq %r14, %r8
pushq %rbx
pushq 0x8a0(%rsp)
pushq 0x8a0(%rsp)
pushq 0x8a0(%rsp)
pushq 0x8a0(%rsp)
pushq 0x8a0(%rsp)
pushq 0x8a0(%rsp)
callq 0x2ee2c
addq $0x40, %rsp
testl %eax, %eax
je 0x2e9ad
leaq 0xb6a6e(%rip), %rsi     # 0xe53f0
leaq 0xbb333(%rip), %rdx     # 0xe9cbc
leaq 0x30(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x2e014
leaq 0xb6692(%rip), %rsi     # 0xe5034
leaq 0xbb308(%rip), %rdx     # 0xe9cb1
pushq $0x52
jmp 0x2e9f0
leaq 0x8(%rsp), %rdi
callq 0x1c593
movl %eax, %ecx
xorl %eax, %eax
cmpl $0x1, %ecx
jne 0x2ea01
leaq 0xb6a29(%rip), %rsi     # 0xe53f0
leaq 0xbb33c(%rip), %rdx     # 0xe9d0a
leaq 0x30(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x2e014
leaq 0xb664d(%rip), %rsi     # 0xe5034
leaq 0xbb2c3(%rip), %rdx     # 0xe9cb1
pushq $0x57
popq %rcx
movq %rbx, %rdi
movq %r14, %r8
xorl %eax, %eax
callq 0x2e014
pushq $0x1
popq %rax
addq $0x830, %rsp            # imm = 0x830
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: background_at_tau
int background_at_tau(
                      struct background *pba,
                      double tau,
                      enum vecback_format return_format,
                      enum interpolation_method inter_mode,
                      int * last_index,
                      double * pvecback /* vector with argument pvecback[index_bg] (must be already allocated with a size compatible with return_format) */
                      ) {

  /** Summary: */

  /** - define local variables */
  double z;

  /** - Get current redshift */
  class_call(background_z_of_tau(pba,tau,&z),
             pba->error_message,
             pba->error_message);

  /** - Get background at corresponding redshift */
  class_call(background_at_z(pba,z,return_format,inter_mode,last_index,pvecback),
             pba->error_message,
             pba->error_message);

  return _SUCCESS_;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x810, %rsp            # imm = 0x810
movq %r8, %r14
movq %rcx, %r15
movl %edx, %ebp
movl %esi, %r12d
movq %rdi, %rbx
leaq 0x8(%rsp), %rsi
callq 0x58962
testl %eax, %eax
je 0x588e7
addq $0xd1c, %rbx            # imm = 0xD1C
leaq 0x8cb3a(%rip), %rsi     # 0xe53f0
leaq 0xa94a4(%rip), %rdx     # 0x101d61
leaq 0x10(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x2e014
leaq 0x8c75e(%rip), %rsi     # 0xe5034
leaq 0xa94a4(%rip), %rdx     # 0x101d81
movq %rbx, %rdi
movl $0xfa, %ecx
jmp 0x58941
movsd 0x8(%rsp), %xmm0
movq %rbx, %rdi
movl %r12d, %esi
movl %ebp, %edx
movq %r15, %rcx
movq %r14, %r8
callq 0x585e8
testl %eax, %eax
je 0x5895e
addq $0xd1c, %rbx            # imm = 0xD1C
leaq 0x8cade(%rip), %rsi     # 0xe53f0
leaq 0xa947a(%rip), %rdx     # 0x101d93
leaq 0x10(%rsp), %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x2e014
leaq 0x8c702(%rip), %rsi     # 0xe5034
leaq 0xa9448(%rip), %rdx     # 0x101d81
movq %rbx, %rdi
movl $0xff, %ecx
movq %r14, %r8
xorl %eax, %eax
callq 0x2e014
pushq $0x1
popq %rax
addq $0x810, %rsp            # imm = 0x810
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
retq
xorl %eax, %eax
jmp 0x5894e
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: V_scf
double V_scf(
             struct background *pba,
             double phi) {
  /** we check first which potential should be considered */
  double result = 0.;
  if(pba->scf_potential == pol_times_exp){
    result =  V_e_scf(pba,phi)*V_p_scf(pba,phi);
  }
  else if(pba->scf_potential == double_exp){
    result = V_double_exp_scf(pba,phi);
  }
  else if(pba->scf_potential == axion){
    result = V_axion_scf(pba,phi);
  }
  else if(pba->scf_potential == phi_2n){
    result = V_phi_2n_scf(pba,phi);
  }
  else if(pba->scf_potential == axionquad){
    result = V_axionquad_scf(pba,phi);
  }
    // printf("result Vf %e\n", result);

  return result;
}

// Assembly:
pushq %rbx
subq $0x20, %rsp
movl 0x160(%rdi), %eax
cmpq $0x4, %rax
ja 0x59be3
movq %rdi, %rbx
leaq 0xa7f2c(%rip), %rcx     # 0x101a60
movslq (%rcx,%rax,4), %rax
addq %rcx, %rax
jmpq *%rax
movq %rbx, %rdi
movsd %xmm0, 0x10(%rsp)
callq 0x618e2
movsd %xmm0, 0x8(%rsp)
movq 0x168(%rbx), %rax
movsd 0x8(%rax), %xmm1
movsd 0x10(%rax), %xmm0
movsd %xmm0, 0x18(%rsp)
movsd 0x10(%rsp), %xmm0
subsd 0x18(%rax), %xmm0
callq 0x8650
addsd 0x18(%rsp), %xmm0
jmp 0x59bce
movsd 0x188(%rbx), %xmm1
mulsd (%rbx), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm0, %xmm0
mulsd %xmm1, %xmm0
mulsd 0x8c110(%rip), %xmm0    # 0xe5cb0
jmp 0x59be7
movq %rbx, %rdi
addq $0x20, %rsp
popq %rbx
jmp 0x61c28
movsd 0x1c0(%rbx), %xmm1
movsd %xmm1, 0x8(%rsp)
movsd 0x210(%rbx), %xmm1
addsd %xmm1, %xmm1
callq 0x8650
mulsd 0x8(%rsp), %xmm0
jmp 0x59be7
movq %rbx, %rdi
addq $0x20, %rsp
popq %rbx
jmp 0x619da
xorpd %xmm0, %xmm0
addq $0x20, %rsp
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: dV_scf
double dV_scf(
              struct background *pba,
	      double phi) {
  /** we check first which potential should be considered */
  double result = 0.;
  if(pba->scf_potential == pol_times_exp){
    result =  dV_e_scf(pba,phi)*V_p_scf(pba,phi) + V_e_scf(pba,phi)*dV_p_scf(pba,phi);
  }
  else if(pba->scf_potential == double_exp){
    result =  dV_double_exp_scf(pba,phi);
  }
  else if(pba->scf_potential == axion){
    result = dV_axion_scf(pba,phi);
  }
  else if(pba->scf_potential == phi_2n){
    result = dV_phi_2n_scf(pba,phi);
  }
  else if(pba->scf_potential == axionquad){
    result = dV_axionquad_scf(pba,phi);
  }

  // printf("result dVf %e\n", result);

  return result;

}

// Assembly:
pushq %rbx
subq $0x30, %rsp
movl 0x160(%rdi), %eax
cmpq $0x4, %rax
ja 0x59d37
movapd %xmm0, %xmm2
movq %rdi, %rbx
leaq 0xa7e64(%rip), %rcx     # 0x101a74
movslq (%rcx,%rax,4), %rax
addq %rcx, %rax
jmpq *%rax
movq %rbx, %rdi
movapd %xmm2, %xmm0
movsd %xmm2, 0x8(%rsp)
callq 0x618fe
movsd %xmm0, 0x28(%rsp)
movq 0x168(%rbx), %rax
movsd 0x8(%rax), %xmm1
movsd 0x10(%rax), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
subsd 0x18(%rax), %xmm0
callq 0x8650
addsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movq %rbx, %rdi
movsd 0x8(%rsp), %xmm0
callq 0x618e2
movsd %xmm0, 0x20(%rsp)
movq 0x168(%rbx), %rax
movsd 0x8(%rax), %xmm2
movsd %xmm2, 0x18(%rsp)
movsd 0x8(%rsp), %xmm0
subsd 0x18(%rax), %xmm0
movsd 0x8bb4b(%rip), %xmm1    # 0xe57e8
addsd %xmm2, %xmm1
callq 0x8650
mulsd 0x18(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
movsd 0x10(%rsp), %xmm1
mulsd 0x28(%rsp), %xmm1
addsd %xmm1, %xmm0
jmp 0x59d3b
movsd 0x188(%rbx), %xmm0
mulsd (%rbx), %xmm0
mulsd %xmm0, %xmm0
mulsd %xmm2, %xmm0
jmp 0x59d3b
movq %rbx, %rdi
movapd %xmm2, %xmm0
addq $0x30, %rsp
popq %rbx
jmp 0x61cb5
movsd 0x1c0(%rbx), %xmm0
movsd 0x210(%rbx), %xmm1
addsd %xmm0, %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, 0x8(%rsp)
addsd %xmm1, %xmm1
addsd 0x8bad3(%rip), %xmm1    # 0xe57e8
movapd %xmm2, %xmm0
callq 0x8650
mulsd 0x8(%rsp), %xmm0
jmp 0x59d3b
movq %rbx, %rdi
movapd %xmm2, %xmm0
addq $0x30, %rsp
popq %rbx
jmp 0x61a82
xorpd %xmm0, %xmm0
addq $0x30, %rsp
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: ddV_scf
double ddV_scf(
               struct background *pba,
               double phi) {
  /** we check first which potential should be considered */
  double result = 0.;

  if(pba->scf_potential == pol_times_exp){
    result =  ddV_e_scf(pba,phi)*V_p_scf(pba,phi) + 2*dV_e_scf(pba,phi)*dV_p_scf(pba,phi) + V_e_scf(pba,phi)*ddV_p_scf(pba,phi);
  }

  else if(pba->scf_potential == double_exp){
    result =  ddV_double_exp_scf(pba,phi);
  }
  else if(pba->scf_potential == axion){
    result = ddV_axion_scf(pba,phi);
  }
  else if(pba->scf_potential == phi_2n){
    result = ddV_phi_2n_scf(pba,phi);
  }
  else if(pba->scf_potential == axionquad){
    result = ddV_axionquad_scf(pba,phi);
  }


  // printf("result ddVf %e\n", result);
  return result;

}

// Assembly:
pushq %rbx
subq $0x30, %rsp
movl 0x160(%rdi), %eax
cmpq $0x4, %rax
ja 0x59e81
movq %rdi, %rbx
leaq 0xa7d28(%rip), %rcx     # 0x101a88
movslq (%rcx,%rax,4), %rax
addq %rcx, %rax
jmpq *%rax
movq %rbx, %rdi
movsd %xmm0, 0x10(%rsp)
callq 0x61927
movsd %xmm0, 0x18(%rsp)
movq 0x168(%rbx), %rax
movsd 0x8(%rax), %xmm1
movsd 0x10(%rax), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
subsd 0x18(%rax), %xmm0
callq 0x8650
addsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movq %rbx, %rdi
movsd 0x10(%rsp), %xmm0
callq 0x618fe
addsd %xmm0, %xmm0
movsd %xmm0, 0x28(%rsp)
movq 0x168(%rbx), %rax
movsd 0x8(%rax), %xmm2
movsd %xmm2, 0x20(%rsp)
movsd 0x10(%rsp), %xmm0
subsd 0x18(%rax), %xmm0
movsd 0x8b9fb(%rip), %xmm1    # 0xe57e8
addsd %xmm2, %xmm1
callq 0x8650
mulsd 0x20(%rsp), %xmm0
mulsd 0x28(%rsp), %xmm0
movsd 0x8(%rsp), %xmm1
mulsd 0x18(%rsp), %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, 0x18(%rsp)
movq %rbx, %rdi
movsd 0x10(%rsp), %xmm0
callq 0x618e2
movsd %xmm0, 0x8(%rsp)
movq %rbx, %rdi
movsd 0x10(%rsp), %xmm0
callq 0x6199f
mulsd 0x8(%rsp), %xmm0
addsd 0x18(%rsp), %xmm0
jmp 0x59e85
movsd 0x188(%rbx), %xmm0
mulsd (%rbx), %xmm0
mulsd %xmm0, %xmm0
jmp 0x59e85
movq %rbx, %rdi
addq $0x30, %rsp
popq %rbx
jmp 0x61d62
movq %rbx, %rdi
addq $0x30, %rsp
popq %rbx
jmp 0x61f51
movq %rbx, %rdi
addq $0x30, %rsp
popq %rbx
jmp 0x61b58
xorpd %xmm0, %xmm0
addq $0x30, %rsp
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../include/background.h
// Function: gsl_sf_gamma
static double gsl_sf_gamma(double z) {
    /* Lanczos coefficients for g = 7 */
    static double p[] = {
        0.99999999999980993227684700473478,
        676.520368121885098567009190444019,
       -1259.13921672240287047156078755283,
        771.3234287776530788486528258894,
       -176.61502916214059906584551354,
        12.507343278686904814458936853,
       -0.13857109526572011689554707,
        9.984369578019570859563e-6,
        1.50563273514931155834e-7
    };

    // if(creal(z) < 0.5)
    //     return M_PI / (sin(M_PI*z)*gsl_sf_gamma(1. - z));

    z -= 1;
    // double complex x = p[0];
    double x = p[0];
    int n;
    for(n = 1; n < 9; n++)
        x += p[n] / (z + 1.*n);
    // double complex t = z + 7.5;
    double t = z + 7.5;
    return sqrt(2*_PI_) * pow(t, z+0.5) * exp(-t) * x;
}

// Assembly:
movapd %xmm0, %xmm1
addsd 0x89b28(%rip), %xmm1    # 0xe57e8
movsd 0xa5c90(%rip), %xmm3    # 0x101958
pushq $0x1
popq %rax
leaq 0xa5dce(%rip), %rcx     # 0x101aa0
cmpq $0x9, %rax
je 0x5bcf2
movsd (%rcx,%rax,8), %xmm0
cvtsi2sd %eax, %xmm2
addsd %xmm1, %xmm2
divsd %xmm2, %xmm0
addsd %xmm0, %xmm3
incq %rax
jmp 0x5bcd2
subq $0x28, %rsp
movsd 0xa5c62(%rip), %xmm0    # 0x101960
addsd %xmm1, %xmm0
movapd %xmm0, 0x10(%rsp)
addsd 0x89fa0(%rip), %xmm1    # 0xe5cb0
movsd %xmm3, 0x8(%rsp)
callq 0x8650
mulsd 0xa5c45(%rip), %xmm0    # 0x101968
movsd %xmm0, (%rsp)
movapd 0x10(%rsp), %xmm0
xorpd 0x894fa(%rip), %xmm0    # 0xe5230
callq 0x86a0
mulsd (%rsp), %xmm0
mulsd 0x8(%rsp), %xmm0
addq $0x28, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: background_free_noinput
int background_free_noinput(
                            struct background *pba
                            ) {

  free(pba->tau_table);
  free(pba->z_table);
  free(pba->loga_table);
  free(pba->d2tau_dz2_table);
  free(pba->d2z_dtau2_table);
  free(pba->background_table);
  free(pba->d2background_dloga2_table);


  return _SUCCESS_;
}

// Assembly:
pushq %rbx
movq %rdi, %rbx
movq 0x450(%rdi), %rdi
callq 0x8430
movq 0x458(%rbx), %rdi
callq 0x8430
movq 0x448(%rbx), %rdi
callq 0x8430
movq 0x468(%rbx), %rdi
callq 0x8430
movq 0x470(%rbx), %rdi
callq 0x8430
movq 0x460(%rbx), %rdi
callq 0x8430
movq 0x478(%rbx), %rdi
callq 0x8430
xorl %eax, %eax
popq %rbx
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: background_free_input
int background_free_input(
                          struct background *pba
                          ) {

  int k;

  if (pba->Omega0_ncdm_tot != 0.) {
    for (k=0; k<pba->N_ncdm; k++) {
      free(pba->q_ncdm[k]);
      free(pba->w_ncdm[k]);
      free(pba->q_ncdm_bg[k]);
      free(pba->w_ncdm_bg[k]);
      free(pba->dlnf0_dlnq_ncdm[k]);
    }

    free(pba->ncdm_quadrature_strategy);
    free(pba->ncdm_input_q_size);
    free(pba->ncdm_qmax);
    free(pba->q_ncdm);
    free(pba->w_ncdm);
    free(pba->q_ncdm_bg);
    free(pba->w_ncdm_bg);
    free(pba->dlnf0_dlnq_ncdm);
    free(pba->q_size_ncdm);
    free(pba->q_size_ncdm_bg);
    free(pba->M_ncdm);
    free(pba->T_ncdm);
    free(pba->ksi_ncdm);
    free(pba->deg_ncdm);
    free(pba->Omega0_ncdm);
    free(pba->m_ncdm_in_eV);
    free(pba->factor_ncdm);
    if (pba->got_files!=NULL)
      free(pba->got_files);
    if (pba->ncdm_psd_files!=NULL)
      free(pba->ncdm_psd_files);
    if (pba->ncdm_psd_parameters!=NULL)
      free(pba->ncdm_psd_parameters);
  }

  if (pba->Omega0_scf != 0.  || pba->log10_fraction_axion_ac > -30. || pba->log10_axion_ac > -30 || pba->m_scf != 0 || pba->f_axion != 0){
    if (pba->scf_parameters != NULL)
      free(pba->scf_parameters);
  }
  return _SUCCESS_;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
movsd 0x288(%rdi), %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0x5e021
jnp 0x5e179
xorl %r14d, %r14d
movslq 0x248(%rbx), %rax
cmpq %rax, %r14
jge 0x5e085
movq 0x4e8(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x8430
movq 0x4f0(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x8430
movq 0x4d8(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x8430
movq 0x4e0(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x8430
movq 0x4f8(%rbx), %rax
movq (%rax,%r14,8), %rdi
callq 0x8430
incq %r14
jmp 0x5e024
movq 0x4d0(%rbx), %rdi
callq 0x8430
movq 0x2c0(%rbx), %rdi
callq 0x8430
movq 0x2c8(%rbx), %rdi
callq 0x8430
movq 0x4e8(%rbx), %rdi
callq 0x8430
movq 0x4f0(%rbx), %rdi
callq 0x8430
movq 0x4d8(%rbx), %rdi
callq 0x8430
movq 0x4e0(%rbx), %rdi
callq 0x8430
movq 0x4f8(%rbx), %rdi
callq 0x8430
movq 0x508(%rbx), %rdi
callq 0x8430
movq 0x500(%rbx), %rdi
callq 0x8430
movq 0x270(%rbx), %rdi
callq 0x8430
movq 0x290(%rbx), %rdi
callq 0x8430
movq 0x2a0(%rbx), %rdi
callq 0x8430
movq 0x2b0(%rbx), %rdi
callq 0x8430
movq 0x280(%rbx), %rdi
callq 0x8430
movq 0x278(%rbx), %rdi
callq 0x8430
movq 0x510(%rbx), %rdi
callq 0x8430
movq 0x258(%rbx), %rdi
callq 0x8430
movq 0x250(%rbx), %rdi
callq 0x8430
movq 0x260(%rbx), %rdi
callq 0x8430
xorpd %xmm1, %xmm1
movsd 0x128(%rbx), %xmm0
ucomisd %xmm1, %xmm0
jne 0x5e1d1
jp 0x5e1d1
movsd 0x1d0(%rbx), %xmm0
ucomisd 0x8b897(%rip), %xmm0    # 0xe9a30
ja 0x5e1d1
movsd 0x1b8(%rbx), %xmm0
ucomisd 0x8b885(%rip), %xmm0    # 0xe9a30
ja 0x5e1d1
movsd 0x188(%rbx), %xmm1
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
jne 0x5e1d1
jp 0x5e1d1
movsd 0x190(%rbx), %xmm1
ucomisd %xmm0, %xmm1
jne 0x5e1d1
jnp 0x5e1dd
movq 0x168(%rbx), %rdi
callq 0x8430
xorl %eax, %eax
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: background_sources
int background_sources(
                       double loga,
                       double * y,
                       double * dy,
                       int index_loga,
                       void * parameters_and_workspace,
                       ErrorMsg error_message
                       ) {

  struct background_parameters_and_workspace * pbpaw;
  struct background * pba;
  double a;
  double * bg_table_row;

  pbpaw = parameters_and_workspace;
  pba =  pbpaw->pba;

  /** - localize the row inside background_table where the current values must be stored */
  bg_table_row = pba->background_table + index_loga*pba->bg_size;

  /** - scale factor a (in fact, given our normalisation conventions, this stands for a/a_0) */
  a = exp(loga);

  /** - corresponding redhsift 1/a-1 */
  pba->z_table[index_loga] = MAX(0.,1./a-1.);

  /** - corresponding conformal time */
  pba->tau_table[index_loga] = y[pba->index_bi_tau];

  /** -> compute all other quantities depending only on a + {B} variables and get them stored
      in one row of background_table
      The value of {B} variables in pData are also copied to pvecback.*/
  class_call(background_functions(pba, a, y, long_info, bg_table_row),
             pba->error_message,
             pba->error_message);

  return _SUCCESS_;

}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x808, %rsp            # imm = 0x808
movq %rdi, %r14
movq (%rcx), %rbx
movslq 0x43c(%rbx), %r15
movslq %edx, %r12
imulq %r12, %r15
shlq $0x3, %r15
addq 0x460(%rbx), %r15
callq 0x86a0
movsd 0x84730(%rip), %xmm1    # 0xe5260
divsd %xmm0, %xmm1
addsd 0x84cac(%rip), %xmm1    # 0xe57e8
xorpd %xmm2, %xmm2
maxsd %xmm2, %xmm1
movq 0x458(%rbx), %rax
movsd %xmm1, (%rax,%r12,8)
movslq 0x4a0(%rbx), %rax
movsd (%r14,%rax,8), %xmm1
movq 0x450(%rbx), %rax
movsd %xmm1, (%rax,%r12,8)
pushq $0x2
popq %rdx
movq %rbx, %rdi
movq %r14, %rsi
movq %r15, %rcx
callq 0x58c54
testl %eax, %eax
je 0x60bca
addq $0xd1c, %rbx            # imm = 0xD1C
leaq 0x84862(%rip), %rsi     # 0xe53f0
leaq 0xa3460(%rip), %rdx     # 0x103ff5
movq %rsp, %r14
movq %r14, %rdi
movq %rbx, %rcx
xorl %eax, %eax
callq 0x2e014
leaq 0x84488(%rip), %rsi     # 0xe5034
leaq 0xa347b(%rip), %rdx     # 0x10402e
movq %rbx, %rdi
movl $0xd1b, %ecx            # imm = 0xD1B
movq %r14, %r8
xorl %eax, %eax
callq 0x2e014
pushq $0x1
popq %rax
jmp 0x60bcc
xorl %eax, %eax
addq $0x808, %rsp            # imm = 0x808
popq %rbx
popq %r12
popq %r14
popq %r15
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: ddV_e_scf
double ddV_e_scf(struct background *pba,
                 double phi
                 ) {
  double scf_lambda = pba->scf_parameters[0];
  //  double scf_alpha  = pba->scf_parameters[1];
  //  double scf_A      = pba->scf_parameters[2];
  //  double scf_B      = pba->scf_parameters[3];

  return pow(-scf_lambda,2)*V_e_scf(pba,phi);
}

// Assembly:
pushq %rax
movq 0x168(%rdi), %rax
movsd (%rax), %xmm1
mulsd %xmm1, %xmm1
movsd %xmm1, (%rsp)
callq 0x618e2
mulsd (%rsp), %xmm0
popq %rax
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: V_p_scf
double V_p_scf(
               struct background *pba,
               double phi) {
  //  double scf_lambda = pba->scf_parameters[0];
  double scf_alpha  = pba->scf_parameters[1];
  double scf_A      = pba->scf_parameters[2];
  double scf_B      = pba->scf_parameters[3];

  return  pow(phi - scf_B,  scf_alpha) +  scf_A;
}

// Assembly:
pushq %rax
movq 0x168(%rdi), %rax
movsd 0x8(%rax), %xmm1
movsd 0x10(%rax), %xmm2
movsd %xmm2, (%rsp)
subsd 0x18(%rax), %xmm0
callq 0x8650
addsd (%rsp), %xmm0
popq %rax
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: dV_p_scf
double dV_p_scf(
                struct background *pba,
                double phi) {

  //  double scf_lambda = pba->scf_parameters[0];
  double scf_alpha  = pba->scf_parameters[1];
  //  double scf_A      = pba->scf_parameters[2];
  double scf_B      = pba->scf_parameters[3];

  return   scf_alpha*pow(phi -  scf_B,  scf_alpha - 1);
}

// Assembly:
pushq %rax
movq 0x168(%rdi), %rax
movsd 0x8(%rax), %xmm2
movsd %xmm2, (%rsp)
subsd 0x18(%rax), %xmm0
movsd 0x83e59(%rip), %xmm1    # 0xe57e8
addsd %xmm2, %xmm1
callq 0x8650
mulsd (%rsp), %xmm0
popq %rax
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: V_axion_scf
double V_axion_scf(
                  struct background *pba,
                  double phi){
    // int n = pba->scf_parameters[0];
    double n = pba->n_axion;
    // double fa = pba->scf_parameters[2];
    double fa = pba->f_axion;
    double m = pba->m_scf*pba->H0;
    double result;
    // printf("n %d fa %e V %e phi/fa %e \n",n,fa,m*m/pow(2,n),phi/fa);
    if(n>1)result = pow(m,2)*pow(fa,2)*pow(1 - cos(phi/fa),n);
    else result = pow(m,2)*pow(fa,2)*(1 - cos(phi/fa));
    // printf("result %e phi %e m^2 %e\n",result,phi,m*m);
    return result;

}

// Assembly:
subq $0x28, %rsp
movsd 0x210(%rdi), %xmm1
movsd %xmm1, 0x8(%rsp)
movsd 0x190(%rdi), %xmm2
movsd (%rdi), %xmm1
movsd %xmm1, 0x10(%rsp)
movsd 0x188(%rdi), %xmm1
movsd %xmm1, 0x20(%rsp)
movsd %xmm2, 0x18(%rsp)
divsd %xmm2, %xmm0
callq 0x86d0
movsd 0x8(%rsp), %xmm1
movapd %xmm0, %xmm2
movsd 0x835e5(%rip), %xmm0    # 0xe5260
ucomisd %xmm0, %xmm1
subsd %xmm2, %xmm0
jbe 0x61c8a
callq 0x8650
movsd 0x18(%rsp), %xmm2
mulsd %xmm2, %xmm2
movsd 0x20(%rsp), %xmm1
mulsd 0x10(%rsp), %xmm1
mulsd %xmm1, %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, %xmm0
addq $0x28, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: dV_axion_scf
double dV_axion_scf(
                  struct background *pba,
                  double phi){
    // int n = pba->scf_parameters[0];
    double n = pba->n_axion;
    // double fa = pba->scf_parameters[2];
    double fa = pba->f_axion;
    double m = pba->m_scf*pba->H0;
    double result;
    if(n>1)result = n*pow(m,2)*fa*pow(1-cos(phi/fa),n-1)*sin(phi/fa);
    else result = pow(m,2)*fa*sin(phi/fa);

    return result;

}

// Assembly:
subq $0x18, %rsp
movapd %xmm0, %xmm1
movsd 0x210(%rdi), %xmm3
movsd 0x188(%rdi), %xmm2
movsd 0x190(%rdi), %xmm0
mulsd (%rdi), %xmm2
ucomisd 0x8357f(%rip), %xmm3    # 0xe5260
mulsd %xmm2, %xmm2
jbe 0x61d42
mulsd %xmm3, %xmm2
mulsd %xmm0, %xmm2
movsd %xmm2, (%rsp)
divsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rsp)
movapd %xmm1, %xmm0
movsd %xmm3, 0x8(%rsp)
callq 0x86d0
movsd 0x8354b(%rip), %xmm2    # 0xe5260
subsd %xmm0, %xmm2
movsd 0x8(%rsp), %xmm1
addsd 0x83ac1(%rip), %xmm1    # 0xe57e8
movapd %xmm2, %xmm0
callq 0x8650
movsd 0x10(%rsp), %xmm1
mulsd (%rsp), %xmm0
movsd %xmm0, (%rsp)
jmp 0x61d4f
mulsd %xmm0, %xmm2
divsd %xmm0, %xmm1
movsd %xmm2, (%rsp)
movapd %xmm1, %xmm0
callq 0x80a0
mulsd (%rsp), %xmm0
addq $0x18, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: ddV_axion_scf
double ddV_axion_scf(
                  struct background *pba,
                  double phi){

     // int n = pba->scf_parameters[0];
     double n = pba->n_axion;
     // double fa = pba->scf_parameters[2];
     double fa = pba->f_axion;
     double m = pba->m_scf*pba->H0;
     double result;
     if(n==1) result = n*pow(m,2)*cos(phi/fa);
     else if (n==2) result =  n*pow(m,2)*(pow(sin(phi/fa),2)+(1-cos(phi/fa))*cos(phi/fa));
     else result = n*pow(m,2)*fa*((n-1)/fa*pow(1-cos(phi/fa),n-2)*pow(sin(phi/fa),2)+pow(1-cos(phi/fa),n-1)/fa*cos(phi/fa)); //this formula bugs sometimes for n=1

     return result;
}

// Assembly:
subq $0x38, %rsp
movsd 0x210(%rdi), %xmm3
movsd 0x188(%rdi), %xmm1
mulsd (%rdi), %xmm1
movsd 0x190(%rdi), %xmm2
ucomisd 0x834d6(%rip), %xmm3    # 0xe5260
jne 0x61da6
jp 0x61da6
mulsd %xmm1, %xmm1
movsd %xmm1, 0x20(%rsp)
divsd %xmm2, %xmm0
callq 0x86d0
jmp 0x61eea
ucomisd 0x83f32(%rip), %xmm3    # 0xe5ce0
mulsd %xmm1, %xmm1
jne 0x61e12
jp 0x61e12
addsd %xmm1, %xmm1
movsd %xmm1, 0x20(%rsp)
divsd %xmm2, %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0x80a0
mulsd %xmm0, %xmm0
movsd %xmm0, 0x18(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x86d0
movsd 0x83474(%rip), %xmm1    # 0xe5260
subsd %xmm0, %xmm1
movsd %xmm1, 0x28(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x86d0
mulsd 0x28(%rsp), %xmm0
addsd 0x18(%rsp), %xmm0
jmp 0x61eea
mulsd %xmm3, %xmm1
mulsd %xmm2, %xmm1
movsd %xmm1, 0x20(%rsp)
movsd 0x839c0(%rip), %xmm1    # 0xe57e8
addsd %xmm3, %xmm1
movsd %xmm1, 0x18(%rsp)
divsd %xmm2, %xmm1
movsd %xmm1, 0x10(%rsp)
divsd %xmm2, %xmm0
movsd %xmm0, 0x8(%rsp)
movsd %xmm2, 0x28(%rsp)
movsd %xmm3, 0x30(%rsp)
callq 0x86d0
movsd 0x83401(%rip), %xmm2    # 0xe5260
subsd %xmm0, %xmm2
movsd 0x30(%rsp), %xmm1
addsd 0x848c7(%rip), %xmm1    # 0xe6738
movapd %xmm2, %xmm0
callq 0x8650
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x80a0
mulsd %xmm0, %xmm0
mulsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x86d0
movsd 0x833ac(%rip), %xmm1    # 0xe5260
subsd %xmm0, %xmm1
movapd %xmm1, %xmm0
movsd 0x18(%rsp), %xmm1
callq 0x8650
divsd 0x28(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
movsd 0x8(%rsp), %xmm0
callq 0x86d0
mulsd 0x18(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
mulsd 0x20(%rsp), %xmm0
addq $0x38, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: dV_phi_2n_scf
double dV_phi_2n_scf(
                  struct background *pba,
                  double phi){
    // int n = pba->scf_parameters[0];
    double n = pba->n_axion;
    double result;

    result = pba->V0_phi2n*2*n*pow((phi),2*n-1);

    return result;

}

// Assembly:
pushq %rax
movsd 0x1c0(%rdi), %xmm2
movsd 0x210(%rdi), %xmm1
addsd %xmm2, %xmm2
mulsd %xmm1, %xmm2
movsd %xmm2, (%rsp)
addsd %xmm1, %xmm1
addsd 0x838a3(%rip), %xmm1    # 0xe57e8
callq 0x8650
mulsd (%rsp), %xmm0
popq %rax
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/background.c
// Function: ddV_phi_2n_scf
double ddV_phi_2n_scf(
                  struct background *pba,
                  double phi){

     // int n = pba->scf_parameters[0];
     double n = pba->n_axion;
     double result;
     if(n==1) result=pba->V0_phi2n*2;
     else result = pba->V0_phi2n*2*n*(2*n-1)*pow((phi),2*n-2);
     return result;
}

// Assembly:
movapd %xmm0, %xmm2
movsd 0x1c0(%rdi), %xmm0
movsd 0x210(%rdi), %xmm1
ucomisd 0x832f3(%rip), %xmm1    # 0xe5260
addsd %xmm0, %xmm0
jne 0x61f76
jp 0x61f76
retq
pushq %rax
mulsd %xmm1, %xmm0
addsd %xmm1, %xmm1
movsd 0x83861(%rip), %xmm3    # 0xe57e8
addsd %xmm1, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, (%rsp)
addsd 0x8479c(%rip), %xmm1    # 0xe6738
movapd %xmm2, %xmm0
callq 0x8650
mulsd (%rsp), %xmm0
addq $0x8, %rsp
retq
// File: /workspace/github_c_makefile_star2_O2/PoulinV[P]AxiCLASS/build/../source/thermodynamics.c
// Function: thermodynamics_checks
int thermodynamics_checks(
                          struct precision * ppr,
                          struct background* pba,
                          struct thermodynamics * pth
                          ) {

  /** Summary: */

  /** - check BBN Y_He fracion */
  class_test((pth->YHe < _YHE_SMALL_)||(pth->YHe > _YHE_BIG_),
             pth->error_message,
             "Y_He=%g out of bounds (%g<Y_He<%g)",pth->YHe,_YHE_SMALL_,_YHE_BIG_);

  /** - tests in order to prevent divisions by zero */
  class_test(pth->YHe == 1.,
             pth->error_message,
             "stop to avoid division by zero");

  /** - test initial condition for recombination */
  class_test(ppr->thermo_z_initial < ppr->recfast_z_He_3,
             pth->error_message,
             "increase z_initial in order to start before HeliumIII recombination");

  return _SUCCESS_;
}

// Assembly:
pushq %r15
pushq %r14
pushq %rbx
subq $0x1000, %rsp           # imm = 0x1000
movq %rdx, %rbx
movsd (%rdx), %xmm0
movsd 0x84048(%rip), %xmm1    # 0xe7f80
ucomisd %xmm0, %xmm1
ja 0x63fb2
ucomisd 0x81d6a(%rip), %xmm0    # 0xe5cb0
ja 0x63fb2
ucomisd 0x81310(%rip), %xmm0    # 0xe5260
jne 0x64030
jp 0x64030
leaq 0x82d9a(%rip), %rsi     # 0xe6cfd
leaq 0x800(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x2e014
leaq 0x81106(%rip), %rsi     # 0xe5082
leaq 0xa14e0(%rip), %rdx     # 0x105463
movq %rsp, %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x2e014
addq $0x1064, %rbx           # imm = 0x1064
leaq 0x81093(%rip), %rsi     # 0xe5034
leaq 0xa14a5(%rip), %rdx     # 0x10544d
movq %rbx, %rdi
movl $0x2f4, %ecx            # imm = 0x2F4
jmp 0x64016
leaq 0xa1440(%rip), %rsi     # 0x1053f9
leaq 0x800(%rsp), %r14
movsd 0x83fb7(%rip), %xmm1    # 0xe7f80
movsd 0x81cdf(%rip), %xmm2    # 0xe5cb0
movq %r14, %rdi
movb $0x3, %al
callq 0x2e014
leaq 0x810a0(%rip), %rsi     # 0xe5082
leaq 0xa1433(%rip), %rdx     # 0x10541c
movq %rsp, %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x2e014
addq $0x1064, %rbx           # imm = 0x1064
leaq 0x8102d(%rip), %rsi     # 0xe5034
leaq 0xa143f(%rip), %rdx     # 0x10544d
movq %rbx, %rdi
movl $0x2ef, %ecx            # imm = 0x2EF
movq %r15, %r8
xorl %eax, %eax
callq 0x2e014
pushq $0x1
popq %rax
addq $0x1000, %rsp           # imm = 0x1000
popq %rbx
popq %r14
popq %r15
retq
movsd 0x278(%rdi), %xmm0
xorl %eax, %eax
ucomisd 0x190(%rdi), %xmm0
jbe 0x64023
leaq 0xa1427(%rip), %rsi     # 0x105472
leaq 0x800(%rsp), %r14
movq %r14, %rdi
xorl %eax, %eax
callq 0x2e014
leaq 0x8101e(%rip), %rsi     # 0xe5082
leaq 0xa144b(%rip), %rdx     # 0x1054b6
movq %rsp, %r15
movq %r15, %rdi
movq %r14, %rcx
xorl %eax, %eax
callq 0x2e014
addq $0x1064, %rbx           # imm = 0x1064
leaq 0x80fab(%rip), %rsi     # 0xe5034
leaq 0xa13bd(%rip), %rdx     # 0x10544d
movq %rbx, %rdi
movl $0x2f9, %ecx            # imm = 0x2F9
jmp 0x64016