Thank you for updating the Docker container! The script above is now working fine. However, this same KSPSolve error comes up in my own code, which I had recently made a separate post about here: [diffusion coupled with boundary equation (mixed-dimensional)].