“…Numerical methods for time-fractional IBVPs with constant or time-independent diffusion parameter have received a huge amount of attention over the last decade. For such problems, several numerical methods have been proposed and analyzed, such as finite difference method [7, 19-21, 27, 36, 38], finite element method [6,32,39,41,43,44,48], discontinuous Galerkin (DG) methods [3, 4, 9-11, 31, 34], spectral method [23], and finite volume method [15,46], etc. The time-fractional IBVPs (1a)-(1b) with time-space dependent diffusivity is indeed very interesting and also practically important, and the numerical solutions of this problems were considered by a few authors only.…”