3,206 research outputs found
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
Numerical methods for calculating poles of the scattering matrix with applications in grating theory
Waveguide and resonant properties of diffractive structures are often
explained through the complex poles of their scattering matrices. Numerical
methods for calculating poles of the scattering matrix with applications in
grating theory are discussed. A new iterative method for computing the matrix
poles is proposed. The method takes account of the scattering matrix form in
the pole vicinity and relies upon solving matrix equations with use of matrix
decompositions. Using the same mathematical approach, we also describe a
Cauchy-integral-based method that allows all the poles in a specified domain to
be calculated. Calculation of the modes of a metal-dielectric diffraction
grating shows that the iterative method proposed has the high rate of
convergence and is numerically stable for large-dimension scattering matrices.
An important advantage of the proposed method is that it usually converges to
the nearest pole.Comment: 9 pages, 2 figures, 4 table
Minimizing Communication for Eigenproblems and the Singular Value Decomposition
Algorithms have two costs: arithmetic and communication. The latter
represents the cost of moving data, either between levels of a memory
hierarchy, or between processors over a network. Communication often dominates
arithmetic and represents a rapidly increasing proportion of the total cost, so
we seek algorithms that minimize communication. In \cite{BDHS10} lower bounds
were presented on the amount of communication required for essentially all
-like algorithms for linear algebra, including eigenvalue problems and
the SVD. Conventional algorithms, including those currently implemented in
(Sca)LAPACK, perform asymptotically more communication than these lower bounds
require. In this paper we present parallel and sequential eigenvalue algorithms
(for pencils, nonsymmetric matrices, and symmetric matrices) and SVD algorithms
that do attain these lower bounds, and analyze their convergence and
communication costs.Comment: 43 pages, 11 figure
- …