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 was adding a fluid described by eqn(1) of https://arxiv.org/pdf/1905.12618.pdf
The code works fine for p = 1. However when I try to decrease p, I get the error in the perturb_init
=>perturb_init(L:784) :error in perturb_solve(ppr, pba, pth, ppt, index_md, index_ic, index_k, pppw[thread]);
=>perturb_solve(L:3045) :error in generic_evolver(perturb_derivs, interval_limit[index_interval], interval_limit[index_interval+1], ppw->pv->y, ppw->pv->used_in_sources, ppw->pv->pt_size, &ppaw, ppr->tol_perturb_integration, ppr->smallest_allowed_variation, perturb_timescale, ppr->perturb_integration_stepsize, ppt->tau_sampling, tau_actual_size, perturb_sources, perhaps_print_variables, ppt->error_message);
=>evolver_ndf15(L:295) :error in new_linearisation(&jac,hinvGak,neq,error_message);
=>new_linearisation(L:976) :condition (funcreturn == FAILURE) is true; Failure in sp_ludcmp. Possibly singular matrix!
I was adding a fluid described by eqn(1) of https://arxiv.org/pdf/1905.12618.pdf The code works fine for p = 1. However when I try to decrease p, I get the error in the perturb_init
=>perturb_init(L:784) :error in perturb_solve(ppr, pba, pth, ppt, index_md, index_ic, index_k, pppw[thread]); =>perturb_solve(L:3045) :error in generic_evolver(perturb_derivs, interval_limit[index_interval], interval_limit[index_interval+1], ppw->pv->y, ppw->pv->used_in_sources, ppw->pv->pt_size, &ppaw, ppr->tol_perturb_integration, ppr->smallest_allowed_variation, perturb_timescale, ppr->perturb_integration_stepsize, ppt->tau_sampling, tau_actual_size, perturb_sources, perhaps_print_variables, ppt->error_message); =>evolver_ndf15(L:295) :error in new_linearisation(&jac,hinvGak,neq,error_message); =>new_linearisation(L:976) :condition (funcreturn == FAILURE) is true; Failure in sp_ludcmp. Possibly singular matrix!
Any hint on what this error is would be helpful!