Stefano Serra Capizzano

2papers

2 Papers

NASep 11, 2012
Optimal Preconditioners for Finite Element Approximations of Convection-Diffusion Equations on structured meshes

Alessandro Russo, Stefano Serra Capizzano, Cristina Tablino Possio

The paper is devoted to the spectral analysis of effective preconditioners for linear systems obtained via a Finite Element approximation to diffusion-dominated convection-diffusion equations. We consider a model setting in which the structured finite element partition is made by equi-lateral triangles. Under such assumptions, if the problem is coercive, and the diffusive and convective coefficients are regular enough, then the proposed preconditioned matrix sequences exhibit a strong clustering at unity, the preconditioning matrix sequence and the original matrix sequence are spectrally equivalent, and the eigenvector matrices have a mild conditioning. The obtained results allow to show the optimality of the related preconditioned Krylov methods. %It is important to stress that The interest of such a study relies on the observation that automatic grid generators tend to construct equi-lateral triangles when the mesh is fine enough. Numerical tests, both on the model setting and in the non-structured case, show the effectiveness of the proposal and the correctness of the theoretical findings.

NADec 20, 2005
The spectral approximation of multiplication operators via asymptotic (structured) linear algebra

Stefano Serra Capizzano

multiplication operator on a Hilbert space may be approximated with finite sections by choosing an orthonormal basis of the Hilbert space. Nonzero multiplication operators on $L^2$ spaces of functions are never compact and then such approximations cannot converge in the norm topology. Instead, we consider how well the spectra of the finite sections approximate the spectrum of the multiplication operator whose expression is simply given by the essential range of the symbol (i.e. the multiplier). We discuss the case of real orthogonal polynomial bases and the relations with the classical Fourier basis whose choice leads to well studied Toeplitz case. The use of circulant approximations leads to constructive algorithms working for the separable multivariate and matrix-valued cases as well.