Algorithms for Large Sparse Matrices Based on Finite Element Method

Many problems arising from engeneering and scientific computing give rise to large, sparce matrices. The aim of this work is to describe a set of highly efficient iterative methods for solving linear systems with large sparse matrices, arising from the analysis of seismic body wave propagation. An "ad hoc" initial boundary value problem is formulated for heterogeneous dissipative media with arbitrary topography. Its numerical implementation is based on Finite Ele- ment Method on non structured mesh. Some results are presented.