Source author record

Ruben Specogna

Ruben Specogna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.AT Computational Engineering, Finance, and Science math.NA Numerical Analysis Computer Vision math.GT

Catalog footprint

What is connected

9works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Explicit geometric construction of sparse inverse mass matrices for arbitrary tetrahedral grids

The geometric reinterpretation of the Finite Element Method (FEM) shows that Raviart Thomas and Nedelec mass matrices map from degrees of freedoms (DoFs) attached to geometric elements of a tetrahedral grid to DoFs attached to the barycentric dual grid. The algebraic inverses of the mass matrices map DoFs attached to the barycentric dual grid back to DoFs attached to the corresponding primal tetrahedral grid, but they are of limited practical use since they are dense. In this paper we present a new geometric construction of sparse inverse mass matrices for arbitrary tetrahedral grids and possibly anisotropic materials, debunking the conventional wisdom that the barycentric dual grid prohibits a sparse representation for inverse mass matrices. In particular, we provide a unified framework for the construction of both edge and face mass matrices and their sparse inverses. Such a unifying principle relies on novel geometric reconstruction formulas, from which, according to a well established design strategy, local mass matrices are constructed as the sum of a consistent and a stabilization term. A major difference with the approaches proposed so far is that the consistent term is defined geometrically and explicitly, that is, without the necessity of computing the inverses of local matrices. This provides a sensible speedup and an easier implementation. We use these new sparse inverse mass matrices to discretize a three dimensional Poisson problem, providing the comparison between the results obtained by various formulations on a benchmark problem with analytical solution.

preprint2021arXiv

Inverting the discrete curl operator: a novel graph algorithm to find a vector potential of a given vector field

We provide a novel framework to compute a discrete vector potential of a given discrete vector field on arbitrary polyhedral meshes. The framework exploits the concept of acyclic matching, a combinatorial tool at the core of discrete Morse theory. We introduce the new concept of complete acyclic matchings and we show that they give the same end result of Gaussian elimination. Basically, instead of doing costly row and column operations on a sparse matrix, we compute equivalent cheap combinatorial operations that preserve the underlying sparsity structure. Currently, the most efficient algorithms proposed in literature to find discrete vector potentials make use of tree-cotree techniques. We show that they compute a special type of complete acyclic matchings. Moreover, we show that the problem of computing them is equivalent to the problem of deciding whether a given mesh has a topological property called collapsibility. This fact gives a topological characterization of well-known termination problems of tree-cotree techniques. We propose a new recursive algorithm to compute discrete vector potentials. It works directly on basis elements of $1$- and $2$-chains by performing elementary Gaussian operations on them associated with acyclic matchings. However, the main novelty is that it can be applied recursively. Indeed, the recursion process allows us to sidetrack termination problems of the standard tree-cotree techniques. We tested the algorithm on pathological triangulations with known topological obstructions. In all tested problems we observe linear computational complexity as a function of mesh size. Moreover, the algorithm is purely graph-based so it is straightforward to implement and does not require specialized external procedures. We believe that our framework could offer new perspectives to sparse matrix computations.

preprint2019arXiv

Efficient computation of Linking number with certification

An efficient numerical algorithm for the computation of linking number is presented. The algorithm keep tracks or rounding error so that it can ensure the correctness of the results.

preprint2016arXiv

Geometric construction of bases of $H_2(\overlineΩ, \partialΩ, \mathbb{Z})$

We present an efficient algorithm for the construction of a basis of $H_2(\overlineΩ,\partialΩ;\mathbb Z)$ via the Poincaré--Lefschetz duality theorem. Denoting by $g$ the first Betti number of $\overline Ω$ the idea is to find, first $g$ different $1$-boundaries of $\overlineΩ$ with supports contained in $\partialΩ$ whose homology classes in $\mathbb R^3 \setminus Ω$ form a basis of $H_1(\mathbb R^3 \setminus Ω;\mathbb Z)$, and then to construct in $\overlineΩ$ a homological Seifert surface of each one of these $1$-boundaries. The Poincaré--Lefschetz duality theorem ensures that the relative homology classes of these homological Seifert surfaces in $\overlineΩ$ modulo $\partialΩ$ form a basis of $H_2(\overlineΩ,\partialΩ;\mathbb Z)$. We devise a simply procedure for the construction of the required set of $1$-boundaries of $\overlineΩ$ that, combined with a fast algorithm for the construction of homological Seifert surfaces, allows the efficient computation of a basis of $H_2(\overlineΩ,\partialΩ;\mathbb Z)$ via this very natural geometrical approach. Some numerical experiments show the efficiency of the method and its performance comparing with other algorithms.

preprint2015arXiv

Modeling of anechoich chambers with equivalent materials and equivalent sources

Numerical simulation of anechoic chambers is a hot topic since it can provide useful data about the performance of the EMC site. However, the mathematical nature of the problem, the physical dimensions of the simulated sites and the frequency ranges pose nontrivial challenges to the simulation. Computational requirements in particular will quickly become unmanageable if adequate techniques are not employed. In this work we describe a novel approach, based on equivalent elements, that enables the simulation of large chambers with modest computational resources. The method is then validated against real measurement results.

preprint2014arXiv

Efficient construction of homological Seifert surfaces

Let $Ω$ be a bounded domain of $\mathbb{R}^3$ whose closure $\overlineΩ$ is polyhedral, and let $\mathcal{T}$ be a triangulation of $\overlineΩ$. Assuming that the boundary of $Ω$ is sufficiently regular, we provide an explicit formula for the computation of homological Seifert surfaces of any $1$-boundary $γ$ of $\mathcal{T}$; namely, $2$-chains of $\mathcal{T}$ whose boundary is $γ$. It is based on the existence of special spanning trees of the complete dual graph of $\mathcal{T}$, and on the computation of certain linking numbers associated with those spanning trees. If the triangulation $\mathcal{T}$ is fine, the explicit formula is too expensive to be used directly. For this reason, making also use of a simple elimination procedure, we devise a fast algorithm for the computation of homological Seifert surfaces. Some numerical experiments illustrate the efficiency of this algorithm.

preprint2014arXiv

Topology preserving thinning for cell complexes

A topology preserving skeleton is a synthetic representation of an object that retains its topology and many of its significant morphological properties. The process of obtaining the skeleton, referred to as skeletonization or thinning, is a very active research area. It plays a central role in reducing the amount of information to be processed during image analysis and visualization, computer-aided diagnosis or by pattern recognition algorithms. This paper introduces a novel topology preserving thinning algorithm which removes \textit{simple cells}---a generalization of simple points---of a given cell complex. The test for simple cells is based on \textit{acyclicity tables} automatically produced in advance with homology computations. Using acyclicity tables render the implementation of thinning algorithms straightforward. Moreover, the fact that tables are automatically filled for all possible configurations allows to rigorously prove the generality of the algorithm and to obtain fool-proof implementations. The novel approach enables, for the first time, according to our knowledge, to thin a general unstructured simplicial complex. Acyclicity tables for cubical and simplicial complexes and an open source implementation of the thinning algorithm are provided as additional material to allow their immediate use in the vast number of practical applications arising in medical imaging and beyond.

preprint2012arXiv

Physics inspired algorithms for (co)homology computation

The issue of computing (co)homology generators of a cell complex is gaining a pivotal role in various branches of science. While this issue can be rigorously solved in polynomial time, it is still overly demanding for large scale problems. Drawing inspiration from low-frequency electrodynamics, this paper presents a physics inspired algorithm for first cohomology group computations on three-dimensional complexes. The algorithm is general and exhibits orders of magnitude speed up with respect to competing ones, allowing to handle problems not addressable before. In particular, when generators are employed in the physical modeling of magneto-quasistatic problems, this algorithm solves one of the most long-lasting problems in low-frequency computational electromagnetics. In this case, the effectiveness of the algorithm and its ease of implementation may be even improved by introducing the novel concept of \textit{lazy cohomology generators}.

preprint2011arXiv

Cohomology in electromagnetic modeling

Electromagnetic modeling provides an interesting context to present a link between physical phenomena and homology and cohomology theories. Over the past twenty-five years, a considerable effort has been invested by the computational electromagnetics community to develop fast and general techniques for potential design. When magneto-quasi-static discrete formulations based on magnetic scalar potential are employed in problems which involve conductive regions with holes, \textit{cuts} are needed to make the boundary value problem well defined. While an intimate connection with homology theory has been quickly recognized, heuristic definitions of cuts are surprisingly still dominant in the literature. The aim of this paper is first to survey several definitions of cuts together with their shortcomings. Then, cuts are defined as generators of the first cohomology group over integers of a finite CW-complex. This provably general definition has also the virtue of providing an automatic, general and efficient algorithm for the computation of cuts. Some counter-examples show that heuristic definitions of cuts should be abandoned. The use of cohomology theory is not an option but the invaluable tool expressly needed to solve this problem.