The integrator produces negative values

[Giappo]

For some parameter choices (e.g. very high numbers, see test-mbd_ml:R147 (pars = c(60, 50, 100, 0.5)). The integration yields negative probabilities. This is best exposed considering the result of the conditioning probability, calculated with "calculate_conditional_prob".

[Giappo]

The function check_pars now has a safety option.
If safety option is enabled then unrealistic parameter values will return a likelihood = -Inf.
The current domains for parameters are:
0 < lambda < 10
0 < mu < 10
0 < nu < 100
1e3 < q < 1 - 1e3