Next Article in Journal
New Developments in Geometric Function Theory II
Previous Article in Journal
Iteration with Bisection to Approximate the Solution of a Boundary Value Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimal Property of B-Bases for the Modified Richardson Method

by
Juan Manuel Peña
Departmento de Matemática Aplicada/IUMA, Universidad de Zaragoza, 50009 Zaragoza, Spain
Submission received: 28 February 2024 / Revised: 21 March 2024 / Accepted: 22 March 2024 / Published: 27 March 2024
(This article belongs to the Section Mathematical Analysis)

Abstract

:
A space with a normalized totally positive basis has a unique normalized B-basis. In computer-aided geometric design, normalized B-bases present optimal shape-preserving properties. More optimal properties of normalized B-bases were proved previously. This paper provides a new optimal property concerning the modified Richardson iterative method when applied to collocation matrices of normalized B-bases. Moreover, a similar optimal property is proved for the tensor product of normalized B-bases.

1. Introduction

Spaces generated by univariate functions often possess totally positive bases, that is, bases whose collocation matrices satisfy the property of being totally positive and, therefore, having all their minors nonnegative (see Section 2). Totally positive matrices are also called totally nonnegative matrices and they belong to the theory of total positivity. This theory has more than one century of history (see [1,2,3,4]). Several applications of totally positive matrices can be seen in the following books or surveys [5,6,7,8,9,10,11]. These applications arise in many different fields. For instance, in approximation theory (e.g., [12]), differential equations (e.g., [11]), chemistry (e.g., [13]), combinatorics (e.g., [14]), mechanics (e.g., [11]), statistics (e.g., [15]), computer-aided geometric design (e.g., [16]), Markov chains (e.g., [17]), numerical analysis (e.g., [18,19]), quantum groups (e.g., [9,20]), accurate computations (e.g., [21]), economics (e.g., [22]) and biomathematics (e.g., [8]). A space with a totally positive basis always has a special totally positive basis that generates all remaining totally positive bases of the space by using totally positive matrices (see [23]). This special basis is called a B-basis.
In the field of computer-aided geometric design, normalized totally positive bases are the bases leading to representations with shape-preserving properties. If a space has a normalized totally positive basis, then it has a unique normalized B-basis. In the case of the space of polynomials of degree less than or equal to n on [ 0 , + ) , the monomial polynomials 1 , t , t 2 , , t n form a B-basis. For the same space of polynomials on a compact interval, the Bernstein polynomials (see Section 2) form the normalized B-basis. A conjecture in the field of computer-aided geometric design proposed by Goodman and Said in [24] claims that the Bernstein basis has optimal shape-preserving properties in the sense that the control polygon is the closest in shape to the generated curve. In [25], Carnicer and Peña proved it, and in [23], they proved that the normalized B-basis is always the normalized totally positive basis with optimal shape-preserving properties. Other optimal properties of B-bases and normalized B-bases are recalled in Section 2, and this paper provides a new optimal property of normalized B-bases.
Another remarkable property of several subclasses of totally positive matrices (e.g., [26,27,28,29,30,31]) comes from the fact that for them, high relative accuratcy algorithms were found to solve fundamental linear algebra computations, such as the computation of their inverses, their eigenvalues, their singular values or the solutions of some systems of linear equations A x = b (where b has alternating signs). However, even for these subclasses of totally positive matrices, there are no accurate methods for solving all linear systems. Hence, the minimal conditioning property of collocation matrices of normalized B-bases is still important, as well as the use of alternative methods, such as iterative methods. Among the classical iterative methods used for totally positive linear systems, we can outline the Richardson method (cf. Example 4.1 of [32,33,34]), which has already played a key role in the progressive iterative approximation (PIA) method. The PIA (see [35]) has been frequently used in computer-aided geometric design and it converges when the bases are normalized totally positive. In this paper, we prove that the nonsingular collocation matrices of normalized B-bases provide the fastest convergence rate of the modified Richardson method among the corresponding collocation matrices of all normalized totally positive bases of the space.
Let us now present the paper’s organization. Section 2 introduces some notations and auxiliary results for totally positive matrices, totally positive bases and B-bases. In particular, some optimal properties of B-bases and normalized B-bases are recalled, as well as many examples of B-bases and normalized B-bases. The modified Richardson iterative method is presented in Section 3, where the mentioned result on its fastest convergence speed for nonsingular collocation matrices of normalized B-bases is proved. In this section, it is also shown that the optimal parameter for the modified Richardson iterative method of a nonsingular collocation matrix of the normalized B-basis of a space U is the minimum among the optimal parameters of the corresponding collocation matrices of all NTP bases of U.
Although the concept of total positivity cannot be extended to multivariate spaces, in Section 4, we extend the result on the fastest convergence speed of the modified Richardson method for nonsingular collocation matrices of normalized B-bases to the case of tensor product spaces. Collocation matrices of tensor product totally positive bases lead to the Kronecker product of their corresponding collocation matrices. We extended the result in spite of the fact that the Kronecker product of totally positive matrices is not a totally positive matrix. Finally, Section 5 summarizes the main conclusions of the paper.

2. B-Bases and Optimal Properties

Let us introduce some matrix notations and basic definitions. A matrix is totally positive (TP) if all its minors are nonnegative (see [5,10]). As already commented in the introduction, TP matrices are also called totally nonnegative matrices in the literature and they arise in many scientific fields, such as approximation theory, differential equations, combinatorics, mechanics, statistics, computer-aided geometric design and biomathematics (e.g., [5,8,10]).
The previous matrix concept leads to a corresponding concept for systems of univariate functions. Let u 0 , , u n be functions of a space U of functions defined on I R . Given a sequence of arbitrary parameters t 1 < < t n + 1 in I, the corresponding collocation matrix at these parameters is given by
M ( t 1 , , t n + 1 ) : = u j 1 ( t i ) 1 i , j n + 1 .
The system ( u 0 , , u n ) of functions defined on I R is totally positive (TP) if all its collocation matrices M ( t 1 , , t n + 1 ) are TP. A TP system of functions on I is called normalized (NTP) if i = 0 n u i ( t ) = 1 for all t I . NTP systems are frequently used in computer-aided geometric design [36] due to their nice shape-preserving properties (see [23,25]) in the sense that the curve imitates the shape of their control polygon. Observe that the collocation matrix of an NTP basis is row stochastic because it is nonnegative and the entries of each row sum to 1.
B-bases are TP bases that generate all TP bases of a space by means of TP matrices. The next characterization of a B-basis is a consequence of Corollary 3.10 of [23] and Proposition 3.11 of [23].
Theorem 1. 
Let ( u 0 , , u n ) be a TP basis of a given space U of univariate functions on I. Then, ( u 0 , , u n ) is a B-basis if and only if for another TP basis ( v 0 , , v n ) of U, the matrix K of the change in basis such that ( v 0 , , v n ) = ( u 0 , , u n ) K is TP.
The following result shows the existence of a B-basis in a space with a TP basis (which was proved in Remark 3.8 of [23]) and the existence and uniqueness of a normalized B-basis if the space has an NTP basis (which was proved in Theorem 4.2 (i) of [23]).
Theorem 2. 
Let U be a space of univariate functions on an interval I. If U has a TP basis, then it has a B-basis, and, if U has an NTP basis, then it has a unique normalized B-basis.
As we shall recall later, among all NTP bases of a given space, the normalized B-basis is the basis satisfying the optimal shape-preserving properties (cf. [23]). As first examples, we can mention that the Bernstein bases and the B-spline bases are the normalized B-bases of their corresponding spaces. As an example of a B-basis that is not normalized, we can mention the monomial basis of the space of polynomials of degree at most n on I = [ 0 , + ) (cf. [25]). At the end of the section, a list of B-bases and normalized B-basis is included.
Now we recall some optimal properties of B-bases from several points of view:
1.
Normalized B-bases present optimal shape-preserving properties in the sense that their control polygon is closer in shape than the control polygon with respect to another NTP basis. In the case of the space of polynomials of degree at most n on [ 0 , 1 ] , this property was first conjectured by Goodman and Said for the Bernstein basis in [24] and later proved by Carnicer and Peña in [25]. In [23], this optimal property was proved for all normalized B-bases.
2.
B-bases are optimally stable with respect to the evaluation among all bases of nonnegative functions. This property was first proved for the Bernstein basis in [37] and later in [38], as well as for many more B-bases in [39].
3.
B-bases are the least supported bases, as shown in [40].
4.
Nonsingular collocation matrices of nomalized B-bases have a minimal condition number for · among the corresponding matrices of all NTP bases of the space, as shown in [41].
In Section 3, we show an additional property of the collocation matrices of normalized B-bases: for them, the modified Richardson method is the fastest among the corresponding matrices of all the NTP bases of the space.
A new field where TP matrices have shown their importance is that of accurate computations. Up to now, algorithms with a high relative accuracy for linear algebra problems have been found only for a few classes of structured matrices. Among them, we have some subclasses of TP matrices (e.g., [26,27,28,29,31]) for which there are high relative accuracy algorithms for computing all their eigenvalues, all their singular values, their inverses and some associated linear systems. In all cases, the starting point is the bidiagonal decomposition of the matrix in order to use the accurate algorithms of [42,43]. The fact that even for these subclasses of TP matrices, we do not have accurate algorithms for solving all linear systems shows the importance of the optimal condition property (mentioned above) of collocation matrices associated with normalized B-bases, as well as the possible convenience of using iterative methods, as shown in the next sections.
We now finish this section by including a list of B-bases and normalized B-bases that illustrate the many types of examples that arise in theoretical and practical applications, as well as the potential applicability of the results of this paper:
1.
In the space P n ( [ a , + ) ) of polynomials of degree less than or equal to n on [ a , + ) , the monomial basis ( 1 , ( t a ) , ( t a ) 2 , , ( t a ) n ) is a B-basis (cf. Section 6 of [25]).
2.
In the space P n ( [ a , b ] ) of polynomials of degree less than or equal to n on a compact interval [ a , b ] , the normalized B-basis is the corresponding Bernstein basis of polynomials ( b 0 n ( t ; a , b ) , , b n n ( t ; a , b ) ) , with
b i n ( t ; a , b ) ( t ) : = n i ( b t ) n i ( t a ) i ( b a ) n , i = 0 , 1 , , n
(cf. [23,25]). The standard Bernstein basis used in computer-aided design is the Bernstein basis on [ 0 , 1 ] :
b i n ( t ) : = n i ( 1 t ) n i ( t ) i , i = 0 , 1 , , n .
3.
Given a sequence ( w i ) 0 i n of positive weights, we can define the following system of rational functions ( r 0 n , , r n n ) on the compact interval [ a , b ] by
r i n ( t ) = w i b i n ( t ; a , b ) j = 0 n w j b j n ( t ; a , b ) , i = 0 , 1 , , n ,
which is the normalized B-basis of the corresponding spanned space of functions. It is the rational Bernstein basis of its space. When the weight w i = 1 for all i, then ( r 0 n , , r n n ) = ( b 0 n ( t ; a , b ) , , b n n ( t ; a , b ) ) is the Bernstein basis on [ a , b ] .
4.
In the space of exponential functions
E : = { i = 0 n c i exp ( λ i t ) | c i , i = 0 , , n } , λ 0 < < λ n , t ,
the basis ( exp ( λ 0 t ) , , exp ( λ n t ) ) is a B-basis.
5.
In the space of Müntz polynomials on [ 0 , + ) , where
M : = { i = 0 n c i t λ i | c i , i = 0 , , n } , λ 0 < < λ n ,
the generalized monomial basis ( t λ 0 , , t λ n ) is a B-basis. The restriction of Müntz polynomials to [ 0 , 1 ] has a B-basis: the generalization of the Bernstein basis given by
p i ( t ) = ( 1 ) n i ( λ i + 1 λ 0 ) ( λ i + 2 λ 0 ) ( λ n λ 0 ) t λ [ λ i , , λ n ] ,
where t λ [ λ i , , λ n ] is the usual notation of divided differences of the function t λ , that is, t λ [ λ i ] : = t λ i and
t λ [ λ i , , λ i + j ] : = t λ [ λ i + 1 , , λ i + j ] t λ [ λ i , , λ i + j 1 ] λ i + j λ i .
6.
The space of trigonometric polynomials
T n = { 1 , cos t , sin t , cos 2 t , sin 2 t , , cos n t , sin n t }
on I = [ A , A ] with A < π 2 has the normalized B-basis ( v 0 , , v m ) , m = 2 n , given by
v i ( t ) = d i sin A + t 2 sin A i sin A t 2 sin A m i , i = 0 , 1 , , m
where
d i = k = 0 [ i / 2 ] m / 2 i k i k k ( 2 cos A ) i 2 k , i = 0 , 1 , , m .
7.
The space of even trigonometric functions given by
C n = s p a n { 1 , cos t , cos 2 t , , cos n t }
on the compact interval [ 0 , π ] has the normalized B-basis ( u 0 n , , u n n ) defined by
u i n ( t ) = n i cos 2 ( n i ) ( t / 2 ) sin 2 i ( t / 2 ) , i = 0 , 1 , , n .
8.
Given a sequence of positive weights ( w i ) 0 i n and a knots vector ( t 0 , , t n + d ) with t i t i + 1 for all i = 0 , 1 , , n + d 1 , we can define the B-spline basis ( N 0 , d , N 1 , d , , N n , d ) over the previous knots vector by
N i , 0 ( t ) = 1 , if t i t < t i + 1 , 0 , otherwise , N i , k ( t ) = t t i t i + k t i N i , k 1 ( t ) + t i + k + 1 t t i + k + 1 t i + 1 N i + 1 , k 1 ( t ) , k = 1 , , d ,
which is the normalized B-basis of the corresponding splines space (see [23]). The basis ( r 0 , , r n ) defined by
r i ( t ) = w i N i , d ( t ) j = 0 n w j N j , d ( t ) , i = 0 , 1 , n ,
is the normalized B-basis of the corresponding NURBS space (see Section 4 of [23]).
9.
The space given by
P 1 = span { 1 , t , cos t , sin t }
on [ 0 , p ] , with p < 2 π , has the normalized B-basis ( Z 0 ( t ) , Z 1 ( t ) , Z 2 ( t ) , Z 3 ( t ) ) . To give the expression of these functions, we use the following cycloidal functions:
sin C ( t ) : = t sin t , cos C ( t ) : = 1 cos t , tan C ( t ) : = sin C ( t ) cos C ( t ) .
Using these cycloidal functions, we now define the functions of the normalized B-basis:
Z 3 ( t ) = sin C ( t ) S , Z 0 ( t ) = Z 3 ( p t ) , Z 2 ( t ) = M cos C ( t ) C Z 3 ( t ) , Z 1 ( t ) = Z 2 ( p t ) ,
where
S = sin C ( p ) , C = cos C ( p ) , T = tan C ( p ) ,
and
M = 1 if p = π , sin p p 2 T otherwise .
Finally, let us mention that for many spaces of the form
P n = span { 1 , t , , t n cos t , sin t } ,
as well as for other mixed spaces generated by algebraic, trigonometric and even hyperbolic polynomials, the normalized B-bases were obtained, as can be seen in papers published by Carnicer, Mainar and Peña.

3. A New Optimal Property of Normalized B-Bases

This section considers Richardson iteration. Before presenting the method, we come back to its connection with the PIA method, which was already announced in the Introduction and for which Richardson iteration plays a key role (cf. [44]). The PIA method was used for curve and spline curve fitting 50 years ago (see [45,46]) and was later extended to curves generated by NTP bases. Other extensions of the method were derived (cf. [47,48]). Nowadays, the PIA is very powerful in computer-aided geometric design for curves and surface fitting.
Let us start by recalling the Richardson iterative method, which is a classical iterative method for solving linear systems of equations. If A is a nonsingular matrix, then the Richardson method for solving the linear system A x = b is given through the recurrence
x m + 1 = x m A x m + b , m = 0 , 1 , 2 ,
If we denote by y the new iteration coming from x, the method can be written in the form
y = ( I A ) x + b ,
where I denotes the identity matrix.
If one wants to accelerate the convergence of the method, then a relaxation of the method is often performed, replacing the x m + 1 by ( 1 w ) x m + w x m + 1 . In the particular case of the Richardson method, this approach derives the modified Richardson method
y = ( I w A ) x + w b ,
The following result is contained in Theorem 2.3 of [44] and analyzes the convergence of the modified Richardson method.
Theorem 3. 
If A is a nonsingular TP row-stochastic matrix, then the modified Richardson iterative method (5) converges to the solution of the system A x = b if and only if w ( 0 , 2 ) .
The following result shows that normalized B-bases satisfy an optimal property for the modified Richardson method in the sense that it converges faster for their nonsingular collocation matrices than for the corresponding collocation matrices of other NTP bases.
Theorem 4. 
Let U be a given space of functions defined on an interval I with an NTP basis, let ( b 0 , , b n ) be a normalized B-basis of U and let w ( 0 , 2 ) . Then, the modified Richardson method for w and for any nonsingular collocation matrix of ( b 0 , , b n ) has the fastest convergence rate among all the corresponding collocation matrices of all NTP bases of U.
Proof. 
The existence of the normalized B-basis ( b 0 , , b n ) is guaranteed by Theorem 1. If B is a nonsingular collocation matrix of ( b 0 , , b n ) , let A be the corresponding collocation matrix at the same points of another NTP basis ( u 0 , , u n ) .
By Corollary 6.6 of [5], a nonsingular TP matrix has all positive eigenvalues. Therefore, let us denote by λ min ( B ) ( > 0 ) and λ min ( A ) ( > 0 ) the minimal eigenvalues of the matrices B and A, respectively. Let us also take into account the fact that since both bases ( b 0 , , b n ) and ( u 0 , , u n ) are NTP, the matrices B and A are both row stochastic. Therefore, 1 is their maximal eigenvalue. Since A and B are nonsingular row-stochastic TP matrices and w ( 0 , 2 ) , the modified Richardson method converges for both matrices by Theorem 3.
Looking at Equation (5), we conclude that the convergence speed of the modified Richardson methods for matrices B and A are related to ρ ( I w B ) and ρ ( I w A ) , respectively. Let us see that ρ ( I w B ) ρ ( I w A ) and the result is proved.
Clearly, the eigenvalues of matrices I w B and I w A are of the form 1 w λ , where λ is an eigenvalue of B or A. Since w > 0 , we have that [ 1 w , 1 w λ min ( B ) ] and [ 1 w , 1 w λ min ( A ) ] are the least intervals containing the spectrum of I w B and the spectrum of I w A , respectively. Since A is a collocation matrix of an NTP basis and B is the corresponding collocation matrix of the normalized B-basis of the space, Corollary 2 of [41] implies that
λ min ( A ) λ min ( B ) .
Then, since w > 0 , we deduce from (6) that
1 w λ min ( A ) 1 w λ min ( B ) .
On the other hand, we also have that
ρ ( I w A ) = max ( | w 1 | , | 1 w λ min ( A ) | ) , ρ ( I w B ) = max ( | w 1 | , | 1 w λ min ( B ) | ) .
Let us assume first that 1 w λ min ( B ) < 0 . Since λ min ( B ) 1 , we can deduce that 1 w λ min ( B ) < 0 implies that w > 1 and | 1 w λ min ( B ) | = w λ min ( B ) 1 . Hence, again taking into account the fact that λ min ( B ) 1 and (8), we can deduce that
| 1 w λ min ( B ) | w 1 ( = ρ ( I w B ) ) max ( | w 1 | , | 1 w λ min ( A ) | ) = ρ ( I w A )
and the result follows in this case.
It remains the case that 1 w λ min ( B ) 0 . Therefore, | 1 w λ min ( B ) | = 1 w λ min ( B ) . By (7), we also have that 1 w λ min ( A ) 0 , and again, | 1 w λ min ( A ) | = 1 w λ min ( A ) . Hence, the result follows from (7) and (8):
ρ ( I w A ) = max ( | w 1 | , | 1 w λ min ( A ) | ) max ( | w 1 | , | 1 w λ min ( B ) | ) = ρ ( I w B ) .
In the second part of Theorem 2.3 of [44], the parameter corresponding to the optimal (fastest) convergence speed for the modified Richardson method applied to a nonsingular totally nonnegative row-stochastic matrix is obtained, as the following result shows.
Theorem 5. 
Let A be a nonsingular TP row-stochastic matrix. Then, the optimal convergence speed for the modified Richardson iterative method with w ( 0 , 2 ) corresponds to ρ opt = ρ ( I w opt A ) and is achieved for
w opt = 2 1 + λ min ( A ) , ρ opt = 1 λ min ( A ) 1 + λ min ( A ) .
We finish this section with an observation derived from Theorem 5.
Remark 1. 
If we denote by w opt ( C ) the parameter corresponding with the optimal convergence speed for the modified Richarson iterative method of a nonsingular totally nonnegative row-stochastic matrix C, by B a nonsingular collocation matrix of the normalized B-basis of a space U and by A the collocation matrix of another NTP basis of U, we can deduce from (6) and (9) that
w opt ( B ) w opt ( A ) ,
that is, the optimal parameter for the modified Richarson iterative method of a nonsingular collocation matrix of the normalized B-basis of a space U is the minimum among the optimal parameters of the corresponding collocation matrices of all NTP bases of U.
In spite of the fact that the concept of total positivity cannot be extended to multivariate spaces, in the following section, we extend Theorem 4 to the case of tensor product spaces.

4. Tensor Product Case

The practical interest of considering the Richardson iteration method for the Kronecker product of matrices can be illustrated again by the PIA method. For instance, in [49], the PIA method, and thus, the Richardson iteration method for the Kronecker product of matrices, was used for the approximation of tensor product surfaces.
We start by introducing the tensor product of systems of functions and the Kronecker product of matrices.
Let us consider two systems u 1 = ( u 0 1 , , u m 1 ) and u 2 = ( u 0 2 , , u n 2 ) of univariate functions defined on [ a , b ] and [ c , d ] , respectively. Then, the system u 1 u 2 : = ( u i 1 ( x ) · u j 2 ( y ) ) i = 0 , , m j = 0 , , n is named a tensor product system and generates a tensor product surface. Now, we introduce the Kronecker product of two square matrices A = ( a i j ) 1 i , j m and B = ( b i j ) 1 i , j n as the matrix A B defined as the m n × m n block matrix
A B : = a 11 B a 1 m B a m 1 B a m m B .
If we consider the collocation matrices B 1 : = ( u j 1 ( x i ) ) 0 i , j m and B 2 : = ( u j 2 ( y i ) ) 0 i , j n of u 1 and u 2 , then B 1 B 2 is the collocation matrix of u 1 u 2 at ( ( x i , y j ) j = 0 , , n ) i = 0 , , m .
The following result shows that the tensor product of normalized B-bases also satisfies an optimal property for the modified Richardson method in the sense that it converges faster for their nonsingular collocation matrices than for the corresponding collocation matrices of other tensor products of NTP bases.
Theorem 6. 
Let U 1 and U 2 be two spaces of univariate functions on [ a , b ] and [ c , d ] , respectively, both with NTP bases. Let b 1 = ( b 0 1 , , b m 1 ) and b 2 = ( b 0 2 , , b n 2 ) be the normalized B-bases of U 1 and U 2 , respectively. Let w ( 0 , 2 ) . Then, the modified Richardson method for w and for any nonsingular collocation matrix of b 1 b 2 has the fastest convergence rate among all the corresponding collocation matrices of the tensor product of NTP bases of U 1 and U 2 .
Proof. 
Let u 1 = ( u 0 1 , , u m 1 ) be any NTP basis on [ a , b ] of U 1 and let u 2 = ( u 0 2 , , u n 2 ) be any NTP basis on [ c , d ] of U 2 . The existence of the normalized B-bases b 1 and b 2 is guaranteed by Theorem 1. If B is a nonsingular collocation matrix of b 1 b 2 , then B = B 1 B 2 , where B 1 and B 2 are collocation matrices of b 1 and b 2 , respectively. Let A be the corresponding collocation matrix at the same points of u 1 u 2 , and thus, A = A 1 A 2 , where A 1 and A 2 are collocation matrices of u 1 and u 2 , respectively.
Since B 1 and B 2 are row-stochastic matrices, it follows from the definition of the Kronecker product that B = B 1 B 2 is also row stochastic and, by the same reason, A = A 1 A 2 is, again, row stochastic.
Given C 1 , which is an n × n matrix, and C 2 , which is an m × m matrix, consider C 1 C 2 . If we have an eigenvalue λ of C 1 and an eigenvalue μ of C 2 , then by Theorem 4.2.12 of [50], λ μ is an eigenvalue of C 1 C 2 and every eigenvalue of C 1 C 2 arises as a product of eigenvalues of C 1 and C 2 . Then, since B = B 1 B 2 is nonsingular, it does not have a zero eigenvalue, and thus, B 1 and B 2 do not have zero eigenvalues either and they are nonsingular. Hence, A 1 and A 2 are nonsingular and the same property holds for A = A 1 A 2 . Taking into account the fact that A 1 , A 2 , B 1 and B 2 are all nonsingular TP matrices, they have all positive eigenvalues and, since they are also row stochastic, their maximal eigenvalue is 1. In conclusion, the eigenvalues of A and B are also positive and their maximal eigenvalue is also 1. In particular, λ min ( A ) > 0 and λ min ( B ) > 0 . The eigenvalues of I w B and I w A are 1 w λ , where λ is an eigenvalue of B or A. Since w > 0 , then [ 1 w , 1 w λ min ( A ) ] is the least interval containing the spectrum of I w A , and the same holds for B. Since w ( 0 , 2 ) , 1 w > 1 , and so ρ ( I w A ) < 1 and ρ ( I w B ) < 1 . Thus, the modified Richardson method converges for A and for B. It remains to compare the convergence speeds. We only need to prove that ρ ( I w B ) ρ ( I w A ) .
Since w ( 0 , 2 ) , [ 1 w , 1 w λ min ( B ) ] and [ 1 w , 1 w λ min ( A ) ] are the least intervals containing the spectrum of I w B and the spectrum of I w A , respectively. Taking into account the fact that by Theorem 1 (ii) of [51], λ min ( A ) λ min ( B ) , and reasoning as in the end of Theorem 4, we can deduce that
ρ ( I w A ) = max ( | w 1 | , | 1 w λ min ( A ) | ) max ( | w 1 | , | 1 w λ min ( B ) | ) = ρ ( I w B ) ,
and the result follows. □
Observe that in the proof of the previous theorem, we guaranteed the convergence of the modified Richardson method for w ( 0 , 2 ) in spite of the fact that the Kronecker product of TP matrices is not necessarily TP. In fact, if we consider the TP matrices
A = 1 0 1 1 , B = 1 1 0 1 ,
then we have that
A B = 1 1 0 0 0 1 0 0 1 1 1 1 0 1 0 1
is not TP because, for instance, its submatrix M formed by rows 3 and 4 and columns 2 and 3:
M = 1 1 1 0
has a negative determinant: det M = 1 < 0 .

5. Conclusions

In addition to several known optimal properties of normalized B-bases of spaces of univariate functions, this paper provides a new one. The modified Richardson method converges faster for their nonsingular collocation matrices than for the corresponding collocation matrices of other normalized totally positive bases of the space. We also prove that the optimal parameter associated with the fastest convergence is smaller for the collocation matrices of the normalized B-basis than the optimal parameter associated with the fastest convergence of the corresponding collocation matrices of the remaining normalized totally positive bases of the space.
The tensor product of normalized B-bases presents, for the generated spaces of bivariate functions, a similar optimal property. The modified Richardson method converges faster for their nonsingular collocation matrices than for the corresponding collocation matrices of other tensor products of normalized totally positive bases. These last collocation matrices are given by the Kronecker product of the collocation matrices of the normalized totally positive bases.

Funding

This research was partially supported by the Spanish research grants PID2022-138569NB-I00 (MCIU/AEI) and RED2022-134176-T (MCI/AEI) and by Gobierno de Aragón (E41_23R).

Data Availability Statement

Data sharing is not applicable.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PIAprogressive iterative approximation
TPtotally positive
NTPnormalized totally positive

References

  1. Fekete, M.; Polya, G. Uber ein Problem von Laguerre. Rend. Circ. Matem. Palermo 1912, 34, 89–120. [Google Scholar] [CrossRef]
  2. Loewner, C. On totally positive matrices. Math. Z. 1955, 63, 338–340. [Google Scholar] [CrossRef]
  3. Schoenberg, I. Über variationsvermindernde lineare Transformationen. Math. Z. 1930, 32, 321–328. [Google Scholar] [CrossRef]
  4. Whitney, A.M. A reduction theorem for totally positive matrices. J. d’Analyse Math. 1952, 2, 88–92. [Google Scholar] [CrossRef]
  5. Ando, T. Totally positive matrices. Linear Algebra Appl. 1987, 90, 165–219. [Google Scholar] [CrossRef]
  6. Fallat, S.M.; Johnson, C.R. Totally Nonnegative Matrices; Princeton Series in Applied Mathematics; Princeton University: Princeton, NJ, USA, 2011. [Google Scholar]
  7. Fomin, S.; Zelevinsky, A. Double Bruhat cells and total positivity. J. Am. Math. Soc. 1999, 12, 335–380. [Google Scholar] [CrossRef]
  8. Gasca, M.; Micchelli, C.A. Total Positivity and Its Applications; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1996. [Google Scholar]
  9. Lusztig, G. A survey of total positivity. Milan J. Math. 2008, 76, 125–134. [Google Scholar] [CrossRef]
  10. Pinkus, A. Totally Positive Matrices; Cambridge Tracts in Mathematics, 181; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  11. Gantmacher, F.P.; Krein, M.G. Oscillation Matrices and Kernels and Small Vibrations of Mechanical Systems; Revised Edition; AMS Chelsea Publishing: Providence, RI, USA, 2002. [Google Scholar]
  12. Schumaker, L.L. Spline Functions: Basic Theory, 3rd ed.; Cambridge Mathematical Library, Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  13. Karlin, S.; McGregor, J.L. Coincidence probabilities of birth and death processes. Pac. J. Math. 1959, 9, 1109–1140. [Google Scholar] [CrossRef]
  14. Brenti, F. Combinatorics and Total positivity. J. Combin. Theory Ser. A 1995, 71, 175–218. [Google Scholar] [CrossRef]
  15. Brown, L.D.; Johnstone, I.M.; MacGibbon, K.B. Variation diminishing transformations: A direct approach to total positivity and its statistical applications. J. Am. Stat. Assoc. 1981, 76, 824–832. [Google Scholar] [CrossRef]
  16. Goodman, T.N.T.; Micchelli, C.A. Corner cutting algorithms for the Bezier representation of free form curves. Linear Alg. Appl. 1988, 99, 225–252. [Google Scholar] [CrossRef]
  17. Frydman, H.; Singer, B. Total positivity and the embedding problem for Markov chains. Math. Proc. Camb. Philos. Soc. 1979, 85, 339–344. [Google Scholar] [CrossRef]
  18. de Boor, C.; Pinkus, A. Backward error analysis for totally positive linear systems. Numer. Math. 1977, 27, 485–490. [Google Scholar] [CrossRef]
  19. Higham, N.J. Accuracy and Stability of Numerical Algorithms, 2nd ed.; SIAM: Philadelphia, PA, USA, 2002. [Google Scholar]
  20. Fomin, S.; Zelevinsky, A. Total positivity: Tests and parameterizations. Math. Intell. 2000, 22, 23–33. [Google Scholar] [CrossRef]
  21. Koev, P. Accurate computations with totally nonnegative matrices. SIAM J. Matrix Anal. Appl. 2007, 29, 731–751. [Google Scholar] [CrossRef]
  22. Karlin, S. Total Positivity; Stanford University Press: Stanford, CA, USA, 1968; Volume 1. [Google Scholar]
  23. Carnicer, J.M.; Peña, J.M. Totally positive bases for shape preserving curve design and optimality of B-splines. Comput. Aided Geom. Des. 1994, 11, 635–656. [Google Scholar] [CrossRef]
  24. Goodman, T.N.T.; Said, H.B. Shape preserving properties of the generalized Ball basis. Comput. Aided Geom. Des. 1991, 8, 115–121. [Google Scholar] [CrossRef]
  25. Carnicer, J.M.; Peña, J.M. Shape preserving representations and optimality of the Bernstein basis. Adv. Comput. Math. 1993, 1, 173–196. [Google Scholar] [CrossRef]
  26. Demmel, J.; Koev, P. The accurate and efficient solution of a totally positive generalized Vandermonde linear system. SIAM J. Matrix Anal. Appl. 2005, 27, 42–52. [Google Scholar] [CrossRef]
  27. Marco, A.; Martínez, J.J. A fast and accurate algorithm for solving Bernstein-Vandermonde linear systems. Linear Algebra Appl. 2007, 422, 616–628. [Google Scholar] [CrossRef]
  28. Marco, A.; Martínez, J.J. Polynomial least squares fitting in the Bernstein basis. Linear Algebra Appl. 2010, 433, 1254–1264. [Google Scholar] [CrossRef]
  29. Marco, A.; Martínez, J.J. Accurate computations with Said-Ball-Vandermonde matrices. Linear Algebra Appl. 2010, 432, 2894–2908. [Google Scholar] [CrossRef]
  30. Marco, A.; Martínez, J.J. Accurate computations with totally positive Bernstein-Vandermonde matrices. Electron. J. Linear Algebra 2013, 26, 357–380. [Google Scholar] [CrossRef]
  31. Marco, A.; Martínez, J.J.; Viaña, R. Accurate bidiagonal decomposition of totally positive h-Bernstein-Vandermonde matrices and applications. Linear Algebra Appl. 2019, 579, 320–335. [Google Scholar] [CrossRef]
  32. Saad, Y. Iterative Methods for Sparse Linear Systems, 2nd ed.; SIAM: Philadelphia, PA, USA, 2003. [Google Scholar]
  33. Trefethen, L.N.; Bau, D. Numerical Linear Algebra; SIAM: Philadelphia, PA, USA, 1997. [Google Scholar]
  34. Watkins, D.S. Fundamentals of Matrix Computations; Wiley: New York, NY, USA, 2002. [Google Scholar]
  35. Delgado, J.; Peña, J.M. Progressive iterative approximation and bases with the fastest convergence rates. Comput. Aided Geom. Des. 2007, 24, 10–18. [Google Scholar] [CrossRef]
  36. Farin, G. Curves and Surfaces for CAGD, 5th ed.; Academic Press: Boston, MA, USA, 2002. [Google Scholar]
  37. Farouki, R.T.; Goodman, T.N.T. On the optimal stability of Bernstein basis. Math. Comp. 1996, 65, 1553–1566. [Google Scholar] [CrossRef]
  38. Peña, J.M. B-splines and optimal stability. Math. Comp. 1997, 66, 1555–1560. [Google Scholar] [CrossRef]
  39. Peña, J.M. On the optimal stability of bases of univariate functions. Numer. Math. 2002, 91, 305–318. [Google Scholar] [CrossRef]
  40. Carnicer, J.M.; Peña, J.M. Total positivity and optimal bases. In Total Positivity and Its Applications; Gasca, M., Micchelli, C.A., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1996; pp. 133–155. [Google Scholar]
  41. Delgado, J.; Peña, J.M. Extremal and optimal properties of B-bases collocation matrices. Numer. Math. 2020, 146, 105–118. [Google Scholar] [CrossRef]
  42. Koev, P. Accurate eigenvalues and SVDs of totally nonnegative matrices. SIAM J. Matrix Anal. Appl. 2005, 27, 1–23. [Google Scholar] [CrossRef]
  43. Marco, A.; Martínez, J.J. Accurate computation of the Moore–Penrose inverse of strictly totally positive matrices. J. Comput. Appl. Math. 2019, 350, 299–308. [Google Scholar] [CrossRef]
  44. Carnicer, J.M.; Delgado, J.; Peña, J.M. Richardson method and totally nonnegative linear systems. Linear Algebra Appl. 2010, 433, 2010–2017. [Google Scholar] [CrossRef]
  45. de Boor, C. How does Agee’s method work? In Proceedings of the 1979 Army Numerical Analysis and Computers Conference, White Sands Missile Range, NM, USA, 14–16 February 1979; ARO Report 79-3. Army Research Office: Research Triangle Park, NC, USA, 1979; pp. 299–302. [Google Scholar]
  46. Qi, D.; Tian, Z.; Zhang, Y.; Zheng, J.B. The method of numeric polish in curve fitting. Acta Math. Sinica (Chin. Ser.) 1975, 18, 173–184. (In Chinese) [Google Scholar]
  47. Maekawa, T.; Matsumoto, Y.; Namiki, K. Interpolation by geometric algorithm. Comput. Aided Des. 2007, 39, 313–323. [Google Scholar] [CrossRef]
  48. Shi, L.; Wang, R. An iterative algorithm of NURBS interpolation and approximation. J. Math. Res. Expo. 2006, 2, 735–743. [Google Scholar]
  49. Carnicer, J.M.; Delgado, J.; Peña, J.M. Richardson’s iterative method for surface interpolation. BIT Numer. Math. 2013, 53, 385–396. [Google Scholar] [CrossRef]
  50. Horn, R.A.; Johnson, C.R. Topics in Matrix Analysis; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  51. Delgado, J.; Orera, H.; Peña, J.M. Optimal properties of tensor product of B-bases. Appl. Math. Lett. 2021, 121, 107473. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peña, J.M. An Optimal Property of B-Bases for the Modified Richardson Method. Axioms 2024, 13, 223. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13040223

AMA Style

Peña JM. An Optimal Property of B-Bases for the Modified Richardson Method. Axioms. 2024; 13(4):223. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13040223

Chicago/Turabian Style

Peña, Juan Manuel. 2024. "An Optimal Property of B-Bases for the Modified Richardson Method" Axioms 13, no. 4: 223. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms13040223

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