Improving Greedy Algorithms for Rational Approximation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Adler, James H., Hu, Xiaozhe, Wang, Xue, Xue, Zhongqin
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866929430713073664
author Adler, James H.
Hu, Xiaozhe
Wang, Xue
Xue, Zhongqin
author_facet Adler, James H.
Hu, Xiaozhe
Wang, Xue
Xue, Zhongqin
contents When developing robust preconditioners for multiphysics problems, fractional functions of the Laplace operator often arise and need to be inverted. Rational approximation in the uniform norm can be used to convert inverting those fractional operators into inverting a series of shifted Laplace operators. Care must be taken in the approximation so that the shifted Laplace operators remain symmetric positive definite, making them better conditioned. In this work, we study two greedy algorithms for finding rational approximations to such fractional operators. The first algorithm improves the orthogonal greedy algorithm discussed in [Li et al., SISC, 2024] by adding one minimization step in the uniform norm to the procedure. The second approach employs the weak Chebyshev greedy algorithm in the uniform norm. Both methods yield non-increasing error. Numerical results confirm the effectiveness of our proposed algorithms, which are also flexible and applicable to other approximation problems. Moreover, with effective rational approximations to the fractional operator, the resulting algorithms show good performance in preconditioning a Darcy-Stokes coupled problem.
format Preprint
id arxiv_https___arxiv_org_abs_2407_15775
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Improving Greedy Algorithms for Rational Approximation
Adler, James H.
Hu, Xiaozhe
Wang, Xue
Xue, Zhongqin
Numerical Analysis
41A20, 65K10, 65N30, 65F08
When developing robust preconditioners for multiphysics problems, fractional functions of the Laplace operator often arise and need to be inverted. Rational approximation in the uniform norm can be used to convert inverting those fractional operators into inverting a series of shifted Laplace operators. Care must be taken in the approximation so that the shifted Laplace operators remain symmetric positive definite, making them better conditioned. In this work, we study two greedy algorithms for finding rational approximations to such fractional operators. The first algorithm improves the orthogonal greedy algorithm discussed in [Li et al., SISC, 2024] by adding one minimization step in the uniform norm to the procedure. The second approach employs the weak Chebyshev greedy algorithm in the uniform norm. Both methods yield non-increasing error. Numerical results confirm the effectiveness of our proposed algorithms, which are also flexible and applicable to other approximation problems. Moreover, with effective rational approximations to the fractional operator, the resulting algorithms show good performance in preconditioning a Darcy-Stokes coupled problem.
title Improving Greedy Algorithms for Rational Approximation
topic Numerical Analysis
41A20, 65K10, 65N30, 65F08
url https://arxiv.org/abs/2407.15775