Add functional form of diffusion

[dimitra-maoutsa]

multiplicative noise

[dimitra-maoutsa]

Current version supports only scalar diffusion.

Find a way to integrate it with the version for matrix valued diffusion.
Preferably without converting the scalar to a matrix.
-Search/create an operator that works both for scalar and matrix mltiplication