Researcher profile

Lyes Rahmouni

Lyes Rahmouni contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
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

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

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.

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.

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.

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.

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.

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.