Domain Embedding Methods for the Stokes Equations.
We consider mixed finite element discretizations of second order elliptic boundary value problems. Emphasis is on the efficient iterative solution by multilevel techniques with respect to an adaptively generated hierarchy of nonuniform triangulations. In particular, we present two multilevel solvers, the first one relying on ideas from domain decomposition and the second one resulting from mixed hybridization. Local refinement of the underlying triangulations is done by efficient and reliable a...
Es wird ein kombinierter Algorithmus zur iterativen Einschlissung der Inversen einer Matrix beschrieben. Es handelt sich dabei um eine intervallmässige Version des Schulz'schen Verfahrens. Es wird bewiesen, dass der Algorithmus genauso effizient ist wie ein hisher bekannter aus [2], dass er aber in Bezug auf den akkumulierten Rundungsfehler dem bisherigen Vorgehen vorzuziehen ist. Ein numerisches Beispiel wird gegeben.
For contractive interval functions we show that results from the iterative process after finitely many iterations if one uses the epsilon-inflated vector as input for instead of the original output vector . Applying Brouwer’s fixed point theorem, zeros of various mathematical problems can be verified in this way.
This paper derives upper and lower bounds for the -condition number of the stiffness matrix resulting from the finite element approximation of a linear, abstract model problem. Sharp estimates in terms of the meshsize h are obtained. The theoretical results are applied to finite element approximations of elliptic PDE's in variational and in mixed form, and to first-order PDE's approximated using the Galerkin–Least Squares technique or by means of a non-standard Galerkin technique in L1(Ω). Numerical...
In this note, we compare some Krylov subspace iterative methods on the MASPAR, a massively parallel computer with 16K processors. In particular, we apply these methods to solve large sparse nonsymmetric linear systems arising from elliptic partial differential equations. The methods under consideration include conjugate gradient type methods, semiiterative methods, and a hybrid variant. Our numerical results show that, on the MASPAR, one should compare iterative methods rather on the basis of total...
An algorithm of the preconditioned conjugate gradient method in which the solution of an auxiliary system is replaced with multiplication by the matrix for suitably chosen is presented.