TY - JOUR A1 - Kleefeld, Andreas A1 - Reißel, Martin T1 - The Levenberg–Marquardt method applied to a parameter estimation problem arising from electrical resistivity tomography JF - Applied Mathematics and Computation Y1 - 2011 SN - 0096-3003 VL - 217 IS - 9 SP - 4490 EP - 4501 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Grajewski, Matthias A1 - Kleefeld, Andreas T1 - Detecting and approximating decision boundaries in low-dimensional spaces JF - Numerical Algorithms N2 - A method for detecting and approximating fault lines or surfaces, respectively, or decision curves in two and three dimensions with guaranteed accuracy is presented. Reformulated as a classification problem, our method starts from a set of scattered points along with the corresponding classification algorithm to construct a representation of a decision curve by points with prescribed maximal distance to the true decision curve. Hereby, our algorithm ensures that the representing point set covers the decision curve in its entire extent and features local refinement based on the geometric properties of the decision curve. We demonstrate applications of our method to problems related to the detection of faults, to multi-criteria decision aid and, in combination with Kirsch’s factorization method, to solving an inverse acoustic scattering problem. In all applications we considered in this work, our method requires significantly less pointwise classifications than previously employed algorithms. KW - MCDA KW - Inverse scattering problem KW - Fault approximation KW - Fault detection Y1 - 2023 SN - 1572-9265 N1 - Corresponding author: Matthias Grajewski VL - 93 IS - 4 PB - Springer Science+Business Media CY - Dordrecht ER - TY - JOUR A1 - Ayala, Rafael Ceja A1 - Harris, Isaac A1 - Kleefeld, Andreas A1 - Pallikarakis, Nikolaos T1 - Analysis of the transmission eigenvalue problem with two conductivity parameters JF - Applicable Analysis N2 - In this paper, we provide an analytical study of the transmission eigenvalue problem with two conductivity parameters. We will assume that the underlying physical model is given by the scattering of a plane wave for an isotropic scatterer. In previous studies, this eigenvalue problem was analyzed with one conductive boundary parameter whereas we will consider the case of two parameters. We prove the existence and discreteness of the transmission eigenvalues as well as study the dependence on the physical parameters. We are able to prove monotonicity of the first transmission eigenvalue with respect to the parameters and consider the limiting procedure as the second boundary parameter vanishes. Lastly, we provide extensive numerical experiments to validate the theoretical work. KW - Transmission Eigenvalues KW - Conductive Boundary Condition KW - Inverse Scattering Y1 - 2023 U6 - https://doi.org/10.1080/00036811.2023.2181167 SN - 0003-6811 PB - Taylor & Francis ER - TY - JOUR A1 - Pieronek, Lukas A1 - Kleefeld, Andreas T1 - On trajectories of complex-valued interior transmission eigenvalues JF - Inverse problems and imaging : IPI N2 - This paper investigates the interior transmission problem for homogeneous media via eigenvalue trajectories parameterized by the magnitude of the refractive index. In the case that the scatterer is the unit disk, we prove that there is a one-to-one correspondence between complex-valued interior transmission eigenvalue trajectories and Dirichlet eigenvalues of the Laplacian which turn out to be exactly the trajectorial limit points as the refractive index tends to infinity. For general simply-connected scatterers in two or three dimensions, a corresponding relation is still open, but further theoretical results and numerical studies indicate a similar connection. KW - Interior transmission problem KW - Eigenvalue trajectories KW - Complex-valued eigenvalues Y1 - 2024 U6 - https://doi.org/10.3934/ipi.2023041 SN - 1930-8337 (Print) SN - 1930-8345 (Online) VL - 18 IS - 2 SP - 480 EP - 516 PB - AIMS CY - Springfield, Mo ER - TY - JOUR A1 - Harris, Isaac A1 - Kleefeld, Andreas T1 - The inverse scattering problem for a conductive boundary condition and transmission eigenvalues JF - Applicable Analysis N2 - In this paper, we consider the inverse scattering problem associated with an inhomogeneous media with a conductive boundary. In particular, we are interested in two problems that arise from this inverse problem: the inverse conductivity problem and the corresponding interior transmission eigenvalue problem. The inverse conductivity problem is to recover the conductive boundary parameter from the measured scattering data. We prove that the measured scatted data uniquely determine the conductivity parameter as well as describe a direct algorithm to recover the conductivity. The interior transmission eigenvalue problem is an eigenvalue problem associated with the inverse scattering of such materials. We investigate the convergence of the eigenvalues as the conductivity parameter tends to zero as well as prove existence and discreteness for the case of an absorbing media. Lastly, several numerical and analytical results support the theory and we show that the inside–outside duality method can be used to reconstruct the interior conductive eigenvalues. KW - Transmission eigenvalues KW - Conductive boundary condition KW - Inverse scattering Y1 - 2018 U6 - https://doi.org/10.1080/00036811.2018.1504028 SN - 1563-504X VL - 99 IS - 3 SP - 508 EP - 529 PB - Taylor & Francis CY - London ER - TY - JOUR A1 - Kleefeld, Andreas A1 - Pieronek, J. T1 - Elastic transmission eigenvalues and their computation via the method of fundamental solutions JF - Applicable Analysis N2 - A stabilized version of the fundamental solution method to catch ill-conditioning effects is investigated with focus on the computation of complex-valued elastic interior transmission eigenvalues in two dimensions for homogeneous and isotropic media. Its algorithm can be implemented very shortly and adopts to many similar partial differential equation-based eigenproblems as long as the underlying fundamental solution function can be easily generated. We develop a corroborative approximation analysis which also implicates new basic results for transmission eigenfunctions and present some numerical examples which together prove successful feasibility of our eigenvalue recovery approach. KW - elastic scattering KW - method of fundamental solutions KW - Interior transmission eigenvalues Y1 - 2020 U6 - https://doi.org/10.1080/00036811.2020.1721473 SN - 1563-504X VL - 100 IS - 16 SP - 3445 EP - 3462 PB - Taylore & Francis CY - London ER - TY - JOUR A1 - Breuß, Michael A1 - Kleefeld, Andreas T1 - Implicit monotone difference methods for scalar conservation laws with source terms JF - Acta Mathematica Vietnamica N2 - In this article, a concept of implicit methods for scalar conservation laws in one or more spatial dimensions allowing also for source terms of various types is presented. This material is a significant extension of previous work of the first author (Breuß SIAM J. Numer. Anal. 43(3), 970–986 2005). Implicit notions are developed that are centered around a monotonicity criterion. We demonstrate a connection between a numerical scheme and a discrete entropy inequality, which is based on a classical approach by Crandall and Majda. Additionally, three implicit methods are investigated using the developed notions. Next, we conduct a convergence proof which is not based on a classical compactness argument. Finally, the theoretical results are confirmed by various numerical tests. KW - Entropy solution KW - Source term KW - Monotone methods KW - Implicit methods KW - Finite difference methods KW - Conservation laws Y1 - 2020 U6 - https://doi.org/10.1007/s40306-019-00354-1 SN - 2315-4144 N1 - Corresponding author: Andreas Kleefeld VL - 45 SP - 709 EP - 738 PB - Springer Singapore CY - Singapore ER - TY - JOUR A1 - Asante-Asamani, E.O. A1 - Kleefeld, Andreas A1 - Wade, B.A. T1 - A second-order exponential time differencing scheme for non-linear reaction-diffusion systems with dimensional splitting JF - Journal of Computational Physics N2 - A second-order L-stable exponential time-differencing (ETD) method is developed by combining an ETD scheme with approximating the matrix exponentials by rational functions having real distinct poles (RDP), together with a dimensional splitting integrating factor technique. A variety of non-linear reaction-diffusion equations in two and three dimensions with either Dirichlet, Neumann, or periodic boundary conditions are solved with this scheme and shown to outperform a variety of other second-order implicit-explicit schemes. An additional performance boost is gained through further use of basic parallelization techniques. KW - Exponential time differencing KW - Real distinct pole KW - Dimensional splitting KW - Reaction-diffusion systems KW - Matrix exponential Y1 - 2020 U6 - https://doi.org/10.1016/j.jcp.2020.109490 SN - 0021-9991 N1 - Corresponding author: Andreas Kleefeld VL - 415 PB - Elsevier CY - Amsterdam ER - TY - CHAP A1 - Kleefeld, Andreas ED - Constanda, Christian T1 - Numerical calculation of interior transmission eigenvalues with mixed boundary conditions T2 - Computational and Analytic Methods in Science and Engineering N2 - Interior transmission eigenvalue problems for the Helmholtz equation play an important role in inverse wave scattering. Some distribution properties of those eigenvalues in the complex plane are reviewed. Further, a new scattering model for the interior transmission eigenvalue problem with mixed boundary conditions is described and an efficient algorithm for computing the interior transmission eigenvalues is proposed. Finally, extensive numerical results for a variety of two-dimensional scatterers are presented to show the validity of the proposed scheme. Y1 - 2020 SN - 978-3-030-48185-8 (Hardcover) U6 - https://doi.org/10.1007/978-3-030-48186-5_9 SP - 173 EP - 195 PB - Birkhäuser CY - Cham ER - TY - JOUR A1 - Martín-Vaquero, J. A1 - Kleefeld, Andreas T1 - Solving nonlinear parabolic PDEs in several dimensions: Parallelized ESERK codes JF - Journal of Computational Physics N2 - There is a very large number of very important situations which can be modeled with nonlinear parabolic partial differential equations (PDEs) in several dimensions. In general, these PDEs can be solved by discretizing in the spatial variables and transforming them into huge systems of ordinary differential equations (ODEs), which are very stiff. Therefore, standard explicit methods require a large number of iterations to solve stiff problems. But implicit schemes are computationally very expensive when solving huge systems of nonlinear ODEs. Several families of Extrapolated Stabilized Explicit Runge-Kutta schemes (ESERK) with different order of accuracy (3 to 6) are derived and analyzed in this work. They are explicit methods, with stability regions extended, along the negative real semi-axis, quadratically with respect to the number of stages s, hence they can be considered to solve stiff problems much faster than traditional explicit schemes. Additionally, they allow the adaptation of the step length easily with a very small cost. Two new families of ESERK schemes (ESERK3 and ESERK6) are derived, and analyzed, in this work. Each family has more than 50 new schemes, with up to 84.000 stages in the case of ESERK6. For the first time, we also parallelized all these new variable step length and variable number of stages algorithms (ESERK3, ESERK4, ESERK5, and ESERK6). These parallelized strategies allow to decrease times significantly, as it is discussed and also shown numerically in two problems. Thus, the new codes provide very good results compared to other well-known ODE solvers. Finally, a new strategy is proposed to increase the efficiency of these schemes, and it is discussed the idea of combining ESERK families in one code, because typically, stiff problems have different zones and according to them and the requested tolerance the optimum order of convergence is different. KW - Multi-dimensional partial differential equations KW - Higher-order codes KW - Nonlinear PDEs Y1 - 2020 U6 - https://doi.org/10.1016/j.jcp.2020.109771 SN - 0021-9991 IS - 423 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Kleefeld, Andreas T1 - The hot spots conjecture can be false: some numerical examples JF - Advances in Computational Mathematics N2 - The hot spots conjecture is only known to be true for special geometries. This paper shows numerically that the hot spots conjecture can fail to be true for easy to construct bounded domains with one hole. The underlying eigenvalue problem for the Laplace equation with Neumann boundary condition is solved with boundary integral equations yielding a non-linear eigenvalue problem. Its discretization via the boundary element collocation method in combination with the algorithm by Beyn yields highly accurate results both for the first non-zero eigenvalue and its corresponding eigenfunction which is due to superconvergence. Additionally, it can be shown numerically that the ratio between the maximal/minimal value inside the domain and its maximal/minimal value on the boundary can be larger than 1 + 10− 3. Finally, numerical examples for easy to construct domains with up to five holes are provided which fail the hot spots conjecture as well. KW - Numerics KW - Boundary integral equations KW - Potential theory KW - Helmholtz equation KW - Interior Neumann eigenvalues Y1 - 2021 U6 - https://doi.org/10.1007/s10444-021-09911-5 SN - 1019-7168 VL - 47 PB - Springer CY - Dordrecht ER - TY - GEN A1 - Burgeth, Bernhard A1 - Kleefeld, Andreas A1 - Naegel, Benoît A1 - Perret, Benjamin T1 - Editorial — Special Issue: ISMM 2019 T2 - Mathematical Morphology - Theory and Applications N2 - This editorial presents the Special Issue dedicated to the conference ISMM 2019 and summarizes the articles published in this Special Issue. Y1 - 2020 U6 - https://doi.org/10.1515/mathm-2020-0200 SN - 2353-3390 VL - 4 IS - 1 SP - 159 EP - 161 PB - De Gruyter CY - Warschau ER - TY - JOUR A1 - Kleefeld, Andreas A1 - Zimmermann, M. ED - Constanda, Christian ED - Bodmann, Bardo E.J. ED - Harris, Paul J. T1 - Computing Elastic Interior Transmission Eigenvalues JF - Integral Methods in Science and Engineering N2 - An alternative method is presented to numerically compute interior elastic transmission eigenvalues for various domains in two dimensions. This is achieved by discretizing the resulting system of boundary integral equations in combination with a nonlinear eigenvalue solver. Numerical results are given to show that this new approach can provide better results than the finite element method when dealing with general domains. Y1 - 2022 SN - 978-3-031-07171-3 U6 - https://doi.org/10.1007/978-3-031-07171-3_10 N1 - Corresponding author: Andreas Kleefeld SP - 139 EP - 155 PB - Birkhäuser CY - Cham ER - TY - CHAP A1 - Burgeth, Bernhard A1 - Kleefeld, Andreas A1 - Zhang, Eugene A1 - Zhang, Yue ED - Baudrier, Étienne ED - Naegel, Benoît ED - Krähenbühl, Adrien ED - Tajine, Mohamed T1 - Towards Topological Analysis of Non-symmetric Tensor Fields via Complexification T2 - Discrete Geometry and Mathematical Morphology N2 - Fields of asymmetric tensors play an important role in many applications such as medical imaging (diffusion tensor magnetic resonance imaging), physics, and civil engineering (for example Cauchy-Green-deformation tensor, strain tensor with local rotations, etc.). However, such asymmetric tensors are usually symmetrized and then further processed. Using this procedure results in a loss of information. A new method for the processing of asymmetric tensor fields is proposed restricting our attention to tensors of second-order given by a 2x2 array or matrix with real entries. This is achieved by a transformation resulting in Hermitian matrices that have an eigendecomposition similar to symmetric matrices. With this new idea numerical results for real-world data arising from a deformation of an object by external forces are given. It is shown that the asymmetric part indeed contains valuable information. Y1 - 2022 SN - 978-3-031-19897-7 U6 - https://doi.org/10.1007/978-3-031-19897-7_5 N1 - Second International Joint Conference, DGMM 2022, Strasbourg, France, October 24–27, 2022 N1 - Corresponding author: Andreas Kleefeld SP - 48 EP - 59 PB - Springer CY - Cham ER - TY - JOUR A1 - Harris, Isaac A1 - Kleefeld, Andreas T1 - Analysis and computation of the transmission eigenvalues with a conductive boundary condition JF - Applicable Analysis N2 - We provide a new analytical and computational study of the transmission eigenvalues with a conductive boundary condition. These eigenvalues are derived from the scalar inverse scattering problem for an inhomogeneous material with a conductive boundary condition. The goal is to study how these eigenvalues depend on the material parameters in order to estimate the refractive index. The analytical questions we study are: deriving Faber–Krahn type lower bounds, the discreteness and limiting behavior of the transmission eigenvalues as the conductivity tends to infinity for a sign changing contrast. We also provide a numerical study of a new boundary integral equation for computing the eigenvalues. Lastly, using the limiting behavior we will numerically estimate the refractive index from the eigenvalues provided the conductivity is sufficiently large but unknown. KW - Boundary integral equations KW - Inverse spectral problem KW - Conductive boundary condition KW - Transmission eigenvalues Y1 - 2020 U6 - https://doi.org/10.1080/00036811.2020.1789598 SN - 1563-504X VL - 101 IS - 6 SP - 1880 EP - 1895 PB - Taylor & Francis CY - London ER - TY - JOUR A1 - Clausnitzer, Julian A1 - Kleefeld, Andreas T1 - A spectral Galerkin exponential Euler time-stepping scheme for parabolic SPDEs on two-dimensional domains with a C² boundary JF - Discrete and Continuous Dynamical Systems - Series B N2 - We consider the numerical approximation of second-order semi-linear parabolic stochastic partial differential equations interpreted in the mild sense which we solve on general two-dimensional domains with a C² boundary with homogeneous Dirichlet boundary conditions. The equations are driven by Gaussian additive noise, and several Lipschitz-like conditions are imposed on the nonlinear function. We discretize in space with a spectral Galerkin method and in time using an explicit Euler-like scheme. For irregular shapes, the necessary Dirichlet eigenvalues and eigenfunctions are obtained from a boundary integral equation method. This yields a nonlinear eigenvalue problem, which is discretized using a boundary element collocation method and is solved with the Beyn contour integral algorithm. We present an error analysis as well as numerical results on an exemplary asymmetric shape, and point out limitations of the approach. KW - Nonlinear eigenvalue problems KW - Boundary integral equations, KW - Exponential Euler scheme, KW - Parabolic SPDEs Y1 - 2024 U6 - https://doi.org/10.3934/dcdsb.2023148 SN - 1531-3492 SN - 1553-524X (eISSN) VL - 29 IS - 4 SP - 1624 EP - 1651 PB - AIMS CY - Springfield ER - TY - JOUR A1 - Ayala, Rafael Ceja A1 - Harris, Isaac A1 - Kleefeld, Andreas T1 - Direct sampling method via Landweber iteration for an absorbing scatterer with a conductive boundary JF - Inverse Problems and Imaging N2 - In this paper, we consider the inverse shape problem of recovering isotropic scatterers with a conductive boundary condition. Here, we assume that the measured far-field data is known at a fixed wave number. Motivated by recent work, we study a new direct sampling indicator based on the Landweber iteration and the factorization method. Therefore, we prove the connection between these reconstruction methods. The method studied here falls under the category of qualitative reconstruction methods where an imaging function is used to recover the absorbing scatterer. We prove stability of our new imaging function as well as derive a discrepancy principle for recovering the regularization parameter. The theoretical results are verified with numerical examples to show how the reconstruction performs by the new Landweber direct sampling method. Y1 - 2024 U6 - https://doi.org/10.3934/ipi.2023051 SN - 1930-8337 SN - 1930-8345 (eISSN) VL - 18 IS - 3 SP - 708 EP - 729 PB - AIMS CY - Springfield ER - TY - CHAP A1 - Kahra, Marvin A1 - Breuß, Michael A1 - Kleefeld, Andreas A1 - Welk, Martin ED - Brunetti, Sara ED - Frosini, Andrea ED - Rinaldi, Simone T1 - An Approach to Colour Morphological Supremum Formation Using the LogSumExp Approximation T2 - Discrete Geometry and Mathematical Morphology N2 - Mathematical morphology is a part of image processing that has proven to be fruitful for numerous applications. Two main operations in mathematical morphology are dilation and erosion. These are based on the construction of a supremum or infimum with respect to an order over the tonal range in a certain section of the image. The tonal ordering can easily be realised in grey-scale morphology, and some morphological methods have been proposed for colour morphology. However, all of these have certain limitations. In this paper we present a novel approach to colour morphology extending upon previous work in the field based on the Loewner order. We propose to consider an approximation of the supremum by means of a log-sum exponentiation introduced by Maslov. We apply this to the embedding of an RGB image in a field of symmetric 2x2 matrices. In this way we obtain nearly isotropic matrices representing colours and the structural advantage of transitivity. In numerical experiments we highlight some remarkable properties of the proposed approach. Y1 - 2024 SN - 978-3-031-57793-2 U6 - https://doi.org/10.1007/978-3-031-57793-2_25 N1 - Third International Joint Conference, DGMM 2024, Florence, Italy, April 15–18, 2024 SP - 325 EP - 337 PB - Springer CY - Cham ER - TY - CHAP A1 - Pieronek, Lukas A1 - Kleefeld, Andreas ED - Constanda, Christian ED - Harris, Paul T1 - The Method of Fundamental Solutions for Computing Interior Transmission Eigenvalues of Inhomogeneous Media T2 - Integral Methods in Science and Engineering: Analytic Treatment and Numerical Approximations N2 - The method of fundamental solutions is applied to the approximate computation of interior transmission eigenvalues for a special class of inhomogeneous media in two dimensions. We give a short approximation analysis accompanied with numerical results that clearly prove practical convenience of our alternative approach. Y1 - 2019 SN - 978-3-030-16077-7 U6 - https://doi.org/10.1007/978-3-030-16077-7_28 SP - 353 EP - 365 PB - Birkhäuser CY - Cham ER - TY - CHAP A1 - Abele, Daniel A1 - Kleefeld, Andreas ED - Constanda, Christian T1 - New Numerical Results for the Optimization of Neumann Eigenvalues T2 - Computational and Analytic Methods in Science and Engineering N2 - We present new numerical results for shape optimization problems of interior Neumann eigenvalues. This field is not well understood from a theoretical standpoint. The existence of shape maximizers is not proven beyond the first two eigenvalues, so we study the problem numerically. We describe a method to compute the eigenvalues for a given shape that combines the boundary element method with an algorithm for nonlinear eigenvalues. As numerical optimization requires many such evaluations, we put a focus on the efficiency of the method and the implemented routine. The method is well suited for parallelization. Using the resulting fast routines and a specialized parametrization of the shapes, we found improved maxima for several eigenvalues. Y1 - 2020 SN - 978-3-030-48185-8 (Print) SN - 978-3-030-48186-5 (Online) U6 - https://doi.org/10.1007/978-3-030-48186-5_1 SP - 1 EP - 20 PB - Birkhäuser CY - Cham ER -