ka=0 leads to exception instead of eliminating absorption in one/two cut models.
Example
y0(0) = 745;
y0(1) = 100;
rate[0] = 1200;
rate[1] = 200;
ka = 0.0;
using model_t = PMXOneCptModel<double, double, double, double>;
model_t model(t0, y0, rate, CL, V2, ka);
auto y = model.solve(ts[0]);
The above model throws exception when check_positive_finite is applied to ka, removing this check will give NaN in output y because ka is used in denominators in the solver.
Description
ka=0leads to exception instead of eliminating absorption in one/two cut models.Example
The above model throws exception when
check_positive_finiteis applied toka, removing this check will giveNaNin outputybecausekais used in denominators in the solver.Expected Output
Solution with no absorption
Current Version:
v0.87.1