BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ruben Specogna

Ruben Specogna contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.ATmath.NANumerical AnalysisComputational Engineering, Finance, and Sciencemath.GT

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access