NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS
|
|
- Audrey Wilkins
- 5 years ago
- Views:
Transcription
1 NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Contents. Introduction. The nonlinear Fourier transform 6 3. Inverse transforms and reconstruction algorithms The low-pass transport matrix method The shortcut method Periodization of the nonlinear inverse transform 4. Computational results Nonlinear Gibbs phenomenon in radial cases Comparison of the two methods in nonsymmetric cases Numerical evidence of the transport matrix efficiency 5 5. Conclusion 7 References 8 Abstract. Two reconstruction methods of Electrical Impedance Tomography (EIT) are numerically compared for nonsmooth conductivities in the plane based on the use of complex geometrical optics (CGO) solutions to D-bar equations involving the global uniqueness proofs for Calderón problem exposed in [Nachman; Annals of Mathematics 43, 996] and [Astala and Päivärinta; Annals of Mathematics 63, 006]: the Astala-Päivärinta theory-based low-pass transport matrix method implemented in [Astala et al.; Inverse Problems and Imaging 5, 0] and the shortcut method which considers ingredients of both theories. The latter method is formally similar to the Nachman theory-based regularized EIT reconstruction algorithm studied in [Knudsen, Lassas, Mueller and Siltanen; Inverse Problems and Imaging 3, 009] and several references from there. New numerical results are presented using parallel computation with size parameters larger than ever, leading mainly to two conclusions as follows. First, both methods can approximate piecewise constant conductivities better and better as the cutoff frequency increases, and there seems to be a Gibbs-like phenomenon producing ringing artifacts. Second, the transport matrix method loses accuracy away from a (freely chosen) pivot point located outside of the object to be studied, whereas the shortcut method produces reconstructions with more uniform quality.
2 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Keywords: Inverse problem, Beltrami equation, Conductivity equation, Inverse conductivity problem, Complex geometrical optics solution, Nonlinear Fourier transform, Scattering transform, Electrical impedance tomography.. Introduction We study a widely applicable nonlinear Fourier transform in dimension two. We perform numerical tests related to the nonlinear Gibbs phenomenon with much larger cutoff frequencies than before. Furthermore, we compare two computational inverse transformations, called low-pass transport matrix method and shortcut method in terms of accuracy. The inverse conductivity problem of Calderón [8] is the main source of applications of the nonlinear Fourier transform we consider. Let Ω R be the unit disc and let σ : Ω (0, ) be an essentially bounded measurable function satisfying σ(x) c > 0 for almost every x Ω. Let u H (Ω) be the unique solution to the following elliptic Dirichlet problem: (.) (.) σ u = 0 in Ω, u Ω = φ H / ( Ω). The inverse conductivity problem consists on recovering σ from the Dirichlet-to- Neumann (DN) map or voltage-to-current map defined by Λ σ : φ σ u. ν Ω Here ν is the unit outer normal to the boundary. Note that the map Λ : σ Λ σ is nonlinear. The inverse conductivity problem is related to many practical applications, including the medical imaging technique called electrical impedance tomography (EIT). There one attaches electrodes to the skin of a patient, feeds electric currents into the body and measures the resulting voltages at the electrodes. Repeating the measurement with several current patterns yields a current-to-voltage data matrix that can be used to compute an approximation Λ δ σ to Λ σ. Since different organs and tissues have different conductivities, recovering σ computationally from Λ δ σ amounts to creating an image of the inner structure of the patient. This technique has some commercial realizations, including breast cancer detection [] and monitoring lung function. See [6, 0] for more information on EIT and its applications. Another application of the nonlinear Fourier transform is the Novikov-Veselov (NV) equation. It is a nonlinear evolution equation generalizing the celebrated Korteweg-de Vries (KdV) equation into dimension (+). There has been significant recent progress in linearizing the NV equation using inverse scattering methods, see [4], [3], [9], [7].
3 NONLINEAR FOURIER ANALYSIS 3 Recovering σ from Λ δ σ is a nonlinear and ill-posed inverse problem, whose computational solution requires regularization. Several categories of solution methods have been suggested and tested in the literature; in this work we focus on so-called D-bar methods based on complex geometrical optics (CGO) solutions. There are three main flavors of D-bar methods for EIT: Schrödinger equation approach for twice differentiable σ. Introduced by Nachman in 996 [8], implemented numerically in [30, 5, 7, 6, 0, 6]. First-order system approach for once differentiable σ. Introduced by Brown and Uhlmann in 997 [7], implemented numerically in [8,, 4]. Beltrami equation approach assuming no smoothness (σ L (Ω)). Introduced by Astala and Päivärinta in 006 [6], implemented numerically in [5, 4, 3, 6]. The assumption σ L (Ω) was the one originally used by Calderón in [8]. Using these approaches, a number of conditional stability results have been studied for the Calderón problem. The most recent results in the plane were obtained by Clop, Faraco and Ruiz in [], where stability in L -norm was proven for conductivities on Lipschitz domains in the fractional Sobolev spaces W α,p with α > 0, < p <, and in [], where the Lipschitz condition on the boundary of the domain was removed. In dimension d 3, conditional stability in Hölder norm for just C +ε conductivities on bounded Lipschitz domains was proved by Caro, García and the third author in [9] using the method presented in [3]. The three aforementioned D-bar methods for EIT in the two-dimensional case are based on the use of nonlinear Fourier transforms specially adapted to the inverse conductivity problem. Schematically, the idea looks like this: Λ σ Data σ(z) Transform τ(k) Inverse transform Nonlinear frequency domain (k-plane) σ(z) Spatial domain (z-plane) The main point above is that the nonlinear Fourier transform can be calculated from the infinite-precision data Λ σ, typically via solving a second-kind Fredholm boundary integral equation for the traces of the CGO solutions on Ω. In practice one is not given the infinite-precision data Λ σ, but rather the noisy and finite-dimensional approximation Λ δ σ. Typically all we know about Λ δ σ is that
4 4 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Λ δ σ Inverse transform Lowpass k < 4 k-plane Inverse transform z-plane σ(z) (a) (b) Figure : Schematic illustration of the nonlinear low-pass filtering approach to regularized EIT. This image corresponds to the so-called shortcut method defined below. The simulated heart-and-lungs phantom σ(z) (bottom left) gives rise to a finite voltage-to-current matrix Λ δ σ (orange square), which can be used to determine the nonlinear Fourier transform. Measurement noise causes numerical instabilities in the transform (irregular white patches), leading to a bad reconstruction (a). However, multiplying the transform by the characteristic function of the disc k < 4 yields a lowpass-filtered transform, which in turn gives a noise-robust approximate reconstruction (b). Λ σ Λ δ σ Y < δ for some (known) noise level δ > 0 measured in an appropriate norm Y. Most CGO-based EIT methods need to be regularized by a truncation k < R in the nonlinear frequency-domain as illustrated in Figure. The regularization step results in a smooth reconstruction. This smoothing property of the nonlinear low-pass filter resembles qualitatively the effect of linear low-pass filtering of images. In particular, the smaller R is, the blurrier the reconstruction becomes. The cut-off frequency R is determined by the noise amplitude δ, and typically R cannot exceed 7 in practical situations. However, it is interesting to understand how the conductivity is represented by its nonlinear Fourier transform. To study this numerically, we compute nonlinear Fourier transforms in large discs such as k < 60 and compute low-pass filtered reconstructions using inverse nonlinear Fourier transform and observe the results.
5 NONLINEAR FOURIER ANALYSIS 5 In the computational experiments we show in this paper, the truncated scattering transform is computed with a very large cutoff frequency through the Beltrami equation solver. In practice, it is not possible to obtain such scattering data from boundary measurements unless we had unrealistically high-precision EIT measurements available. R = 0 R = 40 R = 50 True Figure : Convergence of nonlinearly low-pass filtered conductivity to the true discontinuous conductivity. Note the ringing artefacts reminiscent of the Gibbs effect of linear low-pass filtering. This picture illustrates the so-called shortcut method, and it is the first-ever computation of this kind with cutoff frequencies larger than 0. We compare two reconstruction methods based on the use of the CGO solutions (.): the low-pass transport matrix method implemented in [3], and a shortcut method based on solving a D-bar equation. The latter method does not have rigorous analysis available yet, but it is formally similar to the regularized EIT reconstruction algorithm studied in [8, 30, 5, 7, 6, 9, 0]. Our new computational findings can be roughly summarized by the following two points. First, both methods can approximate piecewise constant conductivities better and better as the cutoff frequency R increases, and there seems to be a Gibbs-like phenomenon producing ringing artifacts. Second, the transport matrix method loses accuracy away from a (freely chosen) pivot point located outside of Ω, whereas the shortcut method produces reconstructions with more uniform quality. Figure shows some of our numerical results via the shortcut method for a nonsymmetric conductivity distribution. The rotationally symmetric examples presented in Section 4. provide numerical evidence of the fact that discontinuous conductivities can be reconstructed with the shortcut method more and more accurately when R tends to infinity.
6 6 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Remark : All the numerical experiments we intend to show in this paper have been computed using either a 9 Gigabyte computer with Linux or parallel computation provided by Techila.. The nonlinear Fourier transform In dimension two Calderón s problem [8] was solved using complex geometrical optics (cgo) solutions by Astala and Päivärinta [6]. In the case of L - conductivities σ the cgo solutions need to be constructed via the Beltrami equation (.) z f µ = µ z f µ, with f µ (z, k) = e ikz ( + ω(z, k)), where µ := σ ( ), ω(z, k) = O as z. + σ z Here z = z + iz C, z = ( / z + i / z )/ and k is a complex parameter. More specifically, for σ L (C) with c σ(z) c almost every z C and σ outside the unit disc, and for κ < with µ(z) κ a.e. z C, Theorem 4. in [6] establishes that for any < p < + /κ and complex parameter k, there exists a unique solution f µ (, k) in the Sobolev space W,p loc (C) such that the Beltrami equation (.) z f µ (z, k) = µ(z) z f µ (z, k), for z C, holds with the asymptotic formula ( )) (.3) f µ (z, k) = e ( ikz + O z as z. In addition, f µ (z, 0). Assuming just σ L (C) bounded away from zero and infinity and writing f µ for the solution to the Beltrami equation with coefficient µ and the same asymptotics (.3), the scattering transform τ of Astala and Päivärinta (AP) is given by the formula (.4) τ(k) = z (ω(z, k) ω (z, k)) dz dz, π C where ω(z, k) := e ikz f µ (z, k), ω (z, k) := e ikz f µ (z, k). Proposition 6.3 in [6] states among other things that the scattering transform τ is bounded. Indeed, τ(k) for all k C. Defining the functions as follows for k, z C: (.5) (.6) h + (z, k) := / (f µ (z, k) + f µ (z, k)), h (z, k) := i/ ( f µ (z, k) f µ (z, k) ), Techila Technologies Ltd is a privately held Finnish provider of High Performance Computing middleware solutions.
7 and NONLINEAR FOURIER ANALYSIS 7 u (z, k) := h + (z, k) i h (z, k) = Ref µ (z, k) + i Imf µ (z, k), u (z, k) := i (h + (z, k) + i h (z, k)) = Imf µ (z, k) + i Ref µ (z, k), it turns out that u, u are the unique solutions to the following elliptic problems: (.7) (.8) σ u = 0, u (z, k) = e ikz( + O ( /z )), σ u = 0, u (z, k) = e ikz( i + O ( /z )), as z. The complex parameter k above can be understood as a nonlinear frequencydomain variable. As well as the function τ(k), defined in (.4), plays the role of the nonlinear Fourier transform in the shortcut method described below, it is played by the following function ν z0 (k) := i h (z 0, k) h + (z 0, k), in the low-pass transport matrix method presented in Section 3., where z 0 is a pivot point outside Ω. As shown in [5, 6], one can determine different functions of k from the knowledge of the infinite-precision Dirichlet-to-Neumann map Λ σ. Furthermore, a practical algorithm was introduced in [3] for recovering σ approximately from a noisy measurement of Λ σ. This algorithm involves a truncation operation (nonlinear low-pass filtering) in the frequency-domain, and has applications in EIT. The scattering transform used in the Schrödinger equation approach by Nachman in [8] is defined as (.9) t(k) = e i(kz+k z) q(z)m(z) dz, R where q = σ / σ / and σ is a strictly positive-real valued function in W,p (R ) for some p > with σ on the unit disc near the boundary and outside the unit disc, and m solves the equation m = g k (qm), where denotes convolution of functions over the plane and g k the Faddeev fundamental solution of the operator 4ik with pole at the origin. For smooth enough σ, one can prove by a straightforward uniqueness argument that the function t and the AP scattering transform τ are related to each other through the following formula: (.0) t(k) = 4πi k τ(k). In the rest of the paper, for non-smooth conductivities t denotes the function defined through formula (.0). Why are t and τ called nonlinear Fourier transforms? There is admittedly some abuse of terminology involved, but the main reason is this. If one replaces m in (.9) with, the result is the linear Fourier
8 8 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN transform of q at ( k, k ) R. But m is only asymptotically close to and actually depends on q. The following theorem holds: Theorem.. Assume σ L (R ) is real-valued, radially symmetric, bounded away from zero and σ outside the unit disc. Then the function t is real-valued and rotationally symmetric. Proof. The rotational symmetry of t follows from the property t(k) = t(ηk) for any η, k C with η =, which is tantamount to τ(k) = η τ(ηk). By uniqueness of the asymptotic boundary value problem (.)-(.3), one can check that f µ (z, ηk) = f µ (ηz, k). Therefore, ω(z, ηk) = ω(ηz, k). For µ the same equalities follow. Using definition (.4) we deduce τ(k) = η τ(ηk). Thanks to t(k) = t(ηk), real-valuedness of t follows from t(k) = t(k) and this fact is equivalent to seing τ(k) = τ( k) due to the identity τ( k ) = η τ(η k ) with η =. Again, by uniqueness of (.)-(.3), it is straightforward to prove f µ (z, k) = f µ (z, k), and ω(z, k) = ω(z, k). deduce τ(k) = τ( k). Attaining the same identity for µ and by (.4) we 3. Inverse transforms and reconstruction algorithms 3.. The low-pass transport matrix method. Given noisy EIT data Λ δ σ, we can evaluate the traces of M µ (, k) = + ω(, k) at the boundary Ω for k < R with some R > 0 depending on the noise level δ. This is achieved by the numerical method described in Section 4. of [3] based on solving the boundary integral equations studied in [5]. Once M µ is known on the boundary, so is f µ (, k) Ω. The fact that µ is supported in Ω implies that f µ (, k) is analytic outside the unit disc. Thus, the coefficients in the Fourier series for the trace of f µ on Ω can also be used to expand f µ as a power series outside Ω. Choose then a point z 0 R \ Ω. As explained above, we know f µ (z 0, k) and f µ (z 0, k) for any k < R. Construct the function (3.) ν (R) z 0 (k) := i h (z 0, k) h + (z 0, k) for k < R, 0 for k R. We next solve the truncated Beltrami equations (3.) (3.3) k α (R) = ν (R) z 0 (k) k α (R), k β (R) = ν (R) z 0 (k) k β (R),
9 with solutions represented in the form (3.4) (3.5) NONLINEAR FOURIER ANALYSIS 9 α (R) (z, z 0, k) = exp(ik(z z 0 ) + ε(k)), β (R) (z, z 0, k) = i exp(ik(z z 0 ) + ε(k)), where ε(k)/k 0 and ε(k)/k 0 as k. Requiring (3.6) α (R) (z, z 0, 0) = and β (R) (z, z 0, 0) = i fixes the solutions uniquely. To compute α (R) we first find solutions η and η to the equation (3.7) k η j = ν (R) z 0 (k) k η j, with asymptotics (3.8) (3.9) η (k) = e ik(z z 0) ( + O(/k) ), η (k) = i e ik(z z 0) ( + O(/k) ), respectively, as k. Such solutions exist and are unique by [6, Theorem 4.]. The solutions η j are complex valued, but pointwise real-linearly independent by [, Theorem 8.4.]. Hence there are constants A, B R such that We now set A η (0) + B η (0) =. α (R) (z, z 0, k) = A η (k) + B η (k). Then α (R) (z, z 0, k) satisfies the equation (3.) and the appropriate asymptotic conditions. Let α (R) µ denote the solution to (3.) with the condition (3.4) so that the coefficient ν z (R) 0 is computed interchanging the roles of µ and µ. Then β (R) µ = i α (R) µ. For computational construction of the frequency-domain cgo solutions η j satisfying ( ), see [3, Section 5.] and [5]. Fix any nonzero k 0 C and choose any point z inside the unit disc. Denote α (R) = a (R) + ia (R) and β (R) = b (R) + ib (R). We can now use the truncated transport matrix (3.0) T (R) = T (R) z,z 0,k 0 := to compute (3.) ( a (R) a (R) b (R) b (R) u (R) (z, k 0 ) = a (R) u (z 0, k 0 ) + a (R) u (z 0, k 0 ), u (R) (z, k 0 ) = b (R) u (z 0, k 0 ) + b (R) u (z 0, k 0 ). In [6] it is proved that for R = the functions u (R) (z, k 0 ), u (R) (z, k 0 ) solve the problems (.7), (.8) with k = k 0, respectively. )
10 0 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN The truncation in (3.) can be interpreted as a nonlinear low-pass filter in the k-plane. This is where the term low-pass transport matrix originates. We know the approximate solutions u (R) j (z, k 0 ) for z Ω and one fixed k 0. We use formulas u (R) = h (R) + ih (R), u (R) = i(h (R) + + ih (R) ) and (.5)-(.6) (with z Ω) to connect u (R), u (R) with f µ (R), f (R) µ. Define (3.) µ (R) (z) = f (R) µ (z, k 0 ). f µ (R) (z, k 0 ) Because of the truncation, the function f (R) µ (, k 0 ) is smooth and there are no difficulties in the numerical differentiation in (3.). Further, in our computations the denominator in (3.) was always strictly positive, so there was no danger of division by zero. Finally we reconstruct the conductivity σ approximatively as (3.3) σ (R) = µ(r) + µ (R). 3.. The shortcut method. Fix R > 0 corresponding to a certain noise level δ in the EIT data. For all k < R, we compute the traces of M µ (, k) and M µ (, k) by solving the boundary integral equation aforementioned in Section 3.. The analyticity of M µ (, k), M µ (, k) outside the unit disc allows to develop these factors as follows M µ (z, k) = + a+ (k) z + a+ (k) +..., M z µ (z, k) = + a (k) + a (k) +..., z z for z >. In Section 5 of [6] it is proved that the scattering transform τ satisfies τ(k) = ( a+ (k) a (k) ), k C. Notice that this formula is consistent with (.4). Let τr δ (k) be a numerical approximation to τ(k) such that τr δ (k) = 0 for k > R. For all z Ω, we solve the D-bar equation (3.4) k m δ R(z, k) = i τr(k)e δ z (k) m δ R (z, k) with m δ R (z, ) Lr (R ) L (R ), for some r >. Finally, we reconstruct σ approximately as σ(z) (m δ R (z, 0)). This method is proven to work for C conductivities σ by computing the approximation τ δ R (k) in a different way, namely the methods described in [8, 30, 5, 7, 6, 9, 0]. Nevertheless, there is no theoretical understanding of what happens when σ is nonsmooth and only essentially bounded.
11 NONLINEAR FOURIER ANALYSIS 3.3. Periodization of the nonlinear inverse transform. The algorithm we use to solve equation (3.4) is presented in [] and is a modification of the method introduced by Vainikko in [3] for the Lippmann-Schwinger equation related to the Helmholtz equation. Such algorithm is based in reducing the equation (3.5) m R (z, k) = + F (z, k ) m R (z, k ) dk, π B R k k to a periodic one. Let us see how this reduction is done. Here F (z, k) := i τ(k)e k (z) and e k (z) := e ikz+ikz. Take ɛ > 0 and set s = R + 3ɛ. Define Q = [ s, s). Choose an infinitely smooth cutoff function η C 0 (R ) verifying η(k) = and 0 η(k) for all k C. { for k < R + ɛ, 0 for k R + ɛ, Definition 3.. We say a function f defined on the plane is s-periodic if f(k) = f(k + s(j + ij )) for any k C, j, j Z. That is to say, we mean that f is s-periodic in each coordinate. Notation. L p (Q) denotes the space of s-periodic functions in L p loc (R ) for p. Note that L (Q) L (R ). For any measurable subset A of R we denote by χ A both the characteristic function of A and the multiplier operator with such function in k. We write E 0 to denote the extension by zero outside the torus Q defined as { f(z, k), if k Q, E 0 f(z, k) = 0, if k C \ Q, for any f(z, ) L (Q). Define a s-periodic approximation to the principal value /(πk) by setting β(k) = η(k)/(πk) for k Q and extending periodically: β(k + j s + ij s) = η(k) for k Q, j, j Z. πk Define a periodic approximate solid Cauchy transform in k by C 0 f(z, k) = ( β f)(z, k) = β(k k )f(z, k ) dk, where denotes convolution on the torus. We shall use another periodization of the Cauchy transform given by Cf(z, k) = Eχ DR+ɛ C0 f(z, k), where E denotes the periodic extension operator defined as Ef(z, k + s(j + ij )) = f(z, k), for any k Q, j, j Z. Q
12 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Define the multiplication operator F R by F R f(z, k) = χ BR (k)f (z, k)f(z, k) if k C and its periodization as F R f(z, k + j s + i j s) = F R f(z, k), for j, j Z, k Q. Let us denote the complex conjugation operator by c, i.e. c(f) = f. We periodize the complex conjugation operator c in the same way. We shall need some boundedness properties of the so-called Cauchy transform. The Cauchy transform, we denote it by, is defined as the following weakly singular integral operator: f(k) = f(k ) π R k k dk, where k = k + ik and dk dk denotes the Lebesgue measure on R. Theorem in [] states that (3.6) φ L q (C) C (q )( q) φ L q (C) for some constant C, < q < and q = q/( q). For any < q < and /q + /q = the space L q (C) L q (C) with the norm L q + L q is a Banach space. Theorem 4.3. in [] claims the following: Theorem 3.. When < q < the Cauchy transform maps continuously the space L q (C) L q (C) into the space of continuous functions on the plane vanishing at infinity with the sup norm and satisfies the estimate φ ( q) /( φ L L (C) q (C) + φ ) L q (C). then It is easy to see that if < q < < q < but not necessarily q + q = (3.7) φ ( ) φ L L q (C) + φ L q (C), (C) and φ is a continuous function assuming the right hand side bounded, where the implicit constant just depends on q, q. We conclude this section with the Theorem 3.. Its proof uses the aforementioned properties of the Cauchy transform and is quite similar to the arguments presented in Sections and 5.4. of the book [6]. Some of the techniques for Theorem 4. in [8] are involved. Theorem 3.. Assume c σ c a.e. in the plane and σ outside the unit disc. Let z C. There exists a unique s-periodic (in k) solution m R (z, k) to (3.8) m R (z, k) = + C F R c m R (z, k), a.e. k Q, with m R (z, ) L r (Q) for some < r <. Furthermore, the solutions of (3.5) and (3.8) agree on the disc of the k-plane with radius R: m R (z, k) = m R (z, k) for k < R.
13 NONLINEAR FOURIER ANALYSIS 3 4. Computational results Once scattering data with big cutoff frequencies are generated, we test the nonlinear inversion procedure consisting of computing the approximate conductivity from the truncated scattering transform by solving (4.) k m R (z, k) = χ BR (k)( i)τ(k)e k (z) m R (z, k), with (4.) m R (z, ) L (R ) L r 0 (R ), for some r 0 (, ). The unique solution m R (z, k) to (4.) with the normalization condition (4.) can be obtained by solving (3.5). Finally, the approximate reconstruction σ R (z) is obtained by doing σ R (z) := m R (z, 0) for z <. 4.. Nonlinear Gibbs phenomenon in radial cases. In this section we test the shortcut method for two radially symmetric conductivities as follows: Firstly, we solve numerically the Beltrami equation and obtain reliable scattering data over a mesh for a real interval of the form [0, R] with R > 0. Secondly, from such scattering data the approximate reconstruction is computed by solving numerically the equation (3.8) on a mesh for the real interval [0, ]. It suffices to confine such computations to real positive numbers since the results on the whole disc can be obtained by simple rotation with respect to the origin by symmetry. The pictures for such reconstructions corresponding to the largest cutoff frequencies exhibit certain artifacts near the jump discontinuities of the actual conductivity quite similar to the Gibbs phenomenon for the truncated linear Fourier transform Example σ. We consider the rotationally symmetric conductivity σ on the plane defined as σ (z) = if z < and σ (z) = if z. Figure 3 shows its representation. Computing the scattering transform The numerical evaluation of the AP scattering transform τ(k) was computed for this example conductivity σ. The used algorithm generates for each fixed k cgo solutions for the Beltrami equation in the z-plane using the Huhtanen- Perämäki preconditioned solver ([5]) by executing standard C-linear GMRES. To do so a z-grid is required. Such a z-grid is determined by a size parameter m z so that mz mz equidistributed points are taken in the z-plane within the square [ s z, s z ) [ s z, s z ) with s z =. and step h z = s z / mz. Let us denote the numerical computation of τ with a z-grid size parameter m z by τ mz. We know that for radial, real-valued conductivities the function t is rotationally symmetric too (Theorem.). Therefore, we can restrict the numerical computation of t to positive values of the k-plane. The chosen k-grid was all the points in the interval [0, 50] with step h = 0.. Thus, the scattering transform is computed farther away from the origin than any previous computation.
14 4 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN The vector τ mz corresponding to the aforementioned k-grid was computed for both m z = 0 and m z =. Next, the function t mz defined as t mz (k) = 4πi k τ mz (k), was evaluated on the same k-grid in [0, 50] for m z = 0, from the AP scattering transform using identity (.0). One can consider the approximate scattering transform τ (with m z = ) reliable since its accuracy is high in a double sense. On one hand, Figure 4 shows the lack of significant differences between t 0 and t for k close to the origin. As a matter of fact, using the notation (4.3) E B := max Im( τ 0(k) τ (k)) k B max Im( τ (k)) k B for a subinterval B of [0, 50], it follows that 00% E [0,0] = %, E [0,40] = 4.058%, E [40,60] = 9.838%. On the other hand, Re( τ ) is very small. By Theorem. and formula (.0), for a rotationally symmetric, real-valued conductivity the scattering transform τ restricted to real numbers has to be pure-imaginary valued. We obtain a nonzero real part Re( τ mz ) for both m z values due to precision errors. In particular, Re( τ ) Note that one can extend t to the disc centered at the origin with radius 50 in the k-plane by simple rotation since t (k ) = t (k ) if k = k. σ ( z ) σ (z) z Figure 3: Rotationally symmetric conductivity example σ. Left: plot of σ as a function of z ranging in the unit disc with jump discontinuities along the curve z =. Right: Profile of σ as a function of z ranging in [0,] with a jump discontinuity at z =.
15 Solving the D-bar equation NONLINEAR FOURIER ANALYSIS 5 We compute a number of reconstructions of σ using the D-bar method truncating the approximate scattering transform τ. This way the computation is optimized by restricting the degrees of freedom to the values of τ (k) at grid points satisfying k < R, for some positive number R called the cutoff frequency. The program constructs the set of evaluation points in the k-plane for which the D-bar equation is solved for every fixed z, as a k-grid of m k m k points within the square [ s k, s k ) [ s k, s k ) with s k =.3 R and step h k = s k / m k. Let σ (z; m k, R) denote the shortcut reconstruction of σ with size parameter m k and truncation radius R. Again we can restrict the computation (of σ (z; m k, R) in this case) to positive values in the z-domain and afterwards extend to the plane by rotation. We evaluated σ (z; m k, R) for each point z in a grid of a certain number N x of equispaced points in the real interval [0,]. The vector σ (z; m k, R) was computed for R = 5, 0, 5, 0, 5, 30, 35, 40, 50,..., 50. For R = 5, 0 we took m k = 8, 9 and for the rest of R values we chose m k =. These are by far the largest nonlinear cut-off frequencies studied numerically so far. In Figure 5 the reader can see three of these reconstruction examples. 30 M = 0 M = k k Figure 4: Comparison between t 0 and t for σ. Left: Profiles of the real parts of t 0 (thin dotted line) and t (thick solid line). Right: Profile of the difference Re( t 0 t ). In both pictures k ranges in [0,50] with step h = 0..
16 6 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Real part; Max =.9 Min = Im 4.67e Real part; Max =.3 Min = Im 3.6e x mk = 8 R = 5 Nx = Real part; Max =.06 Min = x mk = R = 40 Nx = 300 TIME = h 0 m TECHILA Im 4.89e x mk = R = 50 Nx = 50 TIME = 45 m TECHILA Figure 5: Three shortcut method reconstruction examples computed from a scattering transform corresponding to σ. Each picture shows the real part of the approximate reconstruction (thick line) and the original conductivity σ (thin line). In addition, the following information is shown: maximum and minimum values of the real part of the reconstruction, maximum of the absolute value of the imaginary part of the reconstruction (error), taken parameters m k, R, N x and computation time if possible. The examples for m k = and R = 40, 50 were computed using grid computation provided by Techila. Since the shortcut reconstruction for σ is real-valued, the obtained imaginary parts are due to precision errors. Such errors are displayed multiplied by 0 8 in the following tables: m k R error m k R error where error = max z Im( σ (z; m k, R)).
17 NONLINEAR FOURIER ANALYSIS Example σ. Now, we repeat the same experiments for another rotationally symmetric conductivity example σ defined as follows: σ (z) = if z < 0., 0. < z < 0.3 or 0.4 < z < and σ (z) = otherwise. Figure 6 shows its profiles. σ ( z ) σ (z) z Figure 6: Rotationally symmetric conductivity example σ. Left: plot of σ as a function of z ranging in the unit disc with jump discontinuities along the curves z = 0., z = 0., z = 0.3, z = 0.4, z =. Right: Profile of σ as a function of z ranging in [0,] with jump discontinuities at z = 0., z = 0., z = 0.3, z = 0.4, z =. Computing the scattering transform For conductivity σ the scattering transform τ (with m z = ) turns out to be very accurate too. Indeed, t 0 and t corresponding to σ are quite similar for k small in view of Figure 7 and the following facts (see notation (4.3)): E [0,0] = %, E [0,40] = %, E [40,60] = %, for σ, Re( τ ) , for σ. Solving the D-bar equation Likewise the vector σ (z; m k, R) was computed for R = 5, 0, 5, 0, 5, 30, 35, 40, 50,..., 50. For R = 5, 0 we took m k = 8, 9 and for the rest of R values we chose m k =. All these reconstructions were computed using parallel computation provided by Techila. Figure 8 shows some of these numerical results. Below we list the errors error = max z Im( σ (z; m k, R)) : m k R error
18 8 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN m k R error M = 0 M = k k Figure 7: Comparison between t 0 and t for σ. Left: Profiles of the real parts of t 0 (thin dotted line) and t (thick solid line). Right: Profile of the difference Re( t 0 t ). In both pictures k ranges in [0,50] with step h = 0.. Note that the vertical axis scales are different from Figure 4.
19 NONLINEAR FOURIER ANALYSIS 9 Real part; Max =.56 Min = Im.6e 007 Real part; Max =. Min = Im 3.97e x mk = 8 R = 0 Nx = 50 TIME = m TECHILA Real part; Max =.5 Min = 0.93 Im.48e x mk = R = 5 Nx = 50 TIME = h 37 m TECHILA Real part; Max =.4 Min = 0.95 Im.97e x mk = R = 60 Nx = 50 TIME = 30 m TECHILA x mk = R = 50 Nx = 50 TIME = 9 m TECHILA Figure 8: Four shortcut method reconstruction examples computed from a scattering transform corresponding to σ. Each picture shows the real part of the approximate reconstruction (thick line) and the original conductivity σ (thin line). In addition, the following information is shown: maximum and minimum values of the real part of the reconstruction, maximum of the absolute value of the imaginary part of the reconstruction (error), taken parameters m k, R, N x and computation time. Axes scale is the same in the four profiles. These examples were computed using grid computation provided by Techila.
20 0 K. ASTALA, L. PA IVA RINTA, J. M. REYES AND S. SILTANEN 4.. Comparison of the two methods in nonsymmetric cases. We apply both reconstruction methods to some discontinuous non-radial checkerboardstyle conductivities without assuming EIT data. We choose the conductivity examples σ3, σ4 shown in Figure 9. Both take value near the boundary of the unit disc. Notice that the contrasts of these examples σ3, σ4 (defined as the difference between the maximum and the minimum) are.5 and.8, respectively a) b) c) d) Figure 9: Checkerboard style conductivities σ3, σ4 : a), c) show the picture for σ3, σ4, resp., using mesh Matlab function. b), d) show the values for σ3, σ4, resp., over each tile.
21 NONLINEAR FOURIER ANALYSIS Figures, 3 show the approximate reconstructions for both examples and certain cut-off frequencies R using both methods. The relative errors with sup and l norms over the matrix of chosen points in the unit disc are pointed out. We denote such errors as follows: sup = σ σ l σ l 00 %, sqr = σ σ l σ l 00 %, where σ denotes the actual conductivity and σ the approximate reconstruction. Low-pass transport matrix method In order to test the low-pass transport matrix method for each conductivity, we obtain the CGO solutions f µ (z 0, k), f µ (z 0, k) directly from the known conductivities without simulating EIT data. The point z 0 is outside the unit disc and we take the k-grid of m k m k equispaced points in the square [ R, R) with step h k = R/ m k and cut-off frequencies R 50. Next, transportation of f µ (z 0, k), f µ (z 0, k) to the unit disc and the final reconstruction σ (R) are performed following the steps described in Section 3.. To this end, in addition to the k-grid, a grid in the z-variable of mz mz equidistributed points in the square [ s z, s z ), with step h z = s z / mz, is required. For σ 3, σ 4 we choose z 0 = ( , ), z 0 = (0,.656), respectively. For σ 3 we take m z = 7, m k = 8 and for σ 4 m z = 7, m k = 7. The reconstructions of σ 3 could be computed for R 50 but not for σ 4 with R > 0. Shortcut method Concerning the shortcut method, firstly the approximate scattering transform, denote it by τj SC for σ j (j = 3, 4), is computed directly from the known conductivity through the Beltrami equation solver using formula (.4). The code evaluates the approximation over the points inside the disc of radius R (cut-off frequency) from a k-grid of m k m k equispaced points in the square [ R, R]. Another z-grid with mz mz points in [ s z, s z ) (s z > ) is involved here. To compute τ3 SC we take R = 50, m k = 8, m z = 0, and for τ4 SC we choose R = 50, m k = 7, m z = 7. Let us remark that this choice of R enables reconstructions corresponding to cut-off frequencies less than or equal to 50. Secondly, Figures 0, show plots of τj SC and t SC j, where j = 3, 4 and t SC j is computed from τj SC by formula (.0). Finally, we evaluate the approximate reconstruction by solving the D-bar equation with the truncated scattering transform τj SC as a coefficient. The reconstructions are computed over the points of the same kind of z-grid as above belonging to the unit disc with size parameter m z. The code computes the solution to the D-bar equation on every point in the z-grid independently and uses a number of k-grids with size parameter m k and cut-off frequency R.
22 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN These reconstructions were evaluated with m z = 7 and some R values with R 50. In addition, with respect to σ 3 we took m k = 9 and for σ 4, m k = 7. Notice that the case R = 50 was correctly computed for both conductivities σ 3, σ 4 using the shortcut method, but for σ 4 the highest R-value through the low-pass transport matrix method was R = 0. a) τ SC 3 for σ 3 (R = 50, m k = 8, m z = 0). b) τ SC 4 for σ 4 (R = 50, m k = 7, m z = 7). Figure 0: a) Plot of real (left) and imaginary (right) parts of the scattering transform τ3 SC for the conductivity example σ 3. b) Idem for τ4 SC corresponding to σ 4. a) t SC 3 for σ 3 (R = 50, m k = 8, m z = 0). b) t SC 4 for σ 4 (R = 50, m k = 7, m z = 7). Figure : a) Plot of real (left) and imaginary (right) parts of the scattering transform t SC 3 for the conductivity example σ 3. b) Idem for t SC 4 corresponding to σ 4.
23 NONLINEAR FOURIER ANALYSIS 3 z 0 -style shortcut R = 0 R = 40 R = 50 sup(z 0) = 843.7% sqr(z 0) = 4.8% sup(shct) = 49.8% sqr(shct) = 9.6% sup(z 0) = % sqr(z 0) = 73.7% sup(shct) = 49.4% sqr(shct) = 4% sup(z 0) = 667.6% sqr(z 0) = 00% sup(shct) = 49.% sqr(shct) =.5% σ 3 z 0 ( 7/8, 7/8) Figure : Comparison of reconstructions for σ 3 by the low-pass transport matrix method (left) and the shortcut method (right) using different cut-off frequencies R. The chosen point z 0 appears in black. The true conductivity is represented at the bottom row.
24 4 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN z 0 -style shortcut R = 0 sup(z 0) = 9.% sqr(z 0) = 49.% sup(shct) = 65.9% sqr(shct) = 9.6% R = 5 sup(z 0) = 44.6% sqr(z 0) = 50.3% sup(shct) = 67% sqr(shct) = 7.5% R = 0 sup(z 0) = 54.9% sqr(z 0) = 5.% sup(shct) = 69.8% sqr(shct) = 5.8% σ 4 z 0 (0,.66) Figure 3: Comparison of reconstructions for σ 4 by the low-pass transport matrix method (left) and the shortcut method (right) using different cut-off frequencies R. The chosen point z 0 appears in black. The true conductivity is represented at the bottom row.
25 NONLINEAR FOURIER ANALYSIS Numerical evidence of the transport matrix efficiency. In addition, further numerical experiments were made aimed at comparing the actual complex geometric optics solution to the Beltrami equation z f µ (z, k 0 ) = µ(z) z f µ (z, k 0 ), for z Ω, k 0 =, and the transported solution f µ (z, k 0 ) to the unit disc from f µ (z 0, k 0 ) and f µ (z 0, k 0 ) at certain z 0 C \ Ω. Figure 4 shows some pictures for f µ (z, k 0 ) and f µ (z, k 0 ) corresponding to the conductivity example σ 4. Here z <, k 0 and the pivot point outside the unit disc is z 0 (0,.66). To compute f µ the transport matrix was generated with a cut-off frequency R = 0 and we took the size parameters as follows: m k = 7, m z = 7, s z =.5. The relative errors with sup and l norms f µ (, k 0 ) f µ (, k 0 ) f µ (, k 0 ) f µ (, k 0 ) sup := l 00 %, sqr := l 00 % f µ (, k 0 ) l f µ (, k 0 ) l are (4.4) sup = 0.3%, sqr = 0.8%, respectively.
26 6 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Real part Imaginary part a) b) c) Figure 4: Comparison of f µ (z, k 0 ) and f µ (z, k 0 ) for the conductivity σ 4 with R = 0. In the first row we show real and imaginary parts of both functions. In a) the real part of the actual solution f µ (z, k 0 ) is depicted on the left and the real part of the transported f µ (z, k 0 ) on the right (where the point z 0 is represented in black). In b) the same information is showed for the imaginary part of f µ (z, k 0 ) and f µ (z, k 0 ). The last row c) shows the difference f µ (z, k 0 ) f µ (z, k 0 ). The real part is exhibited on the left and the imaginary part on the right. Again, the pivot point z 0 appears in black. All the pictures were generated under the same colormap using Matlab.
27 NONLINEAR FOURIER ANALYSIS 7 5. Conclusion Regarding the radially symmetric examples studied in Section 4., the conclusion is clear: numerical evidence suggests that discontinuous conductivities can be reconstructed with the shortcut method more and more accurately when R tends to infinity. The errors in the reconstructions are very similar to the Gibbs phenomenon observed with truncated linear Fourier transforms. The numerical study of the nonsymmetric examples comparing both reconstruction types leads us to two observations: - For these discontinuous conductivity examples, the shortcut method generates considerably better reconstructions than the transport matrix method. - The approximate solution computed by the transport matrix method is reliable on an area within the unit disc close enough to the selected z 0 point outside the unit disc. Finally, since errors (4.4) are small and plots for both functions f µ (z, k 0 ), f µ (z, k 0 ) are similar in view of Figure 4, we have numerical evidence supporting the fact that the transport matrix method transports well and the worse results by the transport matrix method are explained by the final algebraic steps, including numerical differentiation in (3.), just after the computation of the solutions u (R) (z, k 0 ), u (R) (z, k 0 ) in (3.) through the transport matrix itself. Acknowledgements. This work was supported by the Finnish Centre of Excellence in Inverse Problems Research 0-07 (Academy of Finland CoE-project 505). In addition, KA was supported by Academy of Finland, projects and 34757, and LP was supported by ERC-00 Adv. Grant, InvProb. The third author was also partially supported by Academy of Finland (Decision number 4075) and Ministerio de Ciencia y Tecnología de España (project MTM ).
28 8 K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN References [] M. Assenheimer, O. Laver-Moskovitz, D. Malonek, D. Manor, U. Nahaliel, R. Nitzan, and A. Saad, The T-SCAN technology: electrical impedance as a diagnostic tool for breast cancer detection, Physiological measurement, (00), pp. 8. [] K. Astala, T. Iwaniec, and G. Martin, Elliptic partial differential equations and quasiconformal mappings in the plane, vol. 48 of Princeton Mathematical Series, Princeton University Press, Princeton, NJ, 009. [3] K. Astala, J. Mueller, L. Päivärinta, A. Perämäki, and S. Siltanen, Direct electrical impedance tomography for nonsmooth conductivities, Inverse Problems and Imaging, 5 (0), pp [4] K. Astala, J. Mueller, L. Päivärinta, and S. Siltanen, Numerical computation of complex geometrical optics solutions to the conductivity equation, Applied and Computational Harmonic Analysis, 9 (00), pp [5] K. Astala and L. Päivärinta, A boundary integral equation for Calderón s inverse conductivity problem, in Proc. 7th Internat. Conference on Harmonic Analysis, Collectanea Mathematica, 006. [6] K. Astala and L. Päivärinta, Calderón s inverse conductivity problem in the plane, Annals of Mathematics, 63 (006), pp [7] R. M. Brown and G. Uhlmann, Uniqueness in the inverse conductivity problem for nonsmooth conductivities in two dimensions, Communications in Partial Differential Equations, (997), pp [8] A.-P. Calderón, On an inverse boundary value problem, in Seminar on Numerical Analysis and its Applications to Continuum Physics (Rio de Janeiro, 980), Soc. Brasil. Mat., Rio de Janeiro, 980, pp [9] P. Caro, A. García, and J. M. Reyes, Stability of the Calderón problem for less regular conductivities, Journal of Differential Equations, 54 (03), pp [0] M. Cheney, D. Isaacson, and J. C. Newell, Electrical impedance tomography, SIAM Review, 4 (999), pp [] A. Clop, D. Faraco, and A. Ruiz, Stability of Calderón s inverse conductivity problem in the plane for discontinuous conductivities, Inverse Problems and Imaging, 4 (00), pp [] D. Faraco and K. Rogers, The Sobolev norm of characteristic functions with applications to the Calderón Inverse Problem, The Quarterly Journal of Mathematics, 64 (03), pp [3] B. Haberman and D. Tataru, Uniqueness in Calderón s problem with Lipschitz conductivities, Duke Mathematical Journal, 6 (03), pp [4] S. Hamilton, C. Herrera, J. L. Mueller, and A. Von Herrmann, A direct D- bar reconstruction algorithm for recovering a complex conductivity in -D, Arxiv preprint, (0). [5] M. Huhtanen and A. Perämäki, Numerical solution of the R-linear Beltrami equation, Mathematics of Computation, 8 (0), pp [6] D. Isaacson, J. Mueller, J. Newell, and S. Siltanen, Imaging cardiac activity by the D-bar method for electrical impedance tomography, Physiological Measurement, 7 (006), pp. S43 S50. [7] D. Isaacson, J. L. Mueller, J. C. Newell, and S. Siltanen, Reconstructions of chest phantoms by the D-bar method for electrical impedance tomography, IEEE Transactions on Medical Imaging, 3 (004), pp [8] K. Knudsen, A new direct method for reconstructing isotropic conductivities in the plane, Physiological Measurement, 4 (003), pp
29 NONLINEAR FOURIER ANALYSIS 9 [9] K. Knudsen, M. Lassas, J. Mueller, and S. Siltanen, D-bar method for electrical impedance tomography with discontinuous conductivities, SIAM Journal on Applied Mathematics, 67 (007), p [0] K. Knudsen, M. Lassas, J. Mueller, and S. Siltanen, Regularized D-bar method for the inverse conductivity problem, Inverse Problems and Imaging, 3 (009), pp [] K. Knudsen, J. Mueller, and S. Siltanen, Numerical solution method for the dbarequation in the plane, Journal of Computational Physics, 98 (004), pp [] K. Knudsen and A. Tamasan, Reconstruction of less regular conductivities in the plane, Communications in Partial Differential Equations, 9 (004), pp [3] M. Lassas, J. L. Mueller, S. Siltanen, and A. Stahel, The Novikov-Veselov Equation and the Inverse Scattering Method: II. Computation, Nonlinearity, 5 (0), pp [4], The Novikov-Veselov Equation and the Inverse Scattering Method, Part I: Analysis, Physica D, 4 (0), pp [5] J. Mueller and S. Siltanen, Direct reconstructions of conductivities from boundary measurements, SIAM Journal on Scientific Computing, 4 (003), pp [6], Linear and Nonlinear Inverse Problems with Practical Applications, SIAM, 0. [7] M. Music, P. Perry, and S. Siltanen, Exceptional circles of radial potentials, Inverse Problems, 9 (03). [8] A. I. Nachman, Global uniqueness for a two-dimensional inverse boundary value problem, Annals of Mathematics, 43 (996), pp [9] P. A. Perry, Miura Maps and Inverse Scattering for the Novikov-Veselov Equation, ArXiv e-prints 0.385, (0). [30] S. Siltanen, J. Mueller, and D. Isaacson, An implementation of the reconstruction algorithm of A. Nachman for the -D inverse conductivity problem, Inverse Problems, 6 (000), pp [3] G. Vainikko, Fast solvers of the Lippmann-Schwinger equation, in Direct and inverse problems of mathematical physics (Newark, DE, 997), vol. 5 of Int. Soc. Anal. Appl. Comput., Kluwer Acad. Publ., Dordrecht, 000, pp
arxiv: v2 [math.ap] 17 Jun 2014
NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN arxiv:30.705v [math.ap] 7 Jun 04 Contents. Introduction. The nonlinear
More informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationNumerical computation of complex geometrical optics solutions to the conductivity equation
Numerical computation of complex geometrical optics solutions to the conductivity equation Kari Astala a, Jennifer L Mueller b, Lassi Päivärinta a, Samuli Siltanen a a Department of Mathematics and Statistics,
More informationCalderón s inverse problem in 2D Electrical Impedance Tomography
Calderón s inverse problem in 2D Electrical Impedance Tomography Kari Astala (University of Helsinki) Joint work with: Matti Lassas, Lassi Päivärinta, Samuli Siltanen, Jennifer Mueller and Alan Perämäki
More informationDIRECT ELECTRICAL IMPEDANCE TOMOGRAPHY FOR NONSMOOTH CONDUCTIVITIES
DIRECT ELECTRICAL IMPEDANCE TOMOGRAPHY FOR NONSMOOTH CONDUCTIVITIES K ASTALA, J L MUELLER, A PERÄMÄKI, L PÄIVÄRINTA AND S SILTANEN Abstract. A new reconstruction algorithm is presented for eit in dimension
More informationContributors and Resources
for Introduction to Numerical Methods for d-bar Problems Jennifer Mueller, presented by Andreas Stahel Colorado State University Bern University of Applied Sciences, Switzerland Lexington, Kentucky, May
More informationAALBORG UNIVERSITY. A new direct method for reconstructing isotropic conductivities in the plane. Department of Mathematical Sciences.
AALBORG UNIVERSITY A new direct method for reconstructing isotropic conductivities in the plane by Kim Knudsen January 23 R-23-1 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationA CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA. has a unique solution. The Dirichlet-to-Neumann map
A CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA SUNGWHAN KIM AND ALEXANDRU TAMASAN ABSTRACT. We consider the problem of identifying a complex valued coefficient γ(x, ω) in the conductivity
More informationComputational inversion based on complex geometrical optics solutions
Computational inversion based on complex geometrical optics solutions Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi http://www.siltanen-research.net
More informationarxiv: v1 [math.ap] 25 Sep 2018
arxiv:1809.09272v1 [math.ap] 25 Sep 2018 NACHMAN S RECONSTRUCTION METHO FOR THE CALERÓN PROBLEM WITH ISCONTINUOUS CONUCTIVITIES GEORGE LYTLE, PETER PERRY, AN SAMULI SILTANEN Abstract. We show that Nachman
More informationDownloaded 11/20/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 33, No. 1, pp. 328 341 c 2011 Society for Industrial and Applied Mathematics A SIMPLE NUMERICAL METHOD FOR COMPLEX GEOMETRICAL OPTICS SOLUTIONS TO THE CONDUCTIVITY EQUATION KUI
More informationComputational inversion using complex geometric optics solutions and nonlinear Fourier transforms
Computational inversion using complex geometric optics solutions and nonlinear Fourier transforms Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationSome issues on Electrical Impedance Tomography with complex coefficient
Some issues on Electrical Impedance Tomography with complex coefficient Elisa Francini (Università di Firenze) in collaboration with Elena Beretta and Sergio Vessella E. Francini (Università di Firenze)
More informationMapping Properties of the Nonlinear Fourier Transform in Dimension Two
Communications in Partial Differential Equations, 3: 59 60, 007 Copyright Taylor & Francis Group, LLC ISSN 0360-530 print/53-433 online DOI: 0.080/036053005005304 Mapping Properties of the Nonlinear Fourier
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationStability and instability in inverse problems
Stability and instability in inverse problems Mikhail I. Isaev supervisor: Roman G. Novikov Centre de Mathématiques Appliquées, École Polytechnique November 27, 2013. Plan of the presentation The Gel fand
More informationInverse Electromagnetic Problems
Title: Name: Affil./Addr. 1: Inverse Electromagnetic Problems Gunther Uhlmann, Ting Zhou University of Washington and University of California Irvine /gunther@math.washington.edu Affil./Addr. 2: 340 Rowland
More informationReconstructing conductivities with boundary corrected D-bar method
J. Inverse Ill-Posed Probl., Ahead of Print DOI 10.1515/jip-2013-0042 de Gruyter 2014 Reconstructing conductivities with boundary corrected D-bar method Samuli Siltanen and Janne P. Tamminen Abstract.
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
More informationis measured (n is the unit outward normal) at the boundary, to define the Dirichlet-
SIAM J MATH ANAL Vol 45, No 5, pp 700 709 c 013 Society for Industrial and Applied Mathematics ON A CALDERÓN PROBLEM IN FREQUENCY DIFFERENTIAL ELECTRICAL IMPEDANCE TOMOGRAPHY SUNGWHAN KIM AND ALEXANDRU
More informationResistor Networks and Optimal Grids for Electrical Impedance Tomography with Partial Boundary Measurements
Resistor Networks and Optimal Grids for Electrical Impedance Tomography with Partial Boundary Measurements Alexander Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, Fernando Guevara Vasquez 4 1 Institute
More informationFast shape-reconstruction in electrical impedance tomography
Fast shape-reconstruction in electrical impedance tomography Bastian von Harrach bastian.harrach@uni-wuerzburg.de (joint work with Marcel Ullrich) Institut für Mathematik - IX, Universität Würzburg The
More informationRegularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints
Regularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints Bastian von Harrach (joint work with Mach Nguyet Minh) http://numerical.solutions Institute of
More informationMonotonicity arguments in electrical impedance tomography
Monotonicity arguments in electrical impedance tomography Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany NAM-Kolloquium, Georg-August-Universität
More informationCONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION
CONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 Abstract: For the one-dimensional
More informationInverse parameter identification problems
Inverse parameter identification problems Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany ICP - Institute for Computational
More informationElectrical Impedance Tomography: 3D Reconstructions using Scattering Transforms
Downloaded from orbit.dtu.dk on: Nov 11, 2018 Electrical Impedance Tomography: 3D Reconstructions using Scattering Transforms Delbary, Fabrice; Hansen, Per Christian; Knudsen, Kim Published in: Applicable
More informationarxiv: v2 [math.na] 2 Feb 2016 J. P. TAMMINEN, T. TARVAINEN AND S. SILTANEN
THE D-BAR METHOD FOR DIFFUSE OPTICAL TOMOGRAPHY: A COMPUTATIONAL STUDY arxiv:6.4636v2 [math.na] 2 Feb 26 J. P. TAMMINEN, T. TARVAINEN AND S. SILTANEN Abstract. The D-bar method at negative energy is numerically
More informationLiouville-type theorem for the Lamé system with singular coefficients
Liouville-type theorem for the Lamé system with singular coefficients Blair Davey Ching-Lung Lin Jenn-Nan Wang Abstract In this paper, we study a Liouville-type theorem for the Lamé system with rough coefficients
More information1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0
4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for
More informationWhat are inverse problems?
What are inverse problems? 6th Forum of Polish Mathematicians Warsaw University of Technology Joonas Ilmavirta University of Jyväskylä 8 September 2015 Joonas Ilmavirta (Jyväskylä) Inverse problems 8 September
More informationA note on the MUSIC algorithm for impedance tomography
A note on the MUSIC algorithm for impedance tomography Martin Hanke Institut für Mathematik, Johannes Gutenberg-Universität, 55099 Mainz, Germany E-mail: hanke@math.uni-mainz.de Abstract. We investigate
More informationarxiv: v1 [math.na] 19 Sep 2018
SPECTRAL APPROACHES TO D-BAR PROBLEMS ARISING IN ELECTRICAL IMPEDANCE TOMOGRAPHY arxiv:89.8294v [math.na] 9 Sep 28 CHRISTIAN KLEIN AND NIKOLA STOILOV Abstract. In this work we present spectral approaches
More informationNovel tomography techniques and parameter identification problems
Novel tomography techniques and parameter identification problems Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München Colloquium of the Institute of Biomathematics
More informationAALBORG UNIVERSITY. The Calderón problem with partial data for less smooth conductivities. Kim Knudsen. Department of Mathematical Sciences
AALBORG UNIVERSITY The Calderón problem with partial data for less smooth conductivities by Kim Knudsen R-2005-07 February 2005 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationarxiv: v1 [math.ap] 21 Dec 2014
CONSTRUCTION OF INVISIBLE CONUCTIVITY PERTURBATIONS FOR THE POINT ELECTROE MOEL IN ELECTRICAL IMPEANCE TOMOGRAPHY LUCAS CHESNEL, NUUTTI HYVÖNEN, AN STRATOS STABOULIS Abstract. We explain how to build invisible
More informationAllan Greenleaf, Matti Lassas, and Gunther Uhlmann
Mathematical Research Letters 10, 685 693 (2003) ON NONUNIQUENESS FOR CALDERÓN S INVERSE PROBLEM Allan Greenleaf, Matti Lassas, and Gunther Uhlmann Abstract. We construct anisotropic conductivities with
More informationSurvey of Inverse Problems For Hyperbolic PDEs
Survey of Inverse Problems For Hyperbolic PDEs Rakesh Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA Email: rakesh@math.udel.edu January 14, 2011 1 Problem Formulation
More informationStatistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications
Statistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications Aku Seppänen Inverse Problems Group Department of Applied Physics University
More informationCalderon s inverse conductivity problem in the plane
Calderon s inverse conductivity problem in the plane Kari Astala and Lassi Päivärinta Abstract We show that the Dirichlet to Neumann map for the equation σ u = 0 in a two dimensional domain uniquely determines
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationSensitivity Analysis of 3D Magnetic Induction Tomography (MIT)
Sensitivity Analysis of 3D Magnetic Induction Tomography (MIT) W R B Lionheart 1, M Soleimani 1, A J Peyton 2 1 Department of Mathematics UMIST, Manchester, UK, Email: bill.lionheart@umist.ac.uk, 2 Department
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationON THE ENDPOINT REGULARITY OF DISCRETE MAXIMAL OPERATORS
ON THE ENDPOINT REGULARITY OF DISCRETE MAXIMAL OPERATORS EMANUEL CARNEIRO AND KEVIN HUGHES Abstract. Given a discrete function f : Z d R we consider the maximal operator X Mf n = sup f n m, r 0 Nr m Ω
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationThe Mathematics of Invisibility: Cloaking Devices, Electromagnetic Wormholes, and Inverse Problems. Lectures 1-2
CSU Graduate Workshop on Inverse Problems July 30, 2007 The Mathematics of Invisibility: Cloaking Devices, Electromagnetic Wormholes, and Inverse Problems Lectures 1-2 Allan Greenleaf Department of Mathematics,
More informationSharp estimates for a class of hyperbolic pseudo-differential equations
Results in Math., 41 (2002), 361-368. Sharp estimates for a class of hyperbolic pseudo-differential equations Michael Ruzhansky Abstract In this paper we consider the Cauchy problem for a class of hyperbolic
More informationFollow links Class Use and other Permissions. For more information, send to:
COPYRIGHT NOTICE: Kari Astala, Tadeusz Iwaniec & Gaven Martin: Elliptic Partial Differential Equations and Quasiconformal Mappings in the Plane is published by Princeton University Press and copyrighted,
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationM.P. Ramirez T., IAENG Member, M. C. Robles G., R. A. Hernandez-Becerril
, July 6-8, 2011, London, U.K. On the Construction of Separable-Variables Conductivity Functions, and Their Application for Approaching Numerical Solutions of the Two-Dimensional Electrical Impedance Equation
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationTADAHIRO OH 0, 3 8 (T R), (1.5) The result in [2] is in fact stated for time-periodic functions: 0, 1 3 (T 2 ). (1.4)
PERIODIC L 4 -STRICHARTZ ESTIMATE FOR KDV TADAHIRO OH 1. Introduction In [], Bourgain proved global well-posedness of the periodic KdV in L T): u t + u xxx + uu x 0, x, t) T R. 1.1) The key ingredient
More informationON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics
More informationMapping problems and harmonic univalent mappings
Mapping problems and harmonic univalent mappings Antti Rasila Helsinki University of Technology antti.rasila@tkk.fi (Mainly based on P. Duren s book Harmonic mappings in the plane) Helsinki Analysis Seminar,
More informationarxiv:math/ v1 [math.ap] 6 Feb 2006
On Beltrami equations and Hölder regularity arxiv:math/0602095v [math.ap] 6 Feb 2006 Tonia Ricciardi Dipartimento di Matematica e Applicazioni R. Caccioppoli Università di Napoli Federico II Via Cintia,
More informationTHE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS
THE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS YASUHITO MIYAMOTO Abstract. We prove the hot spots conjecture of J. Rauch in the case that the domain Ω is a planar convex domain satisfying
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationl(y j ) = 0 for all y j (1)
Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that
More informationarxiv: v2 [math.na] 26 Mar 2018
APPROXIMATION OF FULL-BOUNDARY DATA FROM PARTIAL-BOUNDARY ELECTRODE MEASUREMENTS A. HAUPTMANN arxiv:1703.05550v2 [math.na] 26 Mar 2018 Abstract. Measurements on a subset of the boundary are common in electrical
More informationBALANCING GAUSSIAN VECTORS. 1. Introduction
BALANCING GAUSSIAN VECTORS KEVIN P. COSTELLO Abstract. Let x 1,... x n be independent normally distributed vectors on R d. We determine the distribution function of the minimum norm of the 2 n vectors
More informationElectric potentials with localized divergence properties
Electric potentials with localized divergence properties Bastian Gebauer bastian.gebauer@oeaw.ac.at Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences,
More informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
More informationDiscreteness of Transmission Eigenvalues via Upper Triangular Compact Operators
Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationfür Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig A short proof of the self-improving regularity of quasiregular mappings. by Xiao Zhong and Daniel Faraco Preprint no.: 106 2002 A
More informationRecent progress on the factorization method for electrical impedance tomography
Recent progress on the factorization method for electrical impedance tomography Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More informationAn inverse elliptic problem of medical optics with experimental data
An inverse elliptic problem of medical optics with experimental data Jianzhong Su, Michael V. Klibanov, Yueming Liu, Zhijin Lin Natee Pantong and Hanli Liu Abstract A numerical method for an inverse problem
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationSTABILITY OF CALDERÓN S INVERSE CONDUCTIVITY PROBLEM IN THE PLANE FOR DISCONTINUOUS CONDUCTIVITIES. Albert Clop. Daniel Faraco.
Inverse Problems and Imaging Volume 4, No., 00, 49 9 doi:0.3934/ipi.00.4.49 STABILITY OF CALDERÓN S INVERSE CONDUCTIVITY PROBLEM IN THE PLANE FOR DISCONTINUOUS CONDUCTIVITIES Albert Clop Department of
More informationNonlinear aspects of Calderón-Zygmund theory
Ancona, June 7 2011 Overture: The standard CZ theory Consider the model case u = f in R n Overture: The standard CZ theory Consider the model case u = f in R n Then f L q implies D 2 u L q 1 < q < with
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
More informationThe Linear Sampling Method and the MUSIC Algorithm
CODEN:LUTEDX/(TEAT-7089)/1-6/(2000) The Linear Sampling Method and the MUSIC Algorithm Margaret Cheney Department of Electroscience Electromagnetic Theory Lund Institute of Technology Sweden Margaret Cheney
More informationA Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms
Universal Journal of Integral Equations 4 (2016), 54-62 www.papersciences.com A Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms N. H. S. Haidar CRAMS: Center for
More informationDetecting stochastic inclusions in electrical impedance tomography
Detecting stochastic inclusions in electrical impedance tomography Bastian von Harrach harrach@math.uni-frankfurt.de (joint work with A. Barth, N. Hyvönen and L. Mustonen) Institute of Mathematics, Goethe
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationOn Riemann Boundary Value Problem in Hardy Classes with Variable Summability Exponent
Int. Journal of Math. Analysis, Vol. 6, 2012, no. 15, 743-751 On Riemann Boundary Value Problem in Hardy Classes with Variable Summability Exponent Muradov T.R. Institute of Mathematics and Mechanics of
More informationarxiv: v1 [math.ap] 11 Jun 2007
Inverse Conductivity Problem for a Parabolic Equation using a Carleman Estimate with One Observation arxiv:0706.1422v1 [math.ap 11 Jun 2007 November 15, 2018 Patricia Gaitan Laboratoire d Analyse, Topologie,
More informationOn the Optimal Insulation of Conductors 1
JOURNAL OF OPTIMIZATION THEORY AND APPLICATIONS: Vol. 100, No. 2, pp. 253-263, FEBRUARY 1999 On the Optimal Insulation of Conductors 1 S. J. COX, 2 B. KAWOHL, 3 AND P. X. UHLIG 4 Communicated by K. A.
More informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid ge Joint work with P. Caro (U. Helsinki) and D. Dos Santos Ferreira (Paris
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationOn the smoothness of the conjugacy between circle maps with a break
On the smoothness of the conjugacy between circle maps with a break Konstantin Khanin and Saša Kocić 2 Department of Mathematics, University of Toronto, Toronto, ON, Canada M5S 2E4 2 Department of Mathematics,
More informationON THE GAPS BETWEEN ZEROS OF TRIGONOMETRIC POLYNOMIALS
Real Analysis Exchange Vol. 8(), 00/003, pp. 447 454 Gady Kozma, Faculty of Mathematics and Computer Science, Weizmann Institute of Science, POB 6, Rehovot 76100, Israel. e-mail: gadyk@wisdom.weizmann.ac.il,
More informationThe X-ray transform for a non-abelian connection in two dimensions
The X-ray transform for a non-abelian connection in two dimensions David Finch Department of Mathematics Oregon State University Corvallis, OR, 97331, USA Gunther Uhlmann Department of Mathematics University
More informationApproximation of inverse boundary value problems by phase-field methods
Approximation of inverse boundary value problems by phase-field methods Luca RONDI Università degli Studi di Trieste Dipartimento di Matematica e Informatica MSRI Berkeley, 3 December 2010 Luca RONDI (Università
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationL p -boundedness of the Hilbert transform
L p -boundedness of the Hilbert transform Kunal Narayan Chaudhury Abstract The Hilbert transform is essentially the only singular operator in one dimension. This undoubtedly makes it one of the the most
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationD-BAR METHOD AND EXCEPTIONAL POINTS AT POSITIVE ENERGY: A COMPUTATIONAL STUDY
Proceedings of the Project Review, Geo-Mathematical Imaging Grou (Purdue University, West Lafayette IN), Vol. 1 (14). 389-415. D-BAR METHOD AND EXCEPTIONAL POINTS AT POSITIVE ENERGY: A COMPUTATIONAL STUDY
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
More informationCauchy Integrals and Cauchy-Szegö Projections. Elias M. Stein. Conference in honor of Michael Christ Madison, May 2016
1 Cauchy Integrals and Cauchy-Szegö Projections Elias M. Stein Conference in honor of Michael Christ Madison, May 2016 The goal: In C n, to study the L p properties of Cauchy integrals and the Cauchy-Szegö
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationare harmonic functions so by superposition
J. Rauch Applied Complex Analysis The Dirichlet Problem Abstract. We solve, by simple formula, the Dirichlet Problem in a half space with step function boundary data. Uniqueness is proved by complex variable
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More information