“…Moreover, the problem under consideration is a nonlinear time-fractional equation, and classical nonlocal numerical schemes blending with iterative methods would take more expensive computation and have more complexity in the analysis. So it is really competitive and also necessary to construct linearized numerical methods for nonlinear time-fractional problems [12,16,17,18,33]. It is worth to mention that the solution of many time-fractional differential equations typically displays a weak singularity near the initial time [8,9,10,19,20,21,26,29,30], which leads to the loss of time accuracy for many related high-order numerical methods.…”