Search results for "multigrid"
showing 10 items of 30 documents
Comparison between the shifted-Laplacian preconditioning and the controllability methods for computational acoustics
2010
Processes that can be modelled with numerical calculations of acoustic pressure fields include medical and industrial ultrasound, echo sounding, and environmental noise. We present two methods for making these calculations based on Helmholtz equation. The first method is based directly on the complex-valued Helmholtz equation and an algebraic multigrid approximation of the discretized shifted-Laplacian operator; i.e. the damped Helmholtz operator as a preconditioner. The second approach returns to a transient wave equation, and finds the time-periodic solution using a controllability technique. We concentrate on acoustic problems, but our methods can be used for other types of Helmholtz pro…
An algebraic multigrid based shifted-Laplacian preconditioner for the Helmholtz equation
2007
A preconditioner defined by an algebraic multigrid cycle for a damped Helmholtz operator is proposed for the Helmholtz equation. This approach is well suited for acoustic scattering problems in complicated computational domains and with varying material properties. The spectral properties of the preconditioned systems and the convergence of the GMRES method are studied with linear, quadratic, and cubic finite element discretizations. Numerical experiments are performed with two-dimensional problems describing acoustic scattering in a cross-section of a car cabin and in a layered medium. Asymptotically the number of iterations grows linearly with respect to the frequency while for lower freq…
A damping preconditioner for time-harmonic wave equations in fluid and elastic material
2009
A physical damping is considered as a preconditioning technique for acoustic and elastic wave scattering. The earlier preconditioners for the Helmholtz equation are generalized for elastic materials and three-dimensional domains. An algebraic multigrid method is used in approximating the inverse of damped operators. Several numerical experiments demonstrate the behavior of the method in complicated two-dimensional and three-dimensional domains. peerReviewed
Computation of a few smallest eigenvalues of elliptic operators using fast elliptic solvers
2001
The computation of a few smallest eigenvalues of generalized algebraic eigenvalue problems is studied. The considered problems are obtained by discretizing self-adjoint second-order elliptic partial differential eigenvalue problems in two- or three-dimensional domains. The standard Lanczos algorithm with the complete orthogonalization is used to compute some eigenvalues of the inverted eigenvalue problem. Under suitable assumptions, the number of Lanczos iterations is shown to be independent of the problem size. The arising linear problems are solved using some standard fast elliptic solver. Numerical experiments demonstrate that the inverted problem is much easier to solve with the Lanczos…
Qualitative Analysis of Differential, Difference Equations, and Dynamic Equations on Time Scales
2015
and Applied Analysis 3 thank Guest Editors Josef Dibĺik, Alexander Domoshnitsky, Yuriy V. Rogovchenko, Felix Sadyrbaev, and Qi-Ru Wang for their unfailing support with editorial work that ensured timely preparation of this special edition. Tongxing Li Josef Dibĺik Alexander Domoshnitsky Yuriy V. Rogovchenko Felix Sadyrbaev Qi-Ru Wang
(Approximate) Low-Mode Averaging with a new Multigrid Eigensolver
2015
We present a multigrid based eigensolver for computing low-modes of the Hermitian Wilson Dirac operator. For the non-Hermitian case multigrid methods have already replaced conventional Krylov subspace solvers in many lattice QCD computations. Since the $\gamma_5$-preserving aggregation based interpolation used in our multigrid method is valid for both, the Hermitian and the non-Hermitian case, inversions of very ill-conditioned shifted systems with the Hermitian operator become feasible. This enables the use of multigrid within shift-and-invert type eigensolvers. We show numerical results from our MPI-C implementation of a Rayleigh quotient iteration with multigrid. For state-of-the-art lat…
Quasi-continuous-time impurity solver for the dynamical mean-field theory with linear scaling in the inverse temperature
2013
We present an algorithm for solving the self-consistency equations of the dynamical mean-field theory (DMFT) with high precision and efficiency at low temperatures. In each DMFT iteration, the impurity problem is mapped to an auxiliary Hamiltonian, for which the Green function is computed by combining determinantal quantum Monte Carlo (BSS-QMC) calculations with a multigrid extrapolation procedure. The method is numerically exact, i.e., yields results which are free of significant Trotter errors, but retains the BSS advantage, compared to direct QMC impurity solvers, of linear (instead of cubic) scaling with the inverse temperature. The new algorithm is applied to the half-filled Hubbard mo…
Mott transitions in ternary flavor mixtures of ultracold fermions on optical lattices
2009
Ternary flavor mixtures of ultracold fermionic atoms in an optical lattice are studied in the case of equal, repulsive on-site interactions U>0. The corresponding SU(3) invariant Hubbard model is solved numerically exactly within dynamical mean-field theory using multigrid Hirsch-Fye quantum Monte Carlo simulations. We establish Mott transitions close to integer filling at low temperatures and show that the associated signatures in the compressibility and pair occupancy persist to high temperatures, i.e., should be accessible to experiments. In addition, we present spectral functions and discuss the properties of a ``semi-compressible'' state observed for large U near half filling.
A fully adaptive wavelet algorithm for parabolic partial differential equations
2001
We present a fully adaptive numerical scheme for the resolution of parabolic equations. It is based on wavelet approximations of functions and operators. Following the numerical analysis in the case of linear equations, we derive a numerical algorithm essentially based on convolution operators that can be efficiently implemented as soon as a natural condition on the space of approximation is satisfied. The algorithm is extended to semi-linear equations with time dependent (adapted) spaces of approximation. Numerical experiments deal with the heat equation as well as the Burgers equation.
Importance of quantiser design compared to optimal multigrid motion estimation in video coding
2000
Adaptive flow computation and DCT quantisation play complementary roles in motion compensated video coding schemes. Since the introduction of the intuitive entropy-constrained motion estimation of Dufaux et al. (1995), several optimal variable-size block matching algorithms have been proposed. Many of these approaches put forward their intrinsic optimality, but the corresponding visual effect has not been explored. The relative importance of optimal multigrid motion estimation with regard to quantisation is addressed in the context of MPEG-like coding. It is shown that while simpler (suboptimal) motion estimates give subjective results as good as the optimal motion estimates, small enhancem…