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

Francesco P. Andriulli

Francesco P. Andriulli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phmath.NANumerical Analysiseess.IVeess.SPphysics.med-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
6topics
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

14 published item(s)

preprint2022arXiv

A New Refinement-Free Preconditioner for the Symmetric Formulation in Electroencephalography

Widely employed for the accurate solution of the electroencephalography forward problem, the symmetric formulation gives rise to a first kind, ill-conditioned operator ill-suited for complex modelling scenarios. This work presents a novel preconditioning strategy based on an accurate spectral analysis of the operators involved which, differently from other Calderón-based approaches, does not necessitate the barycentric refinement of the primal mesh (i.e., no dual matrix is required). The discretization of the new formulation gives rise to a well-conditioned, symmetric, positive-definite system matrix, which can be efficiently solved via fast iterative techniques. Numerical results for both canonical and realistic head models validate the effectiveness of the proposed formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Brain-Computer Interfaces: Investigating the Transition from Visually Evoked to Purely Imagined Steady-State Potentials

Brain-Computer Interfaces (BCIs) based on Steady State Visually Evoked Potentials (SSVEPs) have proven effective and provide significant accuracy and information-transfer rates. This family of strategies, however, requires external devices that provide the frequency stimuli required by the technique. This limits the scenarios in which they can be applied, especially when compared to other BCI approaches. In this work, we have investigated the possibility of obtaining frequency responses in the EEG output based on the pure visual imagination of SSVEP-eliciting stimuli. Our results show that not only that EEG signals present frequency-specific peaks related to the frequency the user is focusing on, but also that promising classification accuracy can be achieved, paving the way for a robust and reliable visual imagery BCI modality.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Laplacian Filters for Integral Equations: Further Developments and Fast Algorithms

This paper extends the concept of Laplacian filtered quasi-Helmholtz decompositions we have recently introduced, to the basis-free projector-based setting. This extension allows the discrete analyses of electromagnetic integral operators spectra without passing via an explicit Loop-Star decomposition as previously done. We also present a fast scheme for the evaluation of the filters in quasi linear complexity in the total number of unknowns. Together with the fact that only a logarithmic number of these filters are required for solving the h-refinement breakdown of electric field integral equation, this results in an effective preconditioner that rivals Calderón strategies in performance without relying on barycentric refinements. Numerical results confirm the theoretically predicted behavior and the effectiveness of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the Fast Direct Solution of a Preconditioned Electromagnetic Integral Equation

This work presents a fast direct solver strategy for electromagnetic integral equations in the high-frequency regime. The new scheme relies on a suitably preconditioned combined field formulation and results in a single skeleton form plus identity equation. This is obtained after a regularization of the elliptic spectrum through the extraction of a suitably chosen equivalent circulant problem. The inverse of the system matrix is then obtained by leveraging the Woodbury matrix identity, the low-rank representation of the extracted part of the operator, and fast circulant algebra yielding a scheme with a favorable complexity and suitable for the solution of multiple right-hand sides. Theoretical considerations are accompanied by numerical results both of which are confirming and showing the practical relevance of the newly developed scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stabilized Single Current Inverse Source Formulations Based on Steklov-Poincaré Mappings

The inverse source problem in electromagnetics has proved quite relevant for a large class of applications. In antenna diagnostics in particular, Love solutions are often sought at the cost of an increase of the dimension of the linear system to be solved. In this work, instead, we present a reduced-in-size single current formulation of the inverse source problem that obtains one of the Love currents via a stable discretization of the Steklov-Poincaré boundary operator leveraging dual functions. The new approach is enriched by theoretical treatments and by a further low-frequency stabilization of the Steklov-Poincaré operator based on the quasi-Helmholtz projectors that is the first of its kind in this field. The effectiveness and practical relevance of the new schemes are demonstrated via both theoretical and numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Hybrid Volume-Surface-Wire Integral Equation for the Anisotropic Forward Problem in Electroencephalography

Solving the electroencephalography (EEG) forward problem is a fundamental step in a wide range of applications including biomedical imaging techniques based on inverse source localization. State-of-the-art electromagnetic solvers resort to a computationally expensive volumetric discretization of the full head to account for its complex and heterogeneous electric profile. The more efficient, popular in biomedical imaging circles, but unfortunately oversimplifying Boundary Element Method (BEM) relies instead on a piecewise-uniform approximation that severely curbs its application in high resolution EEGs. This contribution lifts the standard BEM contraints by treating the local anisotropies with adequate wire and thin volume integral equations that are tailored to specific structures of the fibrous white matter and the inhomogeneous skull. The proposed hybrid integral equation formulation thereby avoids the full volumetric discretization of the head medium and allows for a realistic and efficient BEM-like solution of the anisotropic EEG forward problem. The accuracy and flexibility of the proposed formulation is demonstrated through numerical experiments involving both canonical and realistic MRI-based head models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A New Preconditioner for the EFIE Based on Primal and Dual Graph Laplacian Spectral Filters

The Electric Field Integral Equation (EFIE) is notorious for its ill-conditioning both in frequency and h-refinement. Several techniques exist for fixing the equation conditioning problems based on hierarchical strategies, Calderon techniques, and related technologies. This work leverages on a new approach, based on the construction of tailored spectral filters for the EFIE components which allow the block renormalization of the EFIE spectrum resulting in a provably constant condition number for the equation. This is achieved without the need for a barycentric refinement and with low computational overhead compared with other schemes. In particular, only sparse matrices are required in addition to the EFIE original matrix. Numerical results will show the robustness of our scheme and its application to the solution of realistic problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Quasi-Helmholtz Projector Stabilized Full Wave Solver Encompassing the Eddy Current Regime

Despite its several qualities, the Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) formulation for simulating scattering by dielectric media suffers from numerical instabilities and severe ill-conditioning at low frequencies. While this drawback has been the object of numerous solution attempts in the standard low-frequency breakdown regime for scattering problems, the eddy-current regime requires a specific treatment. In this contribution, we present an extension of the recently introduced quasi-Helmholtz projectors based preconditioning of the PMCHWT to obtain an equation stable at low frequencies and specifically in three regimes relevant for eddy currents analyses: (i) when the frequency decreases with a constant conductivity, (ii) when the conductivity increases at fixed frequency and (iii) when the frequency decreases while keeping the product of the frequency and the conductivity constant. Being based on quasi-Helmholtz projectors our new strategy does not further degrade the conditioning of the original equation and is compatible with existing fast solvers. The resulting full-wave formulation is capable of smoothly transitioning from simulations at high frequencies to the different low-frequency regimes which encompass, in particular, eddy currents applications. Numerical results demonstrate the validity of our approach in all the regimes with a special emphasis given to eddy currents.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diffusion MRI Consistent Wire Models for Efficient Solutions of the Anisotropic Forward Problem in Electroencephalography

The surface Boundary Element Method (BEM) is one of the most commonly employed formulations to solve the forward problem in electroencephalography, but the applicability of its classical incarnations is lamentably limited to piece-wise homogeneous media. Several head tissues, however, are strongly anisotropic due to their complex underlying micro-structure. This implies that standard boundary integral formulations oversimplify the electrical properties of the head and produce unrealistic solutions, something that drastically limits the suitability and impact of BEM technologies to brain imaging. This contribution addresses this issue by observing that the brain anisotropy in the white matter is due to the presence of neuronal wire-like structures. We then extend the well known wire integral equations used for high frequency problems to the imperfectly conducting quasi-static case and we propose a new hybrid wire/surface/volume integral equation. When applied on multimodal magnetic resonance images combined with tractography, this new approach can flexibly and realistically handle the conductivity anisotropy in any head compartment providing high level of accuracy and efficiency. The beneficial properties of the new formulation together with its impact on brain imaging is demonstrated via numerical results on both canonical and realistic case scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Large Time Step and DC Stable TD-EFIE Discretized with Implicit Runge-Kutta Methods

The Time Domain-Electric Field Integral Equation (TD-EFIE) and its differentiated version are widely used to simulate the transient scattering of a time dependent electromagnetic field by a Perfect Electrical Conductor (PEC). The time discretization of the TD-EFIE can be achieved by a space-time Galerkin approach or, as it is considered in this contribution, by a convolution quadrature using Implicit Runge-Kutta methods. The solution is then computed using the Marching-On-in-Time (MOT) algorithm. The differentiated TD-EFIE has two problems: (i) the system matrix suffers from ill-conditioning when the time step increases (low frequency breakdown) and (ii) it suffers from the DC instability, i.e. the formulation allows for the existence of spurious solenoidal currents that grow slowly in the solution. In this work, we show that (i) and (ii) can be alleviated by leveraging quasi-Helmholtz projectors to separate the Helmholtz components of the induced current and rescale them independently. The efficacy of the approach is demonstrated by numerical examples including benchmarks and real life applications.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On Preconditioning Electromagnetic Integral Equations in the High Frequency Regime via Helmholtz Operators and quasi-Helmholtz Projectors

Fast and accurate resolution of electromagnetic problems via the \ac{BEM} is oftentimes challenged by conditioning issues occurring in three distinct regimes: (i) when the frequency decreases and the discretization density remains constant, (ii) when the frequency is kept constant while the discretization is refined and (iii) when the frequency increases along with the discretization density. While satisfactory remedies to the problems arising in regimes (i) and (ii), respectively based on Helmholtz decompositions and Calderón-like techniques have been presented, the last regime is still challenging. In fact, this last regime is plagued by both spurious resonances and ill-conditioning, the former can be tackled via combined field strategies and is not the topic of this work. In this contribution new symmetric scalar and vectorial electric type formulations that remain well-conditioned in all of the aforementioned regimes and that do not require barycentric discretization of the dense electromagnetic potential operators are presented along with a spherical harmonics analysis illustrating their key properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the modeling of brain fibers in the EEG forward problem via a new family of wire integral equations

Source localization based on electroencephalography (EEG) has become a widely used neuroimagining technique. However its precision has been shown to be very dependent on how accurately the brain, head and scalp can be electrically modeled within the so-called forward problem. The construction of this model is traditionally performed by leveraging Finite Element or Boundary Element Methods (FEM or BEM). Even though the latter is more computationally efficient thanks to the smaller interaction matrices it yields and near-linear solvers, it has traditionally been used on simpler models than the former. Indeed, while FEM models taking into account the different media anisotropies are widely available, BEM models have been limited to isotropic, piecewise homogeneous models. In this work we introduce a new BEM scheme taking into account the anisotropies of the white matter. The boundary nature of the formulation allows for an efficient discretization and modelling of the fibrous nature of the white matter as one-dimensional basis functions, limiting the computational impact of their modelling. We compare our scheme against widely used formulations and establish its correctness in both canonical and realistic cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Spectral Behavior and Normalization of a Resonance-Free and High-Frequency Stable Integral Equation

The CFIE used for solving scattering and radiation problems, although a resonance-free formulation, suffers from an ill-conditioning that strongly depends on the frequency and discretization density, both in the low- and high-frequency regime, resulting in slow convergence rates for iterative solvers. This work presents a new preconditioning scheme for the CFIE that cures the low- and the high-frequency as well as the dense discretization breakdown. The new preconditioner for the CFIE is based on a spherical harmonics analysis and the proper regularization with Helmholtz-type operators. Numerical results have been obtained to prove the effectiveness of this new formulation in real scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Magnetic and Combined Field Integral Equations Based on the Quasi-Helmholtz Projectors

Boundary integral equation methods for analyzing electromagnetic scattering phenomena typically suffer from several of the following problems: (i) ill-conditioning when the frequency is low; (ii) ill-conditioning when the discretization density is high; (iii) ill-conditioning when the structure contains global loops (which are computationally expensive to detect); (iv) incorrect solution at low frequencies due to current cancellations; (v) presence of spurious resonances. In this paper, quasi-Helmholtz projectors are leveraged to obtain a magnetic field integral equation (MFIE) formulation that is immune to drawbacks (i)-(iv). Moreover, when this new MFIE is combined with a regularized electric field integral equation, a new quasi-Helmholtz projector combined field integral equation is obtained that also is immune to (v). Numerical results corroborate the theory and show the practical impact of the newly proposed formulations.

0 citations0 reviewsNo code linkedOpen access