Solve PDE in Implicit Method

[beantsja]

Hi, I am trying to solve simple 1D Heat Equation, but I want the length of the rod to be around 1 micrometer and the total spatial grid points to be around 25. Under this condition, it is forces dt to be extremely small under explicit method, and diffusionPDE seems to use that method.

I noticed there is Implicit Solver in the modules. Based on that, how would I solve PDE using implicit method?

[david-zwicker]

You can simply specify eq.solve(..., solver="implicit") to pick an implicit solver; see documentation. If you want to have more control over the solver, you can also follow this example.

However, it might also be a good idea to increase the spatial discretization to accelerate the procedure.