Our interest lies in the identification of hidden conducting permeable
objects from measurements of the perturbed magnetic field in metal detection
taken over range of low frequencies. The magnetic polarizability tensor (MPT)
provides a characterisation of a conducting permeable object using a small
number of coefficients, has explicit formula for their calculation and a well
understood frequency behaviour, which we call its spectral signature. However,
to compute such signatures, and build a library of them for object
classification, requires repeated solution of a direct (full order) problem,
which is typically accomplished using a finite element discretisation. To
overcome this issue, we propose an efficient reduced order model (ROM) using a
proper orthogonal decomposition (POD) for the rapid computation of MPT spectral
signatures. Our ROM benefits from output certificates, which give bounds on the
accuracy of the predicted outputs with respect to the full order model
solutions. To further increase the efficiency of the...

more |
pdf
| html
None.

mathNAb:
B.A. Wilson, P.D. Ledger : Efficient Computation of the Magnetic Polarizability Tensor Spectral Signature using POD https://t.co/YA97Z8cdhy https://t.co/skInDIG81l

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

A thick-restart Lanczos type algorithm is proposed for Hermitian
$J$-symmetric matrices. Since Hermitian $J$-symmetric matrices possess doubly
degenerate spectra or doubly multiple eigenvalues with a simple relation
between the degenerate eigenvectors, we can improve the convergence of the
Lanczos algorithm by restricting the search space of the Krylov subspace to
that spanned by one of each pair of the degenerate eigenvector pairs. We show
that the Lanczos iteration is compatible with the $J$-symmetry, so that the
subspace can be split into two subspaces that are orthogonal to each other. The
proposed algorithm searches for eigenvectors in one of the two subspaces
without the multiplicity. The other eigenvectors paired to them can be easily
reconstructed with the simple relation from the $J$-symmetry. We test our
algorithm on a matrix originating from a quantum field theory.

more |
pdf
| html
None.

mathNAb:
Ken-Ichi Ishikawa, Tomohiro Sogabe : A thick-restart Lanczos type method for Hermitian $J$-symmetric eigenvalue problems https://t.co/FHTO9B3WHi https://t.co/1BXPEj82RO

LatticePapers:
A thick-restart Lanczos type method for Hermitian $J$-symmetric eigenvalue problems. https://t.co/Q1ZDHQcuYV

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

Automatic differentiation (AD) is an ensemble of techniques that allow to
evaluate accurate numerical derivatives of a mathematical function expressed in
a computer programming language. In this paper we use AD for stating and
solving solid mechanics problems. Given a finite element discretization of the
domain, we evaluate the free energy of the solid as the integral of its strain
energy density, and we make use of AD for directly obtaining the residual force
vector and the tangent stiffness matrix of the problem, as the gradient and the
Hessian of the free energy respectively. The result is a remarkable
simplification in the statement and the solution of complex problems involving
non trivial constraints systems and both geometrical and material non
linearities. Together with the continuum mechanics theoretical basis, and with
a description of the specific AD technique adopted, the paper illustrates the
solution of a number of solid mechanics problems, with the aim of presenting a
convenient numerical implementation approach,...

more |
pdf
| html
None.

mathNAb:
Andrea Vigliotti, Ferdinando Auricchio : Automatic differentiation for solid mechanics https://t.co/jwo8vIHdK7 https://t.co/dmC6rO364D

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

A simple variant of the BDDC preconditioner in which constraints are imposed
on a selected set of subobjects (subdomain subedges, subfaces and vertices
between pairs of subedges) is presented. We are able to show that the condition
number of the preconditioner is bounded by $C \big(1+\log (L/h)\big)^2$, where
$C$ is a constant, and $h$ and $L$ are the characteristic sizes of the mesh and
the subobjects, respectively. As $L$ can be chosen almost freely, the condition
number can theoretically be as small as $O(1)$. We will discuss the pros and
cons of the preconditioner and its application to heterogeneous problems.
Numerical results on supercomputers are provided.

more |
pdf
| html
None.

mathNAb:
Santiago Badia, Alberto F. Martín, Hieu Nguyen : Balancing domain decomposition by constraints associated with subobjects https://t.co/K9phasRkQF https://t.co/kMHptheNqQ

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 0

Unqiue Words: 0

Derksen proved that the spectral norm is multiplicative with respect to
vertical tensor products (also known as tensor Kronecker products). We will use
this result to show that the nuclear norm and other norms of interest are also
multiplicative with respect to vertical tensor products. Finally, we prove an
analogue of Strassen's direct sum conjecture for the nuclear norm.

more |
pdf
| html
None.

mathNAb:
Robert Cochrane : Nuclear Norm Under Tensor Kronecker Products https://t.co/wlWApVoUQH https://t.co/6KXcQE3nJ7

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 0

Unqiue Words: 0

We consider a Lagrangian system $L(q,\dot q) =
\sum_{l=1}^{N}L^{\{l\}}(q,\dot q)$, where the $q$-variable is treated by a
Generalized Additive Runge--Kutta (GARK) method. Applying the technique of
discrete variations, we show how to construct symplectic schemes. Assuming the
diagonal methods for the GARK given, we present some techinques for
constructing the transition matrices. We address the problem of the order of
the methods and discuss some semi-separable and separable problems, showing
some interesting constructions of methods with non-square coefficient matrices.

more |
pdf
| html
None.

mathNAb:
Antonella Zanna : Discrete Variational Methods and Symplectic Generalized Additive Runge--Kutta Methods https://t.co/eEwefdsd1E https://t.co/JsxBSwluSb

mathSGb:
Antonella Zanna : Discrete Variational Methods and Symplectic Generalized Additive Runge--Kutta Methods https://t.co/DFZ92TyjZm https://t.co/MtjZtOJj2z

None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 0

Unqiue Words: 0

In our work, we consider the linear least squares problem for $m\times
n$-systems of linear equations $Ax = b$, $m\geq n$, such that the matrix $A$
and right-hand side vector $b$ can vary within an interval $m\times n$-matrix
and an interval $m$-vector respectively. We have to compute, as sharp as
possible, an interval enclosure of the set of all least squares solutions to
$Ax = b$ when $A$ and $b$ independently vary within their interval bounds. Our
article is devoted to the development of the so-called PPS-methods (based on
Partitioning of the Parameter Set) to solve the above problem. We reduce the
normal equation system, associated with the linear lest squares problem, to a
special extended matrix form and produce a symmetric interval system of linear
equations that is equivalent to the original interval least squares problem. To
solve such symmetric system, we propose a new construction of PPS-methods,
called ILSQ-PPS, which estimates the enclosure of the solution set with
practical efficiency. To demonstrate the capabilities...

more |
pdf
| html
None.

mathNAb:
Sergey P. Shary, Behnam Moradi : Solving interval linear least squares problems by PPS-methods https://t.co/0sRkhpLYue https://t.co/IELasCxpY2

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

Big data analysis has become a crucial part of new emerging technologies such
as Internet of thing (IOT), cyber-physical analysis, deep learning, anomaly
detection etc. Among many other techniques, dimensionality reduction plays a
key role in such analyses and facilitate the procedure of feature selection and
feature extraction. Randomized algorithms are efficient tools for handling big
data tensors. They accelerate decomposing large-scale data tensors by reducing
the computational complexity of deterministic algorithms and also reducing the
communication among different levels of memory hierarchy which is a main
bottleneck in modern computing environments and architectures. In this paper,
we review recent advances in randomization for computation of Tucker
decomposition and Higher Order SVD (HOSVD). We discuss both random projection
and sampling approaches and also single-pass and multi-pass randomized
algorithms and how they can be utilized in computation of Tucker decomposition
and HOSVD. Simulations on real data including...

more |
pdf
| html
None.

mathNAb:
Salman Ahmadi-Asl, Andrzej Cichocki, Anh Huy Phan, Ivan Oseledets, Stanislav Abukhovich, Toshihisa Tanaka : Randomized Algorithms for Computation of Tucker decomposition and Higher Order SVD (HOSVD) https://t.co/eP7J7SAdWz https://t.co/fiRi3pp7PB

None.

None.

Sample Sizes : None.

Authors: 6

Total Words: 0

Unqiue Words: 0

An energy-based model of the ferroelectric polarization process is presented
in the current contribution. In an energy-based setting, dielectric
displacement and strain (or displacement) are the primary independent unknowns.
As an internal variable, the remanent polarization vector is chosen. The model
is then governed by two constitutive functions: the free energy function and
the dissipation function. Choices for both functions are given. As the
dissipation function for rate-independent response is non-differentiable, it is
proposed to regularize the problem. Then, a variational equation can be posed,
which is subsequently discretized using conforming finite elements for each
quantity. We point out which kind of continuity is needed for each field
(displacement, dielectric displacement and remanent polarization) is necessary
to obtain a conforming method, and provide corresponding finite elements. The
elements are chosen such that Gauss' law of zero charges is satisfied exactly.
The discretized variational equations are solved...

more |
pdf
| html
None.

mathNAb:
Astrid S. Pechstein, Martin Meindlhumer, Alexander Humer : High-order mixed finite elements for an energy-based model of the polarization process in ferroelectric materials https://t.co/rbh1Olg1SP https://t.co/fIXNc8CuPN

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 0

Unqiue Words: 0

Adaptive meshing is a fundamental component of adaptive finite element
methods. This includes refining and coarsening meshes locally. In this work, we
are concerned with the red-green-blue refinement strategy and its counterpart -
coarsening. In general, coarsening algorithms are mostly based on an explicitly
given refinement history. In this work, we present a coarsening algorithm on
adaptive red-green-blue meshes without explicitly knowing the refinement
history. To this end, we examine the local structure of these meshes, find an
easy-to-verify criterion to coarsen red-green-blue meshes and prove that this
criterion generates meshes with the desired properties. However, it does not
guarantee that the set of nodes admissible for coarsening is non-empty.
Therefore, we also present an additional algorithm that uses the main ideas of
this criterion to always allow for local coarsening. We present a MATLAB
implementation built on the red-green-blue refinement routine of the
ameshref-package.

more |
pdf
| html
None.

mathNAb:
Stefan A. Funken, Anja Schmidt : A Coarsening Algorithm on Adaptive Red-Green-Blue Refined Meshes https://t.co/eCkQZUljdV https://t.co/0cqvZZZEmh

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 9354

Unqiue Words: 2083

Assert is a website where the best academic papers on arXiv (computer science, math, physics), bioRxiv (biology), BITSS (reproducibility), EarthArXiv (earth science), engrXiv (engineering), LawArXiv (law), PsyArXiv (psychology), SocArXiv (social science), and SportRxiv (sport research) bubble to the top each day.

Papers are scored (in real-time) based on how verifiable they are (as determined by their Github repos) and how interesting they are (based on Twitter).

To see top papers, follow us on twitter @assertpub_ (arXiv), @assert_pub (bioRxiv), and @assertpub_dev (everything else).

To see beautiful figures extracted from papers, follow us on Instagram.

*Tracking 256,578 papers.*

Sort results based on if they are interesting or reproducible.

Interesting

Reproducible