Debora Sesana

Marco Donatelli, Stefano Serra-Capizzano, Debora Sesana

Starting from the spectral analysis of g-circulant matrices, we consider a new multigrid method for circulant and Toeplitz matrices with given generating function. We assume that the size n of the coefficient matrix is divisible by g \geq 2 such that at the lower level the system is reduced to one of size n/g by employing g-circulant based projectors. We perform a rigorous two-grid convergence analysis in the circulant case and we extend experimentally the results to the Toeplitz setting, by employing structure preserving projectors. The optimality of the proposed two-grid method and of the multigrid method is proved, when the number theta \in N of recursive calls is such that 1 < theta < g. The previous analysis is used also to overcome some pathological cases, in which the generating function has zeros located at "mirror points" and the standard two-grid method with g = 2 is not optimal. The numerical experiments show the correctness and applicability of the proposed ideas both for circulant and Toeplitz matrices.

Ernesto Salinelli, Debora Sesana

In this paper we give complete results on the presence of Shift, Slope and Curvature for a correlation model of interest rates, by improving and extending the content of a previous paper on the subject. We get our goal essentially exploiting some properties of Green's matrices and the notion of convexity for eigenvectors.

Eric Ngondiep, Stefano Serra-Capizzano, Debora Sesana

For a given nonnegative integer alpha, a matrix A_{n} of size n is called alpha-Toeplitz if its entries obey the rule A_{n}=[a_{r-alpha*s}]_{r,s=0}^{n-1}. Analogously, a matrix A_{n} again of size n is called alpha-circulant if A_{n}= [a_{(r-alpha*s)mod n}]_{r,s=0}^{n-1}. Such kind of matrices arises in wavelet analysis, subdivision algorithms and more generally when dealing with multigrid/multilevel methods for structured matrices and approximations of boundary value problems. In this paper we study the singular values of alpha-circulants and we provide an asymptotic analysis of the distribution results for the singular values of alpha-Toeplitz sequences in the case where {a_{k}} can be interpreted as the sequence of Fourier coeffcients of an integrable function f over the domain (-pi;pi). Some generalizations to the block, multilevel case, amounting to choose f multivariate and matrix valued, are briefly considered.