Algorithm 886: Padua2D-Lagrange Interpolation at Padua Points on Bivariate Domains

被引：11
作者：
Caliari, Marco [1 ]
De Marchi, Stefano [1 ]
Vianello, Marco [2 ]
机构：
[1] Univ Verona, Dept Comp Sci, I-37134 Verona, Italy
[2] Univ Padua, Dept Pure & Appl Math, I-35121 Padua, Italy
来源：
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE | 2008年 / 35卷 / 03期
关键词：
Algorithms; Bivariate Lagrange interpolation; Padua points; bivariate Chebyshev orthogonal basis; Fortran; 77;
D O I：
10.1145/1391989.1391994
中图分类号：
TP31 [计算机软件];
学科分类号：
081202 ; 0835 ;
摘要：
We present a stable and efficient Fortran implementation of polynomial interpolation at the Padua points on the square [-1, 1](2). These points are unisolvent and their Lebesgue constant has minimal order of growth (log square of the degree). The algorithm is based on the representation of the Lagrange interpolation formula in a suitable orthogonal basis, and takes advantage of a new matrix formulation together with the machine-specific optimized BLAS subroutine DGEMM for the matrix-matrix product. Extension to interpolation on rectangles, triangles and ellipses is also described.
引用
收藏
页数：11
相关论文
共 18 条
[1]  
*AMD, 2006, AMD COR MATH LIB ACM
[2]   Multivariate simultaneous approximation [J].
Bagby, T ;
Bos, L ;
Levenberg, N .
CONSTRUCTIVE APPROXIMATION, 2002, 18 (04) :569-577
[3]   An updated set of Basic Linear Algebra Subprograms (BLAS) [J].
Blackford, LS ;
Demmel, J ;
Dongarra, J ;
Duff, I ;
Hammarling, S ;
Henry, G ;
Heroux, M ;
Kaufman, L ;
Lumsdaine, A ;
Petitet, A ;
Pozo, R ;
Remington, K ;
Whaley, RC .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2002, 28 (02) :135-151
[4]   On polynomial interpolation of two variables [J].
Bojanov, B ;
Xu, Y .
JOURNAL OF APPROXIMATION THEORY, 2003, 120 (02) :267-282
[5]   Bivariate Lagrange interpolation at the Padua points: the ideal theory approach [J].
Bos, Len ;
De Marchi, Stefano ;
Vianello, Marco ;
Xu, Yuan .
NUMERISCHE MATHEMATIK, 2007, 108 (01) :43-57
[6]  
Bos L, 2006, ELECTRON T NUMER ANA, V25, P1
[7]   Bivariate Lagrange interpolation at the Padua points: The generating curve approach [J].
Bos, Len ;
Caliari, Marco ;
De Marchi, Stefano ;
Vianello, Marco ;
Xu, Yuan .
JOURNAL OF APPROXIMATION THEORY, 2006, 143 (01) :15-25
[8]   Bivariate polynomial interpolation on the square at new nodal sets [J].
Caliari, M ;
De Marchi, S ;
Vianello, M .
APPLIED MATHEMATICS AND COMPUTATION, 2005, 165 (02) :261-274
[9]  
CALIARI M, 2007, J COMPUT APPL M 1023
[10]   Hyper2d: A numerical code for hyperinterpolation on rectangles [J].
Caliari, Marco ;
Vianello, Marco ;
De Marchi, Stefano ;
Montagna, Roberto .
APPLIED MATHEMATICS AND COMPUTATION, 2006, 183 (02) :1138-1147
12