From 0404dc297ee01f98dbdc193ae8082eb271f3fabe Mon Sep 17 00:00:00 2001 From: Gabriel Gerlero Date: Tue, 3 Sep 2024 12:09:40 -0300 Subject: [PATCH] Fix capacity handling for FiniteReservoirProblem --- src/finite.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/finite.jl b/src/finite.jl index 7d32ee2..b5d5f0f 100644 --- a/src/finite.jl +++ b/src/finite.jl @@ -309,7 +309,7 @@ function solve(prob::FiniteReservoirProblem{<:DiffusionEquation{1}}, else t = min((prob._rstop / presol.oi)^2, prob._tstop, (prob.capacity / sorptivity(presol))^2) - u[begin] -= sorptivity(presol) * √t + u[begin] = prob.capacity - sorptivity(presol) * √t u[(begin + 1):end] .= presol.(r, t) u[begin + 1] = prob.b end