The team develops constructive, function-theoretic approaches to inverse problems arising in modeling and design, in particular for electro-magnetic systems as well as in the analysis of certain classes of signals.

Data typically consist of measurements or desired behaviors. The general thread is to approximate them by families of solutions to the equations governing the underlying system. This leads us to consider various interpolation and approximation problems in classes of rational and meromorphic functions, harmonic gradients, or solutions to more general elliptic partial differential equations (PDE), in connection with inverse potential problems. A recurring difficulty is to control the singularities of the approximants.

The mathematical tools pertain to complex and harmonic analysis, approximation theory, potential theory, system theory, differential topology, optimization and computer algebra. Targeted applications include:

identification and synthesis of analog microwave devices (filters, amplifiers),

non-destructive control from field measurements in medical engineering (source recovery in magneto/electro-encephalography), and paleomagnetism (determining the magnetization of rock samples).

In each case, the endeavor is to develop algorithms resulting in dedicated software.

Within the extensive field of inverse problems, much of the research by Apics
deals with reconstructing solutions of classical elliptic PDEs from their
boundary behavior. Perhaps the simplest example lies with
harmonic identification of a stable linear dynamical system:
the transfer-function *e.g.* the Cauchy formula.

Practice is not nearly as simple, for *i.e.* to locate the

Step 1 relates to extremal
problems and analytic operator theory, see Section .
Step 2 involves optimization, and some Schur analysis
to parametrize transfer matrices of given Mc-Millan degree
when dealing with systems having several inputs and outputs,
see Section .
It also makes contact with the topology of rational functions, in particular
to count
critical points and to derive bounds, see Section . Step 2 raises
further issues in approximation theory regarding the rate of convergence and
the extent to which singularities of the
approximant (*i.e.* its poles) tend to singularities of the
approximated function; this is where logarithmic potential theory
becomes instrumental, see Section .

Applying a realization procedure to the result of step 2 yields an identification procedure from incomplete frequency data which was first demonstrated in to tune resonant microwave filters. Harmonic identification of nonlinear systems around a stable equilibrium can also be envisaged by combining the previous steps with exact linearization techniques from .

A similar path can be taken to approach design problems in the frequency domain, replacing the measured behavior by some desired behavior. However, describing achievable responses in terms of the design parameters is often cumbersome, and most constructive techniques rely on specific criteria adapted to the physics of the problem. This is especially true of filters, the design of which traditionally appeals to polynomial extremal problems , . Apics contributed to this area the use of Zolotarev-like problems for multi-band synthesis, although we presently favor interpolation techniques in which parameters arise in a more transparent manner, as well as convex relaxation of hyperbolic approximation problems, see Sections and .

The previous example of harmonic identification
quickly suggests a generalization
of itself. Indeed, on identifying *i.e.*, the field) on part of a hypersurface (a curve in 2-D)
encompassing the support of

Inverse potential problems are severely indeterminate because infinitely many
measures within an open set produce the same field outside this set; this phenomenon is called
*balayage* . In the two steps approach
previously described,
we implicitly removed this indeterminacy by requiring in step 1
that the measure
be supported on the boundary (because we seek a function holomorphic
throughout the right half-space), and
by requiring in step 2
that the measure be discrete in the left half-plane (in fact: a sum of
point masses

To recap, the gist of our approach is to approximate boundary data by (boundary traces of) fields arising from potentials of measures with specific support. This differs from standard approaches to inverse problems, where descent algorithms are applied to integration schemes of the direct problem; in such methods, it is the equation which gets approximated (in fact: discretized).

Along these lines, Apics advocates the use of steps 1 and 2 above, along with some singularity analysis, to approach issues of nondestructive control in 2-D and 3-D , , . The team is currently engaged in the generalization to inverse source problems for the Laplace equation in 3-D, to be described further in Section . There, holomorphic functions are replaced by harmonic gradients; applications are to inverse source problems in neurosciences (in particular in EEG/MEG) and inverse magnetization problems in geosciences, see Section .

The approximation-theoretic tools developed by Apics to handle issues mentioned so far are outlined in Section . In Section to come, we describe in more detail which problems are considered and which applications are targeted.

Apics is reaching the end of its 12 years life cycle. A reorganization of the team and of some of its research themes is under way through the project proposal FACTAS.

By standard properties of conjugate differentials, reconstructing Dirichlet-Neumann boundary conditions
for a function harmonic in a plane domain,
when these conditions are already known on a subset

Another application by the team deals with non-constant conductivity
over a doubly connected domain, the set *Tore Supra*
.
The procedure is fast because no numerical integration of
the underlying PDE is needed, as an explicit basis of solutions to the
conjugate Beltrami equation in terms of Bessel functions
was found in this case. Generalizing this approach in a more systematic
manner to free boundary problems of Bernoulli type,
using descent
algorithms based on shape-gradient for such approximation-theoretic
criteria, is an interesting prospect to the team.

The piece of work we just mentioned requires defining and studying Hardy
spaces of the conjugate-Beltrami equation, which is an interesting topic
by itself. For Sobolev-smooth coefficients
of exponent greater than 2,
they were investigated in , .
The case of the critical exponent 2 is treated in , which apparently provides the first example of well-posedness for the Dirichlet problem in the non-strictly elliptic case: the conductivity may be unbounded or zero on sets of zero capacity and, accordingly, solutions need not be locally bounded. More importantly perhaps,
the exponent 2 is also the key to a corresponding theory on very
general (still rectifiable) domains in the plane, as
coefficients of pseudo-holomorphic functions
obtained by conformal transformation onto a disk are merely of

The 3-D version of step 1 in Section is another
subject investigated by Apics: to recover a harmonic function
(up to an additive constant) in a ball or a half-space from partial knowledge of its
gradient. This prototypical inverse problem
(*i.e.* inverse to the Cauchy problem for the Laplace equation)
often recurs in electromagnetism. At present, Apics is involved with
solving instances of this inverse problem arising
in two fields, namely medical imaging
*e.g.* for electroencephalography (EEG)
or magneto-encephalography (MEG), and
paleomagnetism (recovery of rocks magnetization)
, , see Section .
In this connection, we collaborate with two groups of partners:
Athena Inria project-team,
CHU La Timone, and BESA company on the one hand,
Geosciences Lab. at MIT and Cerege CNRS Lab. on the other hand.
The question is considerably more difficult than its 2-D
counterpart, due mainly to the lack of multiplicative structure for harmonic
gradients. Still,
substantial progress has been made over the last years
using methods of harmonic analysis and operator theory.

The team is further concerned with 3-D generalizations and applications to
non-destructive control of step 2 in Section .
A typical problem is here to localize inhomogeneities or defaults such as
cracks, sources or occlusions in a planar or 3-dimensional object,
knowing thermal, electrical, or
magnetic measurements on the boundary.
These defaults can be expressed as a lack of harmonicity
of the solution to the associated Dirichlet-Neumann problem,
thereby posing an inverse potential problem in order to recover them.
In 2-D, finding an optimal discretization of the
potential in Sobolev norm amounts to solve a best rational approximation
problem, and the question arises as to how the location of the
singularities of the approximant (*i.e.* its poles)
reflects the location of the singularities of the potential
(*i.e.* the defaults we seek). This is a fairly deep issue
in approximation theory, to which Apics contributed convergence results
for certain classes of fields
expressed as Cauchy integrals over extremal contours for
the logarithmic potential
, , .
Initial schemes to locate cracks or sources
*via* rational approximation on
planar domains were obtained this way , , . It is remarkable that finite inverse source problems
in 3-D balls, or more general algebraic surfaces,
can be approached using these 2-D techniques upon slicing the
domain into planar sections
, .
More precisely, each section cuts out a planar domain, the boundary of which
carries data which can be proved to match an algebraic function. The
singularities of this algebraic function are not located at the 3-D sources,
but are related to them: the section contains a source if and only if some
function of the singularities in that section meets a relative extremum. Using
bisection it is thus possible to determine an extremal place along all sections
parallel to a given plane direction, up to some threshold which has to be
chosen small enough that one does not miss a source. This way, we reduce the
original source problem in 3-D to a sequence of inverse poles and branchpoints
problems in 2-D.
This bottom line generates a steady research activity
within Apics, and again applications are sought to medical imaging and
geosciences, see Sections ,
and .

Conjectures may be raised on the behavior of optimal potential discretization in 3-D, but answering them is an ambitious program still in its infancy.

Through contacts with CNES (French space agency),
members of the team became involved in identification and tuning
of microwave electromagnetic filters used in space telecommunications,
see Section . The initial problem was
to recover, from band-limited frequency measurements, physical
parameters of the device under examination.
The latter consists of interconnected dual-mode resonant cavities with
negligible loss, hence its scattering matrix is modeled by a

This is where system theory comes into play, through the
so-called *realization* process mapping
a rational transfer function in the frequency domain
to a state-space representation of the underlying system
of linear differential equations in the time domain.
Specifically, realizing the scattering matrix
allows one to construct
a virtual electrical network, equivalent to the filter,
the parameters of which mediate in between the frequency response
and the
geometric characteristics of the cavities (*i.e.* the tuning parameters).

Hardy spaces provide a framework to transform this ill-posed issue into a series of regularized analytic and meromorphic approximation problems. More precisely, the procedure sketched in Section goes as follows:

infer from the pointwise boundary data in the bandwidth
a stable transfer function (*i.e.* one which is holomorphic
in the right half-plane), that may be infinite dimensional
(numerically: of high degree). This is done by solving
a problem analogous to

A stable rational approximation of appropriate degree to the model obtained in the previous step is performed. For this, a descent method on the compact manifold of inner matrices of given size and degree is used, based on an original parametrization of stable transfer functions developed within the team , .

Realizations of this rational approximant are computed. To be useful, they must satisfy certain constraints imposed by the geometry of the device. These constraints typically come from the coupling topology of the equivalent electrical network used to model the filter. This network is composed of resonators, coupled according to some specific graph. This realization step can be recast, under appropriate compatibility conditions , as solving a zero-dimensional multivariate polynomial system. To tackle this problem in practice, we use Gröbner basis techniques and continuation methods which team up in the Dedale-HF software (see Section ).

Let us mention that extensions of classical coupling matrix theory to frequency-dependent (reactive) couplings have been carried-out in recent years for wide-band design applications.

Apics also investigates issues pertaining to
design rather than identification.
Given the topology of the filter,
a basic problem in this connection is to find the optimal response
subject to specifications
that bear on rejection, transmission and group delay of the
scattering parameters.
Generalizing the classical approach based on Chebyshev polynomials
for single band
filters, we recast the problem of multi-band response synthesis
as a generalization of the classical Zolotarev min-max problem
for rational functions .
Thanks to quasi-convexity, the latter
can be solved efficiently using iterative methods relying on linear
programming. These were implemented in the software
easy-FF (see easy-FF). Currently, the team is engaged
in the synthesis of more complex microwave devices
like multiplexers and routers, which connect several
filters through wave guides.
Schur analysis plays an important role here, because
scattering matrices of passive systems are of Schur type
(*i.e.* contractive in the stability region).
The theory originates with the work of I. Schur ,
who devised a recursive test to
check for contractivity of a holomorphic function in the disk.
The so-called Schur parameters of a function
may be viewed as Taylor coefficients for the hyperbolic metric of the disk, and
the fact that Schur functions are contractions for that metric lies at the
root of Schur's test.
Generalizations thereof turn out to be efficient to parametrize
solutions to contractive interpolation problems .
Dwelling on this, Apics contributed
differential parametrizations (atlases of charts) of lossless
matrix functions , , which
are fundamental to our rational approximation
software RARL2 (see Section ).
Schur analysis is also instrumental to approach de-embedding issues,
and provides one with considerable
insight into the so-called matching problem. The latter consists in
maximizing the power a multiport can pass to a given load, and for
reasons of efficiency it
is all-pervasive in microwave and electric network design, *e.g.* of
antennas, multiplexers, wifi cards and more. It can be viewed as a
rational approximation problem in the hyperbolic metric, and the team
presently deals with this hot topic using
contractive interpolation with constraints on boundary peak points,
within the framework of the (defense funded) ANR Cocoram,
see Sections and .

In recent years,
our attention was driven by CNES and UPV (Bilbao)
to questions about stability of high-frequency amplifiers.
Contrary to previously discussed devices, these are *active* components.
The response of an amplifier can be linearized around a
set of primary current and voltages,
and then admittances of the corresponding electrical network
can be computed at various frequencies, using the so-called harmonic
balance method.
The initial goal is to check for stability of the linearized model,
so as to ascertain existence of a well-defined working state.
The network is composed of lumped electrical elements namely
inductors, capacitors, negative *and* positive reactors,
transmission lines, and controlled current sources.
Our research so far has focused on describing the algebraic structure
of admittance functions, so as to set up a function-theoretic framework
where the two-steps approach outlined in Section
can be put to work. The main discovery is that
the unstable part of each partial transfer function is rational and can
be computed by analytic projection,
see Section . We now start investigating the
linearized
harmonic transfer-function around a periodic cycle, to check for stability
under non necessarily small inputs. This generalization
generates both doctoral and postdoctoral work by new students in the team.

To find an analytic function

Here *a priori*
assumptions on
the behavior of the model off

To fix terminology, we refer to *bounded extremal problem*.
As shown in , ,
,
the solution to this convex
infinite-dimensional optimization problem can be obtained
when

(

The case

Various modifications of

The analog of Problem *seek the inner
boundary*, knowing it is a level curve of the solution.
In this case, the Lagrange parameter indicates
how to deform the inner contour in order to improve
data fitting.
Similar topics are discussed in Section for more general equations than the Laplacian, namely
isotropic conductivity equations of the form

Though originally considered in dimension 2,
Problem

When

On the ball, the analog
of Problem

When *Hardy-Hodge* decomposition,
allowing us to express a *i.e.* those generating no field
in the upper half space) .

Just like solving problem

Problem

Companion to problem

Note that

The techniques set forth in this section are used to solve
step 2 in Section and they are instrumental to
approach inverse boundary value problems
for the Poisson equation

We put

A natural generalization of problem

(

Only for

The case where *stable* rational
approximant to *not* be unique.

The former Miaou project (predecessor of Apics) designed a dedicated
steepest-descent algorithm
for the case *local minimum* is
guaranteed; until now it seems to be the only procedure meeting this
property. This gradient algorithm proceeds
recursively with respect to *critical points* of lower degree
(as is done by the RARL2 software, Section ).

In order to establish global convergence results, Apics has undertaken a
deeper study of the number and nature of critical points
(local minima, saddle points, ...), in which
tools from differential topology and
operator theory team up with classical interpolation theory
, .
Based on this work,
uniqueness or asymptotic uniqueness of the approximant
was proved for certain classes of functions like
transfer functions of relaxation
systems (*i.e.*
Markov functions) and more
generally Cauchy integrals over hyperbolic geodesic arcs .
These are the only results of this kind. Research by Apics on this topic
remained dormant for a while by reasons of opportunity,
but revisiting the work in higher dimension is
a worthy and timely endeavor today. Meanwhile,
an analog to AAK theory
was carried out for

A common
feature to the above-mentioned problems
is that critical point equations
yield non-Hermitian orthogonality relations for the denominator
of the approximant. This stresses connections with interpolation,
which is a standard way to build approximants,
and in many respects best or near-best rational approximation
may be regarded as a clever manner to pick interpolation points.
This was exploited in , ,
and is used in an essential manner to assess the
behavior of poles of best approximants to functions with branched
singularities,
which is of particular interest for inverse source problems
(*cf.* Sections
and ).

In higher dimensions, the analog of Problem

Besides,
certain constrained rational approximation problems, of special interest
in identification
and design of passive systems, arise when putting additional
requirements on the approximant, for instance that it should be smaller than 1
in modulus (*i.e.* a Schur function). In particular, Schur interpolation
lately received renewed attention
from the team, in connection with matching problems.
There, interpolation data are subject to
a well-known compatibility condition (positive definiteness of the so-called
Pick matrix), and the main difficulty is to put interpolation
points on the boundary of

Matrix-valued approximation is necessary to handle systems with several
inputs and outputs but it generates additional difficulties
as compared to scalar-valued approximation,
both theoretically and algorithmically. In the matrix case,
the McMillan degree (*i.e.* the degree of a minimal realization in
the System-Theoretic sense) generalizes the usual notion of degree
for rational functions. For instance when poles are simple, the McMillan degree is the sum of the ranks of the residues.

The basic problem that we consider now goes as follows:
*let $\mathcal{F}\in {\left({H}^{2}\right)}^{m\times l}$ and $n$ an
integer; find a rational matrix of size $m\times l$ without
poles in the unit disk and of McMillan degree at most $n$ which is nearest possible
to $\mathcal{F}$ in ${\left({H}^{2}\right)}^{m\times l}$.*
Here the

The scalar approximation algorithm derived in
and mentioned in
Section
generalizes to
the matrix-valued situation . The
first difficulty here is to parametrize
inner matrices (*i.e.* matrix-valued functions
analytic in the unit disk and unitary on the unit circle) of
given McMillan degree degree

Difficulties relative to multiple local minima of course arise in
the matrix-valued case as well, and deriving criteria that
guarantee uniqueness is even
more difficult than in the scalar case. The case of rational functions
of degree

Let us stress that RARL2 seems the only algorithm handling rational approximation in the matrix case that demonstrably converges to a local minimum while meeting stability constraints on the approximant. It is still a working pin of many developments by Apics on frequency optimization and design.

We refer here to the behavior of poles of best
meromorphic approximants, in the

Generally speaking in approximation theory, assessing the
behavior of poles of rational approximants is essential
to obtain error rates as the degree goes large, and to tackle
constructive issues like
uniqueness. However, as explained in Section ,
the original twist by Apics is to consider this issue also as a means
to extract information on
singularities of the solution to a
Dirichlet-Neumann problem.
The general theme is thus: *how do the singularities
of the approximant reflect those of the approximated function?*
This approach to inverse problem for the 2-D Laplacian turns out
to be attractive when singularities
are zero- or one-dimensional (see Section ). It can be used
as a computationally cheap
initial condition for more precise but much heavier
numerical optimizations which often do not even converge
unless properly initialized.
As regards crack detection or source recovery, this approach
boils down to
analyzing the behavior of best meromorphic
approximants of given pole cardinality to a function with branch points, which is the prototype of
a polar singular set.
For piecewise analytic cracks, or in the case of sources, we were able to
prove (, , ),
that the poles of the
approximants accumulate, when the degree goes large,
to some extremal cut of minimum weighted
logarithmic capacity connecting
the singular points of the crack, or the sources
.
Moreover, the asymptotic density
of the poles turns out to be the Green equilibrium distribution
on this cut in

The case of two-dimensional singularities is still an outstanding open problem.

It is remarkable that inverse source problems inside a sphere or an ellipsoid in 3-D can be approached with such 2-D techniques, as applied to planar sections, see Section . The technique is implemented in the software FindSources3D, see Section .

In addition to the above-mentioned research activities, Apics develops and maintains a number of long-term software tools that either implement and illustrate effectiveness of the algorithms theoretically developed by the team or serve as tools to help further research by team members. We present briefly the most important of them.

Scientific Description

Dedale-HF consists in two parts: a database of coupling topologies as well as a dedicated predictor-corrector code. Roughly speaking each reference file of the database contains, for a given coupling topology, the complete solution to the coupling matrix synthesis problem (C.M. problem for short) associated to particular filtering characteristics. The latter is then used as a starting point for a predictor-corrector integration method that computes the solution to the C.M. corresponding to the user-specified filter characteristics. The reference files are computed off-line using Gröbner basis techniques or numerical techniques based on the exploration of a monodromy group. The use of such continuation techniques, combined with an efficient implementation of the integrator, drastically reduces the computational time.

Dedale-HF has been licensed to, and is currently used by TAS-Espana

Functional Description

Dedale-HF is a software dedicated to solve exhaustively the coupling matrix synthesis problem in reasonable time for the filtering community. Given a coupling topology, the coupling matrix synthesis problem consists in finding all possible electromagnetic coupling values between resonators that yield a realization of given filter characteristics. Solving the latter is crucial during the design step of a filter in order to derive its physical dimensions, as well as during the tuning process where coupling values need to be extracted from frequency measurements.

Participant: Fabien Seyfert

Contact: Fabien Seyfert

Keywords: Health - Neuroimaging - Visualization - Compilers - Medical - Image - Processing

FindSources3D is a software program dedicated to the resolution of inverse source problems in electroencephalography (EEG). From pointwise measurements of the electrical potential taken by electrodes on the scalp, FindSources3D estimates pointwise dipolar current sources within the brain in a spherical model.

After a first data transmission “cortical mapping” step, it makes use of best rational approximation on 2-D planar cross-sections and of the software RARL2 in order to locate singularities. From those planar singularities, the 3-D sources are estimated in a last step, see .

The present version of FindSources3D (called FindSources3D-bolis) provides a modular, ergonomic, accessible and interactive platform, with a convenient graphical interface and a tool that can be distributed and used, for EEG medical imaging. Modularity is now granted (using the tools dtk, Qt, with compiled Matlab libraries). It offers a detailed and nice visualization of data and tuning parameters, processing steps, and of the computed results (using VTK).

A new version is being developed that will incorporate a first Singular Value Decomposition (SVD) step in order to be able to handle time dependent data and to find the corresponding principal static components.

Participants: Juliette Leblond, Maureen Clerc (team Athena, Inria Sophia), Jean-Paul Marmorat, Théodore Papadopoulo (team Athena).

Contact: Juliette Leblond

URL: http://

Scientific Description

For the matrix-valued rational approximation step, Presto-HF relies on RARL2. Constrained realizations are computed using the Dedale-HF software. As a toolbox, Presto-HF has a modular structure, which allows one for example to include some building blocks in an already existing software.

The delay compensation algorithm is based on the following assumption: far off the pass-band, one can reasonably expect a good approximation of the rational components of S11 and S22 by the first few terms of their Taylor expansion at infinity, a small degree polynomial in 1/s. Using this idea, a sequence of quadratic convex optimization problems are solved, in order to obtain appropriate compensations. In order to check the previous assumption, one has to measure the filter on a larger band, typically three times the pass band.

This toolbox has been licensed to, and is currently used by Thales Alenia Space in Toulouse and Madrid, Thales airborne systems and Flextronics (two licenses). XLIM (University of Limoges) is a heavy user of Presto-HF among the academic filtering community and some free license agreements have been granted to the microwave department of the University of Erlangen (Germany) and the Royal Military College (Kingston, Canada).

Functional Description

Presto-HF is a toolbox dedicated to low-pass parameter identification for microwave filters. In order to allow the industrial transfer of our methods, a Matlab-based toolbox has been developed, dedicated to the problem of identification of low-pass microwave filter parameters. It allows one to run the following algorithmic steps, either individually or in a single stroke:

• Determination of delay components caused by the access devices (automatic reference plane adjustment),

• Automatic determination of an analytic completion, bounded in modulus for each channel,

• Rational approximation of fixed McMillan degree,

• Determination of a constrained realization.

Participants: Fabien Seyfert, Jean-Paul Marmorat and Martine Olivi

Contact: Fabien Seyfert

Réalisation interne et Approximation Rationnelle L2

Scientific Description

The method is a steepest-descent algorithm. A parametrization of MIMO systems is used, which ensures that the stability constraint on the approximant is met. The implementation, in Matlab, is based on state-space representations.

RARL2 performs the rational approximation step in the software tools PRESTO-HF and FindSources3D. It is distributed under a particular license, allowing unlimited usage for academic research purposes. It was released to the universities of Delft and Maastricht (the Netherlands), Cork (Ireland), Brussels (Belgium), Macao (China) and BITS-Pilani Hyderabad Campus (India).

Functional Description

RARL2 is a software for rational approximation. It computes a stable rational L2-approximation of specified order to a given L2-stable (L2 on the unit circle, analytic in the complement of the unit disk) matrix-valued function. This can be the transfer function of a multivariable discrete-time stable system. RARL2 takes as input either:

• its internal realization,

• its first N Fourier coefficients,

• discretized (uniformly distributed) values on the circle. In this case, a least-square criterion is used instead of the L2 norm.

It thus performs model reduction in the first or the second case, and leans on frequency data identification in the third. For band-limited frequency data, it could be necessary to infer the behavior of the system outside the bandwidth before performing rational approximation.

An appropriate Möbius transformation allows to use the software for continuous-time systems as well.

Participants: Jean-Paul Marmorat and Martine Olivi

Contact: Martine Olivi

Keywords: Numerical algorithm - Supremum norm - Curve plotting - Remez algorithm - Code generator - Proof synthesis

Functional Description

Sollya is an interactive tool where the developers of mathematical floating-point libraries (libm) can experiment before actually developing code. The environment is safe with respect to floating-point errors, i.e. the user precisely knows when rounding errors or approximation errors happen, and rigorous bounds are always provided for these errors.

Among other features, it offers a fast Remez algorithm for computing polynomial approximations of real functions and also an algorithm for finding good polynomial approximants with floating-point coefficients to any real function. As well, it provides algorithms for the certification of numerical codes, such as Taylor Models, interval arithmetic or certified supremum norms.

It is available as a free software under the CeCILL-C license.

Participants: Sylvain Chevillard, Christoph Lauter, Mioara Joldes and Nicolas Jourdan

Partners: CNRS - ENS Lyon - UCBL Lyon 1

Contact: Sylvain Chevillard

Application domains are naturally linked to the problems described in Sections and . By and large, they split into a systems-and-circuits part and an inverse-source-and-boundary-problems part, united under a common umbrella of function-theoretic techniques as described in Section .

Generally speaking, inverse potential problems, similar to the one appearing in Section , occur naturally in connection with systems governed by Maxwell's equation in the quasi-static approximation regime. In particular, they arise in magnetic reconstruction issues. A specific application is to geophysics, which led us to form the Inria Associate Team Impinge (Inverse Magnetization Problems IN GEosciences) together with MIT and Vanderbilt University. A recent collaboration with Cerege (CNRS, Aix-en-Provence), in the framework of the ANR-project MagLune, completes this picture, see Section .

To set up the context, recall that the Earth's geomagnetic field is generated by convection of the liquid metallic core (geodynamo) and that rocks become magnetized by the ambient field as they are formed or after subsequent alteration. Their remanent magnetization provides records of past variations of the geodynamo, which is used to study important processes in Earth sciences like motion of tectonic plates and geomagnetic reversals. Rocks from Mars, the Moon, and asteroids also contain remanent magnetization which indicates the past presence of core dynamos. Magnetization in meteorites may even record fields produced by the young sun and the protoplanetary disk which may have played a key role in solar system formation.

For a long time, paleomagnetic techniques were only capable of analyzing bulk samples and compute their net magnetic moment. The development of SQUID microscopes has recently extended the spatial resolution to sub-millimeter scales, raising new physical and algorithmic challenges. The associate team Impinge aims at tackling them, experimenting with the SQUID microscope set up in the Paleomagnetism Laboratory of the department of Earth, Atmospheric and Planetary Sciences at MIT. Typically, pieces of rock are sanded down to a thin slab, and the magnetization has to be recovered from the field measured on a planar region at small distance from the slab.

Mathematically
speaking, both inverse source problems for EEG from Section and inverse magnetization problems described presently
amount to recover the (3-D valued) quantity

outside the volume

Another timely instance of inverse magnetization problems lies with geomagnetism. Satellites orbiting around the Earth measure the magnetic field at many points, and nowadays it is a challenge to extract global information from those measurements. In collaboration with C. Gerhards from the University of Vienna, Apics has started to work on the problem of separating the magnetic field due to the magnetization of the globe's crust from the magnetic field due to convection in the liquid metallic core. The techniques involved are variants, in a spherical context, from those developed within the Impinge associate team for paleomagnetism, see Section .

This work is conducted in collaboration with Maureen Clerc and Théo Papadopoulo from the team Athena (Inria Sophia).

Solving overdetermined Cauchy problems for the Laplace equation on a
spherical layer (in 3-D) in order to extrapolate
incomplete data (see Section ) is
a necessary
ingredient of the team's approach to inverse source problems, in particular
for applications to EEG, see . Indeed, the latter involves propagating the
initial conditions through several layers of different conductivities,
from the boundary shell
down to the center of the domain where the
singularities (*i.e.* the sources) lie.
Once propagated
to the innermost sphere, it turns out that traces of the
boundary data on 2-D cross sections coincide
with analytic functions with branched singularities
in the slicing plane ,
. The singularities are
related to the actual location of the sources, namely their moduli
reach in turn a
maximum when the plane contains one of the sources. Hence we are
back to the 2-D framework of Section ,
and recovering these singularities
can be performed *via* best rational approximation.
The goal is to produce a fast and sufficiently accurate
initial guess on the number
and location of the sources in order to run heavier
descent algorithms on the direct problem, which are more precise but
computationally costly and often
fail to converge if not properly initialized. Our belief
is that such a localization process can add a geometric, valuable piece of
information to the standard temporal analysis of EEG signal records.

Numerical experiments obtained with our software FindSources3D give
very good results on simulated data and we are now engaged in the process
of handling
real experimental data (see
Sections and ),
in collaboration with the team Athena at Inria Sophia Antipolis,
neuroscience teams in partner-hospitals (Institut de Neurosciences des Systèmes, hospital la Timone, Aix-Marseille Univ., Marseille, http://

This is joint work with Stéphane Bila (XLIM, Limoges).

One of the best training grounds for function-theoretic applications by the team is the identification and design of physical systems whose performance is assessed frequency-wise. This is the case of electromagnetic resonant systems which are of common use in telecommunications.

In space telecommunications (satellite transmissions), constraints specific to on-board technology lead to the use of filters with resonant cavities in the microwave range. These filters serve multiplexing purposes (before or after amplification), and consist of a sequence of cylindrical hollow bodies, magnetically coupled by irises (orthogonal double slits). The electromagnetic wave that traverses the cavities satisfies the Maxwell equations, forcing the tangent electrical field along the body of the cavity to be zero. A deeper study of the Helmholtz equation states that an essentially discrete set of wave vectors is selected. In the considered range of frequency, the electrical field in each cavity can be decomposed along two orthogonal modes, perpendicular to the axis of the cavity (other modes are far off in the frequency domain, and their influence can be neglected).

Near the resonance frequency, a good approximation to the Helmholtz equations is given by a second order differential equation. Thus, one obtains an electrical model of the filter as a sequence of electrically-coupled resonant circuits, each circuit being modeled by two resonators, one per mode, the resonance frequency of which represents the frequency of a mode, and whose resistance accounts for electric losses (surface currents) in the cavities.

This way, the filter can be seen as a quadripole, with two ports, when
plugged onto a resistor at one end and fed with some potential at the other end.
One is now
interested in the power which is transmitted and reflected. This leads
one to define a
scattering matrix

In fact, resonance is not studied via the electrical model,
but via a low-pass
equivalent circuit obtained upon linearizing near the central frequency, which is no
longer
conjugate symmetric (*i.e.* the underlying system may no longer
have real
coefficients) but whose degree is divided by 2 (8 in the example).

In short, the strategy for identification is as follows:

measuring the scattering matrix of the filter near the optimal frequency over twice the pass band (which is 80MHz in the example).

Solving bounded extremal problems for the transmission and the reflection (the modulus of he response being respectively close to 0 and 1 outside the interval measurement, cf. Section ) in order to get a models for the scattering matrix as an analytic matrix-valued function. This provides us with a scattering matrix known to be close to a rational matrix of order roughly 1/4 of the number of data points.

Approximating this scattering matrix by a true rational transfer-function of appropriate degree (8 in this example) via the Endymion or RARL2 software (cf. Section ).

A state space realization of

Finally one builds a realization of the approximant and looks for a change of variables that eliminates non-physical couplings. This is obtained by using algebraic-solvers and continuation algorithms on the group of orthogonal complex matrices (symmetry forces this type of transformation).

The final approximation is of high quality. This can be interpreted as
a confirmation of the linearity assumption on the system:
the relative

The above considerations are valid for a large class of filters. These developments have also been used for the design of non-symmetric filters, which are useful for the synthesis of repeating devices.

The team further investigates problems relative to the design of optimal responses for microwave devices. The resolution of a quasi-convex Zolotarev problems was proposed, in order to derive guaranteed optimal multi-band filter responses subject to modulus constraints . This generalizes the classical single band design techniques based on Chebyshev polynomials and elliptic functions. The approach relies on the fact that the modulus of the scattering parameter

The filtering function appears to be the ratio of two polynomials

The relative simplicity of the derivation of a filter's response, under
modulus constraints, owes much to the possibility of
forgetting about Feldtkeller's equation and express all design constraints
in terms of the filtering function. This no longer the case when
considering the synthesis

Through contacts with CNES (Toulouse) and UPV (Bilbao),
Apics got additionally involved
in the design of amplifiers which, unlike filters, are active devices.
A prominent issue here is stability. A twenty years back, it was not
possible to simulate unstable responses, and only after building a device
could one detect instability. The advent of so-called *harmonic balance*
techniques, which compute steady state responses of linear elements in
the frequency domain and look for a periodic state in the time domain of
a network connecting these linear elements *via*
static non-linearities made it possible to compute the harmonic response
of a (possibly nonlinear and unstable) device .
This has had tremendous impact on
design, and there is a growing demand for software analyzers.
The team is also becoming active in this area.

In this connection, there are two types of stability involved. The first is stability of a fixed
point around which the linearized transfer function
accounts for small signal amplification. The second is stability of a
limit cycle which is reached when the input signal is no longer small
and truly nonlinear amplification is attained
(*e.g.* because of saturation).
Work by the team so far has been concerned with the first type of stability,
and emphasis is put on defining and extracting the “unstable part” of the response, see Section . The stability check for
limit cycles is now under investigation.

This work is carried out in the framework of the Inria Associate Team Impinge, comprising Cauê Borlina, Eduardo Andrade Lima and Benjamin Weiss from the Earth Sciences department at MIT (Boston, USA) and Douglas Hardin, Edward Saff and Cristobal Villalobos from the Mathematics department at Vanderbilt University (Nashville, USA).

The overall goal of Impinge is to determine
magnetic properties of rock
samples (*e.g.* meteorites or stalactites), from weak field measurements
close to the sample that
can nowadays be obtained using SQUIDs (superconducting quantum interference
devices).
Depending on the nature of the rock sample, the magnetization distribution can either be considered to lie in a plane or in a parallelepiped of thickness

We pursued this year our research efforts towards designing algorithms for net moment recovery. The net moment is the integral of the magnetization over its support, and it is a valuable piece of information to physicists which has the advantage of being determined solely by the field: whereas two different magnetizations can generate the same field, the net moment depends only on the field and not on which magnetization produced it. Hence the goal may be described as to build a numerical magnetometer, capable of analyzing data close to the sample. This is in contrast to classical magnetometers which regard the latter as a single dipole, an approximation which is only valid away from the sample and is not suitable to handle weak fields which get quickly blurred by ambient magnetic sources. This research effort was paid in two different, complementary directions.

The first approach consists in using the fact that the integral of

The second approach attempts to generalize the previous expansions in the case when *e.g.*, the level of constraint

We started this year to design an alternate procedure to compute a good linear estimator. It consists in expanding it on a family of piecewise affine functions, with a restricted number of pieces. This still needs to be pushed further in connection with the delicate issue of how dense should the grid of data points be in order to reach a prescribed level of precision. On a related topic, we also derived explicit formulas for the adjoint operator

Concerning full inversion of thin samples, after preliminary experiments on regularization with

Besides, we considered a simplified 2-D setup for magnetizations and magnetic potentials (of which the magnetic field is the gradient). When both the sample and the measurement set are parallel intervals, some best approximation issues related to inverse recovery and relevant BEP problems in Hardy classes of holomorphic functions (see Section ) were solved in . Note that, in the present case, the criterion no longer acts on the boundary of the holomorphy domain (namely, the upper half-plane), but on a strict subset thereof, while the constraint acts on the support of the approximating function. Both involve real parts of functions in the Hilbert Hardy space of the upper half-plane. This is joint work with D. Ponomarev (see Section ). Some extensions are the subject of ongoing work with E. Pozzi (Department of Mathematics and Statistics, St Louis Univ., St Louis, Missouri, USA). They concern more precise approximation criteria, and the development of resolution schemes using the Fourier basis. Meanwhile, BEP in Bergman classes of analytic or generalized analytic functions are under being studied with B. Delgado Lopez (see Sections , ).

For magnetizations supported in a volume

Other types of inverse magnetization problems can be tackled using such techniques, in particular global Geomagnetic issues which arise in spherical geometry. This year, in collaboration with C. Gerhards from the University of Vienna (Austria), we developed a method to separate the crustal component of the Earth's magnetic field from its core component, if an estimate of the field is known on a subregion of the globe . This assumption is not unrealistic: parts of Australia and of northern Europe are considered as fairly well understood from the magnetostatic view point. We look forward to test the algorithm against real data, in collaboration with Geophysicists

The team Apics is a partner of the ANR project MagLune on Lunar magnetism, headed by the Geophysics and Planetology Department of Cerege, CNRS, Aix-en-Provence (see Section ). Recent studies let geoscientists think that the Moon used to have a magnetic dynamo for a while. However, the exact process that triggered and fed this dynamo is still not understood, much less why it stopped. The overall goal of the project is to devise models to explain how this dynamo phenomenon was possible on the Moon.

The geophysicists from Cerege went a couple of times to NASA to perform measurements on a few hundreds of samples brought back from the Moon by Apollo missions. The samples are kept inside bags with a protective atmosphere, and geophysicists are not allowed to open the bags, nor to take out samples from NASA facilities. Moreover, the process must be carried out efficiently as a fee is due to NASA by the time when handling these moon samples. Therefore, measurements were performed with some specific magnetometer designed by our colleagues from Cerege. This device measures the components of the magnetic field produced by the sample, at some discrete set of points located on circles belonging to three cylinders (see Figure ). The objective of Apics is to enhance the numerical efficiency of post-processing data obtained with this magnetometer.

This year, we continued the approach taken in previous years. Under the hypothesis that the field can be well explained by a single magnetic pointwise dipole, and using ideas similar to those underlying the FindSources3D tool (see Sections and ), we try to recover the position and the moment of the dipole using the available measurements.

In a given cylinder, using the associated cylindrical system of coordinates, recovering the position of the dipole boils down to determine its height *e.g.*, by taking the mean value of the argument of

The first strategy consists in observing that

The second strategy directly uses the pole

This work is conducted in collaboration with Maureen Clerc and Théo Papadopoulo, from the team Athena (Inria Sophia).

In 3-D, functional or clinically active regions in the cortex are often modeled by pointwise sources that have to be localized from measurements, taken by electrodes on the scalp, of an electrical potential satisfying a Laplace equation (EEG, electroencephalography). In the works , on the behavior of poles in best rational approximants of fixed degree to functions with branch points, it was shown how to proceed via best rational approximation on a sequence of 2-D disks cut along the inner sphere, for the case where there are finitely many sources (see Section ).

In this connection, a dedicated software FindSources3D (FS3D, see Section )
is being developed, in collaboration with the Inria team Athena and the CMA - Mines ParisTech. In addition to the modular and ergonomic platform version of FS3D,
a new (Matlab) version of the software that automatically performs the estimation of the quantity of sources is being built.
It uses an alignment criterion in addition to other clustering tests for the selection.
It appears that, in the rational approximation step,
*multiple* poles possess a nice behavior with respect to branched
singularities. This is due to the very physical assumptions on the model
(for EEG data that correspond to measurements of the electrical potential, one should consider *triple* poles; for (magnetic) field data however, like in Section or from MEG – magneto-encephalography – data, one should consider poles of order five). Though numerically
observed in , there is no mathematical
justification so far why multiple poles generate such strong accumulation
of the poles of the approximants. This intriguing property, however,
is definitely helping source recovery. It is used in order to automatically estimate the “most plausible”
number of sources (numerically: up to 3, at the moment).
Last but not least, the version of the software currently under development
takes as inputs actual EEG measurements, like time signals, and performs a suitable singular value decomposition in order to separate independent sources.

Magnetic data from MEG recently became available along with EEG data, by our medical partners at the hospital la Timone; indeed, it is now possible to use simultaneously both measurement devices, in order to measure both the electrical potential and a component of the magnetic fields. This should enhance the accuracy of our source recovery algorithms. We will add the treatment of MEG data as another functionality of the software FS3D.

In connection with these and other brain exploration modalities like electrical impedance tomography (EIT), we are now studying conductivity estimation problems. This is the topic of the PhD research work of C. Papageorgakis (co-advised with the Inria team Athena and BESA GmbH, see Section ). In layered models, it concerns the estimation of the conductivity of the skull (an intermediate layer). First, the conductivity of the skull can differ from one individual to another, or for the same person, along the time, and is much smaller than those of the surrounding layers (the brain and the scalp). A preliminary issue in this direction was to estimate a single-valued skull conductivity from one EEG recording. Existence, uniqueness, stability properties and a recovery scheme for this conductivity value were established in the spherical setting when the sources are known, see Part I. When the sources are unknown, we must look for additional data (additional clinical and/or functional EEG, EIT, ...) that could be incorporated in order to recover both the sources locations and the skull conductivity. Second, while the skull essentially consists of a hard bone part, which may be assumed to have constant electrical conductivity, it also contains spongy bone compartments. These two distinct components of the skull actually possess quite different conductivities. The influence of the second on the overall model is also studied in Part II, together with a numerical process allowing to estimate the hard bone conductivity value together with a dipolar source, in realistic geometries.

We also began to consider the inverse problem of recovering the parameters of a skin tumor from thermal measurements, in a 2-D model that takes the form of a static Schrödinger equation. This is joint work with F. Ferranti (IMT Atlantique, Microwave Department) and the topic of the internship of G. Dervaux, see Section .

This is collaborative work with Stéphane Bila (XLIM, Limoges, France), Yohann Sence (XLIM, Limoges, France), Thierry Monediere (XLIM, Limoges, France), Francois Torrès (XLIM, Limoges, France) in the context of the ANR Cocoram (see Section ) as well as with, Fabien Ferrero (LEAT, Sophia-Antipolis, France) Leonardo Lizzi (LEAT, Sophia-Antipolis, France).

Filter synthesis is usually performed under the hypothesis that both ports of the filter are loaded on a constant resistive load (usually 50 Ohm). In complex systems, filters are however cascaded with other devices, and end up being loaded, at least at one port, on a non purely resistive frequency varying load. This is for example the case when synthesizing a multiplexer: each filter is here loaded at one of its ports on a common junction. Thus, the load varies with frequency by construction, and is not purely resistive either. Likewise, in an emitter-receiver, the antenna is followed by a filter. Whereas the antenna can usually be regarded as a resistive load at some frequencies, this is far from being true on the whole pass-band. A mismatch between the antenna and the filter, however, causes irremediable power losses, both in emission and transmission. Our goal is therefore to develop a method for filter synthesis that allows us to match varying loads on specific frequency bands, while enforcing some rejection properties away from the pass-band.

The matching problem of minimizing

When the degree

where

which accounts for the losslessness of the filter. The frequencies

The previous interpolation procedure provides us with a matching/rejecting filtering characteristics at a discrete set of frequencies. It can serve as a
starting point for heavier optimization procedures, where the matching and rejection specifications are expressed uniformly over the bandwidth. Although the practical results thus obtained are quite convincing, we have no proof of their global optimality. This has led us to seek alternative
approaches allowing us to assess, at least in simple cases, global optimality of the obtained response. By optimality of a response we mean, as in classical filtering,
a best match of the response in the uniform norm on a given pass-band,
while meeting given rejection constraints on a stop-band. Following the approach of Fano and Youla, we considered the problem of designing a

The goal is here to help design amplifiers, in particular to detect instability at an early stage of the design. Activity in this area is gaining importance with the coming of a doctoral (S. Fueyo) and a postdoctoral (A. Cooman) student along with planned software developments. Application of our work to oscillator design methodologies started recently with Smain Amari from the Royal Military College of Canada (Kingston, Canada).

As opposed to Filters and Antennas, Amplifiers and Oscillators are active components that intrinsically entail a non-linear functioning. The latter is due to the use of transistors governed by electric laws exhibiting saturation effects, and therefore inducing input/output characteristics that are no longer proportional to the magnitude of the input signal. Hence they typically produce non-linear distortions. A central question arising in the design of amplifiers is to assess stability. The latter may be understood around a functioning point when no input but noise is considered, or else around a periodic trajectory when an input signal at a specified frequency is applied. For oscillators, a precise estimation of their oscillating frequency is crucial during the design process. As regards devices devised to operate at relative low frequencies, time domain simulations, based on the integration of the underlying non-linear dynamical system, answers these questions satisfactorily. For complex microwave amplifiers and oscillators, the situation is however drastically different: the time step necessary to integrate the transmission line's dynamical equations (which behave like simple electrical wire at low frequency) becomes so small that simulations are intractable in reasonable time. In addition to this problem, most linear components of these circuits are known through their frequency response, and require therefore a preliminary, numerically unstable step to obtain their impulse response, prior to any time domain simulation.

For all these reasons it is widely preferred to perform the analysis of such systems in the frequency domain. In the case of stability issues around a functioning point, where only small input signals are considered, the stability of the linearized system obtained by a first order approximation of each non-linear dynamic is considered. This is done by means of the analysis of transfer impedance functions computed at some ports of the circuit. We have shown, that under some realistic hypothesis on the building blocks of the circuit, these transfer functions are meromorphic functions of the frequency variable

When stability is studied around a periodic trajectory, which is determined in practice by Harmonic Balance algorithms, linearization yields a linear time varying dynamical system with periodic coefficients and a periodic trajectory thereof. While in finite dimension the stability of such systems is well understood via the Floquet theory, this is no longer the case in the infinite dimensional setting when delays are considered. Dwelling on the theory of retarded systems, S. Fueyo's PhD work has made remarkable progress on this topic by showing that, for certain simple circuits with properly positioned resistors, the monodromy operator is a compact perturbation of a stable operator, and that only finitely many unstable point of its spectrum can occur. A practical application of this result is to generalize the previously described techniques of stability assessment around a functioning point into a stability assessment technique around periodic trajectories. This can be recast in terms of the finiteness of the number of abscissas of unstable poles of the Harmonic Transfer functions of the circuit. It will be of great importance to generalize such considerations to more complex circuits, whose structure is less well understood at present.

(This is joint work with Qian T. and Dang P. from the university of Macao.)
It was proven in previous year that on a smooth compact hypersurface *i.e.*

This contract (reference Inria: 11282) accompanies the PhD of David Martinez and focuses on the development of efficient techniques for the design of matching network tailored for frequency varying loads. Applications of the latter to the design output multiplexers occurring in space applications will be considered.

This is a research agreement between Inria (Apics and Athena teams) and the German company BESA

This is a scientific consulting activity for the start-up company SKAVENJI. The latter develops an innovative and communicative device to facilitate the production and home consumption of small amounts of energy, produced by one or more local sources of renewable energy. Ongoing work consists in designing a simple controller improving the energy efficiency of the energy production while minimizing the number of charge and discharge cycles of the associated battery. The retained control strategy is based on consumption and production profiles.

Contract Provence Alpes Côte d'Azur (PACA) Region - Inria, BDO (no. 2014-05764) funding the research grant of C. Papageorgakis, see Sections , .

The team participates in the project WIMAG (Wave IMAGing) funded by the Idex

The team co-advises a PhD (G. Bose) with the CMA team of LEAT (http://

The team participates in the transverse action C4PO funded by the Idex

The team also participates in the project ToMaT, “Multiscale Tomography: imaging and modeling ancient materials, technical traditions and transfers”, funded by the Idex

The ANR (Astrid) project Cocoram (Co-design et co-intégration de réseaux d’antennes actives multi-bandes pour systèmes de radionavigation par satellite) started January 2014. We are associated with three other teams from XLIM (Limoges University), geared respectively towards filters, antennas and amplifiers design. The core idea of the project is to realize dual band reception an emission chains by co-conceiving the antenna, the filters, and the amplifier. We are specifically in charge of the theoretical design of the filters, matching the impedance of a bi-polarized dual band antenna. This is a perfect training ground to test, apply and adapt our work on matching problems (see Section ).

The ANR project MagLune (Magnétisme de la Lune) has been approved July 2014. It involves the Cerege (Centre de Recherche et d’Enseignement de Géosciences de l’Environnement, joint laboratory between Université Aix-Marseille, CNRS and IRD), the IPGP (Institut de Physique du Globe de Paris) and ISTerre (Institut des Sciences de la Terre). Associated with Cerege are Inria (Apics team) and Irphe (Institut de Recherche sur les Phénomènes Hors Équilibre, joint laboratory between Université Aix-Marseille, CNRS and École Centrale de Marseille). The goal of this project (led by geologists) is to understand the past magnetic activity of the Moon, especially to answer the question whether it had a dynamo in the past and which mechanisms were at work to generate it. Apics participates in the project by providing mathematical tools and algorithms to recover the remanent magnetization of rock samples from the moon on the basis of measurements of the magnetic field it generates. The techniques described in Section are instrumental for this purpose.

Apics is part of the European Research Network on System Identification (ERNSI) since 1992.

System identification deals with the derivation, estimation and validation of mathematical models of dynamical phenomena from experimental data.

Title: Inverse Magnetization Problems IN GEosciences.

International Partner (Institution - Laboratory - Researcher):

Massachusetts Institute of Technology (United States) - Department of Earth, Atmospheric and Planetary Sciences - Benjamin P. Weiss

Start year: 2016

See also: http://

The associate team Impinge is concerned with the inverse problem of recovering a magnetization distribution from measurements of the magnetic field above rock slabs using a SQUID microscope (developed at MIT). The application domain is to Earth and planetary sciences. Indeed, the remanent magnetization of rocks provides valuable information on their history. This is a renewal of the previous Associate Team Impinge that ended 2015. The US team also involves a group of Mathematicians at Vanderbilt University (see Section ).

**MIT-France seed funding** is a competitive collaborative research
program ran
by the Massachusetts Institute of Technology (Cambridge, Ma, USA). Together with
E. Lima and B. Weiss from the Earth and Planetary Sciences dept. at MIT,
Apics obtained two-years support from the above-mentioned program to run a project entitled:
“Development of Ultra-high Sensitivity Magnetometry for Analyzing Ancient Rock Magnetism”

**NSF Grant** L. Baratchart, S. Chevillard and J. Leblond are
external investigators in the NSF Grant 2015-2018,
“Collaborative Research: Computational
methods for ultra-high sensitivity magnetometry of geological samples”
led by E.B. Saff (Vanderbilt Univ.) and B. Weiss (MIT).

Cauê Borlina (MIT, Boston, Massachusetts, USA, Apr. 24-28).

Nattapong Bosuwan (Mahidol University, Bangkok, Thailand, May-Aug.).

Briceyda Delgado Lopez (PhD student, Cinvestav, Queretaro, Mexico, Jan.-Mar.).

Bernard Hanzon (Univ. Cork, Ireland, Apr.-Jun.).

Douglas Hardin (Vanderbilt University, Nashville, Tennessee, USA, Apr. 24-28).

Eduardo Lima (MIT, Boston, Massachusetts, USA, Apr. 24-28).

Mateusz Rusiniak (BESA GmbH, Gräfelfing, Germany, Dec. 15).

Ioannis Stratis (National and Kapodistrian University of Athens, Greece, Dec. 14-15).

Carsten Wolters (University of Münster, Germany, Dec. 14-15).

Gautier Dervaux (IMT Atlantique, Brest, France, Jul.-Aug.).

L. Baratchart spent the fall semester 2017 at Vanderbilt University, Nashville, Tennessee, teaching a course on inverse problems and pursuing research with with D. Hardin, E.B. Saff and C. Villalobos, as well as E. Lima, all members of the Inria associate team Impinge

Collaboration under contract with Thales Alenia Space (Toulouse, Cannes, and Paris), CNES (Toulouse), XLIM (Limoges), University of Bilbao (Universidad del País Vasco / Euskal Herriko Unibertsitatea, Spain), BESA company (Munich), Flextronics.

Regular contacts with research groups at UST (Villeneuve d'Asq), Universities of Bordeaux-I (Talence), Orléans (MAPMO), Aix-Marseille (CMI-LATP), Nice Sophia Antipolis (Lab. JAD), Grenoble (IJF and LJK), Paris 6 (P. et M. Curie, Lab. JLL), Inria Saclay (Lab. Poems, ENSTA), IMT Atlantique (Institut Mines-Télécom., Brest), Cerege-CNRS (Aix-en-Provence), CWI (the Netherlands), MIT (Boston, USA), Vanderbilt University (Nashville USA), Steklov Institute (Moscow), Michigan State University (East-Lansing, USA), Texas A&M University (College Station USA), Indiana University-Purdue University at Indianapolis, St Louis University (St Louis, Missouri, USA), Cinvestav (Queretaro, Mexico), Politecnico di Milano (Milan, Italy), University of Trieste (Italy), RMC (Kingston, Canada), University of Leeds (UK), of Maastricht (the Netherlands), of Cork (Ireland), Vrije Universiteit Brussel (Belgium), TU-Wien and Universität Wien (Austria), TFH-Berlin (Germany), ENIT (Tunis), KTH (Stockholm), University of Cyprus (Nicosia, Cyprus), University of Macau (Macau, China), SIAE Microelettronica (Milano).

The project is involved in the GDR-project AFHP (CNRS), in the ANR (Astrid program) project Cocoram (with XLIM, Limoges, and DGA), in the ANR (Défis de tous les savoirs program) project MagLune (with Cerege, IPGP, ISTerre, Irphe), in a MIT-France collaborative seed funding, in the Associate Inria Team Impinge (with MIT, Boston), and in a NSF grant (with Vanderbilt University and MIT).

L. Baratchart gave a talk at SEAM 2017 https://

J. Leblond presented oral communications at the Conference “Cerveau.at.UCA”, http://

J.-P. Marmorat presented a poster at the 18th International Symposium on Applied Electromagnetic and Mechanics (ISEM), http://

K. Mavreas presented oral communications at the Conference Tendances des Applications Mathématiques en Tunisie, Algérie, Maroc (TAMTAM),
http://

C. Papageorgakis gave an oral communication at the Conference TAMTAM,
http://

D. Martinez presented a poster at the International Microwave Symposium 2017 https://

G. Bose gave a talk at the International Conference on Electromagnetics in Advanced Applications http://

A. Cooman presented a poster at the Workshop ERNSI 2017, http://

M. Olivi gave a seminar at ONERA, Toulouse, January 12. She presented a poster at the Workshop ERNSI 2017, http://

J. Leblond was the scientific chair of a mini-symposium at the Conference TAMTAM,
http://

K. Mavreas was one of the organizers of the PhD seminar of the Research Center, until Oct.

L. Baratchart was a member of the program committee of “Control and Distributed Parameter Systems”, Bordeaux, 2017.

L. Baratchart is an editor for “Complex Analysis and Operator Theory” (CAOT) and “Constructive Methods and Function Theory” (CMFT).

S. Chevillard was a reviewer for the journal *Transactions on Computers*.

J. Leblond was a reviewer for the journals *Czechoslovak Mathematical Journal, Annals of West University of Timisoara - Mathematics and Computer Science, ESAIM: Control, Optimisation and Calculus of Variations, Adv. in App. Clifford Algebras, ARIMA, Numerical Algorithms*.

F. Seyfert was a reviewer for the journal *IEEE Microwave Theory and Techniques*.

M. Olivi was a reviewer for the journal *Automatica*.

L. Baratchart was an invited speaker at the conference “Complex and Functional Anaysis and teir interactions with Harmonic Analysis”, Bedlewo,
Poland (June)
https://

J. Leblond was an invited speaker at
the 17ème Journée Calcul Scientifique et Modélisation Mathématique d’Amiens, http://

M. Olivi was invited to give a communication at the symposium “Mathematics in Knowledge Engineering and Data Science - Identifying Connections”, June 30, in Maastricht, the Netherlands.

F. Seyfert is a member of the IEEE MTT-8 Technical Committee on Filters and Passive Components

L. Baratchart sits on the committee “Mathématiques et Informatique” of the French Agency for research (ANR).

J. Leblond is an elected member of the “Conseil Scientifique” and of the “Commission Administrative Paritaire” of Inria. Until April, she was a member of the “Conseil Académique” of the Univ. Côte d'Azur (UCA).

**Graduate course**: L. Baratchart gave a graduate course titled “Introduction to Inverse Problems” (40 hours) at Vanderbilt University,
Nashville, Tennessee.

**Colles**: S. Chevillard is giving “Colles” at Centre International de Valbonne (CIV) (2 hours per week).

**Summer School**: M. Olivi gave two lectures at the 38-th summer school of automatic control, Grenoble, September 11-15:
“Optimization-based model reduction” with P. Vuillemin (Onera) and “Applications in electronics: design of frequency filters and stability analysis of amplifiers”.

PhD: C. Papageorgakis, *Conductivity model estimation*, since October 2014 (advisors: J. Leblond, M. Clerc, M. Rusiniak), defended Dec. 15th .

PhD in progress: K. Mavreas, *Inverse source problems in planetary sciences: dipole localization in Moon rocks from sparse magnetic data*, since October 2015 (advisors: S. Chevillard, J. Leblond).

PhD in progress: D. Martinez Martinez, *Méthodologie et Outils de Synthèse pour des Fonctions de Filtrage Chargées par des Impédances complexes*, since October 2015, advisors: L. Baratchart and F. Seyfert.

PhD in progress: G. Bose, Filter Design to Match Antennas, since December 2016, advisors: F. Ferrero and F. Seyfert.

PhD in progress: S. Fueyo, Cycles limites et stabilité dans les circuits,since October 2016, advisors: L. Baratchart and J.-B. Pomet.

L. Baratchart was a referee of the “Mémoire d'habilitation” by Nicolas Brisebarre (ENS de Lyon).

J. Leblond was a member of the “Jury d'admissibilité du concours CR2” of the Inria Research Center Nancy Grand Ouest and of the “Comité de Sélection” for professors at UNSA (Polytech Nice), in March. She was a member of the PhD thesis defense committee of Lobna Merghmi, Aix-Marseille Université, Institut de Math. de Marseille (I2M), Jan.

M. Olivi was a reviewer for the PhD document of Igor Pontes-Duff-Pereira, université de Toulouse, January 11. She was a member of the PhD thesis defense committees of Afrooz Ebadat, KTH, Stockholm, September 8, and Yusuf Bhujwalla, université de Lorraine, Nancy, December 5. She was a member of the “jury d'admission du concours CR” of Inria.

M. Olivi is responsible for Scientific Mediation and president of the Committee MASTIC (Commission d’Animation et de Médiation Scientifique) https://

management of the contract “région PACA: Science Culture Lycée” and organization of about thirty conferences in the high schools of the region (100 students per conference),

co-organization of 10 robotics sessions for 2 classes of middle school students (device “MEDITES” http://

co-organization of the “stage MathC2+”, a four-day internship for 50 high school students (“secondes”, about 16 years old) organized by the Commitee MASTIC and its partners (June 14-17),

co-organization of Inria participation to the event “Le Village des Sciences et de l’Innovation” in Antibes (October 7 & 8, 10000 people),

co-organization of about 10 “cafés scientifiques” (c@fé-in's and cafés Techno, 30 to 80 participants each),

supervision of a (2 months) internship, done by Sabrina Ballauris, for the realization of objects to manipulate, in view of illustrating some mathematical results and scientific principles (Pythagora's puzzles, Galton's board, Galileo's experiment, the fastest toboggan, ...)

She was a member of the scientific committee of the “Forum Mathématiques Vivantes”, a national event (initiated by CFEM, the French Commission for Mathematics Education, http://

M. Olivi animated two half-day workshop sessions “activités débranchées” at “l'ESPE de Nice” for primary school students (March 16 & 24), 200 students each session). She gave three presentations for high school students in Gap and Marseille. With K. Mavreas, she participates to the event “Le Village des Sciences et de l’Innovation” in Antibes (October 7 & 8, 10000 people).

A. Cooman gave a presentation at the “stage MathC2+”, a four-day internship for high-school students organized by the Committee MASTIC and its partners (June 14-17).