“…Even under the lowering of regularity have enormous difficulties occurred in the existence of weak solutions. Among many different ap-proaches and attempts in this direction, e.g., the Γ-limit method in Bellettini & Fusco [3], the Young measure solutions in Chen & Zhang [4], and numerical scheme analyses in Esedoglu [9] and Esedoglu & Greer [10], to our best knowledge, Zhang [30] was the first to successfully prove that, for n = 1, there are infinitely many (Lipschitz) weak solutions to (1.1) for any given smooth initial data u 0 . His method uses the variational technique of differential inclusion together with the so called in-approximation method or convex integration; this new method can also deal with other ill-posed forward-backward diffusion problems (see, e.g., the pioneering work of Höllig [16] and its recent generalization by Zhang [31].)…”