# Recent Progress on the Factorization Method for Electrical Impedance Tomography

Authors: Bastian Harrach
Journal: Computational and Mathematical Methods in Medicine
 (2013)

---

## Abstract

The Factorization Method is a noniterative method to detect the shape and position of conductivity anomalies inside an object. The
method was introduced by Kirsch for inverse scattering problems and extended to electrical impedance tomography (EIT) by Brühl and Hanke. Since these pioneering works, substantial progress has been made on the theoretical foundations of the method. The necessary assumptions have been weakened, and the proofs have been considerably simplified. In this work, we aim to summarize this progress and present a state-of-the-art formulation of the Factorization Method for EIT with continuous data. In particular, we formulate the method for general piecewise analytic conductivities and give short and self-contained proofs.


---

## Body

## 1. Introduction

Electrical impedance tomography (EIT) aims to reconstruct the spatial conductivity distribution inside an imaging subjectΩ⊆ℝn from current-voltage measurements on a part of its surface Σ⊆∂Ω. Mathematically, this leads to the problem of recovering the coefficient σ(x) in the elliptic partial differential equation
(1)∇·σ∇uσg=0inΩ,σ∂νuσg|∂Ω={gonΣ,0else
from knowledge of the corresponding Neumann-to-Dirichlet operator (NtD)
(2)Λ(σ):g⟼uσg|Σ,
where uσg is the solution of  (1). We describe the precise mathematical setting in Section 2.1.In several applications, EIT is used to determine the position of conductivity changes. This includes anomaly detection problems, whereΛ(σ) is compared to a reference NtD Λ(σ0) in order to determine, if and where σ differs from a known background conductivity σ0. This problem also appears in time-difference EIT, where measurements at different times are compared to monitor temporal conductivity changes. These applications lead to the shape reconstruction problem of determining the support of σ-σ0 from Λ(σ) and Λ(σ0).A prominent noniterative shape reconstruction method is the Factorization Method. It was introduced by Kirsch [1] for inverse scattering problems and extended to EIT by Brühl and Hanke [2–4]. In its original form (cf. [4]), the method assumes that
(3)σ0(x)=1,σ(x)=1+κ(x)χD(x),
where D⊆Ω is a union of separated, smoothly bounded, and simply connected domains, on which there is a conductivity jump of at least ϵ>0; that is,
(4)κ(x)≥ϵ∀x∈D,orκ(x)≤-ϵ∀x∈D.
The method then characterizes the unknown shape D by a range criterion. For all unit vectors d∈ℝn, ∥d∥=1,
(5)z∈DiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)|1/2),
where Φz,d is the so-called dipole function, that is, the solution of
(6)ΔΦz,d=d·∇δzinΩ,∂νΦz,d|∂Ω=0.
The range criterion  (5) can be implemented numerically, so that each point z∈Ω can be tested whether it belongs to the unknown inclusion or not.Substantial progress has been made on the Factorization Method since the original works of Kirsch, Brühl, and Hanke. In the following, we restrict ourselves to progress in the context of EIT. Overviews on the FM for EIT have been given by Hanke and Brühl [5], in the book of Kirsch and Grinberg [6], and in a recent chapter of Hanke and Kirsch in Scherzer's Handbook of Mathematical Methods in Imaging [7]. The FM for EIT has been treated as a special case of more general elliptic problems by Kirsch [8], the author [9], and by Nachman, Päivärinta and Teirilä [10]. A half-space setting has been considered by Hanke and Schappel [11]. Electrode models have been covered in the works of Brühl, Lechleiter, Hakula, Hanke, Hyvönen and Pursiainen [5, 12–15]. The FM has been extended to the complex conductivity case arising in frequency-difference EIT by Seo, Woo and the author [16, 17]. A priori separated indefinite inclusions have been treated by Schmitt [18], and Schmitt and Kirsch discussed the determination of the contrast level in the context of the FM for EIT in [19]. Hyvönen and the author have removed the assumptions on the inclusion contrast and boundary regularity in [20], and [21] discusses the relation of the FM to localized potentials.In this work, we aim to summarize the theoretical progress and present a state-of-the-art formulation of the Factorization Method for EIT with continuous data. In particular, we will formulate the method for general piecewise analytic conductivities and give short and self-contained proofs.


## 2. Setting and Auxiliary Results

### 2.1. The Setting

We start by making the mathematical setting precise. LetΩ⊆ℝn, n≥2, denote a bounded domain with smooth boundary ∂Ω and outer normal vector ν. Let Σ⊆∂Ω be an open part of the boundary. L+∞(Ω) denotes the subspace of L∞(Ω)-functions with positive essential infima. H⋄1(Ω) and L⋄2(Σ) denote the spaces of H1- and L2-functions with vanishing integral mean on ∂Ω (resp., Σ).Forσ∈L+∞(Ω) and g∈L⋄2(Σ), there exists a unique solution uσg∈H⋄1(Ω) of the elliptic partial differential equation
(7)∇·σ∇uσg=0inΩ,σ∂νuσg|∂Ω={gonΣ,0else,
so that we can define the Neumann-to-Dirichlet operator (NtD)
(8)Λ(σ):L⋄2(Σ)⟶L⋄2(Σ),g⟼uσg|Σ,
where uσg∈H⋄1(Ω) solves  (7). Λ(σ) is a self-adjoint, compact linear operator.Letσ0∈L+∞(Ω) be piecewise analytic. For each point z∈Ω that has a neighborhood in which σ0 is analytic, and each unit vector d∈ℝn, ∥d∥=1, let Φz,d be the solution of
(9)∇·σ0∇Φz,d=d·∇δzinΩ,σ0∂νΦz,d|∂Ω=0.Φz,d is called a dipole function.


### 2.2. Auxiliary Results

Our presentation of the Factorization Method in the next section relies on the following four lemmas. The first lemma is frequently called amonotony lemma since it shows that a larger conductivity leads to a smaller NtD. More precisely, it shows a relation between the difference of two NtDs and the difference of the corresponding conductivities and the interior energy of an electric potential. The second lemma shows that this energy term is the image of the adjoint of an auxiliary virtual measurement operator that is defined on a subregion of Ω. The third lemma is a functional analytic relation between the norm of an image of an operator and the range of its adjoint. Together with the first two lemmas, it implies that the range of the auxiliary virtual measurement operator can be calculated from the NtDs. Finally, using the previous dipole functions, the last lemma shows that the range of the auxiliary virtual measurement operator determines the region on which they are defined.We start with the monotony lemma.Lemma 1.
Letσ1,σ0∈L+∞(Ω). Then, for all g∈L⋄2(∂Ω), (10)∫Ω‍(σ0-σ1)|∇u0|2dx≤∫Σ‍g(Λ1-Λ0)gds≤∫Ω‍σ0σ1(σ0-σ1)|∇u0|2dx,where we abbreviated Λj:=Λ(σj), j=0,1, and u0:=uσ0g.Proof.
The lemma seems to go back to Ikehata, Kang, Seo, and Sheen [22, 23], cf. also the similar arguments in Kirsch [8], Ide et al. [24], and in the works of Seo and the author [16, 25]. For the sake of completeness, we copy the short proof from [25]. For all g∈L⋄2(∂Ω), we have that
(11)∫Ω‍σ1∇u1·∇u0dx=∫Σ‍gu0ds=∫Ω‍σ0∇u0·∇u0dx=∫Σ‍gΛ0gds.
Hence, from
(12)∫Ω‍σ1∇(u1-u0)·∇(u1-u0)dx=∫Ω‍σ1|∇u1|2dx-∫Ω‍σ0|∇u0|2dx+∫Ω‍(σ1-σ0)|∇u0|2dx,
we obtain that
(13)∫Σ‍g(Λ1-Λ0)gds=∫Ω‍(σ0-σ1)|∇u0|2dx+∫Ω‍σ1|∇(u1-u0)|2dx,
which already yields the first asserted inequality.
By interchangingσ1 and σ0, we conclude that
(14)∫Σ‍g(Λ0-Λ1)gds=∫Ω‍(σ1-σ0)|∇u1|2dx+∫Ω‍σ0|∇(u0-u1)|2dx=∫Ω‍(σ1|∇u1-σ0σ1∇u0|2+(σ0-σ02σ1)|∇u0|2)dx
and, hence, obtain the second inequality.Given a reference conductivityσ0∈L+∞(Ω) and a measurable subset D⊆Ω, we define the virtual measurement operatorLD by
(15)LD:L2(D)n⟶L⋄2(Σ),F⟼v|Σ,
where v∈H⋄1(Ω) solves
(16)∫Ωσ0∇v·∇wdx=∫DF·∇wdx∀w∈H⋄1(Ω).
The energy term |∇u0|2 in Lemma 1 can be identified with the norm of the adjoint of this virtual measurement operator.Lemma 2.
The adjoint operator ofLD is given by
(17)LD*:L⋄2(Σ)⟶L2(D)n,g⟼∇u0|D,
where u0∈H⋄1(Ω) solves
(18)∇·σ0∇u0=0inΩ,σ0∂νu0|∂Ω={gonΣ,0else.Proof.
For allg∈L⋄2(Σ) and F∈L2(D)n, we have that
(19)∫D(LD*g)·Fdx=∫Σg(LDF)dx=∫Σgv|Σdx=∫Ωσ0∇u0·∇vdx=∫D∇u0·Fdx,
which shows the assertion.The following functional analytic lemma uses bounds on the image of an operator to characterize the range of its dual operator.Lemma 3.
LetX and Y be real Hilbert spaces with inner products (·,·)X and (·,·)Y, respectively. Let A∈ℒ(X;Y) and x′∈X. Then,(20)x′∈ℛ(A*)iff∃C>0:|(x′,x)X|≤C∥Ax∥∀x∈X.In particular, if X, Y1, and Y2 are three real Hilbert spaces, Ai∈ℒ(Yi,X), i=1,2, and if there exists C>0 with
(21)∥A1*x∥≤C∥A2*x∥∀x∈X,
then ℛ(A1)⊆ℛ(A2).Proof.
The assertion can be generalized to Banach spaces, and, in that context, it is called the “14th important property of Banach spaces” in Bourbaki [26]. For the sake of completeness, we rewrite the proof from [27] to Hilbert spaces.
Ifx′∈ℛ(A*), then there exists y′∈Y such that x′=A*y′. Hence,
(22)|(x′,x)X|=|(A*y′,x)X|=|(y′,Ax)Y|≤∥y′∥∥Ax∥∀x∈X,
so that the assertion holds with C=∥y′∥.
Now letx′∈X be such that there exists C>0 with |(x′,x)X|≤C∥Ax∥ for all x∈X. We define
(23)f(y):=(x′,x)Xforeveryy=Ax∈ℛ(A).
Then, f is a well-defined, continuous linear functional on ℛ(A). By setting it to zero on ℛ(A)⊥, we can extend f to a continuous linear functional on Y. Using the Riesz theorem, it follows that there exists y′∈Y with
(24)(y′,y)Y=f(y)∀y∈ℛ(A).
Hence, for all x∈X, we have
(25)(A*y′,x)X=(y′,Ax)Y=f(Ax)=(x′,x)X,
so that x′=A*y'∈ℛ(A*).The last lemma shows that the range of the virtual measurement operatorLD determines the region D on which it is defined. We state the lemma for a simple special case, a generalized version of the lemma will be formulated in Section 3.2.Lemma 4.
Letσ0=1, D⊆Ω be open, and D¯⊆Ω have a connected complement Ω∖D¯.
Then, for all unit vectorsd∈ℝn, ∥d∥=1, and every point z∈Ω∖∂D, it holds that
(26)z∈DiffΦz,d|Σ∈ℛ(LD).Proof.
The proof is similar to the one of [21, Lemma 2.9].
First, letz∈D and ϵ>0 be such that Bϵ(z)¯⊆D. We choose
(27)f1∈H1(Bϵ(z))withf1|∂Bϵ(z)=Φz,d|∂Bϵ(z),f2∈H1(Bϵ(z))withΔf2=0,∂νf2|∂Bϵ(z)=∂νΦz,d|∂Bϵ(z),
and let F∈L2(D)n be the zero continuation of ∇(f1-f2) to D.
Then, the function(28)v:={Φz,dinΩ∖Bϵ(z)¯,f1inBϵ(z)
fulfills v∈H⋄1(Ω), and, for all w∈H⋄1(Ω),
(29)∫Ω‍∇v·∇wdx=∫Ω∖Bϵ(z)¯‍∇Φz,d·∇wdx+∫Bϵ(z)‍∇f1·∇wdx=-∫∂Bϵ(z)‍∂νΦz,dw|∂Bϵ(z)ds+∫Bϵ(z)‍∇f1·∇wdx=∫Bϵ(z)‍∇(f1-f2)·∇wdx=∫D‍F·∇wdx.
This shows that Φz,d|Σ=v|Σ=LD(F)∈ℛ(LD).
Now letΦz,d|Σ∈ℛ(LD). Let v∈H⋄1(Ω) be the function from the definition of LD. Then,
(30)v|Σ=Φz,d|Σ,∂νv|Σ=0=∂νΦz,d|Σ,
so that it follows by unique continuation that v=Φz,d in the connected set Ω∖(D¯∪{z}).
Ifz∉D-, then d·∇δz∉H-2(Ω∖D¯), and thus Φz,d∉L2(Ω∖D¯), which contradicts that v=Φz,d in Ω∖(D¯∪{z}). Hence, z∈D¯.


## 2.1. The Setting

We start by making the mathematical setting precise. LetΩ⊆ℝn, n≥2, denote a bounded domain with smooth boundary ∂Ω and outer normal vector ν. Let Σ⊆∂Ω be an open part of the boundary. L+∞(Ω) denotes the subspace of L∞(Ω)-functions with positive essential infima. H⋄1(Ω) and L⋄2(Σ) denote the spaces of H1- and L2-functions with vanishing integral mean on ∂Ω (resp., Σ).Forσ∈L+∞(Ω) and g∈L⋄2(Σ), there exists a unique solution uσg∈H⋄1(Ω) of the elliptic partial differential equation
(7)∇·σ∇uσg=0inΩ,σ∂νuσg|∂Ω={gonΣ,0else,
so that we can define the Neumann-to-Dirichlet operator (NtD)
(8)Λ(σ):L⋄2(Σ)⟶L⋄2(Σ),g⟼uσg|Σ,
where uσg∈H⋄1(Ω) solves  (7). Λ(σ) is a self-adjoint, compact linear operator.Letσ0∈L+∞(Ω) be piecewise analytic. For each point z∈Ω that has a neighborhood in which σ0 is analytic, and each unit vector d∈ℝn, ∥d∥=1, let Φz,d be the solution of
(9)∇·σ0∇Φz,d=d·∇δzinΩ,σ0∂νΦz,d|∂Ω=0.Φz,d is called a dipole function.


## 2.2. Auxiliary Results

Our presentation of the Factorization Method in the next section relies on the following four lemmas. The first lemma is frequently called amonotony lemma since it shows that a larger conductivity leads to a smaller NtD. More precisely, it shows a relation between the difference of two NtDs and the difference of the corresponding conductivities and the interior energy of an electric potential. The second lemma shows that this energy term is the image of the adjoint of an auxiliary virtual measurement operator that is defined on a subregion of Ω. The third lemma is a functional analytic relation between the norm of an image of an operator and the range of its adjoint. Together with the first two lemmas, it implies that the range of the auxiliary virtual measurement operator can be calculated from the NtDs. Finally, using the previous dipole functions, the last lemma shows that the range of the auxiliary virtual measurement operator determines the region on which they are defined.We start with the monotony lemma.Lemma 1.
Letσ1,σ0∈L+∞(Ω). Then, for all g∈L⋄2(∂Ω), (10)∫Ω‍(σ0-σ1)|∇u0|2dx≤∫Σ‍g(Λ1-Λ0)gds≤∫Ω‍σ0σ1(σ0-σ1)|∇u0|2dx,where we abbreviated Λj:=Λ(σj), j=0,1, and u0:=uσ0g.Proof.
The lemma seems to go back to Ikehata, Kang, Seo, and Sheen [22, 23], cf. also the similar arguments in Kirsch [8], Ide et al. [24], and in the works of Seo and the author [16, 25]. For the sake of completeness, we copy the short proof from [25]. For all g∈L⋄2(∂Ω), we have that
(11)∫Ω‍σ1∇u1·∇u0dx=∫Σ‍gu0ds=∫Ω‍σ0∇u0·∇u0dx=∫Σ‍gΛ0gds.
Hence, from
(12)∫Ω‍σ1∇(u1-u0)·∇(u1-u0)dx=∫Ω‍σ1|∇u1|2dx-∫Ω‍σ0|∇u0|2dx+∫Ω‍(σ1-σ0)|∇u0|2dx,
we obtain that
(13)∫Σ‍g(Λ1-Λ0)gds=∫Ω‍(σ0-σ1)|∇u0|2dx+∫Ω‍σ1|∇(u1-u0)|2dx,
which already yields the first asserted inequality.
By interchangingσ1 and σ0, we conclude that
(14)∫Σ‍g(Λ0-Λ1)gds=∫Ω‍(σ1-σ0)|∇u1|2dx+∫Ω‍σ0|∇(u0-u1)|2dx=∫Ω‍(σ1|∇u1-σ0σ1∇u0|2+(σ0-σ02σ1)|∇u0|2)dx
and, hence, obtain the second inequality.Given a reference conductivityσ0∈L+∞(Ω) and a measurable subset D⊆Ω, we define the virtual measurement operatorLD by
(15)LD:L2(D)n⟶L⋄2(Σ),F⟼v|Σ,
where v∈H⋄1(Ω) solves
(16)∫Ωσ0∇v·∇wdx=∫DF·∇wdx∀w∈H⋄1(Ω).
The energy term |∇u0|2 in Lemma 1 can be identified with the norm of the adjoint of this virtual measurement operator.Lemma 2.
The adjoint operator ofLD is given by
(17)LD*:L⋄2(Σ)⟶L2(D)n,g⟼∇u0|D,
where u0∈H⋄1(Ω) solves
(18)∇·σ0∇u0=0inΩ,σ0∂νu0|∂Ω={gonΣ,0else.Proof.
For allg∈L⋄2(Σ) and F∈L2(D)n, we have that
(19)∫D(LD*g)·Fdx=∫Σg(LDF)dx=∫Σgv|Σdx=∫Ωσ0∇u0·∇vdx=∫D∇u0·Fdx,
which shows the assertion.The following functional analytic lemma uses bounds on the image of an operator to characterize the range of its dual operator.Lemma 3.
LetX and Y be real Hilbert spaces with inner products (·,·)X and (·,·)Y, respectively. Let A∈ℒ(X;Y) and x′∈X. Then,(20)x′∈ℛ(A*)iff∃C>0:|(x′,x)X|≤C∥Ax∥∀x∈X.In particular, if X, Y1, and Y2 are three real Hilbert spaces, Ai∈ℒ(Yi,X), i=1,2, and if there exists C>0 with
(21)∥A1*x∥≤C∥A2*x∥∀x∈X,
then ℛ(A1)⊆ℛ(A2).Proof.
The assertion can be generalized to Banach spaces, and, in that context, it is called the “14th important property of Banach spaces” in Bourbaki [26]. For the sake of completeness, we rewrite the proof from [27] to Hilbert spaces.
Ifx′∈ℛ(A*), then there exists y′∈Y such that x′=A*y′. Hence,
(22)|(x′,x)X|=|(A*y′,x)X|=|(y′,Ax)Y|≤∥y′∥∥Ax∥∀x∈X,
so that the assertion holds with C=∥y′∥.
Now letx′∈X be such that there exists C>0 with |(x′,x)X|≤C∥Ax∥ for all x∈X. We define
(23)f(y):=(x′,x)Xforeveryy=Ax∈ℛ(A).
Then, f is a well-defined, continuous linear functional on ℛ(A). By setting it to zero on ℛ(A)⊥, we can extend f to a continuous linear functional on Y. Using the Riesz theorem, it follows that there exists y′∈Y with
(24)(y′,y)Y=f(y)∀y∈ℛ(A).
Hence, for all x∈X, we have
(25)(A*y′,x)X=(y′,Ax)Y=f(Ax)=(x′,x)X,
so that x′=A*y'∈ℛ(A*).The last lemma shows that the range of the virtual measurement operatorLD determines the region D on which it is defined. We state the lemma for a simple special case, a generalized version of the lemma will be formulated in Section 3.2.Lemma 4.
Letσ0=1, D⊆Ω be open, and D¯⊆Ω have a connected complement Ω∖D¯.
Then, for all unit vectorsd∈ℝn, ∥d∥=1, and every point z∈Ω∖∂D, it holds that
(26)z∈DiffΦz,d|Σ∈ℛ(LD).Proof.
The proof is similar to the one of [21, Lemma 2.9].
First, letz∈D and ϵ>0 be such that Bϵ(z)¯⊆D. We choose
(27)f1∈H1(Bϵ(z))withf1|∂Bϵ(z)=Φz,d|∂Bϵ(z),f2∈H1(Bϵ(z))withΔf2=0,∂νf2|∂Bϵ(z)=∂νΦz,d|∂Bϵ(z),
and let F∈L2(D)n be the zero continuation of ∇(f1-f2) to D.
Then, the function(28)v:={Φz,dinΩ∖Bϵ(z)¯,f1inBϵ(z)
fulfills v∈H⋄1(Ω), and, for all w∈H⋄1(Ω),
(29)∫Ω‍∇v·∇wdx=∫Ω∖Bϵ(z)¯‍∇Φz,d·∇wdx+∫Bϵ(z)‍∇f1·∇wdx=-∫∂Bϵ(z)‍∂νΦz,dw|∂Bϵ(z)ds+∫Bϵ(z)‍∇f1·∇wdx=∫Bϵ(z)‍∇(f1-f2)·∇wdx=∫D‍F·∇wdx.
This shows that Φz,d|Σ=v|Σ=LD(F)∈ℛ(LD).
Now letΦz,d|Σ∈ℛ(LD). Let v∈H⋄1(Ω) be the function from the definition of LD. Then,
(30)v|Σ=Φz,d|Σ,∂νv|Σ=0=∂νΦz,d|Σ,
so that it follows by unique continuation that v=Φz,d in the connected set Ω∖(D¯∪{z}).
Ifz∉D-, then d·∇δz∉H-2(Ω∖D¯), and thus Φz,d∉L2(Ω∖D¯), which contradicts that v=Φz,d in Ω∖(D¯∪{z}). Hence, z∈D¯.


## 3. The Factorization Method

Now we will formulate the Factorization Method and characterize a region where a conductivityσ differs from a reference conductivity σ0 by a range criterion. Before we turn to a new general formulation of the method, we first state it for a special case that is similar to the one that was treated in the original works of Brühl and Hanke [3, 4].


### 3.1. The Factorization Method for a Simple Special Case

Theorem 5.
Letσ0=1 and σ=1+χD, where D⊆Ω is an open set so that D¯⊆Ω has a connected complement Ω∖D¯. Then, for all z∈Ω, z∉∂D, and all dipole directions d∈ℝn, ∥d∥=1,
(31)z∈DiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)|1/2).Proof.
The monotony Lemma1 yields that for all g∈L⋄2(∂Ω),
(32)∫D‍|∇u0|2dx≥∫Σ‍g(Λ(1)-Λ(σ))gds≥∫D‍12|∇u0|2dx.
Hence, |Λ(σ)-Λ(1)|=Λ(1)-Λ(σ), and, using Lemma 2, we can restate this in the form
(33)∥LD*g∥2≥∥|Λ(σ)-Λ(1)|1/2g∥2≥12∥LD*g∥2.
Using the functional analytic Lemma 3, this implies that
(34)ℛ(LD)=ℛ(|Λ(σ)-Λ(1)|1/2),
and thus the assertion follows from the relation between D and ℛ(LD) in Lemma 4.Obviously, the same arguments can be used to treat the caseσ(x)=1+κ(x)χD(x), when there exists a conductivity jump ϵ>0 so that either
(35)κ(x)≥ϵ∀x∈Dorκ(x)≤-ϵ∀x∈D.


### 3.2. The Factorization Method for the General Piecewise Analytic Case

Now we drop the assumptions that the background is constant, that there is a clear conductivity jump, and that the complement of the inclusions is connected. We will merely assume that the reference conductivityσ0 is a piecewise analytic function and that either σ-σ0≥0 or σ-σ0≤0. Roughly speaking, under this general assumption, the Factorization Method then characterizes the support of σ-σ0 up to holes in the support that have no connections to Σ. For a precise formulation, we use the concept of the inner and outer support from [28] that has been inspired by the use of the infinity support of Kusiak and Sylvester [29]; see also [25, 30].Definition 6.
A relatively open setU⊆Ω¯ is called connected to Σ if U∩Ω is connected and U∩Σ≠∅.
For a measurable functionκ:Ω→ℝ, we define (a)
thesupportsupp(κ) as the complement (in Ω¯) of the union of those relatively open U⊆Ω¯, for which κ|U≡0,(b)
theinner supportinnsuppκ as the union of those open sets U⊆Ω, for which essinfx∈U|κ(x)|>0,(c)
theouter supportoutΣsuppκ as the complement (in Ω¯) of the union of those relatively open U⊆Ω¯ that are connected to Σ and for which κ|U≡0.The interior of a set M⊆Ω is denoted by intM and its closure (with respect to ℝn) by M¯. If M is measurable, we also define(d)
outΣM=outΣsuppχM. It is easily checked that outΣ(suppκ)=outΣsuppκ.
With this concept, we can extend the range characterization in Lemma4 to a general setting (see also Remark 9 later).Lemma 7.
Letσ0∈L+∞(Ω) be piecewise analytic. Let D⊆Ω be measurable.
Then, for all unit vectorsd∈ℝn, ∥d∥=1, and every point z∈Ω that has a neighborhood in which σ0 is analytic,
(36)z∈intDimpliesΦz,d|Σ∈ℛ(LD),
and
(37)Φz,d|Σ∈ℛ(LD)impliesz∈outΣD.Proof.
Ifz∈intD, then there exists a small ball Bϵ(z)¯⊆D, and the first assertion follows as in the proof of Lemma 4.
To show the second assertion, letΦz,d|Σ∈ℛ(LD), and let v∈H⋄1(Ω) be the function from the definition of LD, so that (as in the proof of Lemma 4)
(38)v|Σ=Φz,d|Σ,∂νv|Σ=0=∂νΦz,d|Σ.
Assume that z∉outΣD. Then, there exists a relatively open U⊆Ω¯ that is connected to Σ and contains z. Hence, by unique continuation, it follows that v|U=Φz,d|U, and we obtain the same contradiction as in the proof of Lemma 4.Now, we can formulate and prove the Factorization Method for general piecewise analytic conductivities.Theorem 8.
Letσ∈L+∞(Ω), and let σ0∈L+∞(Ω) be a piecewise analytic function. Let either
(39)σ(x)≥σ0(x)∀x∈Ωorσ(x)≤σ0(x)∀x∈Ω.
Then, for all z∈Ω that have a neighborhood in which σ0 is analytic, as well as all unit vectors d∈ℝn, ∥d∥=1,
(40)z∈innsupp(σ-σ0)impliesΦz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2),
and
(41)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)impliesz∈outΣsupp(σ-σ0).Proof.
Letz∈Ω have a neighborhood in which σ0 is analytic, and let d∈ℝn be a unit vector with ∥d∥=1. We only prove the assertions for σ≥σ0. The other case is completely analogous.
First, letz∈innsupp(σ-σ0). Then there exists a small ball Bϵ(z) and δ>0 so that σ-σ0≥δ on Bϵ(z). Using the monotony Lemma 1, it follows that, for all g∈L⋄2(∂Ω),
(42)∫Σ‍g(Λ(σ0)-Λ(σ))gds≥∫Ω‍σ0σ(σ-σ0)|∇u0|2dx≥δ∥σ0∥L∞(Ω)∥1σ∥L∞(Ω)∫Bϵ(z)‍|∇u0|2dx≥δ∥σ0∥L∞(Ω)∥1σ∥L∞(Ω)∥LBϵ(z)*g∥2.
Using the functional analytic Lemma 3, we obtain that
(43)ℛ(LBϵ(z))⊆ℛ(|Λ(σ)-Λ(σ0)|1/2),
and Lemma 7 yields that
(44)Φz,d|Σ∈ℛ(LBϵ(z))⊆ℛ(|Λ(σ)-Λ(σ0)|1/2).
On the other hand, with D:=supp(σ-σ0), the monotony Lemma 1 shows that for all g∈L⋄2(∂Ω),
(45)∫Σ‍g(Λ(σ0)-Λ(σ1))gds≤∫Ω‍(σ1-σ0)|∇u0|2dx≤∥σ1-σ0∥L∞(Ω)∫D‍|∇u0|2dx≤∥σ1-σ0∥L∞(Ω)∥LD*g∥2,
so that we obtain from the functional analytic Lemma 3(46)ℛ(|Λ(σ)-Λ(σ0)|1/2)⊆ℛ(LD).
Hence, Lemma 7 yields that
(47)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)⊆ℛ(LD)impliesz∈outΣ(supp(σ-σ0))=outΣsupp(σ-σ0).Remark 9.
Theorem8 shows that the Factorization Method is able to detect the support of a conductivity difference up to the difference between the outer and the inner support, that is, roughly speaking, up to holes in the support that have no connections to the boundary. It leaves open whether points in such holes will fulfill the range criterion of the Factorization Method or not.
A result of Hyvönen and the author [30, Lemma 2.5] shows that for every smooth domain D with D¯⊂Ω and every unit vector d∈ℝn, ∥d∥=1,
(48)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)∀z∈∂D
implies
(49)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)∀z∈D.
In that sense, we can expect that holes in the support will be filled up and that the set detected by the Factorization Method is essentially the outer support of the conductivity difference.


### 3.3. The Factorization Method for the Indefinite Case

It is a long standing open theoretical problem whether the range criterion of the Factorization Method holds true without thedefiniteness assumption that σ≥σ0 on Ω or σ≤σ0 on Ω. However, Grinberg, Kirsch, and Schmitt [18, 31] showed how to exclude a region E⊆Ω from Ω, in such a way that the Factorization Method only requires the definiteness assumption on Ω∖E. In this subsection, we show how their idea can be incorporated into our formulation of the method.To point out the main idea, we first formulate the result for a simple special case. Let us stress that, forσ=1+χD+-(1/2)χD-, it is not known whether
(50)z∈D+∪D-iffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)|1/2).
However, we can still use the Factorization Method if we have some a priori knowledge that separates D+ and D-. More precisely, if we know a subset E that contains D- without intersecting D+, then we can use the Factorization Method to find D+∪E (and thus D+).Theorem 10.
Letσ0=1 and σ=1+χD+-(1/2)χD-, where D+,D-⊆Ω are open. Let E⊆Ω be an open set.(a)
IfD+⊆E and D-∪E¯⊆Ω has a connected complement, then for all z∈Ω, z∉∂(D-∪E), and all dipole directions d∈ℝn, ∥d∥=1,
(51)z∈D-∪EiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)+LELE*|1/2).(b)
IfD-⊆E and D+∪E¯⊆Ω has a connected complement, then for all z∈Ω, z∉∂(D+∪E), and all dipole directions d∈ℝn, ∥d∥=1,
(52)z∈D+∪EiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)-LELE*|1/2).Proof.
The monotony Lemma1 yields that for all g∈L⋄2(∂Ω),
(53)12∫D-‍|∇u0|2dx-∫D+‍|∇u0|2dx≤∫Σ‍g(Λ(σ)-Λ(1))gds≤∫D-‍|∇u0|2dx-12∫D+‍|∇u0|2dx.
Since (cf. Lemma 2)
(54)∫Σ‍g(LELE*)gds=∥LE*g∥2=∫E‍|∇u0|2dx,
it follows for case (a) that
(55)12∫D-∪E‍|∇u0|2dx≤∫Σ‍g(Λ(σ)-Λ(1)+2LELE*)gds≤∫D-∪E‍|∇u0|2dx.
Using the functional analytic Lemma 3, this implies that
(56)ℛ(LD-∪E)=ℛ(|Λ(σ)-Λ(1)+2LELE*|1/2),
so that the assertion (a) follows from Lemma 4.
In case (b), we obtain that(57)32∫D+∪E‍|∇u0|2dx≥∫Σ‍g(Λ(1)-Λ(σ)+2LELE*)gds≥12∫D+∪E‍|∇u0|2dx,
and the same arguments as above yield the assertion.We can also extend these ideas to the general setting of Section3.2.Theorem 11.
Letσ∈L+∞(Ω) and let σ0∈L+∞(Ω) be a piecewise analytic function. Let E⊆Ω be a measurable set.
Chooseα,β∈ℝ such that
(58)α>∥σ-σ0∥L∞(Ω),β>∥σ0σ(σ0-σ)∥L∞(Ω).(a)
Ifσ≤σ0 on Ω∖E, then for all z∈Ω that have a neighborhood in which σ0 is analytic, as well as all unit vectors d∈ℝn, ∥d∥=1,
(59)z∈innsupp(σ-σ0)∪EimpliesΦz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)+αLELE*|1/2),and(60)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)+αLELE*|1/2)impliesz∈outΣ(supp(σ-σ0)∪E).(b)
Ifσ≥σ0 on Ω∖E, then for all z∈Ω that have a neighborhood in which σ0 is analytic, as well as all unit vectors d∈ℝn, ∥d∥=1,
(61)z∈innsupp(σ-σ0)∪E,impliesΦz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)-βLELE*|1/2)and(62)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)-βLELE*|1/2)impliesz∈outΣ(supp(σ-σ0)∪E).Proof.
For everyz∈innsupp(σ-σ0) with z∉E, there exists a small ball Bϵ(z) and δ>0 so that σ0-σ≥δ on Bϵ(z). Using the monotony Lemma 1, it follows that, for all g∈L⋄2(∂Ω),
(63)∫Σ‍g(Λ(σ)-Λ(σ0)+αLELE*)gds≥∫Ω‍(σ0-σ1)|∇u0|2dx+α∫E‍|∇u0|2dx≥δ∫Bϵ(z)‍|∇u0|2dx+(α-∥σ-σ0∥L∞(Ω))∫E‍|∇u0|2dx.
As in the previous proofs, we obtain from Lemmas 2 and 3 that
(64)ℛ(|Λ(σ)-Λ(σ0)+αLELE*|1/2)⊇ℛ(LBϵ(z)∪E),
so that the first implication of (a) follows from Lemma 7.
The monotony Lemma1 also implies that
(65)∫Σ‍g(Λ(σ)-Λ(σ0)+αLELE*)gds≤∫Ω‍σ0σ(σ0-σ)|∇u0|2dx+α∫E‍|∇u0|2dx≤∥σ0σ(σ0-σ)∥L∞(Ω)∫supp(σ-σ0)‍|∇u0|2dx+α∫E‍|∇u0|2dx,
so that the second implication of (a) follows from Lemmas 2, 3, and 7. Assertion (b) can be proven analogously.Remark9 also applies to this case.


## 3.1. The Factorization Method for a Simple Special Case

Theorem 5.
Letσ0=1 and σ=1+χD, where D⊆Ω is an open set so that D¯⊆Ω has a connected complement Ω∖D¯. Then, for all z∈Ω, z∉∂D, and all dipole directions d∈ℝn, ∥d∥=1,
(31)z∈DiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)|1/2).Proof.
The monotony Lemma1 yields that for all g∈L⋄2(∂Ω),
(32)∫D‍|∇u0|2dx≥∫Σ‍g(Λ(1)-Λ(σ))gds≥∫D‍12|∇u0|2dx.
Hence, |Λ(σ)-Λ(1)|=Λ(1)-Λ(σ), and, using Lemma 2, we can restate this in the form
(33)∥LD*g∥2≥∥|Λ(σ)-Λ(1)|1/2g∥2≥12∥LD*g∥2.
Using the functional analytic Lemma 3, this implies that
(34)ℛ(LD)=ℛ(|Λ(σ)-Λ(1)|1/2),
and thus the assertion follows from the relation between D and ℛ(LD) in Lemma 4.Obviously, the same arguments can be used to treat the caseσ(x)=1+κ(x)χD(x), when there exists a conductivity jump ϵ>0 so that either
(35)κ(x)≥ϵ∀x∈Dorκ(x)≤-ϵ∀x∈D.


## 3.2. The Factorization Method for the General Piecewise Analytic Case

Now we drop the assumptions that the background is constant, that there is a clear conductivity jump, and that the complement of the inclusions is connected. We will merely assume that the reference conductivityσ0 is a piecewise analytic function and that either σ-σ0≥0 or σ-σ0≤0. Roughly speaking, under this general assumption, the Factorization Method then characterizes the support of σ-σ0 up to holes in the support that have no connections to Σ. For a precise formulation, we use the concept of the inner and outer support from [28] that has been inspired by the use of the infinity support of Kusiak and Sylvester [29]; see also [25, 30].Definition 6.
A relatively open setU⊆Ω¯ is called connected to Σ if U∩Ω is connected and U∩Σ≠∅.
For a measurable functionκ:Ω→ℝ, we define (a)
thesupportsupp(κ) as the complement (in Ω¯) of the union of those relatively open U⊆Ω¯, for which κ|U≡0,(b)
theinner supportinnsuppκ as the union of those open sets U⊆Ω, for which essinfx∈U|κ(x)|>0,(c)
theouter supportoutΣsuppκ as the complement (in Ω¯) of the union of those relatively open U⊆Ω¯ that are connected to Σ and for which κ|U≡0.The interior of a set M⊆Ω is denoted by intM and its closure (with respect to ℝn) by M¯. If M is measurable, we also define(d)
outΣM=outΣsuppχM. It is easily checked that outΣ(suppκ)=outΣsuppκ.
With this concept, we can extend the range characterization in Lemma4 to a general setting (see also Remark 9 later).Lemma 7.
Letσ0∈L+∞(Ω) be piecewise analytic. Let D⊆Ω be measurable.
Then, for all unit vectorsd∈ℝn, ∥d∥=1, and every point z∈Ω that has a neighborhood in which σ0 is analytic,
(36)z∈intDimpliesΦz,d|Σ∈ℛ(LD),
and
(37)Φz,d|Σ∈ℛ(LD)impliesz∈outΣD.Proof.
Ifz∈intD, then there exists a small ball Bϵ(z)¯⊆D, and the first assertion follows as in the proof of Lemma 4.
To show the second assertion, letΦz,d|Σ∈ℛ(LD), and let v∈H⋄1(Ω) be the function from the definition of LD, so that (as in the proof of Lemma 4)
(38)v|Σ=Φz,d|Σ,∂νv|Σ=0=∂νΦz,d|Σ.
Assume that z∉outΣD. Then, there exists a relatively open U⊆Ω¯ that is connected to Σ and contains z. Hence, by unique continuation, it follows that v|U=Φz,d|U, and we obtain the same contradiction as in the proof of Lemma 4.Now, we can formulate and prove the Factorization Method for general piecewise analytic conductivities.Theorem 8.
Letσ∈L+∞(Ω), and let σ0∈L+∞(Ω) be a piecewise analytic function. Let either
(39)σ(x)≥σ0(x)∀x∈Ωorσ(x)≤σ0(x)∀x∈Ω.
Then, for all z∈Ω that have a neighborhood in which σ0 is analytic, as well as all unit vectors d∈ℝn, ∥d∥=1,
(40)z∈innsupp(σ-σ0)impliesΦz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2),
and
(41)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)impliesz∈outΣsupp(σ-σ0).Proof.
Letz∈Ω have a neighborhood in which σ0 is analytic, and let d∈ℝn be a unit vector with ∥d∥=1. We only prove the assertions for σ≥σ0. The other case is completely analogous.
First, letz∈innsupp(σ-σ0). Then there exists a small ball Bϵ(z) and δ>0 so that σ-σ0≥δ on Bϵ(z). Using the monotony Lemma 1, it follows that, for all g∈L⋄2(∂Ω),
(42)∫Σ‍g(Λ(σ0)-Λ(σ))gds≥∫Ω‍σ0σ(σ-σ0)|∇u0|2dx≥δ∥σ0∥L∞(Ω)∥1σ∥L∞(Ω)∫Bϵ(z)‍|∇u0|2dx≥δ∥σ0∥L∞(Ω)∥1σ∥L∞(Ω)∥LBϵ(z)*g∥2.
Using the functional analytic Lemma 3, we obtain that
(43)ℛ(LBϵ(z))⊆ℛ(|Λ(σ)-Λ(σ0)|1/2),
and Lemma 7 yields that
(44)Φz,d|Σ∈ℛ(LBϵ(z))⊆ℛ(|Λ(σ)-Λ(σ0)|1/2).
On the other hand, with D:=supp(σ-σ0), the monotony Lemma 1 shows that for all g∈L⋄2(∂Ω),
(45)∫Σ‍g(Λ(σ0)-Λ(σ1))gds≤∫Ω‍(σ1-σ0)|∇u0|2dx≤∥σ1-σ0∥L∞(Ω)∫D‍|∇u0|2dx≤∥σ1-σ0∥L∞(Ω)∥LD*g∥2,
so that we obtain from the functional analytic Lemma 3(46)ℛ(|Λ(σ)-Λ(σ0)|1/2)⊆ℛ(LD).
Hence, Lemma 7 yields that
(47)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)⊆ℛ(LD)impliesz∈outΣ(supp(σ-σ0))=outΣsupp(σ-σ0).Remark 9.
Theorem8 shows that the Factorization Method is able to detect the support of a conductivity difference up to the difference between the outer and the inner support, that is, roughly speaking, up to holes in the support that have no connections to the boundary. It leaves open whether points in such holes will fulfill the range criterion of the Factorization Method or not.
A result of Hyvönen and the author [30, Lemma 2.5] shows that for every smooth domain D with D¯⊂Ω and every unit vector d∈ℝn, ∥d∥=1,
(48)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)∀z∈∂D
implies
(49)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)|1/2)∀z∈D.
In that sense, we can expect that holes in the support will be filled up and that the set detected by the Factorization Method is essentially the outer support of the conductivity difference.


## 3.3. The Factorization Method for the Indefinite Case

It is a long standing open theoretical problem whether the range criterion of the Factorization Method holds true without thedefiniteness assumption that σ≥σ0 on Ω or σ≤σ0 on Ω. However, Grinberg, Kirsch, and Schmitt [18, 31] showed how to exclude a region E⊆Ω from Ω, in such a way that the Factorization Method only requires the definiteness assumption on Ω∖E. In this subsection, we show how their idea can be incorporated into our formulation of the method.To point out the main idea, we first formulate the result for a simple special case. Let us stress that, forσ=1+χD+-(1/2)χD-, it is not known whether
(50)z∈D+∪D-iffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)|1/2).
However, we can still use the Factorization Method if we have some a priori knowledge that separates D+ and D-. More precisely, if we know a subset E that contains D- without intersecting D+, then we can use the Factorization Method to find D+∪E (and thus D+).Theorem 10.
Letσ0=1 and σ=1+χD+-(1/2)χD-, where D+,D-⊆Ω are open. Let E⊆Ω be an open set.(a)
IfD+⊆E and D-∪E¯⊆Ω has a connected complement, then for all z∈Ω, z∉∂(D-∪E), and all dipole directions d∈ℝn, ∥d∥=1,
(51)z∈D-∪EiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)+LELE*|1/2).(b)
IfD-⊆E and D+∪E¯⊆Ω has a connected complement, then for all z∈Ω, z∉∂(D+∪E), and all dipole directions d∈ℝn, ∥d∥=1,
(52)z∈D+∪EiffΦz,d|Σ∈ℛ(|Λ(σ)-Λ(1)-LELE*|1/2).Proof.
The monotony Lemma1 yields that for all g∈L⋄2(∂Ω),
(53)12∫D-‍|∇u0|2dx-∫D+‍|∇u0|2dx≤∫Σ‍g(Λ(σ)-Λ(1))gds≤∫D-‍|∇u0|2dx-12∫D+‍|∇u0|2dx.
Since (cf. Lemma 2)
(54)∫Σ‍g(LELE*)gds=∥LE*g∥2=∫E‍|∇u0|2dx,
it follows for case (a) that
(55)12∫D-∪E‍|∇u0|2dx≤∫Σ‍g(Λ(σ)-Λ(1)+2LELE*)gds≤∫D-∪E‍|∇u0|2dx.
Using the functional analytic Lemma 3, this implies that
(56)ℛ(LD-∪E)=ℛ(|Λ(σ)-Λ(1)+2LELE*|1/2),
so that the assertion (a) follows from Lemma 4.
In case (b), we obtain that(57)32∫D+∪E‍|∇u0|2dx≥∫Σ‍g(Λ(1)-Λ(σ)+2LELE*)gds≥12∫D+∪E‍|∇u0|2dx,
and the same arguments as above yield the assertion.We can also extend these ideas to the general setting of Section3.2.Theorem 11.
Letσ∈L+∞(Ω) and let σ0∈L+∞(Ω) be a piecewise analytic function. Let E⊆Ω be a measurable set.
Chooseα,β∈ℝ such that
(58)α>∥σ-σ0∥L∞(Ω),β>∥σ0σ(σ0-σ)∥L∞(Ω).(a)
Ifσ≤σ0 on Ω∖E, then for all z∈Ω that have a neighborhood in which σ0 is analytic, as well as all unit vectors d∈ℝn, ∥d∥=1,
(59)z∈innsupp(σ-σ0)∪EimpliesΦz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)+αLELE*|1/2),and(60)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)+αLELE*|1/2)impliesz∈outΣ(supp(σ-σ0)∪E).(b)
Ifσ≥σ0 on Ω∖E, then for all z∈Ω that have a neighborhood in which σ0 is analytic, as well as all unit vectors d∈ℝn, ∥d∥=1,
(61)z∈innsupp(σ-σ0)∪E,impliesΦz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)-βLELE*|1/2)and(62)Φz,d|Σ∈ℛ(|Λ(σ)-Λ(σ0)-βLELE*|1/2)impliesz∈outΣ(supp(σ-σ0)∪E).Proof.
For everyz∈innsupp(σ-σ0) with z∉E, there exists a small ball Bϵ(z) and δ>0 so that σ0-σ≥δ on Bϵ(z). Using the monotony Lemma 1, it follows that, for all g∈L⋄2(∂Ω),
(63)∫Σ‍g(Λ(σ)-Λ(σ0)+αLELE*)gds≥∫Ω‍(σ0-σ1)|∇u0|2dx+α∫E‍|∇u0|2dx≥δ∫Bϵ(z)‍|∇u0|2dx+(α-∥σ-σ0∥L∞(Ω))∫E‍|∇u0|2dx.
As in the previous proofs, we obtain from Lemmas 2 and 3 that
(64)ℛ(|Λ(σ)-Λ(σ0)+αLELE*|1/2)⊇ℛ(LBϵ(z)∪E),
so that the first implication of (a) follows from Lemma 7.
The monotony Lemma1 also implies that
(65)∫Σ‍g(Λ(σ)-Λ(σ0)+αLELE*)gds≤∫Ω‍σ0σ(σ0-σ)|∇u0|2dx+α∫E‍|∇u0|2dx≤∥σ0σ(σ0-σ)∥L∞(Ω)∫supp(σ-σ0)‍|∇u0|2dx+α∫E‍|∇u0|2dx,
so that the second implication of (a) follows from Lemmas 2, 3, and 7. Assertion (b) can be proven analogously.Remark9 also applies to this case.


## 4. Conclusions and Remarks

The Factorization Method can be used to detect regions in which a conductivity differs from a known reference conductivity. In this work, we summarized the progress on the method's theoretical foundation. We formulated the method for general piecewise analytic conductivities and gave comparatively simple and self-contained proofs. We also showed how the idea of excluding a part of the imaging region can be incorporated into this formulation.The regularity assumptions can be weakened even further. Our proofs only require unique continuation arguments for the reference conductivityσ0 and the existence of the dipole functions.Two major open theoretical questions still exist in the context of the Factorization Method. The theoretical justification of the method requires a definiteness condition (on the whole domain or after excluding an a priori known part of the domain). It is unknown whether the method's range criterion holds without such a definiteness condition. The second open question concerns the numerical stability of the method's range criterion. So far, there are no rigorous convergence results for numerical implementations of this range criterion (see, however, Lechleiter [32] for a first step in this direction). As a promising approach to overcome both problems, we would like to point out the recent work on monotony-based methods [28].



---

*Source: 425184-2013-08-27.xml*