“…Its main idea is to select a coarse-grid space to produce a rough approximation for the solution of nonlinear problems, and then use it as the initial guess on the fine-grid space. Later, two-grid methods combined with all kinds of spatial discretization methods are successfully extended to semilinear elliptic interface problems [25,7], semilinear elliptic optimal control problems [27,15], eigenvalue problems [31,32,29], and further developed for other applications, such as nonlinear Sobolev equations [30], Cahn-Hilliard equations [16], Stokes-Darcy problems [14,20], Darcy-Brinkman fracture models [6], nonlinear time-fractional parabolic equations [33], semilinear parabolic integrodifferential equations [17,4,26], nonlinear hyperbolic integro-differential equation [21], time-dependent Schrödinger equation [3,24,23], etc.…”