“…Due to the nonlocality of the fractional operator, the calculation of the L1 method and its variants is costly compared with the integer-order derivative operators. Generally speaking, we need O(N ) storage and O(N 2 ) arithmetic operations for calculating the L1 method δ α t u n , 1 ≤ n ≤ N (see [11,9,22]), while for the integer-order derivative, i.e., α = 1, the L1 method reduces to the firstorder backward formula, which needs O(1) storage and O(N ) arithmetic operations. The fast L1 method has been extensively studied in the literature (see, e.g., [11]).…”