In trying to find a smoother that would work with Algorithm 2, our idea is to address the problem described in §3.3 by overcoming the singularity associated with the FP method. Motivated by the Euler’s Elastica model [9], we decided to introduce two positive parameters a and b in the function g. Specifically our proposal for the diffusion coefficient G is