Fast algorithms for convolution quadrature of Riemann-Liouville fractional derivative
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2019
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866916533231419392 |
|---|---|
| author | Sun, Jing Nie, Daxin Deng, Weihua |
| author_facet | Sun, Jing Nie, Daxin Deng, Weihua |
| contents | Recently, the numerical schemes of the Fokker-Planck equations describing anomalous diffusion with two internal states have been proposed in [Nie, Sun and Deng, arXiv: 1811.04723], which use convolution quadrature to approximate the Riemann-Liouville fractional derivative; and the schemes need huge storage and computational cost because of the non-locality of fractional derivative and the large scale of the system. This paper first provides the fast algorithms for computing the Riemann-Liouville derivative based on convolution quadrature with the generating function given by the backward Euler and second-order backward difference methods; the algorithms don't require the assumption of the regularity of the solution in time, while the computation time and the total memory requirement are greatly reduced. Then we apply the fast algorithms to solve the homogeneous fractional Fokker-Planck equations with two internal states for nonsmooth data and get the first- and second-order accuracy in time. Lastly, numerical examples are presented to verify the convergence and the effectiveness of the fast algorithms. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_1901_05128 |
| institution | arXiv |
| publishDate | 2019 |
| record_format | arxiv |
| spellingShingle | Fast algorithms for convolution quadrature of Riemann-Liouville fractional derivative Sun, Jing Nie, Daxin Deng, Weihua Numerical Analysis 26A33, 44A35, 65M06 Recently, the numerical schemes of the Fokker-Planck equations describing anomalous diffusion with two internal states have been proposed in [Nie, Sun and Deng, arXiv: 1811.04723], which use convolution quadrature to approximate the Riemann-Liouville fractional derivative; and the schemes need huge storage and computational cost because of the non-locality of fractional derivative and the large scale of the system. This paper first provides the fast algorithms for computing the Riemann-Liouville derivative based on convolution quadrature with the generating function given by the backward Euler and second-order backward difference methods; the algorithms don't require the assumption of the regularity of the solution in time, while the computation time and the total memory requirement are greatly reduced. Then we apply the fast algorithms to solve the homogeneous fractional Fokker-Planck equations with two internal states for nonsmooth data and get the first- and second-order accuracy in time. Lastly, numerical examples are presented to verify the convergence and the effectiveness of the fast algorithms. |
| title | Fast algorithms for convolution quadrature of Riemann-Liouville fractional derivative |
| topic | Numerical Analysis 26A33, 44A35, 65M06 |
| url | https://arxiv.org/abs/1901.05128 |