Public repository of the Cosmic Linear Anisotropy Solving System (master for the most recent version of the standard code; GW_CLASS to include Cosmic Gravitational Wave Background anisotropies; classnet branch for acceleration with neutral networks; ExoCLASS branch for exotic energy injection; class_matter branch for FFTlog)
I am running CLASS v 2.9.4. I was making some changes in primordial module. As it was given in the comment section of line 1039 for the shape of potential. I wrote the potential as mentioned and when I ran the class from ./class .ini file commenting the shape of potential as my case. But it was running the case of polynomial one. Then I got to know that it can run only the polynomial and higgs inflation one . So , I made the chage in the polynomial case in place of V0+phi*V1+phi*2V2 + ...... one , I changed it to my desired potential form. And I made the .ini as follows :
modes = 's,t'
But I am getting the following error:
shanno@shanno:~/Desktop/new_class/class_public-2.9.4$ ./class inflation.ini
Reading input parameters
-> matched budget equations by adjusting Omega_Lambda = 6.878622e-01
Running CLASS version v2.9.4
Computing background
-> age = 13.795359 Gyr
-> conformal age = 14165.045412 Mpc
-> radiation/matter equality at z = 3402.964903
corresponding to conformal time = 112.800022 Mpc
Computing thermodynamics with Y_He=0.2453
-> recombination at z = 1089.184869 (max of visibility function)
corresponding to conformal time = 280.591903 Mpc
with comoving sound horizon = 144.703150 Mpc
angular diameter distance = 12.735871 Mpc
and sound horizon angle 100thetas = 1.042196
Thomson optical depth crosses one at z = 1082.099523
giving an angle 100theta_ = 1.046775
-> baryon drag stops at z = 1059.089246
corresponding to conformal time = 286.504945 Mpc
with comoving sound horizon rs = 147.384053 Mpc
-> reionization with optical depth = 0.092664
corresponding to conformal time = 4255.316282 Mpc
Computing sources
Computing primordial spectra (simulating inflation)
Error in primordial_init
=>primordial_init(L:365) :error in primordial_inflation_solve_inflation(ppt,ppm,ppr);
=>primordial_inflation_solve_inflation(L:1458) :error in primordial_inflation_spectra(ppt, ppm, ppr, y_ini);
=>primordial_inflation_spectra(L:1675) :error in primordial_inflation_one_wavenumber(ppt,ppm,ppr,y_ini,index_k);
=>primordial_inflation_one_wavenumber(L:1759) :error in primordial_inflation_one_k(ppm, ppr, k, y, dy, &curvature, &tensors);
=>primordial_inflation_one_k(L:1887) :error in generic_integrator(primordial_inflation_derivs, tau_start, tau_end, y, &pipaw, ppr->primordial_inflation_tol_integration, ppr->smallest_allowed_variation, &gi);
=>generic_integrator(L:103) :error in rkqs(&x, h, eps, &hdid, &hnext, derivs, parameters_and_workspace_for_derivs, pgi);
=>rkqs(L:154) :condition (xnew == *x) is true; stepsize underflow at x=0.000000e+00
Can somebody suggest me what could be the problem. This would be very helpful.
Hi everyone ,
I am running CLASS v 2.9.4. I was making some changes in primordial module. As it was given in the comment section of line 1039 for the shape of potential. I wrote the potential as mentioned and when I ran the class from ./class .ini file commenting the shape of potential as my case. But it was running the case of polynomial one. Then I got to know that it can run only the polynomial and higgs inflation one . So , I made the chage in the polynomial case in place of V0+phi*V1+phi*2V2 + ...... one , I changed it to my desired potential form. And I made the .ini as follows : modes = 's,t'
mode = logarithmic
lensing = yes output = tCl,pCl,lCl,mPk P_k_ini type = inflation_V k_pivot = 0.05 potential = polynomial V_0=1.e+8 V_1=-2 V_2=1 V_3=3.e-4 V_4=0. full_potential = polynomial l_max_scalars = 2500 l_max_tensors = 500 P_k_max_h/Mpc = 1. z_pk = 0., 1.2, 3.5 root = output/inflation write primordial = yes write parameters = yeap write warnings = y input_verbose = 1 background_verbose = 1 thermodynamics_verbose = 1 perturbations_verbose = 1 transfer_verbose = 1 primordial_verbose = 1 lensing_verbose = 1 output_verbose = 1
But I am getting the following error: shanno@shanno:~/Desktop/new_class/class_public-2.9.4$ ./class inflation.ini Reading input parameters -> matched budget equations by adjusting Omega_Lambda = 6.878622e-01 Running CLASS version v2.9.4 Computing background -> age = 13.795359 Gyr -> conformal age = 14165.045412 Mpc -> radiation/matter equality at z = 3402.964903 corresponding to conformal time = 112.800022 Mpc Computing thermodynamics with Y_He=0.2453 -> recombination at z = 1089.184869 (max of visibility function) corresponding to conformal time = 280.591903 Mpc with comoving sound horizon = 144.703150 Mpc angular diameter distance = 12.735871 Mpc and sound horizon angle 100thetas = 1.042196 Thomson optical depth crosses one at z = 1082.099523 giving an angle 100theta_ = 1.046775 -> baryon drag stops at z = 1059.089246 corresponding to conformal time = 286.504945 Mpc with comoving sound horizon rs = 147.384053 Mpc -> reionization with optical depth = 0.092664 corresponding to conformal time = 4255.316282 Mpc Computing sources Computing primordial spectra (simulating inflation)
Error in primordial_init =>primordial_init(L:365) :error in primordial_inflation_solve_inflation(ppt,ppm,ppr); =>primordial_inflation_solve_inflation(L:1458) :error in primordial_inflation_spectra(ppt, ppm, ppr, y_ini); =>primordial_inflation_spectra(L:1675) :error in primordial_inflation_one_wavenumber(ppt,ppm,ppr,y_ini,index_k); =>primordial_inflation_one_wavenumber(L:1759) :error in primordial_inflation_one_k(ppm, ppr, k, y, dy, &curvature, &tensors); =>primordial_inflation_one_k(L:1887) :error in generic_integrator(primordial_inflation_derivs, tau_start, tau_end, y, &pipaw, ppr->primordial_inflation_tol_integration, ppr->smallest_allowed_variation, &gi); =>generic_integrator(L:103) :error in rkqs(&x, h, eps, &hdid, &hnext, derivs, parameters_and_workspace_for_derivs, pgi); =>rkqs(L:154) :condition (xnew == *x) is true; stepsize underflow at x=0.000000e+00
Can somebody suggest me what could be the problem. This would be very helpful.
regards, shan.