Add a generic implicit solver for non-linear diffusion

[benbovy]

The solver could be reused for implementing marine transport / erosion (see https://doi.org/10.1016/j.epsl.2019.115728) or other processes.