This is a quick fix for a regridding issue that was present. I think we need to be careful about how and when solvers and operators (and temporary vectors) are destroyed. If they are destroyed after a regridding operation, then I think any scratch data allocated by the operator will be lost.
This is a quick fix for the issue. We just destroy all the operators and solvers after a time step, and recreate them at the begging of a time step.
This is a quick fix for a regridding issue that was present. I think we need to be careful about how and when solvers and operators (and temporary vectors) are destroyed. If they are destroyed after a regridding operation, then I think any scratch data allocated by the operator will be lost.
This is a quick fix for the issue. We just destroy all the operators and solvers after a time step, and recreate them at the begging of a time step.