The Scheduled Relaxation Jacobi (SRJ) method is a viable candidate as a high performance linear solver for elliptic partial differential equations (PDEs). The method greatly improves the convergence of the standard Jacobi iteration by applying a sequence of $M$ well-chosen overrelaxation and underrelaxation factors in each cycle of the algorithm to effectively attenuate the solution error. In previous work, optimal SRJ schemes (sets of relaxation factors) have been derived to accelerate convergence for specific discretizations of elliptic PDEs. In this work, we develop a family of SRJ schemes which can be applied to solve elliptic PDEs regardless of the specific discretization employed. To achieve favorable convergence, we train an algorithm to select which scheme in this family to apply at each cycle of the linear solve process, based on convergence data collected from applying these schemes to the one-dimensional Poisson equation. The automatic selection heuristic that is developed based on this limited data is found to provide good convergence for a wide range of problems.
翻译:Jacobi(SRJ)是作为高性能线性线性求解器的可行选择方程式,用于对椭圆部分差异方程式(PDEs)进行高性能分解。该方法通过在算法的每个周期中应用一个序列,即$Mo$的优选过度放宽和低放宽系数,以有效减轻溶解错误。在以往的工作中,为加速对椭圆部分方程式具体离散的趋同,制定了最佳SRJ(放松因子集)办法。在这项工作中,我们开发了一套可适用于解决离散性部分方程式的SRJ方案,无论采用何种特定的离散化,都可用于解决离异性方程式的趋同。为了实现有利的趋同,我们培训了一种算法,根据从将这些计划应用于单维Poisson方程式所收集的趋同数据,在线性解法进程的每一周期中选择适用哪种办法。根据这一有限数据开发的自动选择超常法,可以为广泛的问题提供良好的趋同。