FEM Heat Solvers
For nonstiff problems it's recommended you use
FEMDiffEqHeatEuler, while for stiff problems it's recommended that you use
Full List of Methods
Additionally, for linear solves, one can choose the method by which the linear solve takes place via the
method keyword argument.
sol = solve(prob,FEMDiffEqHeatCrankNicholson(),solver=:CG)