Comparison of Numerically Stable Methods for Implementation of a Double Porosity Model with First-Order Reaction Terms