Results 1  10
of
83
On the multilevel splitting of finite element spaces
 Numer. Math
, 1986
"... Federal Republic of Germany ..."
Stability of Multiscale Transformations
 J. Fourier Anal. Appl
, 1996
"... After briefly reviewing the interrelation between Rieszbases, biorthogonality and a certain stability notion for multiscale basis transformations we establish a basic stability criterion for a general Hilbert space setting. An important tool in this context is a strengthened Cauchy inequality. It i ..."
Abstract

Cited by 92 (22 self)
 Add to MetaCart
After briefly reviewing the interrelation between Rieszbases, biorthogonality and a certain stability notion for multiscale basis transformations we establish a basic stability criterion for a general Hilbert space setting. An important tool in this context is a strengthened Cauchy inequality. It is based on direct and inverse estimates for a certain scale of spaces induced by the underlying multiresolution sequence. Furthermore, we highlight some properties of these spaces pertaining to duality, interpolation, and applications to norm equivalences for Sobolev spaces. AMS Subject Classification: 41A17, 41A65, 46A35, 46B70, 46E35 Key Words: Riesz bases, biorthogonality, stability, projectors, interpolation theory, Kmethod, duality, Jackson, Bernstein inequalities 1 Background and Motivation A standard framework for approximately recapturing a function v in some infinite dimensional separable Hilbert space V , say, either from explicitly given data or as a solution of an operator equ...
Hierarchical Bases and the Finite Element Method
, 1997
"... CONTENTS 1 Introduction 1 2 Preliminaries 3 3 Fundamental TwoLevel Estimates 7 4 A Posteriori Error Estimates 16 5 TwoLevel Iterative Methods 23 6 Multilevel Cauchy Inequalities 30 7 Multilevel Iterative Methods 34 References 41 1. Introduction In this work we present a brief introduction to hie ..."
Abstract

Cited by 68 (4 self)
 Add to MetaCart
CONTENTS 1 Introduction 1 2 Preliminaries 3 3 Fundamental TwoLevel Estimates 7 4 A Posteriori Error Estimates 16 5 TwoLevel Iterative Methods 23 6 Multilevel Cauchy Inequalities 30 7 Multilevel Iterative Methods 34 References 41 1. Introduction In this work we present a brief introduction to hierarchical bases, and the important part they play in contemporary finite element calculations. In particular, we examine their role in a posteriori error estimation, and in the Department of Mathematics, University of California at San Diego, La Jolla, CA 92093. The work of this author was supported by the Office of Naval Research under contract N0001489J1440. 2 Randolph E. Bank formulation of iterative methods for solving the large sparse sets of linear equations arising from the finite element discretization. Our goal is that the development should be largely selfcontained, but at the same time accessible and interest
Overlapping Schwarz Methods On Unstructured Meshes Using NonMatching Coarse Grids
 Numer. Math
, 1996
"... . We consider two level overlapping Schwarz domain decomposition methods for solving the finite element problems that arise from discretizations of elliptic problems on general unstructured meshes in two and three dimensions. Standard finite element interpolation from the coarse to the fine grid may ..."
Abstract

Cited by 52 (17 self)
 Add to MetaCart
. We consider two level overlapping Schwarz domain decomposition methods for solving the finite element problems that arise from discretizations of elliptic problems on general unstructured meshes in two and three dimensions. Standard finite element interpolation from the coarse to the fine grid may be used. Our theory requires no assumption on the substructures which constitute the whole domain, so each substructure can be of arbitrary shape and of different size. The global coarse mesh is allowed to be nonnested to the fine grid on which the discrete problem is to be solved and both the coarse meshes and the fine meshes need not be quasiuniform. In addition, the domains defined by the fine and coarse grid need not be identical. The one important constraint is that the closure of the coarse grid must cover any portion of the fine grid boundary for which Neumann boundary conditions are given. In this general setting, our algorithms have the same optimal convergence rate of the usual ...
Domain decomposition algorithms for indefinite elliptic problems
 SIAM J. Sci. Stat. Comput
, 1992
"... Iterative methods for linear systems of algebraic equations arising from the finite element discretization of nonsymmetric and indefinite elliptic problems are considered. Methods previously known to work well for positive definite, symmetric problems are extended to certain nonsymmetric problems, w ..."
Abstract

Cited by 51 (16 self)
 Add to MetaCart
Iterative methods for linear systems of algebraic equations arising from the finite element discretization of nonsymmetric and indefinite elliptic problems are considered. Methods previously known to work well for positive definite, symmetric problems are extended to certain nonsymmetric problems, which can also have some eigenvalues in the left half plane. We first consider an additive Schwarz method applied to linear, second order, symmetric or nonsymmetric, indefinite elliptic boundary value problems in two and three dimensions. An alternative linear system, which has the same solution as the original problem, is derived and this system is then solved by using GMRES, an iterative method of conjugate gradient type. In each iteration step, a coarse mesh finite element problem and a number of local problems are solved on small, overlapping subregions into which the original region is subdivided. We show that the rate of convergence is independent of the number of degrees of freedom and the number of local problems if the coarse mesh is fine enough. The performance of the method in two dimensions is illustrated by results of several numerical experiments. We also consider two other iterative method for solving the same class of elliptic problems in two dimensions. Using an observation of Dryja and Widlund, we show that the rate of convergence of certain iterative substructuring methods deteriorates only quite slowly when the local problems increase in size. A similar result is established for Yserentant’s hierarchical basis method.
An Algorithm for Coarsening Unstructured Meshes
 Numer. Math
, 1996
"... . We develop and analyze a procedure for creating a hierarchical basis of continuous piecewise linear polynomials on an arbitrary, unstructured, nonuniform triangular mesh. Using these hierarchical basis functions, we are able to define and analyze corresponding iterative methods for solving the lin ..."
Abstract

Cited by 49 (5 self)
 Add to MetaCart
. We develop and analyze a procedure for creating a hierarchical basis of continuous piecewise linear polynomials on an arbitrary, unstructured, nonuniform triangular mesh. Using these hierarchical basis functions, we are able to define and analyze corresponding iterative methods for solving the linear systems arising from finite element discretizations of elliptic partial differential equations. We show that such iterative methods perform as well as those developed for the usual case of structured, locally refined meshes. In particular, we show that the generalized condition numbers for such iterative methods are of order J 2 , where J is the number of hierarchical basis levels. Key words. Finite element, hierarchical basis, multigrid, unstructured mesh. AMS subject classifications. 65F10, 65N20 1. Introduction. Iterative methods using the hierarchical basis decomposition have proved to be among the most robust for solving broad classes of elliptic partial differential equations, ...
Some Nonoverlapping Domain Decomposition Methods
, 1998
"... . The purpose of this paper is to give a unified investigation of a class of nonoverlapping domain decomposition methods for solving secondorder elliptic problems in two and three dimensions. The methods under scrutiny fall into two major categories: the substructuringtype methods and the Neumann ..."
Abstract

Cited by 47 (8 self)
 Add to MetaCart
. The purpose of this paper is to give a unified investigation of a class of nonoverlapping domain decomposition methods for solving secondorder elliptic problems in two and three dimensions. The methods under scrutiny fall into two major categories: the substructuringtype methods and the NeumannNeumanntype methods. The basic framework used for analysis is the parallel subspace correction method or additive Schwarz method, and other technical tools include localglobal and globallocal techniques. The analyses for both two and threedimensional cases are carried out simultaneously. Some internal relationships between various algorithms are observed and several new variants of the algorithms are also derived. Key words. nonoverlapping domain decomposition, Schur complement, localglobal and globallocal techniques, jumps in coe#cients, substructuring, NeumannNeumann, balancing methods AMS subject classifications. 65N30, 65N55, 65F10 PII. S0036144596306800 1. Introduction. T...
A nonoverlapping domain decomposition method for Maxwell’s equations in three dimensions
 SIAM J. Numer. Anal
"... Abstract. We propose a substructuring preconditioner for solving threedimensional elliptic equations with strongly discontinuous coefficients. The new preconditioner can be viewed as a variant of the classical substructuring preconditioner proposed by Bramble, Pasiack and Schatz (1989), but with muc ..."
Abstract

Cited by 44 (16 self)
 Add to MetaCart
Abstract. We propose a substructuring preconditioner for solving threedimensional elliptic equations with strongly discontinuous coefficients. The new preconditioner can be viewed as a variant of the classical substructuring preconditioner proposed by Bramble, Pasiack and Schatz (1989), but with much simpler coarse solvers. Though the condition number of the preconditioned system may not have a good bound, we are able to show that the convergence rate of the PCG method with such substructuring preconditioner is nearly optimal, and also robust with respect to the (possibly large) jumps of the coefficient in the elliptic equation. 1.
An EnergyMinimizing Interpolation For Robust Multigrid Methods
 SIAM J. SCI. COMPUT
, 1998
"... We propose a robust interpolation for multigrid based on the concepts of energy minimization and approximation. The formulation is general; it can be applied to any dimensions. The analysis for one dimension proves that the convergence rate of the resulting multigrid method is independent of the coe ..."
Abstract

Cited by 43 (6 self)
 Add to MetaCart
We propose a robust interpolation for multigrid based on the concepts of energy minimization and approximation. The formulation is general; it can be applied to any dimensions. The analysis for one dimension proves that the convergence rate of the resulting multigrid method is independent of the coefficient of the underlying PDE, in addition to being independent of the mesh size. We demonstrate numerically the effectiveness of the multigrid method in two dimensions by applying it to a discontinuous coefficient problem and an oscillatory coefficient problem. We also show using a onedimensional Helmholtz problem that the energy minimization principle can be applied to solving elliptic problems that are not positive definite.
The method of alternating projections and the method of subspace corrections in Hilbert space
 J. Am. Math. Soc
"... Abstract. A new identity is given in this paper for estimating the norm of the product of nonexpansive operators in Hilbert space. This identity can be applied for the design and analysis of the method of alternating projections and the method of subspace corrections. The method of alternating proje ..."
Abstract

Cited by 42 (7 self)
 Add to MetaCart
Abstract. A new identity is given in this paper for estimating the norm of the product of nonexpansive operators in Hilbert space. This identity can be applied for the design and analysis of the method of alternating projections and the method of subspace corrections. The method of alternating projections is an iterative algorithm for determining the best approximation to any given point in a Hilbert space from the intersection of a nite number of subspaces by alternatively computing the best approximations from the individual subspaces which make up the intersection. The method of subspace corrections is an iterative algorithm for nding the solution of a linear equation in a Hilbert space by approximately solving equations restricted on a number of closed subspaces which make up the entire space. The new identity given in the paper provides a sharpest possible estimate for the rate of convergence of these algorithms. It is also proved in the paper that the method of alternating projections is essentially equivalent to the method of subspace corrections. Some simple examples of multigrid and domain decomposition methods are given to illustrate the application of the new identity. 1.