Save time-evolution results of PDE

[mistriotis-lab]

Hello,

First of all, kudos to the pypde team for this amazing package!

I am solving the Brusselator model at pypde and I would like to track and save the entire time-evolution of u,v. Perhaps save it in the 'sol._data' instead of just saving the final snapshot. What should I change in the tracker object to achieve that?

############
from pde import PDE, FieldCollection, PlotTracker, ScalarField, UnitGrid

define the PDE

a, b = 1, 3
d0, d1 = 1, 0.1
eq = PDE(
{
"u": f"{d0} * laplace(u) + {a} - ({b} + 1) * u + u2 * v",
"v": f"{d1} * laplace(v) + {b} * u - u2 * v",
}
)

initialize state

grid = UnitGrid([64, 64])
u = ScalarField(grid, a, label="Field $u$")
v = b / a + 0.1 * ScalarField.random_normal(grid, label="Field $v$")
state = FieldCollection([u, v])

simulate the pde

tracker = PlotTracker(interval=1, plot_args={"vmin": 0, "vmax": 5})
sol = eq.solve(state, t_range=20, dt=1e-3, tracker=tracker)
#############

[david-zwicker]

You can save intermediate data at regular intervals using the storage module; see this example or this example.

[mistriotis-lab]

Thank you for your response!