Results 1  10
of
1,544
An optimal control approach to a posteriori error estimation in finite element methods
, 2001
"... ..."
Surface Parameterization: a Tutorial and Survey
 In Advances in Multiresolution for Geometric Modelling, Mathematics and Visualization
, 2005
"... Summary. This paper provides a tutorial and survey of methods for parameterizing surfaces with a view to applications in geometric modelling and computer graphics. We gather various concepts from differential geometry which are relevant to surface mapping and use them to understand the strengths and ..."
Abstract

Cited by 239 (7 self)
 Add to MetaCart
(Show Context)
Summary. This paper provides a tutorial and survey of methods for parameterizing surfaces with a view to applications in geometric modelling and computer graphics. We gather various concepts from differential geometry which are relevant to surface mapping and use them to understand the strengths and weaknesses of the many methods for parameterizing piecewise linear surfaces and their relationship to one another. 1
GALERKIN FINITE ELEMENT APPROXIMATIONS OF STOCHASTIC ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS
, 2004
"... We describe and analyze two numerical methods for a linear elliptic problem with stochastic coefficients and homogeneous Dirichlet boundary conditions. Here the aim of the computations is to approximate statistical moments of the solution, and, in particular, we give a priori error estimates for the ..."
Abstract

Cited by 193 (11 self)
 Add to MetaCart
We describe and analyze two numerical methods for a linear elliptic problem with stochastic coefficients and homogeneous Dirichlet boundary conditions. Here the aim of the computations is to approximate statistical moments of the solution, and, in particular, we give a priori error estimates for the computation of the expected value of the solution. The first method generates independent identically distributed approximations of the solution by sampling the coefficients of the equation and using a standard Galerkin finite element variational formulation. The Monte Carlo method then uses these approximations to compute corresponding sample averages. The second method is based on a finite dimensional approximation of the stochastic coefficients, turning the original stochastic problem into a deterministic parametric elliptic problem. A Galerkin finite element method, of either the h or pversion, then approximates the corresponding deterministic solution, yielding approximations of the desired statistics. We present a priori error estimates and include a comparison of the computational work required by each numerical approximation to achieve a given accuracy. This comparison suggests intuitive conditions for an optimal selection of the numerical approximation.
Finite element exterior calculus, homological techniques, and applications
 ACTA NUMERICA
, 2006
"... ..."
Adaptive wavelet methods for elliptic operator equations— convergence rates
 Math. Comput
, 2001
"... Abstract. This paper is concerned with the construction and analysis of waveletbased adaptive algorithms for the numerical solution of elliptic equations. These algorithms approximate the solution u of the equation by a linear combination of N wavelets. Therefore, a benchmark for their performance ..."
Abstract

Cited by 171 (33 self)
 Add to MetaCart
(Show Context)
Abstract. This paper is concerned with the construction and analysis of waveletbased adaptive algorithms for the numerical solution of elliptic equations. These algorithms approximate the solution u of the equation by a linear combination of N wavelets. Therefore, a benchmark for their performance is provided by the rate of best approximation to u by an arbitrary linear combination of N wavelets (so called Nterm approximation), which would be obtained by keeping the N largest wavelet coefficients of the real solution (which of course is unknown). The main result of the paper is the construction of an adaptive scheme which produces an approximation to u with error O(N −s)in the energy norm, whenever such a rate is possible by Nterm approximation. The range of s>0 for which this holds is only limited by the approximation properties of the wavelets together with their ability to compress the elliptic operator. Moreover, it is shown that the number of arithmetic operations needed to compute the approximate solution stays proportional to N. The adaptive algorithm applies to a wide class of elliptic problems and wavelet bases. The analysis in this paper puts forward new techniques for treating elliptic problems as well as the linear systems of equations that arise from the wavelet discretization. 1.
A FeedBack Approach to Error Control in Finite Element Methods: Basic Analysis and Examples
 EastWest J. Numer. Math
, 1996
"... this paper. ..."
A restricted additive Schwarz preconditioner for general sparse linear systems
 SIAM J. Sci. Comput
, 1999
"... Abstract. We introduce some cheaper and faster variants of the classical additive Schwarz preconditioner (AS) for general sparse linear systems and show, by numerical examples, that the new methods are superior to AS in terms of both iteration counts and CPU time, as well as the communication cost w ..."
Abstract

Cited by 129 (24 self)
 Add to MetaCart
(Show Context)
Abstract. We introduce some cheaper and faster variants of the classical additive Schwarz preconditioner (AS) for general sparse linear systems and show, by numerical examples, that the new methods are superior to AS in terms of both iteration counts and CPU time, as well as the communication cost when implemented on distributed memory computers. This is especially true for harder problems such as indefinite complex linear systems and systems of convectiondiffusion equations from threedimensional compressible flows. Both sequential and parallel results are reported. Key words. Overlapping domain decomposition, preconditioner, iterative method, sparse matrix AMS(MOS) subject classifications. 65N30, 65F10
Adaptive Discontinuous Galerkin Finite Element Methods for Compressible Fluid Flows
 SIAM J. Sci. Comput
"... this paper is to discuss the a posteriori error analysis and adaptive mesh design for discontinuous Galerkin finite element approximations to systems of conservation laws. In Section 2, we introduce the model problem and formulate its discontinuous Galerkin finite element approximation. Section 3 is ..."
Abstract

Cited by 122 (17 self)
 Add to MetaCart
this paper is to discuss the a posteriori error analysis and adaptive mesh design for discontinuous Galerkin finite element approximations to systems of conservation laws. In Section 2, we introduce the model problem and formulate its discontinuous Galerkin finite element approximation. Section 3 is devoted to the derivation of weighted a posteriori error bounds for linear functionals of the solution. Finally, in Section 4 we present some numerical examples to demonstrate the performance of the resulting adaptive finite element algorithm. 2 Model problem and discretisation Given an open bounded polyhedral domain fl in lI n, n _> 1, with boundary 0fl, we consider the following problem: find u: fl > lI m, m _> 1, such that div(u) = 0 in , (2.1) where, ,: m __> mxn is continuously differentiable. We assume that the system of conservation laws (2.1) may be supplemented by appropriate initial/boundary conditions. For example, assuming that B(u, y) := EiL1 biVu'(u) has m real eigenvalues and a complete set of linearly independent eigenvectors for all y = (yl,, Yn) C n; then at inflow/outflow boundaries, we require that B(u, n)(u g) = 0, where n denotes the unit outward normal vector to 0fl, B(u, n) is the negative part of B(u, n) and g is a (given) realvalued function. To formulate the discontinuous Galerkin finite element method (DGFEM, for short) for (2.1), we first introduce some notation. Let 7 = {n} be an admissible subdivision of fl into open element domains n; here h is a piecewise constant mesh function with h(x) = diam(n) 2 Houston e al. when x is in element n. For p Iq0, we define the following finite element space n,  {v [L()]": vl [%(n)] " Vn }, where Pp(n) denotes the set of polynomials of degree at most p over n. Given that v [Hi(n)] m for each n...
An Explicit Link between Gaussian Fields and Gaussian Markov random fields: the SPDE approach
 PREPRINTS IN MATHEMATICAL SCIENCES
, 2010
"... Continuously indexed Gaussian fields (GFs) is the most important ingredient in spatial statistical modelling and geostatistics. The specification through the covariance function gives an intuitive interpretation of its properties. On the computational side, GFs are hampered with the bign problem, ..."
Abstract

Cited by 115 (17 self)
 Add to MetaCart
Continuously indexed Gaussian fields (GFs) is the most important ingredient in spatial statistical modelling and geostatistics. The specification through the covariance function gives an intuitive interpretation of its properties. On the computational side, GFs are hampered with the bign problem, since the cost of factorising dense matrices is cubic in the dimension. Although the computational power today is alltimehigh, this fact seems still to be a computational bottleneck in applications. Along with GFs, there is the class of Gaussian Markov random fields (GMRFs) which are discretely indexed. The Markov property makes the involved precision matrix sparse which enables the use of numerical algorithms for sparse matrices, that for fields in R 2 only use the squareroot of the time required by general algorithms. The specification of a GMRF is through its full conditional distributions but its marginal properties are not transparent in such a parametrisation. In this paper, we show that using an approximate stochastic weak solution to (linear) stochastic partial differential equations (SPDEs), we can, for some GFs in the Matérn class, provide an explicit link, for any triangulation of R d, between GFs and GMRFs. The consequence is that we can take the best from the two worlds and do the modelling using GFs but do the computations using GMRFs. Perhaps more importantly,