2,167 research outputs found
Tensor-based multiscale method for diffusion problems in quasi-periodic heterogeneous media
This paper proposes to address the issue of complexity reduction for the
numerical simulation of multiscale media in a quasi-periodic setting. We
consider a stationary elliptic diffusion equation defined on a domain such
that is the union of cells and we
introduce a two-scale representation by identifying any function defined
on with a bi-variate function , where relates to the
index of the cell containing the point and relates to a local
coordinate in a reference cell . We introduce a weak formulation of the
problem in a broken Sobolev space using a discontinuous Galerkin
framework. The problem is then interpreted as a tensor-structured equation by
identifying with a tensor product space of
functions defined over the product set . Tensor numerical methods
are then used in order to exploit approximability properties of quasi-periodic
solutions by low-rank tensors.Comment: Changed the choice of test spaces V(D) and X (with regard to
regularity) and the argumentation thereof. Corrected proof of proposition 3.
Corrected wrong multiplicative factor in proposition 4 and its proof (was 2
instead of 1). Added remark 6 at the end of section 2. Extended remark 7.
Added references. Some minor improvements (typos, typesetting
Block Circulant and Toeplitz Structures in the Linearized Hartree–Fock Equation on Finite Lattices: Tensor Approach
This paper introduces and analyses the new grid-based tensor approach to
approximate solution of the elliptic eigenvalue problem for the 3D
lattice-structured systems. We consider the linearized Hartree-Fock equation
over a spatial lattice for both periodic and
non-periodic problem setting, discretized in the localized Gaussian-type
orbitals basis. In the periodic case, the Galerkin system matrix obeys a
three-level block-circulant structure that allows the FFT-based
diagonalization, while for the finite extended systems in a box (Dirichlet
boundary conditions) we arrive at the perturbed block-Toeplitz representation
providing fast matrix-vector multiplication and low storage size. The proposed
grid-based tensor techniques manifest the twofold benefits: (a) the entries of
the Fock matrix are computed by 1D operations using low-rank tensors
represented on a 3D grid, (b) in the periodic case the low-rank tensor
structure in the diagonal blocks of the Fock matrix in the Fourier space
reduces the conventional 3D FFT to the product of 1D FFTs. Lattice type systems
in a box with Dirichlet boundary conditions are treated numerically by our
previous tensor solver for single molecules, which makes possible calculations
on rather large lattices due to reduced numerical
cost for 3D problems. The numerical simulations for both box-type and periodic
lattice chain in a 3D rectangular "tube" with up to
several hundred confirm the theoretical complexity bounds for the
block-structured eigenvalue solvers in the limit of large .Comment: 30 pages, 12 figures. arXiv admin note: substantial text overlap with
arXiv:1408.383
A literature survey of low-rank tensor approximation techniques
During the last years, low-rank tensor approximation has been established as
a new tool in scientific computing to address large-scale linear and
multilinear algebra problems, which would be intractable by classical
techniques. This survey attempts to give a literature overview of current
developments in this area, with an emphasis on function-related tensors
Application of Hierarchical Matrix Techniques To The Homogenization of Composite Materials
In this paper, we study numerical homogenization methods based on integral
equations. Our work is motivated by materials such as concrete, modeled as
composites structured as randomly distributed inclusions imbedded in a matrix.
We investigate two integral reformulations of the corrector problem to be
solved, namely the equivalent inclusion method based on the Lippmann-Schwinger
equation, and a method based on boundary integral equations. The fully
populated matrices obtained by the discretization of the integral operators are
successfully dealt with using the H-matrix format
- …