“…Therefore, it is vital and necessary to design numerical methods to preserve discrete energy dissipation laws. Many efforts have been devoted to the development of numerical methods for energy stability in this active research field, which include, but are not limited to, the convex splitting method [4, 23, 28-30, 33, 45, 54, 56], the average vector field method [6,55], exponential time differencing (ETD) method [10,11,13,22,36,37,44,51] and the invariant energy quadratization (IEQ) method [34,[66][67][68]70]. In addition, the scalar auxiliary variable (SAV) method [1,14,15,27,57,58,62,69] has been successfully developed, inspired by a similar idea of the IEQ method.…”