×

张量积结构线性系统的Krylov子空间方法。 (英语) Zbl 1208.65044号

应用低张量秩近似,作者发展了求解具有张量积结构的线性系统的有效Krylov子空间方法。
例如,这种结构使用张量化有限元基将泊松方程的变分公式离散化为多维超立方体。
这些系统可以写成矩阵的(d)Kronecker积之和,例如,在(d=2)的情况下
\[(A_1\otimes I_{n_2}+I_{n_1}\otimesA_{n_2})x=b_1\otIMesb_2,\]
其中,(A_s\in\mathbb{R}^{n_s\times n_s},b_s\in\fathbb{R}^{ns},s=1,2,\)和(I_{ns}\)表示(ns\times n_s)恒等式矩阵。在重新计算后,给出了控制理论中遇到的著名Sylvester方程。虽然已经为该等式开发了Krylov子空间算法的几个变体,但需要对情况(d>2)进行改进,因为指数缩放的计算成本和内存需求对于更高的维度来说是不利的。高维的低张量秩近似已经由引入L.格拉瑟迪克[计算,72,第3-4期,247-265(2004;Zbl 1058.65036号)]但对于大规模矩阵,该算法中标量倍数的矩阵指数的计算变得昂贵。本文提出的方法仅依赖于矩阵-向量乘法与(A_s)的乘法,计算成本与(d)呈线性增长,而不是指数增长。
在得到张量算法、分解和低张量秩近似的一些结果之后,描述了新提出的带有Arnoldi分解的Krylov子空间方法的基本算法以及系数为上Hessenberg形式的压缩线性系统的证明,并给出了一些实现细节。分析了对称正定和非对称正定情形的收敛性。此外,还提出了一种扩展的张量Krylov子空间方法,该方法还允许矩阵-向量乘积与((A_s)^{-1})。
这些方法是使用MATLAB软件张量工具箱和进行了学术示例验证。应用标准的有限差分离散化方法,以具有由均匀分布的伪随机数组成的可分离右侧的一维泊松方程为对称示例,以对流扩散方程为非对称示例。

MSC公司:

65层10 线性系统的迭代数值方法
15A69号 多线性代数,张量演算
65号06 含偏微分方程边值问题的有限差分方法
65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
35J05型 拉普拉斯算子、亥姆霍兹方程(约化波动方程)、泊松方程
35年25日 二阶椭圆方程的边值问题
65楼30 其他矩阵算法(MSC2010)
15A24号 矩阵方程和恒等式
15甲16 矩阵的指数函数和相似函数
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接