BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Lyes Rahmouni

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

ResearcherUnclaimed source record
physics.comp-pheess.IVmath.NANumerical Analysisphysics.med-ph

Catalog footprint

What is connected

8works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

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
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

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
preprint2016arXiv

Conforming Discretizations of Boundary Element Solutions of the Electroencephalography Forward Problem

In this paper we present a new discretization strategy for the boundary element formulation of the Electroencephalography (EEG) forward problem. Boundary integral formulations, classically solved with the Boundary Element Method (BEM), are widely used in high resolution EEG imaging because of their recognized advantages in several real case scenarios. Unfortunately however, it is widely reported that the accuracy of standard BEM schemes is limited, especially when the current source density is dipolar and its location approaches one of the brain boundary surfaces. This is a particularly limiting problem given that during an high-resolution EEG imaging procedure, several EEG forward problem solutions are required for which the source currents are near or on top of a boundary surface. This work will first present an analysis of standardly discretized EEG forward problems, reporting on a theoretical issue of some of the formulations that have been used so far in the community. We report on the fact that several standardly used discretizations requires the expansion term to be a square integrable function. Instead, those techniques are not consistent when a more appropriate mapping is considered. Such a mapping allows the expansion function term to be a less regular function, thus sensibly reducing the need for mesh refinements and low-precisions handling strategies that are currently required. These mappings, however, require a different and conforming discretization which must be suitably adapted to them. To do this we adopt a mixed discretization based on dual boundary elements residing on a suitably defined dual mesh. Finally we show how the resulting EEG problems has favorable properties with respect to previously proposed schemes and we show their applicability to real case modeling scenarios obtained from MRI data.

0 citations0 reviewsNo code linkedOpen access