Next Article in Journal
The Weitzenböck Type Curvature Operator for Singular Distributions
Next Article in Special Issue
On an Energy-Dependent Quantum System with Solutions in Terms of a Class of Hypergeometric Para-Orthogonal Polynomials on the Unit Circle
Previous Article in Journal
Collaborative Learning Based on Harry Potter for Learning Geometric Figures in the Subject of Mathematics
Previous Article in Special Issue
On Differential Equations Associated with Perturbations of Orthogonal Polynomials on the Unit Circle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

On Row Sequences of Hermite–Padé Approximation and Its Generalizations

by
Nattapong Bosuwan
1,2
1
Department of Mathematics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
2
Centre of Excellence in Mathematics, Commission on Higher Education, Ministry of Education, Si Ayutthaya Road, Bangkok 10400, Thailand
Mathematics 2020, 8(3), 366; https://doi.org/10.3390/math8030366
Submission received: 20 January 2020 / Revised: 1 March 2020 / Accepted: 3 March 2020 / Published: 6 March 2020

Abstract

:
Hermite–Padé approximation has been a mainstay of approximation theory since the concept was introduced by Charles Hermite in his proof of the transcendence of e in 1873. This subject occupies a large place in the literature and it has applications in different subjects. Most of the studies of Hermite–Padé approximation have mainly concentrated on diagonal sequences. Recently, there were some significant contributions in the direction of row sequences of Type II Hermite–Padé approximation. Moreover, various generalizations of Type II Hermite–Padé approximation were introduced and studied on row sequences. The purpose of this paper is to reflect the current state of the study of Type II Hermite–Padé approximation and its generalizations on row sequences. In particular, we focus on the relationship between the convergence of zeros of the common denominators of such approximants and singularities of the vector of approximated functions. Some conjectures concerning these studies are posed.

1. Introduction

Charles Hermite [1] was the first who introduced the idea of Hermite–Padé approximation. Particularly, he used it for systems of exponential functions to prove that e is transcendental. A formal study of Hermite–Padé approximation for general systems of functions was initiated by Mahler [2] (see also the papers by his students, Coates and Jager [3,4], for important results in this regard). There are basically two different types of Hermite–Padé approximation, namely Types I and II. However, the combination of both (called mixed type) is possible. In this paper, we primarily consider Type II.
Hermite–Padé approximation and its relatives have applications in various areas, for example, in number theory (see [1,5,6,7,8,9]), numerical analysis (see [10,11,12,13,14,15,16,17,18]), multiple orthogonal polynomials (see [18,19,20,21]), linear algebraic equations (see [22]), nonlinear dynamical systems (see [23]), Brownian motion (see [24]), in random matrices (see [19,25]), Gibbs phenomenon (see [26]), and Lie algebra solution of differential equations (see [27]). In addition to the proof of the transcendence of e , Hermite–Padé approximation was used in various irrationality and transcendence proofs of important numbers (see, e.g., [1,5,6,7,8]). Moreover, one can say that the theory of multiple orthogonal polynomials originated from Hermite–Padé approximation (see, e.g., Section 2.2 of [28] for the explanation).
An optimal choice of the coefficients of the denominators and numerators of Hermite–Padé approximants makes it an important tool to study analytic continuation of functions and localization of their singularities. We can make this point more clearly by considering classical Padé approximation (scalar Hermite–Padé approximation) stated as follows.
Given a formal Taylor series at the origin
F ( z ) = k = 0 f k z k ,
for any integers n , m 0 , we can find polynomials P P n and polynomials Q P m , Q 0 , such that
( Q F P ) ( z ) = O ( z n + m + 1 ) , a s z 0
( P n is the set of all polynomials of degree at most n). The rational function
R n , m : = P Q = P n , m Q n , m
is uniquely defined and is called the ( n , m ) classical Padé approximant of F. The polynomials P n , m and Q n , m in Equation (2) are selected so that Q n , m is monic and g c d ( P n , m , Q n , m ) = 1 . For a function F as in Equation (1), we denote by R 0 ( F ) the radius of the largest open disk at the origin to which F can be extended analytically and by R m ( F ) the radius of the largest open disk at the origin to which F can be extended so that F has at most m poles counting multiplicities. Set
B R : = { z C : | z | < R } .
By Q m ( F ) , we denote the monic polynomial whose zeros are the poles of F in B R m ( F ) counting multiplicities. The set of all distinct zeros of Q m ( F ) is denoted by P m ( F ) .
Analytic continuation and locations of poles including their multiplicities of F ( z ) = k = 0 f k z k in B R m ( F ) when F has exactly m poles in B R m ( F ) can be completely described by the convergences of { R n , m } n 0 and { Q n , m } n 0 in the following theorem (see [29,30]).
Theorem 1
(Montessus de Ballore–Gonchar’s Theorem). Let F be defined as in Equation (1) and fix m N . Then, the following statements are equivalent:
(a) 
R 0 ( F ) > 0 and F has exactly m poles in B R m ( F ) counting multiplicities.
(b) 
There is a polynomial Q m of degree m , Q m ( 0 ) 0 , such that the sequence of { Q n , m } n 0 satisfies
lim sup n Q m Q n , m 1 / n = θ < 1 ,
where · denotes the coefficient norm in the space of all polynomials.
Moreover, if either (a) or (b) holds, then Q m = Q m ( F ) ,
θ = max { | λ | : λ P m ( F ) } R m ( F ) ,
and
lim sup n F R n , m K 1 / n = z K R m ( F ) ,
where K is any compact subset of B R m ( F ) \ P m ( F ) and · K denotes the sup-norm on K .
In [29], Montessus de Ballore proved that (a) implies (b) with Q m = Q m ( F ) and the inequalities with the sign “≤” instead of “=” in Equations (3) and (4). This part of the above theorem is commonly known as Montessus de Ballore’s theorem. The implication (b)⇒(a) with the inequalities with the sign “≥” instead of “=” in Equations (3) and (4) was proved by Gonchar in [30]. In this current paper, we refer to this part as the inverse statement of Montessus de Ballore–Gonchar’s Theorem. We also note that the sequence { R n , m ( z ) } n 0 diverges at every point z , | z | > R m ( F ) (see Section 7 of [31]).
Later, Gonchar [31] studied the rate of attraction of each individual pole of F in B R m ( F ) to zeros of Q n , m . He introduced several indicators describing the asymptotic behavior of the zeros of Q n , m to a point λ C \ { 0 } . Set
| z w | 1 : = min { 1 , | z w | } , z , w C .
Fix m N and a function F defined as in Equation (1). Let
P n , m : = { λ n , 1 , λ n , 2 , , λ n , m n } , m n m , n 0 ,
denote the collection of zeros of Q n , m repeated according to their multiplicity and enumerated in nondecreasing distance to the point λ . The first indicator is defined by
Δ ( λ ) : = lim sup n j = 1 m n | λ n , j λ | 1 1 / n = lim sup n | λ n , j λ | < 1 | λ n , j λ | 1 / n .
It is easy to check that 0 Δ ( λ ) 1 under the convention when m n = 0 or | λ n , j λ | 1 for all j = 1 , 2 , , m n , the product is taken to be 1. The second and third indicators, nonnegative integers σ ( λ ) and γ ( λ ) , are defined as follows. We say that γ ( λ ) : = ν if
lim n | λ n , ν λ | = 0 a n d lim sup n | λ n , ν + 1 λ | > 0
(for ν > m n by convention | λ n , ν λ | : = 1 , and when lim sup n | λ n , 1 λ | > 0 , we take γ ( λ ) = 0 ). Similarly, σ ( λ ) : = ν if
lim sup n | λ n , ν λ | 1 / n < 1 a n d lim sup n | λ n , ν + 1 λ | 1 / n 1 .
Moreover, we define
δ j ( λ ) : = lim sup n | λ n , j λ | 1 1 / n .
Clearly, γ ( λ ) σ ( λ ) , and the statements Δ ( λ ) < 1 and σ ( λ ) 1 are equivalent.
The following theorem (see [31], Theorem 1) asserts that a pole of F of order ν in B R m ( F ) attracts with geometric rate exactly ν zeros of the polynomials Q n , m .
Theorem 2
(Gonchar’s theorem). Let F be defined as in Equation (1), m N , and let 0 λ C . The following two assertions are equivalent:
(a) 
λ B R m ( F ) and F has a pole at λ .
(b) 
Δ ( λ ) < 1 .
If either (a) or (b) takes place, then
Δ ( λ ) = | λ | R m ( F ) a n d σ ( λ ) = ν ,
where ν is the order of the pole at λ . Moreover, if we assume further that lim inf n | λ λ n , ν + 1 | > 0 , then
δ 1 ( λ ) = δ 2 ( λ ) = = δ ν ( λ ) = | λ | R m ( F ) 1 / ν .
Naively, Gonchar asked what we can say about λ if λ attracts a certain numbers of zeros of the polynomials Q n , m without such geometric rate in the above theorem. In [31], Gonchar also proposed the following conjecture.
Theorem 3
(Gonchar’s conjecture). Fix m N . Let F be defined as in Equation (1). Assume that γ ( λ ) 1 and λ 0 . Then, this series defines a function which is holomorphic at z = 0 , R m 1 ( F ) | λ | , and λ is a singularity of F .
Gonchar’s conjecture remains open. Special cases of the conjecture were proved by Vavilov, López, Prokhorov, and Suetin (see [32,33,34,35]). In the final progress, the following weaker version of Gonchar’s conjecture was proved by Suetin [35].
Theorem 4
(Suetin’s theorem). Assume that the formal power series F ( z ) = k = 0 f k z k has coefficients such that for fixed m N and sufficiently large n N the approximants R n , m have precisely m finite poles λ n , 1 , , λ n , m , which are convergent:
lim n λ n , j = λ j 0 , j = 1 , , m .
Then,
(i) 
The power series defines a holomorphic function F in the disk B R min , where R min : = min 1 j m | λ j | .
(ii) 
R m 1 ( F ) = max 1 j m | λ j | .
(iii) 
All points λ 1 , , λ m are singularities of F , the ones lying in the disk B R m 1 ( F ) are poles, and F has no other poles in this disk.
It is easy to check that if the Taylor coefficients of F ( z ) = k = 0 f k z k satisfy f k 0 and f k + 1 0 , then f k / f k + 1 is a zero of Q n , 1 . Therefore, when m = 1 , Gonchar’s conjecture and Suetin’s theorem reduce to the following Fabry ratio theorem (see [36]).
Theorem 5
(Fabry ratio theorem). Suppose that the coefficients of a power series F ( z ) = k = 0 f k z k are such that the limit
lim k f k f k + 1 = λ 0
exists. Then, the series converges uniformly on each compact subset of the disk B | λ | and λ is a singular point of the function F .
Recently, extensions of Montessus de Ballore–Gonchar’s theorem to Type II Hermite–Padé approximation and its generalizations were proved. As time progresses, there were many results on Type II Hermite–Padé approximation in the direction of Gonchar’s conjecture. The purpose of this survey paper is to review all those results and collect open problems in this respect.

2. Hermite–Padé Approximation

2.1. Definition and Notation

Type II Hermite–Padé approximation involves the approximation of a vector of functions by a vector of rational functions with the same denominator. Let F = ( F 1 , F 2 , , F d ) be a system of d formal Taylor expansions at the origin; that is, for each = 1 , 2 , , d , we have
F ( z ) = k = 0 f k , z k , f k , C .
In what follows, N : = { 1 , 2 , 3 , } , P n is the set of all polynomials of degree at most n , and for a given multi-index m = ( m 1 , m 2 , , m d ) N d , we define
| m | : = m 1 + m 2 + + m d .
Definition 1.
Let F = ( F 1 , F 2 , , F d ) be a system of d formal Taylor expansions as in Equation (5) and m = ( m 1 , m 2 , , m d ) N d be a fixed multi-index. Then, for each n max { m 1 , m 2 , , m d } , there exist Q P | m | and P P n m for all = 1 , 2 , , d such that Q 0 and
Q ( z ) F ( z ) P ( z ) = O ( z n + 1 ) , as   z 0 .
The vector of rational functions R n , m = ( P 1 / Q , P 2 / Q , , P d / Q ) is called an ( n , m ) Type II Hermite–Padé (HP) approximant of F .
Traditionally, the numbers of interpolation conditions of F , for = 1 , 2 , , d , in Equation (6) at 0 are selected to be the same, which is n + 1 . Since m = 0 does not provide any advantage in locating singularities of F , we can restrict our multi-indices m N d as stated in Definition 1. For any fixed ( n , m ) N × N d , in general, R n , m may not be unique and we assume that, given ( n , m ) , one particular solution is assigned. We write
R n , m = ( R n , m , 1 , R n , m , 2 , , R n , m , d ) = ( P n , m , 1 , P n , m , 2 , , P n , m , d ) / Q n , m ,
where P n , m , , = 1 , 2 , , d , and Q n , m are chosen so that Q n , m is a monic polynomial that has no common zero with all the P n , m , . When m remains fixed, we call the sequences { R n , m } n | m | m row sequences. When n and m have the relation,
m 1 = = m d = m , n = ( d + 1 ) m , m N
(or nearby configurations of multi-indices), the sequences { R n , m } are called diagonal sequences. Another construction called Type I HP approximants has very close algebraic relation to Type II HP approximants. For Type I HP approximation, one wants to approximate polynomial combination of the vector ( F 1 , F 2 , , F d ) by a polynomial. However, in Type II HP approximation, one wants to approximate F separately by rational functions with the same denominator. For interested readers, we refer to Chapter 4 of [37] for more details about Type I HP approximation. Since we only provide the survey of the studies of Type II HP approximants, we omit the word “Type II” when we refer to Type II HP approximants. Moreover, we would like to emphasize that, for d = 1 , the HP approximant R n , m reduces to the ( n m , m ) classical Padé approximant. However, because we are interested in the cases when m is fixed and n , all theorems in Section 1 hold for the scalar HP approximation.

2.2. Results and Conjectures

In the direction of row sequences, the paper [38] by Graves-Morris and Saff was a pioneering result in the sense that it initiated a convergence theory for HP approximants. In [38], they proved a Montessus de Ballore type theorem for HP approximants under the concept of polewise independence. Later, Cacoq, de la Calle, and López [39] improved the results [38] in several ways; namely, improving the estimates on the rates of convergences of { R n , m } n | m | and { Q n , m } n | m | and weakening the assumption of polewise independence. Note that in the Montessus de Ballore type theorem for HP approximants in [39], they found the exact rate of convergence of { R n , m } n | m | but not of { Q n , m } n | m | . A significant contribution in the direction of row sequences is due to Cacoq, de la Calle, and López [40] where they gave necessary and sufficient conditions for the convergence with geometric rate of { Q n , m } n | m | and calculated the exact rate of convergence of { Q n , m } n | m | . To explain the results in [40], we need to state some definitions.
Definition 2.
Let Ω : = ( Ω 1 , Ω 2 , , Ω d ) be a system of domains such that, for each = 1 , 2 , , d , F is meromorphic in Ω . We say that λ C is a pole of F in Ω of order τ if there exists an index { 1 , 2 , , d } such that λ Ω is a pole of F of order τ , and for β either λ is a pole of F β of order less than or equal to τ or λ Ω β . When Ω = ( Ω , Ω , , Ω ) , we say that λ is a pole of F in Ω .
Let R 0 ( F ) be the radius of the largest disk B R 0 ( F ) to which all the expansions F , = 1 , 2 , , d can be extended analytically. If R 0 ( F ) = 0 , we take B R m ( F ) = , m 0 ; otherwise, R m ( F ) is the radius of the largest disk B R m ( F ) centered at the origin to which all the analytic elements ( F , B R 0 ( F ) ) can be extended so that F has at most m poles counting multiplicities.
To prove an analog of Montessus de Ballore–Gonchar’s theorem for HP approximants, we have to decide what actually is the limit of { Q n , m } n | m | . We found that the convergence of { Q n , m } n | m | for the vector case is more complicated than the one for the scalar case, as the following example shows.
Example 1.
Let F = ( F 1 , F 2 ) , where
F 1 ( z ) = 1 z 2 + k = 0 z 3 2 k + 1 z 4 , F 2 ( z ) = 1 z 2 + k = 0 z 3 2 k
and fix m = ( 1 , 1 ) . Define F ^ : = ( F ^ 1 , F ^ 2 ) , where
F ^ 1 ( z ) = F 1 ( z ) F 2 ( z ) = 1 z 4 , F ^ 2 ( z ) = F 2 ( z ) = 1 z 2 + k = 0 z 3 2 k
It is not difficult to check that HP approximants of ( F , m ) and ( F ^ , m ) have the same Q n , m . Applying Theorem 4.4 in [39] for F ^ , 2 and 4 which are poles of F attract zeros of Q n , m with geometric rate as n . To be precise, we obtain
lim sup n Q n , m ( z 2 ) ( z 4 ) 1 / n < 1 .
However, since both F 1 and F 2 are meromorphic up to B 3 , 4 is not a pole of F .
Based on the above example, the idea of pole in Definition 2 is not suitable when we study a system of functions. The authors of [40] proposed a new definition of pole (called “system pole”) below.
Definition 3.
Given F = ( F 1 , F 2 , , F d ) in as Equation (5) and m = ( m 1 , m 2 , , m d ) N d , we say that λ C \ { 0 } is a system pole of order τ of F with respect to m if τ is the largest positive integer such that for each t = 1 , 2 , , τ , there exists at least one polynomial combination of the form
= 1 d v F , deg v < m , = 1 , 2 , , d ,
which is holomorphic on a neighborhood of B | λ | ¯ except for a pole at z = λ of exact order t .
Notice that in Definition 3, instead of examining poles of each function F separately, we examine poles of the polynomial combinations of the vector ( F 1 , F 2 , , F d ) . When d = 1 , the statement F has a system pole of order τ with respect to m = m simply means that F has a pole of order τ in B R m ( F ) . From Example 1, 2 and 4 are system poles of order 1 of F with respect to m but 4 is not a pole of F . Conversely, the following example shows that a pole of F may not be a system pole.
Example 2.
Let F = ( F 1 , F 2 ) , where
F 1 ( z ) = 1 z 1 + 1 z 2 , F 2 ( z ) = 1 z 3
and fix m = ( 1 , 1 ) . Clearly, 2 is a pole of F . Since 1 and 3 are only system poles of F , 2 is not a system pole of F .
In conclusion, a system pole may not be a pole of F or vice versa. Note that a system F cannot have more than | m | system poles with respect to m counting their order (see Lemma 3.5 of [40]).
To state the main result (see Theorem 6 below) in [40], we need a generalization of the notion R m ( F ) for a system of functions F . For each system pole λ of F with respect to m , we want to define a corresponding characteristic number R λ ( F , m ) as follows. Let τ be the order of λ as a system pole of F . For each t = 1 , , τ , denote by r λ , t ( F , m ) the largest of all the numbers R t ( g ) (the radius of the largest disk containing at most t poles of g), where g is a polynomial combination of type Equation (8) that is analytic on a neighborhood of B | λ | ¯ except for a pole at z = λ of order t. Then, we define
R λ ( F , m ) : = min k = 1 , 2 , , τ r λ , k ( F , m ) .
As in the scalar case, we denote by Q m ( F ) the monic polynomial whose zeros are the system poles of F with respect to m taking account of their order. The set of distinct zeros of Q m ( F ) is denoted by P ( F , m ) .
An analog of Montessus de Ballore–Gonchar’s theorem for HP approximation (see Theorem 1.4 of [40]) is stated below.
Theorem 6.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . Then, the following two assertions are equivalent:
(a) 
R 0 ( F ) > 0 and F has exactly | m | system poles with respect to m counting multiplicities.
(b) 
The denominators Q n , m , n | m | , of the HP approximants of F are uniquely determined for all sufficiently large n , and there exists a polynomial Q | m | of degree | m | , Q | m | ( 0 ) 0 , such that
lim sup n Q | m | Q n , m 1 / n = θ < 1 ,
where · denotes the coefficient norm in the space of polynomials.
Moreover, if either (a) or (b) takes place, then Q | m | Q m ( F ) and
θ = max | λ | R λ ( F , m ) : P ( F , m ) .
An exact expression for the rate and region of convergence of { R n , m } n | m | to F is given in Theorem 3.7 of [40]. Theorem 6 is considered to be a cornerstone in the study of HP approximation on row sequences in the sense that it contains the first inverse type result.
In the spirit of Gonchar’s theorem and Gonchar’s conjecture, López and Zaldivar [41] proposed the two conjectures below (Conjectures 1 and 2). Note that the notations Δ , σ , γ , and δ j in Conjectures 1–3 and Theorem 7 are defined as in Section 1 taking
P n , m : = { λ n , 1 , λ n , 2 , , λ n , u n } , u n | m | , n N ,
to be the collection of zeros of Q n , m .
Conjecture 1.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . If the denominators Q n , m are uniquely determined for all sufficiently large n and σ ( λ ) 1 , then λ is a system pole of F with respect to m of order τ = σ ( λ ) .
The generalization of “singularity” when we consider a system of functions is stated below.
Definition 4.
Given F = ( F 1 , F 2 , , F d ) as in Equation (5) and m = ( m 1 , m 2 , , m d ) N d , we say that λ C \ { 0 } is a system singularity of F with respect to m if there exists at least one polynomial combination G of the form in Equation (8) analytic on B | λ | and λ is a singular point of G.
Conjecture 2.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . If the denominators Q n , m are uniquely determined for all sufficiently large n and γ ( λ ) 1 , then λ is a system singularity of F with respect to m .
In ([41], p. 155), López and Zaldivar gave an example to support their conjectures. That example shows that a point in C \ { 0 } can be a system pole or a system singularity depending how it attracts zeros of Q n , m . For the self-contained purpose, we show their example here.
Example 3.
Consider F = ( F 1 , F 2 ) , where
F 1 ( z ) : = 1 z 1 + e z and F 2 ( z ) : = log ( z 1 )
and m : = ( 1 , 1 ) . Clear, 1 is a system pole of order 1 and a system singularity of F with respect to m . Their experiment states that deg ( Q n , m ) = 2 for n sufficiently large and
lim sup n | λ n , 1 1 | 1 / n = 0 and | λ n , 2 1 | 1 n , n ,
where Q n , m ( z ) = ( z λ n , 1 ) ( z λ n , 2 ) and | λ n , 1 1 | | λ n , 2 1 | .
Conjecture 1 for the scalar case ( d = 1 ) is the part “(b)⇒(a)” in Gonchar’s theorem. However, Conjecture 2 for the scalar case is Gonchar’s conjecture which remains open. The converse statement of Conjecture 1 proved in Theorem 2.1 of [41] is stated as follows.
Theorem 7.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . Assume that λ is a system pole of order τ of F with respect to m . Then,
Δ ( λ ) | λ | R λ ( F , m ) a n d σ ( λ ) τ .
According to Gonchar’s theorem, a natural conjecture is
Conjecture 3.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . Assume that λ is a system pole of order τ of F with respect to m . Then, the inequalities in Equation (10) are equalities. Moreover, if we assume further that lim inf n | λ λ n , τ + 1 | > 0 , then
δ 1 ( λ ) = δ 2 ( λ ) = = δ τ ( λ ) = | λ | R λ ( F , m ) 1 / τ .
In the past few years, some progress (see [41,42,43]) on Conjecture 2 was made. In those papers, the authors studied Conjecture 2 with an additional assumption, namely,
lim n Q n , m = Q | m | , deg ( Q | m | ) = | m | , Q | m | ( 0 ) 0 .
Without loss of generality, let λ 1 , λ 2 , , λ | m | be the zeros of Q | m | in such a way that
0 < | λ 1 | | λ 2 | | λ | m | |
and λ n , 1 , λ n , 2 , , λ n , | m | be the zeros of Q n , m indexed so that
lim n λ n , k = λ k , k = 1 , 2 , , | m | .
The latest progress on Conjecture 2 with the condition in Equation (11) is the following theorem (see Theorem 4.4 of [42]).
Theorem 8.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . Assume that Equations (11) and (12) hold. Assume further that, if λ j , λ j + 1 , , λ j + N 1 , where N > 1 , are the zeros of Q | m | lying on the same circle (namely, they have the same modulus), then
lim sup n | λ n , k λ k | 1 / n < 1 , k = j , j + 1 , , j + N 1 .
Then, each λ k , k = 1 , 2 , , | m | , is a system singularity of F with respect to m . Moreover, if λ k is a zero of multiplicity τ k of Q | m | which verifies Equation (13), then it is a system pole of F with respect to m of order τ k .
Considering the coefficients of Q n , m as the coefficients of a certain recurrence relation, López and Gerpe [42] employed Buslaev’s results (extensions of Poincaré’s theorem on recursion relations) in Theorems 1 and 2 of [44] to prove the above theorem.
Note that the following special case of Conjecture 2 (an analog of Suetin’s theorem) remains unsolved.
Conjecture 4.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . Assume that Q n , m is unique for all sufficiently large n , Equation (11) takes place, and Q | m | ( λ ) = 0 . Then, λ is a system singularity of F with respect to m .

2.3. Some Remarks

  • The statement (b)⇒(a) in Theorem 6 was the first inverse-type result on the study of HP approximation on row sequences. Its proof is very constructive. It suggests how to find a polynomial combination in Equation (8) verifying that all zeros of Q | m | in Equation (9) are the system poles of F with respect to m .
  • The study of zeros of P n , m , , = 1 , 2 , , d , in Equation (7) is irrelevant to our interest in this paper. However, it is worth mentioning the paper [45] where la Calle Ysern and Mínguez Ceniceros studied the distribution of zeros of P n , m , , = 1 , 2 , , d , as m is fixed and n .

3. Generalized Hermite–Padé Approximations

3.1. Definitions and Notation

After an analog of Montessus de Ballore–Gonchar’s theorem for HP approximation (Theorem 6) was proved, Theorem 6 was extended for various generalizations of HP approximation. These generalizations are formulated to approximate a vector of functions holomorphic on a neighborhood of the following sets E .
Let E be an infinite compact subset of the complex plane C such that C ¯ \ E is simply connected. Whenever we consider set E, E is described as above, unless we specifically say otherwise. Denote by H ( E ) the space of all functions holomorphic in some neighborhood of E . We define
H ( E ) d : = { ( F 1 , F 2 , , F d ) : F H ( E ) , = 1 , 2 , , d } .
In Section 3, we are mainly interested in approximating a vector in the space H ( E ) d .

3.1.1. Orthogonal Hermite–Padé Approximations

The first two approximations are constructed from orthogonal polynomials on E. Let μ be a finite positive Borel measure with infinite support supp ( μ ) contained in E. We write μ M ( E ) and define the associated inner product,
g , h μ : = g ( ζ ) h ( ζ ) ¯ d μ ( ζ ) , g , h L 2 ( μ ) .
Using the Gram–Schmidt process, we can generate the sequence of the orthonormal polynomials with positive leading coefficients
p n ( z ) : = κ n z n + , κ n > 0 , n = 0 , 1 , 2 , ,
corresponding to the inner product in Equation (14). It is well-known that, for each n 0 , such orthonormal polynomial p n is unique. Combining the concepts of HP approximation and orthogonal polynomials, we define two types of orthogonal Hermite–Padé approximations. The first one is a natural extension of HP approximation.
Definition 5.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d and μ M ( E ) . Fix a multi-index m = ( m 1 , m 2 , , m d ) N d . Then, for each n max { m 1 , m 2 , , m d } , there exists Q n , m μ , S P | m | such that Q n , m μ , S 0 and
Q n , m μ , S F , p k μ = 0 for   all   k = n m + 1 , n m + 2 , , n   and   = 1 , 2 , , d .
The corresponding vector of rational functions
R n , m μ , S : = ( R n , m , 1 μ , S , R n , m , 2 μ , S , , R n , m , d μ , S )
= j = 0 n m 1 Q n , m μ , S F 1 , p j μ p j Q n , m μ , S , j = 0 n m 2 Q n , m μ , S F 2 , p j μ p j Q n , m μ , S , , j = 0 n m d Q n , m μ , S F d , p j μ p j Q n , m μ , S
are called an ( n , m ) standard orthogonal Hermite–Padé (SOHP) approximant of F with respect to μ .
When E = { z C : | z | 1 } and μ is the normalized arc length on the unit circle, p n ( z ) = z n for all n 0 and the system of linear equations in Equation (15) reduces to
1 2 π i | z | = 1 Q n , m μ , S ( z ) F ( z ) z k + 1 d z = 0 for   all   k = n m + 1 , n m + 2 , , n   and   = 1 , 2 , , d .
Then, the above polynomial Q n , m μ , S coincides with the polynomial Q in Equation (6), which further implies that, for each = 1 , 2 , , d , the corresponding polynomial j = 0 n m Q n , m μ , S F , p j μ p j coincides with the polynomial P in Equation (6). Moreover, when d = 1 , the rational function R n , m μ , S is the usual ( n m , m ) orthogonal Padé approximant defined in [46]. Therefore, the approximation in Definition 5 is a natural generalization of orthogonal Padé approximation to the vector case.
The definition of the other orthogonal Hermite–Padé approximation was recently introduced by Bosuwan and López [47] in the problem concerning localization of system poles of F H ( E ) d around the set E.
Definition 6.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d and μ M ( E ) . Fix a multi-index m = ( m 1 , m 2 , , m d ) N d and n N . Then, there exists Q n , m μ , M P | m | such that Q n , m μ , M 0 and
z k Q n , m μ , M F , p n μ = 0 for   all   k = 0 , 1 , , m 1   and   = 1 , 2 , , d .
The corresponding vector of rational functions
R n , m μ , M : = ( R n , m , 1 μ , M , R n , m , 2 μ , M , , R n , m , d μ , M )
= j = 0 n 1 Q n , m μ , M F 1 , p j μ p j Q n , m μ , M , j = 0 n 1 Q n , m μ , M F 2 , p j μ p j Q n , m μ , M , , j = 0 n 1 Q n , m μ , M F d , p j μ p j Q n , m μ , M
are called an ( n , m ) modified orthogonal Hermite–Padé (MOHP) approximant of F with respect to μ .
The motivation of the above definition is in the paragraph that contains Equation (39).

3.1.2. Faber–Hermite–Padé Approximations

Now, we want to combine the ideas of HP approximation and Faber polynomials on E. Let Φ be the unique Riemann mapping function from C ¯ \ E to the exterior of the closed unit disk verifying Φ ( ) = and Φ ( ) > 0 . For each ρ > 1 , the level curve of index ρ and the canonical domain of index ρ are defined by
Γ ρ : = { z C : | Φ ( z ) | = ρ } and D ρ : = E { z C : | Φ ( z ) | < ρ } ,
respectively. Let F H ( E ) d . Denote by ρ m ( F ) the index ρ > 1 of the largest canonical domain D ρ to which F has at most m poles counting multiplicities.
The Faber polynomial of E of degree n is
Φ n ( z ) : = 1 2 π i Γ ρ Φ n ( t ) t z d t , z D ρ , n = 0 , 1 , 2 , .
One can also define Φ n as the polynomial part of the Laurent expansion of Φ n at infinity. The nth Faber coefficient of F H ( E ) with respect to Φ n is defined by the formula
[ F ] n : = 1 2 π i Γ ρ F ( t ) Φ ( t ) Φ n + 1 ( t ) d t ,
where ρ ( 1 , ρ 0 ( F ) ) .
As with the SOHP and MOHP approximations, we have two ways to define Faber–Hermite–Padé approximations.
Definition 7.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d . Fix m = ( m 1 , m 2 , , m d ) N d . Then, for each n max { m 1 , m 2 , , m d } , there exists Q n , m E , S P | m | such that Q n , m E , S 0 and
[ Q n , m E , S F ] k = 0 for   all   k = n m + 1 , n m + 2 , , n   and   = 1 , 2 , , d .
The corresponding vector of rational functions
R n , m E , S : = ( R n , m , 1 E , S , R n , m , 2 E , S , , R n , m , d E , S )
= j = 0 n m 1 [ Q n , m E , S F 1 ] j Φ j Q n , m E , S , j = 0 n m 2 [ Q n , m E , S F 2 ] j Φ j Q n , m E , S , , j = 0 n m d [ Q n , m E , S F d ] j Φ j Q n , m E , S
are called an ( n , m ) standard Faber–Hermite–Padé (SFHP) approximant of F with respect to E .
One can easily check that, when E = { z C : | z | 1 } , the SFHP approximant R n , m E , S is the same as a HP approximant R n , m . Furthermore, when d = 1 , the approximant R n , m E , S reduces to the usual ( n m , m ) Padé-Faber approximant (see its definition in Section 3 of [48]).
Definition 8.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d . Fix m = ( m 1 , m 2 , , m d ) N d and n N . Then, there exists Q n , m E , M P | m | such that Q n , m E , M 0 and
[ z k Q n , m E , M F ] n = 0 for   all   l l k = 0 , 1 , , m 1   and   = 1 , 2 , , d .
The corresponding vector of rational functions
R n , m E , M : = ( R n , m , 1 E , M , R n , m , 2 E , M , , R n , m , d E , M )
= j = 0 n 1 [ Q n , m E , M F 1 ] j Φ j Q n , m E , M , j = 0 n 1 [ Q n , m E , M F 2 ] j Φ j Q n , m E , M , , j = 0 n 1 [ Q n , m E , M F d ] j Φ j Q n , m E , M
are called an ( n , m ) modified Faber–Hermite–Padé (MFHP) approximant of F with respect to E .
The above MFHP approximation was recently introduced by Bosuwan and López [49]. In [49], MFHP approximation was proved to be an effective tool to locate system poles of F H ( E ) d around the set E.

3.1.3. Multipiont Hermite–Padé Approximations

Let α E be a table of points; more precisely, α = { α n , k } , k = 1 , , n , n = 1 , 2 , . The definition of multipoint Hermite–Padé approximation is
Definition 9.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d . Set
a n ( z ) : = k = 1 n ( z α n , k )
and fix a multi-index m = ( m 1 , m 2 , , m d ) N d . Then, for each n max { m 1 , m 2 , , m d } , there exist Q n , m α P | m | and P n , m , α P n m for all = 1 , 2 , , d such that Q n , m α 0 and
( Q n , m α F P n , m , α ) / a n + 1 H ( E ) , for   all   = 1 , 2 , , d .
The corresponding vector of rational functions
R n , m α = R n , m , 1 α , R n , m , 2 α , , R n , m , d α = P n , m , 1 α Q n , m α , P n , m , 2 α Q n , m α , , P n , m , d α Q n , m α
is called an ( n , m ) multipoint Hermite–Padé (MHP) approximant of F with respect α.
This approximation was introduced in [2] by Mahler long before other generalizations. Note that, if all the interpolation points are 0 , then the MHP approximants reduce to HP approximants.

3.1.4. Some Remarks

  • Finding Q n , m μ , S , Q n , m μ , M , Q n , m E , S , o r Q n , m E , M is equivalent to solving | m | + 1 unknowns from linear system of | m | equations in Equations (15)–(18), respectively. Moreover, finding P n , m , α , = 1 , 2 , , d and Q n , m α is equivalent to solving ( n + 1 ) d + 1 unknowns from linear system of ( n + 1 ) d equations in Equation (19). Therefore, these polynomials, Q n , m μ , S , Q n , m μ , M , Q n , m E , S , Q n , m E , M , and Q n , m α , always exist but may not be unique. Since such polynomials are not the zero function, we normalize them to be “monic” polynomials. Moreover, Q n , m α in Equation (20) is chosen so that it does not have a common zero with all the P n , m , α . We would like to emphasize that for any ( n , m ) N × N d , R n , m μ , S , R n , m μ , M , R n , m E , S , R n , m E , M , and R n , m α may not be unique.
  • Extensions of generalized HP approximations in Definitions 5–9 to a compact set E whose complement is connected are possible. However, the results in this survey paper are restricted to the case when E is a compact subset of the complex plane with simply connected complement in the extended complex plane. This is because, for the sets E which are disconnected, the zeros of the corresponding orthonormal polynomials p n , second type functions s n (defined in Equation (21)), or Faber polynomials Φ n may lie in C \ E which may be the locations of system poles.

3.1.5. Classes of Measures in M ( E )

For the studies of SOHP and MOHP approximations, we need to impose asymptotic properties of the sequences of the orthonormal polynomials { p n } n 0 and the corresponding second type functions { s n } n 0 defined below
s n ( z ) : = p n ( ζ ) ¯ z ζ d μ ( ζ ) , z C ¯ \ s u p p ( μ ) .
We keep classes of measures in M ( E ) here so that they do not disturb the flow of our paper. The readers who are not interested in SOHP or MOHP approximations may skip Section 3.1.5.
Definition 10.
Let μ M ( E ) .
(a) 
μ Reg 1 ( E ) if and only if
lim n | p n ( z ) | 1 / n = | Φ ( z ) | .
(b) 
μ Reg 2 ( E ) if and only if
lim n | s n ( z ) | 1 / n = | Φ ( z ) | 1 .
(c) 
μ Reg 1 , 2 ( E ) if and only if μ Reg 1 ( E ) Reg 2 ( E ) .
(d) 
μ Reg 1 , 2 m ( E ) if and only if μ Reg 1 , 2 ( E ) and there exists a positive constant c such that
κ n m κ n c , n n 0 .
(e) 
μ Rat 1 ( E ) if and only if
lim n p n ( z ) p n + 1 ( z ) = 1 Φ ( z ) .
(f) 
μ Rat 2 ( E ) if and only if
lim n s n + 1 ( z ) s n ( z ) = 1 Φ ( z ) .
(g) 
μ Rat 1 , 2 ( E ) if and only if μ Rat 1 ( E ) Rat 2 ( E ) .
(h) 
μ S ( E ) if and only if lim n c n c n + 1 = 1 and
lim n p n ( z ) c n Φ n ( z ) = S ( z ) ,
where the c n s are positive constants and S is a non-vanishing holomorphic function on C ¯ \ E .
The limits in Equations (22) and (23) are assumed to hold uniformly on compact subsets of C \ E and the ones in Equations (24)–(26) are assumed to hold uniformly on compact subsets of C ¯ \ E .
The classes Reg 1 ( E ) and Reg 2 ( E ) are more or less the same in some cases. In particular, if E is convex, then Reg 1 ( E ) = Reg 2 ( E ) and these two classes coincide with the regular class in the usual sense (see Definition 3.1.2 of [50] for the definition of the regular class in the usual sense). Clearly, S ( E ) Rat 1 ( E ) Reg 1 ( E ) and Rat 2 ( E ) Reg 2 ( E ) .
Denote by K 1 the collection of all compact sets E (stated at the beginning of Section 3) satisfying the condition that the inverse of the exterior conformal function Φ 1 can be extended continuously to C ¯ \ { w C : | w | < 1 } . If we assume that E K 1 and μ Rat 1 ( E ) , then
lim n p n ( z ) s n ( z ) = Φ ( z ) Φ ( z ) ,
uniformly on compact subsets of C ¯ \ E (see Lemma 3.1 of [46]). This statement implies that, if E K 1 , then S ( E ) Rat 1 ( E ) Rat 2 ( E ) . Moreover, it is well-known that, if μ Rat 1 ( E ) , then
lim n κ n κ n + 1 = c a p ( E )
(see, e.g., Lemma 2 of [51]). Therefore, for each m N , Rat 1 ( E ) Reg 1 m ( E ) and Rat 1 , 2 ( E ) Reg 1 , 2 m ( E ) . For a general compact set E , the class of measures Reg 1 ( E ) has been well studied and characterized in terms of the analytic properties of the measure or of the corresponding sequence of leading coefficients (see, e.g., [50], Theorem 3.1.1). The situation is not quite the same for other classes of measures. To discuss them would take us too far from our main direction; rather we refer the reader to Pages 532–533 in [52] for more details and references.

3.1.6. Classes of Tables of Interpolation Points

When we state results on MHP approximation, we refer to two classes of tables of interpolation points.
Definition 11.
Let α E be a table of interpolation points, namely
α = { α n , k } , k = 1 , , n , n = 1 , 2 , .
Set
a n ( z ) : = k = 1 n ( z α n , k ) .
(a) 
α Strong ( E ) if and only if the corresponding polynomials a n satisfy the following strong asymptotics:
lim n a n ( z ) ( c a p ( E ) ) n Φ n ( z ) = G ( z ) 0 ,
uniformly on compact subsets of C ¯ \ E .
(b) 
α Root ( E ) if and only if the corresponding polynomials a n satisfy the following nth root asymptotics:
lim n | a n ( z ) | 1 / n = c a p ( E ) | Φ ( z ) | ,
uniformly on compact subsets of C ¯ \ E ,
Clearly, Strong ( E ) Root ( E ) . Moreover, it is well known that Strong ( E ) and Root ( E ) (see Chapters 8 and 9 of [53] for more details about both classes).

3.2. Results and Conjectures

3.2.1. The Scalar Case

We add a section dedicated to the scalar case of all generalizations (defined in Section 3.1) here. The purpose is not only to help the readers consolidate their understanding of the vector case but also to discuss some interesting theorems and conjectures corresponding to such scalar case. When d = 1 , we have m = m and write
R n , m μ , S = R n , m μ , S , R n , m μ , M = R n , m μ , M , R n , m E , S = R n , m E , S , R n , m E , M = R n , m E , M , R n , m α = R n , m α .
The first theorem simply says that one can use the poles of the mth row sequences { R n , m μ , S } n m , { R n , m μ , M } n m , { R n , m E , S } n m , { R n , m E , M } n m , or { R n , m α } n m to detect the m poles of F H ( E ) nearest E (with respect to the level curves).
Theorem 9.
Let F H ( E ) and m N be fixed. Then, the following statements are equivalent.
(a) 
F has exactly m poles counting multiplicities in D ρ m ( F ) .
(b) 
For all n sufficiently large, Q n , m μ , S corresponding to μ Rat 1 , 2 ( E ) has degree m and there exists a polynomial Q ^ m of degree m such that
lim sup n Q n , m μ , S Q ^ m 1 / n = θ ^ < 1 .
(For μ Rat 1 , 2 ( E ) , the statement of the theorem does not include this assertion.)
(c) 
For all n sufficiently large, Q n , m μ , M corresponding to μ Reg 1 , 2 m ( E ) has degree m and there exists a polynomial Q ˇ m of degree m such that
lim sup n Q n , m μ , M Q ˇ m 1 / n = θ ˇ < 1 .
(For μ Reg 1 , 2 m ( E ) , the statement of the theorem does not include this assertion.)
(d) 
For all n sufficiently large, Q n , m E , S has degree m and there exists a polynomial Q ¨ m of degree m such that
lim sup n Q n , m E , S Q ¨ m 1 / n = θ ¨ < 1 .
(e) 
For all n sufficiently large, Q n , m E , M has degree m and there exists a polynomial Q ˜ m of degree m such that
lim sup n Q n , m E , M Q ˜ m 1 / n = θ ˜ < 1 .
(f) 
For all n sufficiently large, Q n , m α corresponding to α Root ( E ) has degree m and there exists a polynomial Q ˘ m of degree m such that
lim sup n Q n , m α Q ˘ m 1 / n = θ ˘ < 1 .
(For α Root ( E ) , the statement of the theorem does not include this assertion.)
Moreover, if one of Assertions (a)–(f) takes place, then
(i) 
Q ^ m = Q ˇ m = Q ¨ m = Q ˜ m = Q ˘ m = Q m ( F ) ,
where Q m ( F ) is the monic polynomial whose zeros are the poles of F in D ρ m ( F ) taking account of their order.
(ii) 
θ ^ = θ ˇ = θ ¨ = θ ˜ = θ ˘ = max { | Φ ( λ ) | : λ P m ( F ) } ρ m ( F ) ,
where P m ( F ) is the set of the distinct zeros of Q m ( F ) .
(iii) 
For any compact subset K of D ρ m ( F ) \ P m ( F ) ,
lim sup n F R n , m μ , S K 1 / n Φ K ρ m ( F ) ,
lim sup n F R n , m μ , M K 1 / n Φ K ρ m ( F ) ,
lim sup n F R n , m E , S K 1 / n Φ K ρ m ( F ) ,
lim sup n F R n , m E , M K 1 / n Φ K ρ m ( F ) ,
lim sup n F R n , m α K 1 / n Φ K ρ m ( F ) ,
under the convention that if K E , then Φ K is replaced by 1 .
Remark 1.
The inequality in Equation (31) is equality when α is a Newton type.
The proofs concerning (a)⇔(b), (a)⇔(c), (a)⇔(d), (a)⇔(e), and (a)⇔(f) are in Corollary 1 of [52], Theorem 1.2 of [47], Theorem 2 of [54], Corollary 1.6 of [49], and Theorem 1.3 of [55], respectively.
We strongly believe that for this scalar case, the exact rates of convergences of { R n , m μ , S } n m , { R n , m μ , M } n m , { R n , m E , S } n m , { R n , m E , M } n m , and { R n , m α } n m , should be as indicated in the following conjecture.
Conjecture 5.
Let F H ( E ) and m N . If F has exactly m poles in D ρ m ( F ) , then the inequalities in Equations (27)–(31) are equalities.
Clearly, if E = { z C : | z | 1 } and d μ = d θ / 2 π on the boundary of E , then R n , m μ , S , R n , m μ , M , R n , m E , S and R n , m E , M are the ( n m , m ) classical Padé approximant and Conjecture 5 holds true because of Equation (4) in Montessus de Ballore–Gonchar’s theorem. Another supporting evidence for the validity of Conjecture 5 is provided for the case when K = E , namely if F has exactly m poles in D ρ m ( F ) , then
lim sup n F R n , m μ , S E 1 / n = lim sup n F R n , m μ , M E 1 / n = lim sup n F R n , m E , S E 1 / n = lim sup n F R n , m E , M E 1 / n
= lim sup n F R n , m α E 1 / n = lim sup n σ n , m 1 / n = 1 ρ m ( F ) ,
where
σ n , m : = inf r F r E ,
and the infimum is taken over the class of all rational functions of type ( n , m )
r ( z ) = a n z n + a n 1 z n 1 + + a 0 b m z m + b m 1 z m 1 + + b 0 .
We refer the reader to [56] for the proof of the fifth equality in Equation (32).
For Gonchar’s conjecture, we made good progress in the case when the denominators Q n , m μ , S , Q n , m μ , M , Q n , m E , S , Q n , m E , M , and Q n , m α converge to a polynomial of degree m as n . An analog of Suetin’s theorem for the scalar case of all generalized HP approximations (see Theorem 4 of [44], Theorem 2 of [57], Theorem 2.7 of [46], and Theorems 2.5 and 2.6 of [58]) is the following.
Theorem 10.
Let F H ( E ) and m N be fixed. If one of the following holds:
(a) 
E K 1 , μ S ( E ) , for all n sufficiently large, Q n , m μ , S has degree m, and lim n Q n , m μ , S = Q m ;
(b) 
E K 1 , μ S ( E ) , for all n sufficiently large, Q n , m μ , M has degree m, and lim n Q n , m μ , M = Q m ;
(c) 
for all n sufficiently large, Q n , m E , S has degree m and lim n Q n , m E , S = Q m ;
(d) 
for all n sufficiently large, Q n , m E , M hhas degree m and lim n Q n , m E , M = Q m ; or
(e) 
α Strong ( E ) , for all n sufficiently large, Q n , m α has degree m, and lim n Q n , m α = Q m ,
then all of the following hold:
(i) 
ρ 0 ( F ) = min { | Φ ( λ ) | : Q m ( λ ) = 0 } ;
(ii) 
ρ m 1 ( F ) = max { | Φ ( λ ) | : Q m ( λ ) = 0 } ;
and
(iii) 
all zeros of Q m are singularities of F ; those lying in D ρ m 1 ( F ) are poles (counting multiplicities), and F has no other poles in D ρ m 1 ( F ) .
The proof of the above theorem relies on deep results on refinements of Poincaré’s theorem on recurrence relations developed by Buslaev in Theorems 5 and 6 of [57]. These Buslaev’s results connect Suetin’s theorem (for classical Padé approximants) and analogous ones for generalized Padé approximants (see Theorem 3 of [57] and Theorems 2.7 and 2.8 of [58]).
An analog of Gonchar’s theorem for Q n . m α when α is a Newton type was proved by himself in Theorem 2 of [31]. For other generalizations, we have proofs only for an analog of (a)⇒(b) in Gonchar’s theorem (see Corollary 5.6 of [51], Corollary 1 of [59], Theorems 2.1 and 2.3 of [58]).
A direct analog of the Fabry ratio theorem for orthogonal and Faber polynomial expansions is
Theorem 11.
Let F H ( E ) . If one of the following holds:
(a) 
lim n [ F ] n [ F ] n + 1 = ξ ;
or
(b) 
E K 1 , μ S ( E ) , and
lim n F , p n μ F , p n + 1 μ = ξ ,
then Φ 1 ( ξ ) is a singularity of F and ρ 0 ( F ) = | ξ | .
The part of the above theorem concerning the limit of the ratio of Faber coefficients follows from a straightforward change of variables:
[ F ] n = 1 2 π i Γ ρ F ( t ) Φ ( t ) Φ n + 1 ( t ) d t = 1 2 π i | z | = ρ F ( Φ 1 ( z ) ) z n + 1 d z , where   t = Φ 1 ( z )
and the use of the Fabry ratio theorem for the analytic part of F Φ 1 . The part of the above theorem concerning the limit of the ratio of Fourier coefficients is in Theorem 2.3 of [46].
We end this section by stating an application of Theorem 11 to orthogonal polynomial theory on the unit circle. Firstly, we recall some facts. If E = B 1 ¯ and the measure μ supported on the unit circle satisfies the Szegő condition,
0 2 π log w ( θ ) d θ > ,
where d μ ( θ ) = w ( θ ) d θ / 2 π + d μ s ( θ ) is the Radon–Nikodym decomposition of μ , then it is well known that μ S ( E ) , the leading coefficients of the orthonormal polynomials p n satisfy
lim n κ n = κ : = exp 1 4 π 0 2 π log w ( θ ) d θ ,
and
1 S i n t ( z ) = 1 κ k = 0 p k ( 0 ) ¯ p k ( z ) , uniformly   on   compact   subsets   of   B ,
where
S int ( z ) : = exp 1 4 π 0 2 π log w ( θ ) e i θ + z e i θ z d θ , z B ,
denotes the interior Szegő function (see ([60], pp. 19–20) for the proof). Therefore, Theorem 11 enables us to locate the singularity nearest the origin of the reciprocal of the interior Szegő function 1 / S i n t in terms of the Verblunsky coefficients
α n : = p n ( 0 ) ¯ / κ n .
Theorem 12.
Let μ satisfy the Szegő condition in Equation (33) and assume that 1 / S i n t H ( B ¯ ) . Suppose that the Verblunsky coefficients α n corresponding to μ verify
lim n α n α n + 1 = λ .
Then, λ is a singularity of 1 / S i n t and 1 / S i n t is holomorphic on B | λ | .
The above result was stated in Corollary 2.4 of [46].

3.2.2. The Vector Case

Graves-Morris and Saff were the first to prove an extension of Montessus de Ballore’s theorem to MHP approximation in Theorem 3 of [38]. Their theorem relies on the concept of polewise independence, which influenced the current author to define the following adapted polewise independence.
Definition 12.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d be a vector of functions meromorphic in some canonical domain D ρ and let m = ( m 1 , m 2 , , m d ) N d be the multi-index. Then, the function F is said to be polewise independent with respect to the multi-index m in D ρ if and only if there do not exist polynomials v 1 , v 2 , , v d , at least one of which being non-null, satisfying:
(i) 
deg v m 1 , = 1 , 2 , , d .
(ii) 
= 1 d ( v Φ ) · F H ( D ρ \ E ) .
Note that, if we replace Φ in (ii) of Definition 12 by the identity mapping, then the above definition reduces the definition of polewise independence in [38].
Putting together Theorem 2.3 in [51] and Theorem 1 in [61], we obtain the following analog of Montessus de Ballore’s theorem for SOHP and SFHP approximations under the concept of polewise independence in Definition 12.
By Q | m | F , we denote the monic polynomial whose zeros are the poles of F in D ρ | m | ( F ) taking account of their order. The set of distinct zeros of Q | m | F is denoted by P | m | ( F ) .
Theorem 13.
Let m N d be a fixed multi-index, F = ( F 1 , F 2 , , F d ) H ( E ) d , and μ Rat 1 , 2 ( E ) . Suppose that F is polewise independent (in the sense of Definition 12) with respect to the multi-index m in D ρ | m | ( F ) . Then, R n , m μ , S and R n , m E , S are uniquely determined for all sufficiently large n , and for each = 1 , 2 , , d and for any compact set K D ρ | m | ( F ) \ P | m | ( F ) ,
lim sup n F R n , m , μ , S K 1 / n Φ K ρ | m | ( F ) and lim sup n F R n , m , E , S K 1 / n Φ K ρ | m | ( F ) .
Moreover,
lim sup n Q n , m μ , S Q | m | F 1 / n max λ P | m | ( F ) | Φ ( λ ) | ρ | m | ( F ) ,
lim sup n Q n , m E , S Q | m | F 1 / n max λ P | m | ( F ) | Φ ( λ ) | ρ | m | ( F ) .
In fact, Cacoq and López [62] were the first who proved the above theorem but in the context of SOHP approximation for the case when E is the closed unit disk and μ is supported on the unit circle.
Making use of incomplete orthogonal Padé and incomplete Padé–Faber approximations (see Definition 5.1 of [51] and Definition 5 of [59]), Bosuwan proved another Montessus de Ballore type theorem for SOHP and SFHP approximations in Theorem 2.4 of [51] and Theorem 1 of [59]. The idea of incomplete approximants allows us to analyze each function F individually. Let us define some more notation about the region of convergence.
Given a system F = ( F 1 , F 2 , , F d ) and a multi-index m = ( m 1 , m 2 , , m d ) N d \ { 0 } , we define
D m ( F ) : = ( D ρ m 1 ( F 1 ) , D ρ m 2 ( F 2 ) , , D ρ m d ( F d ) ) .
By Q m F , we denote the monic polynomial whose zeros are the poles of F in D m ( F ) counting multiplicities. This set of poles is denoted by P m ( F ) . For = 1 , 2 , , d , set P m , ( F ) : = P m ( F ) D ρ m ( F ) .
To each pole λ of F in this system of domains
D m ( F ) : = ( D ρ m 1 ( F 1 ) , D ρ m 2 ( F 2 ) , , D ρ m d ( F d ) ) ,
we associate an index ( λ ) { 1 , 2 , , d } as follows. The index ( λ ) verifies that λ D ρ m ( λ ) ( F ( λ ) ) and λ is a pole of F ( λ ) of the same order as is a pole of F in D m ( F ) . If there are several indices satisfying this condition, then we choose one among those with greatest ρ m ( F ) .
Theorem 14.
Let m N d be a fixed multi-index, F = ( F 1 , F 2 , , F d ) H ( E ) d , P m ( F ) = { λ 1 , λ 2 , , λ q } , and μ Rat 1 , 2 ( E ) . Suppose that F H ( E ) d has exactly | m | poles in D m ( F ) . Then, R n , m μ , S and R n , m E , S is uniquely determined for all sufficiently large n and for each = 1 , 2 , , d , R n , m , μ , S and R n , m , E , S converge uniformly to F on compact subsets of D ρ m ( F ) \ P m , ( F ) . Moreover, for all = 1 , 2 , , d ,
lim sup n F R n , m , μ , S K 1 / n Φ K ρ m ( F ) and lim sup n F R n , m , E , S K 1 / n Φ K ρ m ( F ) ,
where K is any compact subset of D ρ m ( F ) \ P m , ( F ) and if K E , then Φ K is replaced by 1 . Additionally, we have
lim sup n Q m F Q n , m μ , S 1 / n max j = 1 , 2 , , q | Φ ( λ j ) | ρ m ( λ j ) ( F ( λ j ) ) ,
lim sup n Q m F Q n , m E , S 1 / n max j = 1 , 2 , , q | Φ ( λ j ) | ρ m ( λ j ) ( F ( λ j ) ) ,
Note that Theorems 13 and 14 have their own values. Section 3 in [51] and Section 2 in [59] gave examples showing when Theorem 13 is applicable but Theorem 14 is not applicable or vice versa.
During my visit to Universidad Carlos III de Madrid in the summer of 2016, López and I initiated a project emphasizing the study of analogs of the inverse statement of Montessus de Ballore–Gonchar’s Theorem for SOHP and SFHP approximations. As guided by the principle of the definition of system pole in [40] (see Definition 3), we define
Definition 13.
Given F = ( F 1 , F 2 , , F d ) H ( E ) d and m = ( m 1 , m 2 , , m d ) N d , we say that λ C is a system pole of order τ of F with respect to m if τ is the largest positive integer such that, for each t = 1 , 2 , , τ , there exists at least one polynomial combination of the form
= 1 d v F , deg ( v ) < m , = 1 , 2 , , d ,
which is holomorphic in a neighborhood of D ¯ | Φ ( λ ) | except for a pole at z = λ of exact order t .
We had an application in my mind that we wanted to use the zeros of { Q n , m μ , S } n | m | and { Q n , m E , S } n | m | to detect | m | system poles of F H ( E ) d nearest the set E . Although we are not able to solve our initial aim (Conjecture 6 below), we observed that the system poles of F with respect to m and their orders are the same as the system poles of
F ¯ : = ( F 1 , , z m 1 1 F 1 , F 2 , , z m d 1 F d )
with respect to
m ¯ = ( 1 , 1 , , 1 ) , w i t h | m | = | m ¯ |
and their orders. The creation of the vector in Equation (39) motivated us to define MOHP and MFHP approximations. Importantly, the zeros of { Q n , m μ , M } n | m | and { Q n , m E , M } n | m | enable us to detect | m | system poles of F H ( E ) d nearest the set E (see [47,49]).
To state main results in [47,49,55], we need a generalization of R λ ( F , m ) . For each system pole λ of F with respect to m , we define a characteristic index ρ ¯ λ ( F , m ) as follows. Let τ be the order of λ as a system pole of F . For each t = 1 , 2 , , τ , denote by ρ λ , t ( F , m ) the largest of all the numbers ρ t ( G ) (the index of the largest canonical domain containing at most t poles of G), where G is a polynomial combination of the type in Equation (38) that is holomorphic in a neighborhood of D ¯ | Φ ( λ ) | except for a pole at z = λ of order t . There is only a finite number of such possible values so the maximum is indeed attained. Then, we define
ρ ¯ λ ( F , m ) : = min t = 1 , 2 , , τ ρ λ , t ( F , m ) .
Combining Theorem 1.2 in [47], Corollary 1.6 in [49], and Theorem 1.3 in [55], we arrive at the following theorem which contains analogs of Montessus de Ballore–Gonchar’s Theorem for MOHP, MFHP, and MHP approximations.
Theorem 15.
Let F = ( F 1 , F 2 , , F d ) H ( E ) d , m N d be a fixed multi-index, μ Reg 1 , 2 | m | ( E ) , and α Strong ( E ) . Denote by Q m ( F ) the monic polynomial whose zeros are the system poles of F with respect to m taking account of their order and by P ( F , m ) the set of all zeros of Q m ( F ) . Then, the following assertions are equivalent:
(a) 
F has exactly | m | system poles with respect to m counting multiplicities.
(b) 
The polynomials Q n , m μ , M are uniquely determined for all sufficiently large n , and there exists a polynomial Q ^ | m | of degree | m | such that
lim sup n Q n , m μ , M Q ^ | m | 1 / n = θ ^ < 1 .
(c) 
The polynomials Q n , m E , M are uniquely determined for all sufficiently large n , and there exists a polynomial Q ˜ | m | of degree | m | such that
lim sup n Q n , m E , M Q ˜ | m | 1 / n = θ ˜ < 1 .
(d) 
The polynomials Q n , m α are uniquely determined for all sufficiently large n , and there exists a polynomial Q ˇ | m | of degree | m | such that
lim sup n Q n , m α Q ˇ | m | 1 / n = θ ˇ < 1 .
Moreover, if one of Assertions (a)–(d) takes place, then Q ^ | m | = Q ˜ | m | = Q ˇ | m | = Q m ( F ) and
θ ^ = θ ˜ = θ ˇ = max | Φ ( λ ) | ρ ¯ λ ( F , m ) : λ P ( F , m ) .
Remark 2.
(i) 
To prove that (a) implies (b) with Equation (40) replaced by
lim sup n Q n , m μ , M Q m ( F ) 1 / n < 1 ,
we only need to impose that μ Reg 2 ( E ) .
(ii) 
To prove that (a) implies (d) with Equation (41) replaced by
lim sup n Q n , m α Q m ( F ) 1 / n < 1 ,
we only need to impose that α Root ( E ) .
(iii) 
For the estimates on the convergences of { R n , m μ , M } n | m | , { R n , m E , M } n | m | , and { R n , m α } n | m | , we refer the readers to Theorem 1.2 of [47], Theorem 1.4 of [49], and Theorem 1.3 of [55] to avoid introducing a complicated collection of notations.
As Conjecture 5 (the scalar case) stands open, the exact rate of convergences { R n , m μ , M } n | m | , { R n , m E , M } n | m | , and { R n , m α } n | m | are still unknown.
Our expected result for the study SOHP and SFHP approximations is stated below.
Conjecture 6.
Let F H ( E ) d , m N d be a fixed multi-index, and μ Rat 1 , 2 ( E ) . Then, the following assertions are equivalent:
(a) 
F has exactly | m | system poles with respect to m counting multiplicities.
(b) 
The polynomials Q n , m μ , S of F are uniquely determined for all sufficiently large n , and there exists a polynomial Q ^ | m | of degree | m | such that
lim sup n Q n , m μ , S Q ^ | m | 1 / n = θ ^ < 1 .
(c) 
The polynomials Q n , m E , S of F are uniquely determined for all sufficiently large n , and there exists a polynomial Q ˜ | m | of degree | m | such that
lim sup n Q n , m E , S Q ˜ | m | 1 / n = θ ˜ < 1 .
Moreover, if one of Assertions (a)–(c) takes place, then Q ^ | m | = Q ˜ | m | = Q m ( F ) and
θ ^ = θ ˜ = θ ˇ = max | Φ ( λ ) | ρ ¯ λ ( F , m ) : λ P ( F , m ) .
Concerning Conjectures 1–4 replacing Q n , m by Q n , m μ , S , Q n , m μ , M , Q n , m E , S , Q n , m E , M , and Q n , m α , much less is known except the scalar case (see Theorem 10 above). Especially, the structures of R n , m μ , S , R n , m μ , M , R n , m E , S , and R n , m E , M are much more complicated than HP approximants. There is no nice equality (similar to Equality (3.3) in [41]) which allows us to say something about the convergence and divergence of incomplete Padé approximation appear in SOHP, MOHP, SFHP, and MFHP approximations. Such equality is also a main ingredient for the proofs of Gonchar’s theorem and Suetin’s theorem for classical Padé approximation stated in Section 1 and an analog of Gonchar’s theorem for the scalar MHP approximation when α is a Newton type in Theorem 2 of [31].

4. Conclusions

The study of classical Padé approximants on row sequences roots from the work of Montessus de Ballore [29] on the uniform convergence of row sequences of the approximants. The subject received renewed interest by Gonchar [30], who proved the converse statement of Montessus de Ballore’s theorem. Soon after that, he studied an attraction of an individual pole of F in B R m ( F ) to the zeros of Q n , m when m is fixed as n in [31]. In the same paper, Gonchar proposed his conjecture (which is commonly called Gonchar’s conjecture and remains unsolved). Since its introduction, several positive answers have supported the conjecture (see, e.g., [32,33,34,35]).
Graves-Morris and Saff [38] extended Montessus de Ballore’s theorem to the vector case (which we call HP approximation). Recently, Cacoq, de la Calle Ysern, and López [40] defined the notion of system pole and proved the inverse result of Graves-Morris and Saff’s result. In the spirit of Gonchar’s conjecture, López and Gerpe proposed the following conjectures.
Conjecture 1.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . If the denominators Q n , m are uniquely determined for all sufficiently large n and σ ( λ ) 1 , then λ is a system pole of F with respect to m of order τ = σ ( λ ) .
Conjecture 2.
Let F be a vector of formal Taylor expansions at the origin and fix a multi-index m N d . If the denominators Q n , m are uniquely determined for all sufficiently large n and γ ( λ ) 1 , then λ is a system singularity of F with respect to m .
Some progress on Conjecture 2 was made in [41,42,43]. However, both conjectures are still open. Concurrently, generalizations of HP approximation were introduced and studied on row sequences [47,49,51,55,61]. Both conjectures can be asked for all generalizations of HP approximation. Although SOHP and SFHP approximations are natural, the proofs of analogs of Montessus de Ballore–Gonchar’s theorem for such approximations are not available (see Conjecture 6 for our expectation).
Analogs of Montessus de Ballore–Gonchar’s theorem and Suetin’s theorem for the scalar case of all generalizations were completely proved. However, the exact rates of convergences of { R n , m μ , S } n m , { R n , m μ , M } n m , { R n , m E , S } n m , { R n , m E , M } n m , and { R n , m α } n m (when α is not a Newton type) are unknown (see Conjecture 5 for our expectation).

Funding

The research of N. Bosuwan was supported by the Strengthen Research Grant for New Lecturer from the Thailand Research Fund and the Office of the Higher Education Commission (MRG6080133) and the Faculty of Science, Mahidol University.

Acknowledgments

I wish to express my gratitude toward the anonymous referees and the editor for helpful comments and suggestions leading to improvements of this work. I also want to thank Guillermo López Lagomasino and Edward Barry Saff for introducing me to the area of Padé approximation theory.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
HP approximationHermite–Padé approximation
SOHP approximationstandard orthogonal Hermite–Padé approximation
MOHP approximationmodified orthogonal Hermite–Padé approximation
SFHP approximationstandard Faber–Hermite–Padé approximation
MFHP approximationmodified Faber–Hermite–Padé approximation
MHP approximationmultipoint Hermite–Padé approximation

References

  1. Hermite, C. Sur la fonction exponentielle. C R. Acad. Sci. Paris 1873, 77, 18–24, 74–79, 226–233, 285–293. [Google Scholar]
  2. Mahler, K. Perfect systems. Compos. Math. 1968, 19, 95–166. [Google Scholar]
  3. Coates, J. On the algebraic approximation of functions. I, II, III. Indag. Math. 1966, 28, 421–461. [Google Scholar] [CrossRef] [Green Version]
  4. Jager, H. A simultaneous generalization of the Padé table. I–VI. Indag. Math. 1964, 26, 193–249. [Google Scholar] [CrossRef] [Green Version]
  5. Zudilin, W. Arithmetic of linear forms involving odd zeta values. J. Théorie Nombres Bordx. 2004, 16, 251–291. [Google Scholar] [CrossRef]
  6. Ball, K.; Rivoal, T. Irrationalité d’une infinit’e de valeurs de la fonction zêta aux entiers impairs. Invent. Math. 2001, 146, 193–207. [Google Scholar] [CrossRef]
  7. Lindemann, F. Über die Zahl π. Math. Ann. 1882, 20, 213–225. [Google Scholar] [CrossRef]
  8. Apéry, R. Irrationalité de ζ(2) et ζ(3). Astérisque 1979, 61, 11–13. [Google Scholar]
  9. Van Assche, W. Analytic number theory and rational approximation. In Coimbra Lecture Notes on Orthogonal Polynomials; Branquinho, A., Foulquié, A., Eds.; Nova Science Pub.: New York, NY, USA, 2008; pp. 197–229. [Google Scholar]
  10. Beckermann, B.; Labahn, G. A uniform approach for Hermite Padé and simultaneous Padé approximants and their matrix-type generalizations. Numer. Algorithms 1992, 3, 45–54. [Google Scholar] [CrossRef]
  11. Beckermann, B.; Labahn, G. A uniform approach for the fast computation of matrix-type Padé approximants. SIAM J. Matrix Anal. Appl. 1994, 15, 804–823. [Google Scholar] [CrossRef]
  12. Beckermann, B.; Labahn, G. Fraction-free computation of matrix rational interpolants and matrix GCDs. SIAM J. Matrix Anal. Appl. 2000, 22, 114–144. [Google Scholar] [CrossRef]
  13. Borges, C.F. On a class of Gauss-like quadrature rules. Numer. Math. 1994, 67, 271–288. [Google Scholar] [CrossRef]
  14. Cabay, S.; Jones, A.R.; Labahn, G. Computation of numerical Padé-Hermite and simultaneous Padé systems II: A weakly stable algorithm. SIAM J. Matrix Anal. Appl. 1996, 17, 268–297. [Google Scholar] [CrossRef] [Green Version]
  15. Cabay, S.; Labahn, G. A superfast algorithm for multi-dimensional Padé systems. Numer. Algorithms 1992, 2, 201–224. [Google Scholar] [CrossRef]
  16. Fidalgo Prieto, U.; Illán, J.; López Lagomasino, G. Hermite-Padé approximation and simultaneous quadrature formulas. J. Approx. Theory 2004, 126, 171–197. [Google Scholar] [CrossRef] [Green Version]
  17. Lindman, E.L. Free-space boundary conditions for the time dependent wave equation. J. Comput. Phys. 1975, 18, 66–78. [Google Scholar] [CrossRef]
  18. Coussement, J.; Van Assche, W. Gaussian quadrature for multiple orthogonal polynomials. J. Comput. Appl. Math. 2005, 178, 131–145. [Google Scholar] [CrossRef] [Green Version]
  19. Kuijlaars, A.B.J. Multiple orthogonal polynomial ensembles. Recent trends in orthogonal polynomials and approximation theory. Contemp. Math. 2010, 507, 155–176. [Google Scholar]
  20. Aptekarev, A.I. Asymptotics of simultaneously orthogonal polynomials in the Angelesco case. Math. USSR Sb. 1989, 64, 57–84. [Google Scholar] [CrossRef]
  21. Aptekarev, A.I. Strong asymptotics of multiply orthogonal polynomials for Nikishin systems. Sb. Math. 1999, 190, 631–669. [Google Scholar] [CrossRef]
  22. Martín, P.; Baker, G.A., Jr. Two-point quasifractional approximant in physics. Truncation error. J. Math. Phys. 1991, 32, 1470–1477. [Google Scholar] [CrossRef]
  23. Aptekarev, A.; Kaliaguine, V.; Iseghem, J.V. The genetic sums’ representation for the moments of a system of Stieltjes functions and its application. Constr. Approx. 2000, 16, 487–524. [Google Scholar] [CrossRef]
  24. Daems, E.; Kuijlaars, A.B.J. Multiple orthogonal polynomials of mixed type and non-intersecting Brownian motions. J. Approx. Theory 2007, 146, 91–114. [Google Scholar] [CrossRef] [Green Version]
  25. Bleher, P.M.; Kuijlaars, A.B.J. Random matrices with external source and multiple orthogonal polynomials. Int. Math. Res. Not. 2004, 2004, 109–129. [Google Scholar] [CrossRef]
  26. Beckermann, B.; Kalyagin, V.; Matos, A.; Wielonsky, F. How well does the Hermite-Padé approximation smooth the Gibbs phenomenon? Math. Comp. 2011, 80, 931–958. [Google Scholar] [CrossRef] [Green Version]
  27. Shang, Y. Analytical solution for an in-host viral infection model with time-inhomogeneous rates. Acta Phys. Polon. B 2015, 46, 1567–1577. [Google Scholar] [CrossRef] [Green Version]
  28. Van Assche, W. Padé and Hermite-Padé approximation and orthogonality. Surv. Approx. Theory 2006, 2, 61–91. [Google Scholar]
  29. De Montessus de Ballore, R. Sur les fractions continues algébrique. Bull. Soc. Math. Fr. 1902, 30, 28–36. [Google Scholar] [CrossRef]
  30. Gonchar, A.A. On convergence of Padé approximants for some classes of meromorphic functions. Sb. Math. 1975, 26, 555–575. [Google Scholar] [CrossRef]
  31. Gonchar, A.A. Poles of rows of the Padé table and meromorphic continuation of functions. Sb. Math. 1981, 43, 527–546. [Google Scholar] [CrossRef]
  32. Vavilov, V.V. On the singular points of a meromorphic function given by Its Taylor series. Dokl. Akad. Nauk SSSR 1976, 231, 1281–1284. [Google Scholar]
  33. Vavilov, V.V.; López Lagomasino, G.; Prokhorov, V.A. On an Inverse Problem for the Rows of a Padé Table. Mat. Sb. 1979, 110, 117–129. [Google Scholar] [CrossRef]
  34. Vavilov, V.V.; Prokhorov, V.A.; Suetin, S.P. The poles of the mth row of the Padé table and the singular points of a function. Mat. Sb. 1983, 122, 475–480. [Google Scholar]
  35. Suetin, S.P. On an inverse problem for the mth row of the Padé table. Sb. Math. 1985, 52, 231–244. [Google Scholar] [CrossRef]
  36. Fabry, E. Sur les points singuliers d’une fonction données par son développement de Taylor. Ann. École Norm. Sup. Paris 1896, 13, 367–399. [Google Scholar] [CrossRef] [Green Version]
  37. Nikishin, E.M. Rational Approximations and Orthogonality; Amer. Math. Soc.: Providence, RI, USA, 1991. [Google Scholar]
  38. Graves-Morris, P.R.; Saff, E.B. A de Montessus theorem for vector valued rational interpolants. In Rational Approximation and Interpolation; Springer: Berlin/Heidelberg, Germany, 1984; pp. 227–242. [Google Scholar]
  39. Cacoq, J.; de la Calle Ysern, B.; López Lagomasino, G. Incomplete Padé approximation and convergence of row sequences of Hermite-Padé approximants. J. Approx. Theory 2013, 170, 59–77. [Google Scholar] [CrossRef]
  40. Cacoq, J.; de la Calle Ysern, B.; López Lagomasino, G. Direct and inverse results on row sequences of Hermite-Padé approximation. Constr. Approx. 2013, 38, 133–160. [Google Scholar] [CrossRef] [Green Version]
  41. Lagomasino, G.L.; Gerpe, Y.Z. Inverse results on row sequences of Hermite-Padé approximation. Proc. Steklov Inst. Math. 2017, 298, 152–169. [Google Scholar] [CrossRef]
  42. López, G.L.; Gerpe, Y.Z. Higher order recurrences and row sequences of Hermite-Padé approximation. J. Differ. Equ. Appl. 2018, 24, 1830–1845. [Google Scholar] [CrossRef]
  43. Lagomasino, G.L. On row sequences of Padé and Hermite-Padé approximation. In Proceedings of the Modern Trends in Constructive Function Theory: Conference in Honor of Ed Saff’s 70th Birthday: Constructive Functions 2014, Vanderbilt University, Nashville, TN, USA, 26–30 May 2014; Volume 661. [Google Scholar]
  44. Buslaev, V.I. Relations for the coefficients, and singular points of a function. Math. USSR-Sb. 1988, 59, 349–377. [Google Scholar] [CrossRef]
  45. De la Calle Ysern, B.; Mínguez Ceniceros, J. Zero distribution of incomplete Padé and Hermite-Padé approximations. J. Approx. Theory 2016, 201, 13–29. [Google Scholar] [CrossRef]
  46. Bosuwan, N.; López Lagomasino, G.; Saff, E.B. Determining singularities using row sequences of Padé-orthogonal approximants. Jaen J. Approx. 2013, 5, 179–208. [Google Scholar]
  47. Bosuwan, N.; López Lagomasino, G. Determining system poles using row sequences of orthogonal Hermite-Padé approximants. J. Approx. Theory 2018, 231, 15–40. [Google Scholar] [CrossRef] [Green Version]
  48. Suetin, S.P. On the convergence of rational approximations to polynomial expansions in domains of meromorphy of a given function. Math. USSR Sb. 1978, 34, 367–381. [Google Scholar] [CrossRef]
  49. Bosuwan, N.; López Lagomasino, G. Direct and inverse results on row sequences of simultaneous Padé-Faber approximants. Mediterr. J. Math. 2019, 16, 36. [Google Scholar] [CrossRef] [Green Version]
  50. Stahl, H.; Totik, V. General Orthogonal Polynomials; Cambridge University Press: Cambridge, UK, 1992; Volume 43. [Google Scholar]
  51. Bosuwan, N. Convergence of row sequences of simultaneous Padé-orthogonal approximants. Comput. Methods Funct. Theory 2017, 17, 525–556. [Google Scholar] [CrossRef]
  52. Bosuwan, N.; López Lagomasino, G. Inverse theorem on row sequences of linear Padé-orthogonal approximants. Comput. Methods Funct. Theory 2015, 15, 529–554. [Google Scholar] [CrossRef] [Green Version]
  53. Walsh, J.L. Interpolation and Approximation by Rational Functions in the Complex Domain, 5th ed.; Colloquium Publications, American Mathematical Society: Providence, RI, USA, 1969. [Google Scholar]
  54. Suetin, S.P. Inverse theorems on generalized Padé approximants. Math. USSR Sb. 1980, 37, 581–597. [Google Scholar] [CrossRef]
  55. Bosuwan, N.; Lagomasino, G.L.; Gerpe, Y.Z. Direct and inverse results for multipoint Hermite-Padé approximants. Anal. Math. Phys. 2019. [Google Scholar] [CrossRef] [Green Version]
  56. Saff, E.B. Regions of meromorphy determined by the degree of best rational approximation. Proc. Am. Math. Soc. 1971, 29, 30–38. [Google Scholar] [CrossRef]
  57. Buslaev, V.I. An analogue of Fabry’s theorem for generalized Padé approximants. Math. Sb. 2009, 200, 39–106. [Google Scholar] [CrossRef]
  58. Bosuwan, N. Direct and inverse results on row sequences of generalized Padé approximants to polynomial expansions. Acta Math. Hung. 2019, 157, 191–219. [Google Scholar] [CrossRef]
  59. Bosuwan, N. On Montessus de Ballore’s theorem for simultaneous Padé-Faber approximants. Demonstr. Math. 2018, 51, 45–61. [Google Scholar] [CrossRef]
  60. Geronimus, L.Y. Orthogonal Polynomials on a Circle and Interval; Pergamon Press: Oxford, UK, 1960. [Google Scholar]
  61. Bosuwan, N. Convergence of row sequences of simultaneous Padé-Faber approximants. Math. Notes 2018, 103, 643–656. [Google Scholar] [CrossRef] [Green Version]
  62. Cacoq, J.; López Lagomasino, G. Convergence of row sequences of simultaneous Fourier-Padé approximation. Jaen J. Approx. 2012, 4, 101–120. [Google Scholar]

Share and Cite

MDPI and ACS Style

Bosuwan, N. On Row Sequences of Hermite–Padé Approximation and Its Generalizations. Mathematics 2020, 8, 366. https://doi.org/10.3390/math8030366

AMA Style

Bosuwan N. On Row Sequences of Hermite–Padé Approximation and Its Generalizations. Mathematics. 2020; 8(3):366. https://doi.org/10.3390/math8030366

Chicago/Turabian Style

Bosuwan, Nattapong. 2020. "On Row Sequences of Hermite–Padé Approximation and Its Generalizations" Mathematics 8, no. 3: 366. https://doi.org/10.3390/math8030366

APA Style

Bosuwan, N. (2020). On Row Sequences of Hermite–Padé Approximation and Its Generalizations. Mathematics, 8(3), 366. https://doi.org/10.3390/math8030366

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop