The numerical solution of spectral fractional diffusion problems in the form ${\mathcal A}^\alpha u = f$ is studied, where $\mathcal A$ is a selfadjoint elliptic operator in a bounded domain $\Omega\subset {\mathbb R}^d$, and $\alpha \in (0,1]$. The finite difference approximation of the problem leads to the system ${\mathbb A}^\alpha {\mathbf u} = {\mathbf f}$, where ${\mathbb A}$ is a sparse, symmetric and positive definite (SPD) matrix, and ${\mathbb A}^\alpha$ is defined by its spectral decomposition. In the case of finite element approximation, ${\mathbb A}$ is SPD with respect to the dot product associated with the mass matrix. The BURA method is introduced by the best uniform rational approximation of degree $k$ of $t^{\alpha}$ in $[0,1]$, denoted by $r_{\alpha,k}$. Then the approximation ${\bf u}_k\approx {\bf u}$ has the form ${\bf u}_k = c_0 {\mathbf f} +\sum_{i=1}^k c_i({\mathbb A} - {\widetilde{d}}_i {\mathbb I})^{-1}{\mathbf f}$, ${\widetilde{d}}_i<0$, thus requiring the solving of $k$ auxiliary linear systems with sparse SPD matrices. The BURA method has almost optimal computational complexity, assuming that an optimal PCG iterative solution method is applied to the involved auxiliary linear systems. The presented analysis shows that the absolute values of first %${\widetilde{d}}_i$ $\left\{{\widetilde{d}}_i\right\}_{i=1}^{k'}$ can be extremely large. In such a case the condition number of ${\mathbb A} - {\widetilde{d}}_i {\mathbb I}$ is practically equal to one. Obviously, such systems do not need preconditioning. The next question is if we can replace their solution by directly multiplying ${\mathbf f}$ with $-c_i/{\widetilde{d}}_i$. Comparative analysis of numerical results is presented as a proof-of-concept for the proposed RS-BURA method.


翻译:=============================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================

0
下载
关闭预览

相关内容

Linux导论,Introduction to Linux,96页ppt
专知会员服务
77+阅读 · 2020年7月26日
Stabilizing Transformers for Reinforcement Learning
专知会员服务
58+阅读 · 2019年10月17日
学习自然语言处理路线图
专知会员服务
137+阅读 · 2019年9月24日
LibRec 精选:AutoML for Contextual Bandits
LibRec智能推荐
7+阅读 · 2019年9月19日
《自然》(20190221出版)一周论文导读
科学网
6+阅读 · 2019年2月23日
Unsupervised Learning via Meta-Learning
CreateAMind
42+阅读 · 2019年1月3日
Ray RLlib: Scalable 降龙十八掌
CreateAMind
9+阅读 · 2018年12月28日
disentangled-representation-papers
CreateAMind
26+阅读 · 2018年9月12日
数据分析师应该知道的16种回归技术:Lasso回归
数萃大数据
16+阅读 · 2018年8月13日
BAT机器学习面试题1000题(316~320题)
七月在线实验室
14+阅读 · 2018年1月18日
【推荐】决策树/随机森林深入解析
机器学习研究会
5+阅读 · 2017年9月21日
【学习】Hierarchical Softmax
机器学习研究会
4+阅读 · 2017年8月6日
Auto-Encoding GAN
CreateAMind
7+阅读 · 2017年8月4日
Arxiv
0+阅读 · 2021年7月8日
Arxiv
0+阅读 · 2021年7月8日
VIP会员
相关资讯
LibRec 精选:AutoML for Contextual Bandits
LibRec智能推荐
7+阅读 · 2019年9月19日
《自然》(20190221出版)一周论文导读
科学网
6+阅读 · 2019年2月23日
Unsupervised Learning via Meta-Learning
CreateAMind
42+阅读 · 2019年1月3日
Ray RLlib: Scalable 降龙十八掌
CreateAMind
9+阅读 · 2018年12月28日
disentangled-representation-papers
CreateAMind
26+阅读 · 2018年9月12日
数据分析师应该知道的16种回归技术:Lasso回归
数萃大数据
16+阅读 · 2018年8月13日
BAT机器学习面试题1000题(316~320题)
七月在线实验室
14+阅读 · 2018年1月18日
【推荐】决策树/随机森林深入解析
机器学习研究会
5+阅读 · 2017年9月21日
【学习】Hierarchical Softmax
机器学习研究会
4+阅读 · 2017年8月6日
Auto-Encoding GAN
CreateAMind
7+阅读 · 2017年8月4日
Top
微信扫码咨询专知VIP会员