### Top 10 Arxiv Papers Today in Numerical Analysis

##### #1. Efficient Computation of the Magnetic Polarizability Tensor Spectral Signature using POD
###### B. A. Wilson, P. D. Ledger
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #2. A thick-restart Lanczos type method for Hermitian $J$-symmetric eigenvalue problems
###### Ken-Ichi Ishikawa, Tomohiro Sogabe
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #3. Automatic differentiation for solid mechanics
###### Andrea Vigliotti, Ferdinando Auricchio
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.
###### Tweets
mathNAb: Andrea Vigliotti, Ferdinando Auricchio : Automatic differentiation for solid mechanics https://t.co/jwo8vIHdK7 https://t.co/dmC6rO364D
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #4. Balancing domain decomposition by constraints associated with subobjects
###### Santiago Badia, Alberto F. Martín, Hieu Nguyen
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 0
Unqiue Words: 0

##### #5. Nuclear Norm Under Tensor Kronecker Products
###### Robert Cochrane
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.
###### Tweets
mathNAb: Robert Cochrane : Nuclear Norm Under Tensor Kronecker Products https://t.co/wlWApVoUQH https://t.co/6KXcQE3nJ7
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 0
Unqiue Words: 0

##### #6. Discrete Variational Methods and Symplectic Generalized Additive Runge--Kutta Methods
###### Antonella Zanna
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 1
Total Words: 0
Unqiue Words: 0

##### #7. Solving interval linear least squares problems by PPS-methods
###### Sergey P. Shary, Behnam Moradi
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #8. Randomized Algorithms for Computation of Tucker decomposition and Higher Order SVD (HOSVD)
###### Salman Ahmadi-Asl, Andrzej Cichocki, Anh Huy Phan, Ivan Oseledets, Stanislav Abukhovich, Toshihisa Tanaka
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 6
Total Words: 0
Unqiue Words: 0

##### #9. High-order mixed finite elements for an energy-based model of the polarization process in ferroelectric materials
###### Astrid S. Pechstein, Martin Meindlhumer, Alexander Humer
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.
###### Tweets
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.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 0
Unqiue Words: 0

##### #10. A Coarsening Algorithm on Adaptive Red-Green-Blue Refined Meshes
###### Stefan A. Funken, Anja Schmidt
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.
###### Tweets
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.
###### Other stats
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.

###### Search
Sort results based on if they are interesting or reproducible.
Interesting
Reproducible
Online
###### Stats
Tracking 256,578 papers.