pith. sign in

arxiv: 2606.25520 · v1 · pith:V2JMS46Unew · submitted 2026-06-24 · 🧮 math.NA · cs.NA

Efficient Krylov solvers for inverse source problem in 2D space-time fractional diffusion equation

Asim Ilyas , Stefano Serra-Capizzano This is my paper

Pith reviewed 2026-06-25 20:53 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords inverse source problemspace-time fractional diffusionGLT theorymultilevel Toeplitz matricespreconditionersGMRESquasi-boundary regularizationspectral analysis
0
0 comments X

The pith

GLT-derived preconditioners for multilevel Toeplitz systems from regularized 2D fractional inverse problems produce eigenvalue clustering around one.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper addresses the inverse source problem for a two-dimensional space-time fractional diffusion equation with variable coefficients by first applying quasi-boundary regularization to handle ill-posedness. Finite-difference discretization yields large linear systems possessing a multilevel Toeplitz-like block structure. Generalized Locally Toeplitz theory supplies a spectral analysis that guides the construction of structure-preserving preconditioners. These preconditioners induce eigenvalue clustering around one in the preconditioned matrices, which in turn accelerates GMRES convergence when recovering the source term, as confirmed by the reported numerical tests.

Core claim

The proposed preconditioners based on GLT spectral analysis preserve the multilevel structure of the discretization matrices arising from the quasi-boundary regularized inverse source problem and produce a general eigenvalue clustering around one for the preconditioned sequence, thereby accelerating the convergence of GMRES in source reconstruction.

What carries the argument

GLT-based preconditioners constructed from the spectral symbol of the multilevel Toeplitz-like discretization matrices of the regularized problem.

If this is right

  • GMRES convergence for source reconstruction becomes independent of mesh size once the preconditioner is applied.
  • The same preconditioning strategy remains effective when the diffusion coefficient varies in space.
  • The multilevel block structure is retained after preconditioning, allowing reuse of fast matrix-vector products.
  • The approach extends the applicability of GLT analysis from direct to inverse problems for fractional operators.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same clustering property may hold for other regularization strategies that produce similar Toeplitz-like patterns.
  • Higher-dimensional or nonlinear extensions could exploit the same GLT symbol to design cheap approximate inverses.
  • If the variable coefficient is smooth, the clustering radius may be bounded independently of the fractional orders.

Load-bearing premise

The GLT spectral analysis accurately captures the eigenvalue distribution of the preconditioned matrices arising from the variable-coefficient discretization of the quasi-boundary regularized problem.

What would settle it

A numerical test with variable coefficients in which the eigenvalues of the preconditioned matrix fail to cluster around one, or in which GMRES iteration counts do not decrease substantially compared with the unpreconditioned case, would falsify the claim.

read the original abstract

In this work, we consider a two-dimensional time-space fractional diffusion equation with a variable coefficient and investigate the inverse source problem of reconstructing the source term f(x,y) , after regularizing the problem using the quasi-boundary value method to mitigate ill-posedness. A finite difference discretization results in a large-scale linear system with a multilevel Toeplitz-like block structure. We perform a spectral analysis of the associated matrix sequences, employing tools from Generalized Locally Toeplitz (GLT) theory, and construct efficient preconditioners based on the GLT analysis. The proposed preconditioners preserve the multilevel structure of the discretization matrices and leads to a general eigenvalue clustering around one for the preconditioned sequence. Numerical experiments validate the theoretical findings and demonstrate that the proposed approach significantly accelerates the convergence of the GMRES method in reconstructing the source term in the two-dimensional space-time fractional diffusion equation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper addresses the inverse source problem for a 2D space-time fractional diffusion equation with variable coefficient by applying quasi-boundary regularization, finite-difference discretization to multilevel Toeplitz-like systems, GLT-based spectral analysis, and construction of structure-preserving preconditioners claimed to produce eigenvalue clustering at 1, with numerical experiments showing accelerated GMRES convergence for source reconstruction.

Significance. If the GLT analysis and clustering hold for variable coefficients, the work supplies a theoretically grounded preconditioning strategy that exploits the multilevel structure for efficient solution of large ill-posed fractional inverse problems; this is potentially useful in applications requiring repeated solves of such systems. The explicit use of GLT theory to guide preconditioner design is a methodological strength.

major comments (2)
  1. [Spectral analysis section (GLT symbol derivation)] Spectral analysis section: the claim that the preconditioned matrix sequence exhibits 'general eigenvalue clustering around one' for variable c(x,y) requires an explicit derivation showing that the GLT symbol of the preconditioner exactly cancels the pointwise factor c(x,y) (or its regularized analogue) so that the essential range collapses to the singleton {1}. If the symbol instead retains an interval whose length scales with the oscillation of c, the clustering guarantee does not hold for general variable coefficients and the central theoretical claim is weakened.
  2. [Numerical experiments section] Numerical experiments section: the reported GMRES iteration counts and eigenvalue plots must be accompanied by tests on coefficients c(x,y) with controlled oscillation amplitude (e.g., c ranging over [1,10] versus nearly constant); without such tests it is unclear whether the observed clustering and acceleration are artifacts of mild variation rather than a general property.
minor comments (2)
  1. [Abstract] Abstract: subject-verb agreement error ('leads' should be 'lead').
  2. [Discretization section] Notation for the multilevel block structure should be introduced with an explicit matrix pattern or small example before the GLT analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments identify two areas where the presentation of the GLT analysis and the supporting numerical evidence can be strengthened. We address each point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: Spectral analysis section: the claim that the preconditioned matrix sequence exhibits 'general eigenvalue clustering around one' for variable c(x,y) requires an explicit derivation showing that the GLT symbol of the preconditioner exactly cancels the pointwise factor c(x,y) (or its regularized analogue) so that the essential range collapses to the singleton {1}. If the symbol instead retains an interval whose length scales with the oscillation of c, the clustering guarantee does not hold for general variable coefficients and the central theoretical claim is weakened.

    Authors: We agree that an explicit cancellation step is needed for full rigor. In the current Section 3 the GLT symbol of the discretization operator is derived as a product of the multilevel Toeplitz symbol and the diagonal factor c(x,y). The preconditioner is defined via the reciprocal of that symbol (including the pointwise division by c). We will insert a new lemma that writes the symbol of the preconditioned operator explicitly and shows that, for any continuous positive c bounded away from zero, the essential range reduces to the singleton {1}. The revised proof will also address the regularized quasi-boundary term. revision: yes

  2. Referee: Numerical experiments section: the reported GMRES iteration counts and eigenvalue plots must be accompanied by tests on coefficients c(x,y) with controlled oscillation amplitude (e.g., c ranging over [1,10] versus nearly constant); without such tests it is unclear whether the observed clustering and acceleration are artifacts of mild variation rather than a general property.

    Authors: We accept the suggestion. The present experiments use a single moderately varying coefficient. In the revision we will add a new subsection (or table) that repeats the eigenvalue plots and GMRES counts for three families: nearly constant (range [0.95,1.05]), moderate ([1,5]), and strongly varying ([1,10] and [1,50]). All other parameters (grid size, fractional orders, regularization parameter) will be held fixed so that the dependence on oscillation amplitude is isolated. revision: yes

Circularity Check

0 steps flagged

No significant circularity; standard GLT application to discretization

full rationale

The derivation applies established GLT theory to the multilevel Toeplitz structure arising from finite-difference discretization of the regularized inverse source problem. The eigenvalue clustering claim follows directly from the GLT symbol of the preconditioned operator rather than any fitted parameter, self-definition, or reduction to the numerical test data. No load-bearing step reduces to a self-citation whose validity depends on the present paper, and the central spectral result is independent of the specific source reconstruction experiments.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The abstract introduces no free parameters, additional axioms, or invented entities beyond standard finite-difference discretization and GLT theory.

pith-pipeline@v0.9.1-grok · 5681 in / 1072 out tokens · 27753 ms · 2026-06-25T20:53:01.808148+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

51 extracted references · 36 canonical work pages

  1. [1]

    Aliu, E.E.M

    M., Abdulgaffar, J.N. Aliu, E.E.M. Micah, I. Taiwo, A.U. Abasido, and A.K. AdesugbaFractional Stochastic Calculus and Multifractional Processes: Applications in Financial Modeling, Polac. Int. J. Econs MGT Sci, 10 (2024): 79-86

    2024

  2. [2]

    Barbarino,A systematic approach to reduced GLT, BIT 62 (3) (2022), PP

    G. Barbarino,A systematic approach to reduced GLT, BIT 62 (3) (2022), PP. 681–743

    2022

  3. [3]

    Barbarino, C

    G. Barbarino, C. Garoni, and S. Serra-Capizzano,Block generalized locally Toeplitz sequences: theory and applications in the unidimensional case, Electron. Trans. Numer. Anal. 53 (2020), pp. 28–112. https://doi.org/10.1553/etna_vol53s28

    work page doi:10.1553/etna_vol53s28 2020

  4. [4]

    Barbarino, C

    G. Barbarino, C. Garoni, and S. Serra-Capizzano,Block generalized locally Toeplitz sequences: theory and applications in the multidimensional case, Electron. Trans. Numer. Anal. 53 (2020), pp. 113–216. https://doi.org/10.1553/etna_vol53s113

    work page doi:10.1553/etna_vol53s113 2020

  5. [5]

    Bini.Matrix structures in parallel matrix computations, Calcolo 25 (1-2) (1988), pp

    D. Bini.Matrix structures in parallel matrix computations, Calcolo 25 (1-2) (1988), pp. 37–51

    1988

  6. [6]

    W. Bu, Y. Tang, and J. Yang,Galerkin finite element method for two-dimensional Riesz space frac- tional diffusion equations, J. Comput. Phys. 276 (2014), pp. 26-38.https://doi.org/10.1016/j. jcp.2014.07.023

    work page doi:10.1016/j 2014

  7. [7]

    Chen and W

    M. Chen and W. Deng,Fourth order accurate scheme for the space fractional diffusion equations, SIAM J. Numer. Anal. 52 (3) (2014), pp. 1418-1438.https://doi.org/10.1137/130933447

    work page doi:10.1137/130933447 2014

  8. [8]

    Denche and K

    M. Denche and K. Bessila,A modified quasi-boundary value method for ill-posed problems, J. Math. Anal. Appl. 301 (2) (2005), pp. 419-426.https://doi.org/10.1016/j.jmaa.2004.08.001

    work page doi:10.1016/j.jmaa.2004.08.001 2005

  9. [9]

    Dhatt, E

    G. Dhatt, E. Lefrançois, and G. Touzot,Finite element method, John Wiley & Sons, 2012

    2012

  10. [10]

    Diaz-Fuentes, S

    R. Diaz-Fuentes, S. Serra-Capizzano, R. Sormani.Aω-circulant regularization for subdivision scheme linear systems arising in curve representation and generation, Algorithms 16 (2023), paper 328; online 8-7-23.https://doi.org/10.3390/a16070328

    work page doi:10.3390/a16070328 2023

  11. [11]

    Donatelli, M

    M. Donatelli, M. Mazza, and S. Serra-Capizzano,Spectral analysis and structure preserving pre- conditioners for fractional diffusion equations, J. Comput. Phys. 307 (2016), pp. 262-279.https: //doi.org/10.1016/j.jcp.2015.11.061

    work page doi:10.1016/j.jcp.2015.11.061 2016

  12. [12]

    Eymard, T

    R. Eymard, T. Gallouët, and R. Herbin,Finite volume methods, Handb. Numer. Anal., 7 (2000), pp. 713-1018

    2000

  13. [13]

    Ghafoor, M

    A. Ghafoor, M. Fiaz, M. Hussain, A. Ullah, E.A.A. Ismail, and F.A. Awwad,Dynamics of the time- fractional reaction–diffusion coupled equations in biological and chemical processes, Sci. Rep., 14(1) (2024): 7549.https://doi.org/10.1038/s41598-024-58073-z

    work page doi:10.1038/s41598-024-58073-z 2024

  14. [14]

    Garoni and S

    C. Garoni and S. Serra-Capizzano,Generalized locally Toeplitz sequences: theory and applications. Vol. I, Springer, Cham, 2017.https://doi.org/10.1007/978-3-319-53679-8

    work page doi:10.1007/978-3-319-53679-8 2017

  15. [15]

    Garoni and S

    C. Garoni and S. Serra-Capizzano,Generalized locally Toeplitz sequences: theory and applications. Vol. II, Springer, Cham, 2018.https://doi.org/10.1007/978-3-030-02233-4

    work page doi:10.1007/978-3-030-02233-4 2018

  16. [16]

    Garoni and S

    C. Garoni and S. Serra-Capizzano,Generalized locally Toeplitz sequences: a spectral analysis tool for discretized differential equations, Splines and PDEs: from approximation theory to numerical linear algebra, 161–236, Lecture Notes in Math., 2219, Fond. CIME/CIME Found. Subser., Springer, Cham, 2018.https://doi.org/10.1007/978-3-319-94911-6_3 22

    work page doi:10.1007/978-3-319-94911-6_3 2018

  17. [18]

    Golinskii and S

    L. Golinskii and S. Serra-Capizzano,The asymptotic properties of the spectrum of nonsymmetrically perturbed Jacobi matrix sequences, J. Approx. Theory 144 (1) (2007), pp. 84-102.https://doi.org/ 10.1016/j.jat.2006.05.002

    work page doi:10.1016/j.jat.2006.05.002 2007

  18. [19]

    Gomez-Aguilar and A

    J.F. Gomez-Aguilar and A. Atangana, eds.Applications of fractional calculus to modeling in dynamics and chaos. CRC Press, 2022

    2022

  19. [20]

    Gottlieb and S.A

    D. Gottlieb and S.A. Orszag,Numerical analysis of spectral methods: theory and applications, Society for Industrial and Applied Mathematics, Philadelphia, 1977

    1977

  20. [21]

    Z. Hao, Z. Zhang, and R. Du,Fractional centered difference scheme for high-dimensional integral fractional Laplacian, J. Comput. Phys. 424 (2021), 109851.https://doi.org/10.1016/j.jcp.2020. 109851

    work page doi:10.1016/j.jcp.2020 2021

  21. [22]

    Hào, N.V

    D.N. Hào, N.V. Duc, and H. Sahli,A non-local boundary value problem method for parabolic equations backward in time, J. Math. Anal. Appl. 345 (2) (2008), 805-815.https://doi.org/10.1016/j.jmaa. 2008.04.064

    work page doi:10.1016/j.jmaa 2008

  22. [23]

    S.Y. Hon, C. Li, R.L. Sormani, R. Krause, and S. Serra-Capizzano,Symbol-based multilevel blockτ preconditioners for multilevel block Toeplitz systems: GLT-based analysis and applications, SIAM J. Matrix Anal. Appl. 46 (4) (2025), 2331–2359.https://doi.org/10.1137/24M1702088

    work page doi:10.1137/24m1702088 2025

  23. [24]

    Ilyas, R.A

    A. Ilyas, R.A. Khalid, and S.A. Malik,Identifying temperature distribution and source term for gen- eralized diffusion equation with arbitrary memory kernel, Math. Methods Appl. Sci., 47 (7) (2024), 5894-5915.https://doi.org/10.1002/mma.9896

    work page doi:10.1002/mma.9896 2024

  24. [25]

    Ilyas and S.A

    A. Ilyas and S.A. Malik,Direct and some inverse problems for a generalized diffusion equation with variable coefficients, Comput. Appl. Math. 43 (6) (2024), 364.https://doi.org/10.1007/ s40314-024-02869-2

    2024

  25. [26]

    Ilyas and S

    A. Ilyas and S. Serra-Capizzano,Inverse source problems for identifying time and space-dependent coefficients in a 2D generalized diffusion equation, Appl. Math. Comput. 507 (2025), 129597.https: //doi.org/10.1016/j.amc.2025.129597

    work page doi:10.1016/j.amc.2025.129597 2025

  26. [27]

    Ilyas, M.F

    A. Ilyas, M.F. Khan, R.L. Sormani, G. Tento, and S. Serra-Capizzano.Determining the space de- pendent coefficients in space-time fractional diffusion equations via Krylov preconditioning, Linear Algebra Appl. (2026).https://doi.org/10.1016/j.laa.2026.02.024

    work page doi:10.1016/j.laa.2026.02.024 2026

  27. [28]

    Jin and W

    B. Jin and W. Rundell,A tutorial on inverse problems for anomalous diffusion processes, Inverse prob. 31 (3) (2015), 035003.https://doi.org/10.1088/0266-5611/31/3/035003

    work page doi:10.1088/0266-5611/31/3/035003 2015

  28. [29]

    J. Jia, H. Wang, and X. Zheng,A fast algorithm for time-fractional diffusion equation with space- time-dependent variable order, Numer. Algorithms, 94(4) (2023): 1705-1730.https://doi.org/10. 1007/s11075-023-01552-7

    2023

  29. [30]

    Kilbas, H.M

    A. Kilbas, H.M. Srivastava, and J.J. Trujillo,Theory and applications of fractional differential equa- tions, vol. 204, Elsevier, Amsterdam, 2006

    2006

  30. [31]

    Lei and H.-W

    S.-L. Lei and H.-W. Sun,A circulant preconditioner for fractional diffusion equations, J. Comput. Phys. 242 (2013), 715-725.https://doi.org/10.1016/j.jcp.2013.02.025. 23

    work page doi:10.1016/j.jcp.2013.02.025 2013

  31. [32]

    Li, X.-M

    M. Li, X.-M. Gu, C. Huang, M. Fei, and G. Zhang,A fast linearized conservative finite element method for the strongly coupled nonlinear fractional Schrödinger equations, J. Comput. Phys. 358 (2018), 256-282.https://doi.org/10.1016/j.jcp.2017.12.044

    work page doi:10.1016/j.jcp.2017.12.044 2018

  32. [33]

    Metzler, J.H

    R. Metzler, J.H. Jeon, A.G. Cherstvy, and E. Barkai,Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking, Phys. Chem. Chem. Phys., 16(44) (2014), 24128-24164.https://doi.org/10.1039/C4CP03465A

    work page doi:10.1039/c4cp03465a 2014

  33. [34]

    Owolabi, E

    K.M. Owolabi, E. Pindza, B. Karaagac, and G. Oguz,Laplace transform-homotopy perturbation method for fractional time diffusive predator–prey models in ecology, Partial Differ. Equ. Appl. Math. 9 (2024), 100607.https://doi.org/10.1016/j.padiff.2023.100607

    work page doi:10.1016/j.padiff.2023.100607 2024

  34. [35]

    H.-K.Pang, H.-H.Qin, H.-W.Sun, andT.-T.Ma,Circulant-based approximate inverse preconditioners for a class of fractional diffusion equations, Comput. Math. Appl. 85 (2021), 18-29.https://doi. org/10.1016/j.camwa.2021.01.007

    work page doi:10.1016/j.camwa.2021.01.007 2021

  35. [36]

    Rajyaguru, R

    A. Rajyaguru, R. Metzler, A.G. Cherstvy, and B. Berkowitz,Quantifying anomalous chemical diffu- sion through disordered porous rock materials, Phys. Chem. Chem. Phys., 27(17) (2025): 9056-9067. https://doi.org/10.1039/D5CP00654F

    work page doi:10.1039/d5cp00654f 2025

  36. [37]

    Saad and M.H

    Y. Saad and M.H. Schultz,GMRES: A generalized minimal residual algorithm for solving nonsym- metric linear systems, SIAM J. Sci. Stat. Comput. 7 (3) (1986), 856-869.https://doi.org/10.1137/ 0907058

    1986

  37. [38]

    Sandev, I.M

    T. Sandev, I.M. Sokolov, R. Metzler, and A. Chechkin,Beyond monofractional kinetics, Chaos Soli- tons Fract. 102 (2017), 210–217.https://doi.org/10.1016/j.chaos.2017.05.001

    work page doi:10.1016/j.chaos.2017.05.001 2017

  38. [39]

    Scalas, R

    E. Scalas, R. Gorenflo, and F. Mainardi,Fractional calculus and continuous-time finance, Physica A Stat. Mech. Appl. 284 (1-4) (2000), 376-384.https://doi.org/10.1016/S0378-4371(00)00255-7

    work page doi:10.1016/s0378-4371(00)00255-7 2000

  39. [40]

    Serra-Capizzano, S.:The rate of convergence of Toeplitz based PCG methods for second order non- linear boundary value problems. Numer. Math. 81 (3) (1999), 461–495

    1999

  40. [41]

    337 (2001), 37–78.https://doi.org/10.1016/S0024-3795(01)00335-4

    Serra-Capizzano, S.:Spectral behavior of matrix sequences and discretized boundary value problems, Linear Algebra Appl. 337 (2001), 37–78.https://doi.org/10.1016/S0024-3795(01)00335-4

    work page doi:10.1016/s0024-3795(01)00335-4 2001

  41. [42]

    Serra-Capizzano,Matrix algebra preconditioners for multilevel Toeplitz matrices are not super- linear, Linear Algebra Appl

    S. Serra-Capizzano,Matrix algebra preconditioners for multilevel Toeplitz matrices are not super- linear, Linear Algebra Appl. 343–344 (2002), 303-319.https://doi.org/10.1016/S0024-3795(01) 00361-5

    work page doi:10.1016/s0024-3795(01 2002

  42. [43]

    S. Serra-Capizzano,Generalized locally Toeplitz sequences: spectral analysis and applications to dis- cretized partial differential equations, Linear Algebra and its Applications, 366 (2003), pp. 371–402. https://doi.org/10.1016/S0024-3795(02)00504-9

    work page doi:10.1016/s0024-3795(02)00504-9 2003

  43. [44]

    Serra-Capizzano,The GLT class as a generalized Fourier analysis and applications, Linear Algebra and its Applications, 419 (2006), pp

    S. Serra-Capizzano,The GLT class as a generalized Fourier analysis and applications, Linear Algebra and its Applications, 419 (2006), pp. 180–233.https://doi.org/10.1016/j.laa.2006.04.012

    work page doi:10.1016/j.laa.2006.04.012 2006

  44. [45]

    Simmons, Q

    A. Simmons, Q. Yang, and T. Moroney,A finite volume method for two-sided fractional diffusion equations on non-uniform meshes, J. Comput. Phys. 335 (2017), 747-759.https://doi.org/10. 1016/j.jcp.2017.01.061

    2017

  45. [46]

    Suhaib, A

    K. Suhaib, A. Ilyas, and S.A. Malik,On the inverse problems for a family of integro-differential equations, Math. Model. Anal. 28 (2) (2023), 255-270.https://doi.org/10.3846/mma.2023.16139

    work page doi:10.3846/mma.2023.16139 2023

  46. [47]

    Sun and G

    Z. Sun and G. Gao,Fractional Differential Equations: Finite Difference Methods, De Gruyter, Berlin, 2020. 24

    2020

  47. [48]

    Thomas,Numerical Partial Differential Equations: Finite Difference Methods, vol

    J.W. Thomas,Numerical Partial Differential Equations: Finite Difference Methods, vol. 22, Springer, Berlin, 2013

    2013

  48. [49]

    Tuan, L.D

    N.H. Tuan, L.D. Long, V.T. Nguyen, and T. Tran,On a final value problem for the time-fractional diffusion equation with inhomogeneous source, Inverse Probl. Sci. Eng. 25 (9) (2017), 1367-1395. https://doi.org/10.1080/17415977.2016.1259316

    work page doi:10.1080/17415977.2016.1259316 2017

  49. [50]

    Z. Wang, X. Huang, and G. Shi,Analysis of nonlinear dynamics and chaos in a fractional order financial system with time delay, Comput. Math. Appl. 62 (3) (2011), 1531-1539.https://doi.org/ 10.1016/j.camwa.2011.04.057

    work page doi:10.1016/j.camwa.2011.04.057 2011

  50. [51]

    L. Zhao, C. Li, and F. Zhao,Efficient difference schemes for the Caputo-tempered fractional diffusion equations based on polynomial interpolation, Comm. Appl. Math. Comput. 3, no. 1 (2021), 1-40. ttps://doi.org/10.1007/s42967-020-00067-5

    work page doi:10.1007/s42967-020-00067-5 2021

  51. [52]

    Zhou, L.Y

    F.X. Zhou, L.Y. Wang, Z.Y. Liu, and W.C. Zhao.A viscoelastic-viscoplastic mechanical model of time-dependent materials based on variable-order fractional derivative, Mech. Time-Depend. Mater. 26, (3) (2022), 699-717.https://doi.org/10.1007/s11043-021-09508-x. 25

    work page doi:10.1007/s11043-021-09508-x 2022