This avoids type instability in rate and rate_derivative due to the small type union.
Using a single type is probably better in the long run, but need to keep an eye on rate_derivative for higher-order (order >> 2)reactions.
Coverage decreased (-0.8%) to 72.289% when pulling 743785577cfde91c1e81e5e8f494e12e96f5275e on develop into a2904f296a893930196df66cf4da0d9f98685438 on master.
This avoids type instability in
rate
andrate_derivative
due to the small type union. Using a single type is probably better in the long run, but need to keep an eye onrate_derivative
for higher-order (order >> 2)reactions.