Page 1 of 1

### Adding a perturbation discontinuity in CAMB

Posted: March 15 2022
Hi everyone,
I am not an expert on CAMB, and currently I'm trying to introduce a model like the one on https://arxiv.org/pdf/2006.06686.pdf (in specific eqs. 85) where at an $a_{critic}$ some Early Dark Energy perturbations change immediatly from a value $0$ to some value $\delta^*$, $\theta^*$.
The problem is that when I manually try to introduce on camb ay(w_ix)$=\delta^*$, ay(w_ix+1)$=\theta^*$ when $\text{abs}(a-a_{critic})<\varepsilon$, I get a compilation time extremely large or directly I get the following error

Code: Select all

 Dverk error -3: the subroutine was unable  to  satisfy  the  error requirement  with a particular step-size that is less than or * equal to hmin, which may mean that tol is too small--- but most likely you've messed up the y array indexing; compiling with bounds checking may (or may not) help find the problem.
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG
ERROR STOP Invalid parameter

Error termination. Backtrace:
#0  0x10353d147
#1  0x10353ddbf
#2  0x10353f0ab
#3  0x10303034f
#4  0x103053e07
Is there some hack or any other way to define this perturbation discontinuities? Thanks in advance for your help