基于张量乘积的快速谱元算法

摘要

针对椭圆型方程的谱元离散系统构造了一种基于张量乘积的快速直接解法。 分析显示, 新算法的计算量仅相当于迭代方法迭代Kx+Ky次的计算量(这里Kx,Ky分别为x,y方向的区域剖分数), 特别适合那些网格不多但多项式阶数较高的谱元离散.我们还将张量乘积方法推广到具有Neumann边界条件的奇异泊松问题的求解, 给出了具体的实现方法.最后, 利用张量乘积构造了变形区域上椭圆型方程的预条件子, 数值结果显示预条件系统的条件数与多项式阶数无关