115 research outputs found
Block-adaptive Cross Approximation of Discrete Integral Operators
In this article we extend the adaptive cross approximation (ACA) method known
for the efficient approximation of discretisations of integral operators to a
block-adaptive version. While ACA is usually employed to assemble hierarchical
matrix approximations having the same prescribed accuracy on all blocks of the
partition, for the solution of linear systems it may be more efficient to adapt
the accuracy of each block to the actual error of the solution as some blocks
may be more important for the solution error than others. To this end, error
estimation techniques known from adaptive mesh refinement are applied to
automatically improve the block-wise matrix approximation. This allows to
interlace the assembling of the coefficient matrix with the iterative solution
Numerics of boundary-domain integral and integro-differential equations for BVP with variable coefficient in 3D
This is the post-print version of the article. The official published version can be accessed from the links below - Copyright @ 2013 Springer-VerlagA numerical implementation of the direct boundary-domain integral and integro-differential equations, BDIDEs, for treatment of the Dirichlet problem for a scalar elliptic PDE with variable coefficient in a three-dimensional domain is discussed. The mesh-based discretisation of the BDIEs with tetrahedron domain elements in conjunction with collocation method leads to a system of linear algebraic equations (discretised BDIE). The involved fully populated matrices are approximated by means of the H-Matrix/adaptive cross approximation technique. Convergence of the method is investigated.This study is partially supported by the EPSRC grant EP/H020497/1:"Mathematical Analysis of Localised-Boundary-Domain Integral Equations for Variable-Coefficients
Boundary Value Problems"
Comparison of some Reduced Representation Approximations
In the field of numerical approximation, specialists considering highly
complex problems have recently proposed various ways to simplify their
underlying problems. In this field, depending on the problem they were tackling
and the community that are at work, different approaches have been developed
with some success and have even gained some maturity, the applications can now
be applied to information analysis or for numerical simulation of PDE's. At
this point, a crossed analysis and effort for understanding the similarities
and the differences between these approaches that found their starting points
in different backgrounds is of interest. It is the purpose of this paper to
contribute to this effort by comparing some constructive reduced
representations of complex functions. We present here in full details the
Adaptive Cross Approximation (ACA) and the Empirical Interpolation Method (EIM)
together with other approaches that enter in the same category
Complexity Analysis of a Fast Directional Matrix-Vector Multiplication
We consider a fast, data-sparse directional method to realize matrix-vector
products related to point evaluations of the Helmholtz kernel. The method is
based on a hierarchical partitioning of the point sets and the matrix. The
considered directional multi-level approximation of the Helmholtz kernel can be
applied even on high-frequency levels efficiently. We provide a detailed
analysis of the almost linear asymptotic complexity of the presented method.
Our numerical experiments are in good agreement with the provided theory.Comment: 20 pages, 2 figures, 1 tabl
On thin plate spline interpolation
We present a simple, PDE-based proof of the result [M. Johnson, 2001] that
the error estimates of [J. Duchon, 1978] for thin plate spline interpolation
can be improved by . We illustrate that -matrix
techniques can successfully be employed to solve very large thin plate spline
interpolation problem
Application of hierarchical matrices for computing the Karhunen-Loève expansion
Realistic mathematical models of physical processes contain uncertainties. These models are often described by stochastic differential equations (SDEs) or stochastic partial differential equations (SPDEs) with multiplicative noise. The uncertainties in the right-hand side or the coefficients are represented as random fields. To solve a given SPDE numerically one has to discretise the deterministic operator as well as the stochastic fields. The total dimension of the SPDE is the product of the dimensions of the deterministic part and the stochastic part. To approximate random fields with as few random variables as possible, but still retaining the essential information, the Karhunen-Lo`eve expansion (KLE) becomes important. The KLE of a random field requires the solution of a large eigenvalue problem. Usually it is solved by a Krylov subspace method with a sparse matrix approximation. We demonstrate the use of sparse hierarchical matrix techniques for this. A log-linear computational cost of the matrix-vector product and a log-linear storage requirement yield an efficient and fast discretisation of the random fields presented
Wavelet boundary element methods – Adaptivity and goal-oriented error estimation
This article is dedicated to the adaptive wavelet boundary element method. It computes an approximation to the unknown solution of the boundary integral equation under consideration with a rate , whenever the solution can be approximated with this rate in the setting determined by the underlying wavelet basis. The computational cost scale linearly in the number of degrees of freedom. Goal-oriented error estimation for evaluating linear output functionals of the solution is also considered. An algorithm is proposed that approximately evaluates a linear output functional with a rate , whenever the primal solution can be approximated with a rate and the dual solution can be approximated with a rate , while the cost still scale linearly in . Numerical results for an acoustic scattering problem and for the point evaluation of the potential in case of the Laplace equation are reported to validate and quantify the approach
- …