Furthermore they, together with Kang, introduced a higher order variational model [9] based on the Euler’s elastica which connects the level lines by using Euler elastica curves [19] instead of using straight lines as the TV model does. Unfortunately for the latter two models there appear to exist no fast methods to find the numerical solution. The aim of this paper is to develop a fast multigrid algorithm for the CDD model.