In the present paper, we consider a preconditioning strategy for Finite Element (FE) matrix sequences {A_(n) (a)}_(n) discretizing the elliptic problem [GRAPHICS] with a(x, y) being a uniformly positive function and nu denoting the unit outward normal direction. More precisely, in connection with preconditioned conjugate gradient (PCG) like methods, we define the preconditioning sequence: {P_n(a)}_(n), P_n(a) := (D) over tilde ^(1/2)_(n) (a)A_(n)(1)(D) over tilde ^(1/2)_(n)(a), where (D) over tilde _(n)(a) is the suitable scaled main diagonal of A_(n)(a). In fact, under the mild assumption of Lebesgue integrability of a(x), the weak clustering at the unity of the corresponding preconditioned sequence is proved. Moreover, if a(x, y) is regular enough and if a uniform triangulation is considered, then the preconditioned sequence shows a strong clustering at the unity so that the sequence {P_n(a)}_(n) turns out to be a superlinear preconditioning sequence for {A_(n)(a)}_(n). The computational interest is due to the fact that the computation with A_(n)(a) is reduced to computations involving diagonals and two- level Toeplitz structures {A_(n)(1)}_(n) with banded pattern. Some numerical experimentations confirm the efficiency of the discussed proposal
Serra Capizzano, S., Tablino Possio, C. (2001). Finite element matrix sequences: the case of rectangular domains. NUMERICAL ALGORITHMS, 28(1-4), 309-327 [10.1023/A:1014088011253].
Finite element matrix sequences: the case of rectangular domains
Tablino Possio, C.
2001
Abstract
In the present paper, we consider a preconditioning strategy for Finite Element (FE) matrix sequences {A_(n) (a)}_(n) discretizing the elliptic problem [GRAPHICS] with a(x, y) being a uniformly positive function and nu denoting the unit outward normal direction. More precisely, in connection with preconditioned conjugate gradient (PCG) like methods, we define the preconditioning sequence: {P_n(a)}_(n), P_n(a) := (D) over tilde ^(1/2)_(n) (a)A_(n)(1)(D) over tilde ^(1/2)_(n)(a), where (D) over tilde _(n)(a) is the suitable scaled main diagonal of A_(n)(a). In fact, under the mild assumption of Lebesgue integrability of a(x), the weak clustering at the unity of the corresponding preconditioned sequence is proved. Moreover, if a(x, y) is regular enough and if a uniform triangulation is considered, then the preconditioned sequence shows a strong clustering at the unity so that the sequence {P_n(a)}_(n) turns out to be a superlinear preconditioning sequence for {A_(n)(a)}_(n). The computational interest is due to the fact that the computation with A_(n)(a) is reduced to computations involving diagonals and two- level Toeplitz structures {A_(n)(1)}_(n) with banded pattern. Some numerical experimentations confirm the efficiency of the discussed proposalI documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.