@misc{WeiserFreytagErdmannetal.2018, author = {Weiser, Martin and Freytag, Yvonne and Erdmann, Bodo and Hubig, Michael and Mall, Gita}, title = {Optimal Design of Experiments for Estimating the Time of Death in Forensic Medicine}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-67247}, year = {2018}, abstract = {Estimation of time of death based on a single measurement of body core temperature is a standard procedure in forensic medicine. Mechanistic models using simulation of heat transport promise higher accuracy than established phenomenological models in particular in nonstandard situations, but involve many not exactly known physical parameters. Identifying both time of death and physical parameters from multiple temperature measurements is one possibility to reduce the uncertainty significantly. In this paper, we consider the inverse problem in a Bayesian setting and perform both local and sampling-based uncertainty quantification, where proper orthogonal decomposition is used as model reduction for fast solution of the forward model. Based on the local uncertainty quantification, optimal design of experiments is performed in order to minimize the uncertainty in the time of death estimate for a given number of measurements. For reasons of practicability, temperature acquisition points are selected from a set of candidates in different spatial and temporal locations. Applied to a real corpse model, a significant accuracy improvement is obtained already with a small number of measurements.}, language = {en} } @misc{Binkowski2017, type = {Master Thesis}, author = {Binkowski, Felix}, title = {On the convergence behavior of spectral deferred correction methods for convection-diffusion equations}, pages = {65}, year = {2017}, abstract = {Spectral deferred correction (SDC) Methoden, vorgestellt von Dutt, Greengard und Rokhlin in [1], sind iterative Verfahren zur numerischen L{\"o}sung von Anfangswertproblemen f{\"u}r gew{\"o}hnliche Differentialgleichungen. Wenn diese Methoden konvergieren, dann wird unter Verwendung von Zeitschrittverfahren niedriger Ordnung eine Kollokationsl{\"o}sung berechnet. Die L{\"o}sung von steifen Anfangswertproblemen ist eine relevante Problemstellung in der numerischen Mathematik. SDC-Methoden, speziell f{\"u}r steife Probleme, werden von Martin Weiser in [2] konstruiert. Die Theorie und die Experimente beziehen sich dabei auf Probleme, die aus r{\"a}umlich semidiskretisierten Reaktions-Diffusions-Gleichungen entstehen. In dieser Arbeit werden die Ans{\"a}tze aus [2] auf Konvektions-Diffusions-Gleichungen angewendet und das resultierende Konvergenzverhalten von SDC-Methoden untersucht. Basierend auf einem einfachen Konvektions-Diffusions-Operator, dessen spektrale Eigenschaften umfassend studiert werden, wird ein Schema zur Verbesserung dieses Verhaltens entwickelt. Numerische Experimente zeigen, dass eine Verbesserung der in [1] eingef{\"u}hrten SDC-Methoden m{\"o}glich ist. Die Untersuchungen ergeben weiterhin, dass das auch f{\"u}r komplexere Konvektions-Diffusions-Probleme gilt. [1] Alok Dutt, Leslie Greengard, und Vladimir Rokhlin. "Spectral deferred correction methods for ordinary differential equations." In: BIT 40.2 (2000), pp. 241-266. [2] Martin Weiser. "Faster SDC convergence on non-equidistant grids by DIRK sweeps." In: BIT 55.4 (2015), pp. 1219-1241.}, language = {en} } @article{MoldenhauerWeiserZachow2017, author = {Moldenhauer, Marian and Weiser, Martin and Zachow, Stefan}, title = {Adaptive Algorithms for Optimal Hip Implant Positioning}, volume = {17}, journal = {PAMM}, number = {1}, doi = {10.1002/pamm.201710071}, pages = {203 -- 204}, year = {2017}, abstract = {In an aging society where the number of joint replacements rises, it is important to also increase the longevity of implants. In particular hip implants have a lifetime of at most 15 years. This derives primarily from pain due to implant migration, wear, inflammation, and dislocation, which is affected by the positioning of the implant during the surgery. Current joint replacement practice uses 2D software tools and relies on the experience of surgeons. Especially the 2D tools fail to take the patients' natural range of motion as well as stress distribution in the 3D joint induced by different daily motions into account. Optimizing the hip joint implant position for all possible parametrized motions under the constraint of a contact problem is prohibitively expensive as there are too many motions and every position change demands a recalculation of the contact problem. For the reduction of the computational effort, we use adaptive refinement on the parameter domain coupled with the interpolation method of Kriging. A coarse initial grid is to be locally refined using goal-oriented error estimation, reducing locally high variances. This approach will be combined with multi-grid optimization such that numerical errors are reduced.}, language = {en} } @misc{GoetschelMinion2018, author = {G{\"o}tschel, Sebastian and Minion, Michael L.}, title = {Parallel-in-Time for Parabolic Optimal Control Problems Using PFASST}, journal = {Domain Decomposition Methods in Science and Engineering XXIV}, publisher = {Springer}, address = {Cham}, doi = {10.1007/978-3-319-93873-8_34}, pages = {363 -- 371}, year = {2018}, abstract = {In gradient-based methods for parabolic optimal control problems, it is necessary to solve both the state equation and a backward-in-time adjoint equation in each iteration of the optimization method. In order to facilitate fully parallel gradient-type and nonlinear conjugate gradient methods for the solution of such optimal control problems, we discuss the application of the parallel-in-time method PFASST to adjoint gradient computation. In addition to enabling time parallelism, PFASST provides high flexibility for handling nonlinear equations, as well as potential extra computational savings from reusing previous solutions in the optimization loop. The approach is demonstrated here for a model reaction-diffusion optimal control problem.}, language = {en} } @article{WeiserErdmannSchenkletal.2018, author = {Weiser, Martin and Erdmann, Bodo and Schenkl, Sebastian and Muggenthaler, Holger and Hubig, Michael and Mall, Gita and Zachow, Stefan}, title = {Uncertainty in Temperature-Based Determination of Time of Death}, volume = {54}, journal = {Heat and Mass Transfer}, number = {9}, publisher = {Springer}, doi = {10.1007/s00231-018-2324-4}, pages = {2815 -- 2826}, year = {2018}, abstract = {Temperature-based estimation of time of death (ToD) can be per- formed either with the help of simple phenomenological models of corpse cooling or with detailed mechanistic (thermodynamic) heat transfer mod- els. The latter are much more complex, but allow a higher accuracy of ToD estimation as in principle all relevant cooling mechanisms can be taken into account. The potentially higher accuracy depends on the accuracy of tissue and environmental parameters as well as on the geometric resolution. We in- vestigate the impact of parameter variations and geometry representation on the estimated ToD based on a highly detailed 3D corpse model, that has been segmented and geometrically reconstructed from a computed to- mography (CT) data set, differentiating various organs and tissue types.}, language = {en} } @misc{GoetschelWeiser2019, author = {G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Lossy Compression for Large Scale PDE Problems}, issn = {1438-0064}, doi = {10.1101/506378}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-73817}, year = {2019}, abstract = {Solvers for partial differential equations (PDE) are one of the cornerstones of computational science. For large problems, they involve huge amounts of data that needs to be stored and transmitted on all levels of the memory hierarchy. Often, bandwidth is the limiting factor due to relatively small arithmetic intensity, and increasingly so due to the growing disparity between computing power and bandwidth. Consequently, data compression techniques have been investigated and tailored towards the specific requirements of PDE solvers during the last decades. This paper surveys data compression challenges and corresponding solution approaches for PDE problems, covering all levels of the memory hierarchy from mass storage up to main memory. Exemplarily, we illustrate concepts at particular methods, and give references to alternatives.}, language = {en} } @misc{GoetschelSchielaWeiser2019, author = {G{\"o}tschel, Sebastian and Schiela, Anton and Weiser, Martin}, title = {Kaskade 7 -- a Flexible Finite Element Toolbox}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-74616}, year = {2019}, abstract = {Kaskade 7 is a finite element toolbox for the solution of stationary or transient systems of partial differential equations, aimed at supporting application-oriented research in numerical analysis and scientific computing. The library is written in C++ and is based on the Dune interface. The code is independent of spatial dimension and works with different grid managers. An important feature is the mix-and-match approach to discretizing systems of PDEs with different ansatz and test spaces for all variables. We describe the mathematical concepts behind the library as well as its structure, illustrating its use at several examples on the way.}, language = {en} } @misc{Schiela2008, author = {Schiela, Anton}, title = {State constrained optimal control problems with states of low regularity}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-509}, number = {08-24}, year = {2008}, abstract = {We consider first order optimality conditions for state constrained optimal control problems. In particular we study the case where the state equation has not enough regularity to admit existence of a Slater point in function space. We overcome this difficulty by a special transformation. Under a density condition we show existence of Lagrange multipliers, which have a representation via measures and additional regularity properties.}, language = {en} } @misc{ErdmannLangSeebass1997, author = {Erdmann, Bodo and Lang, Jens and Seebass, Martin}, title = {Adaptive Solutions of Nonlinear Parabolic Equations with Application to Hyperthermia Treatments}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3130}, number = {SC-97-44}, year = {1997}, abstract = {We present a self-adaptive finite element method to solve nonlinear evolution problems in 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. The proposed method is applied to hyperthermia treatments to demonstrate its potential for the solving of complicated problems.}, language = {en} } @misc{Weiser2003, author = {Weiser, Martin}, title = {Interior Point Methods in Function Space}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7578}, number = {03-35}, year = {2003}, abstract = {A primal-dual interior point method for optimal control problems is considered. The algorithm is directly applied to the infinite dimensional problem. Existence and convergence of the central path are analyzed, and linear convergence of a short step pathfollowing method is established.}, language = {en} } @misc{WeiserSchielaDeuflhard2003, author = {Weiser, Martin and Schiela, Anton and Deuflhard, Peter}, title = {Asymptotic Mesh Independence of Newton's Method Revisited}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7352}, number = {03-13}, year = {2003}, abstract = {The paper presents a new affine invariant theory on asymptotic mesh independence of Newton's method in nonlinear PDEs. Compared to earlier attempts, the new approach is both much simpler and more natural from the algorithmic point of view. The theory is exemplified at collocation methods for ODE boundary value problems and at finite element methods for elliptic PDE problems.}, language = {en} } @misc{Deuflhard2002, author = {Deuflhard, Peter}, title = {A Comparison of Related Concepts in Computational Chemistry and Mathematics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7138}, number = {02-46}, year = {2002}, abstract = {This article studies the relation of the two scientific languages Chemistry and Mathematics via three selected comparisons: (a) QSSA versus dynamic ILDM in reaction kinetics, (b) lumping versus discrete Galerkin methods in polymer chemistry, and (c) geometrical conformations versus metastable conformations in drug design. The common clear message from these comparisons is that chemical intuition may pave the way for mathematical concepts just as chemical concepts may gain from mathematical precisioning. Along this line, significant improvements in chemical research and engineering have already been possible -- and can be further expected in the future from the dialogue between the two scientific languages.}, language = {en} } @misc{Deuflhard2002, author = {Deuflhard, Peter}, title = {From Molecular Dynamics to Conformational Dynamics in Drug Design}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6878}, number = {02-20}, year = {2002}, abstract = {Computational drug design studies molecular recognition in the {\em virtual lab}. The arising Hamiltonian dynamics is known to be chaotic and ill-conditioned already after picoseconds, whereas times are \$msec\$ up to \$min\$. Classical molecular dynamics with long term trajectory computation gives, at best, information about time and statistical ensemble averages. The present paper surveys a recent new modeling approach called {\em conformational dynamics}, which is due to the author and Ch. Sch{\"u}tte. This approach achieves information about the dy time scales by telescoping a short term deterministic model with a statistical model. Examples of small biomolecules are included.}, language = {en} } @misc{DeuflhardHochmuth2002, author = {Deuflhard, Peter and Hochmuth, Reinhard}, title = {Multiscale Analysis of Thermoregulation in the Human Microvascular System}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6985}, number = {02-31}, year = {2002}, abstract = {The bio-heat transfer equation is a macroscopic model for describing the heat transfer in microvascular tissue. So far the deduction of the Helmholtz term in the bio-heat transfer equation is not co role. In view of a future numerical application of this new mathematical model to treatment planning in hyperthermia we derive asymptotic estimates for first and second order correctors.}, language = {en} } @misc{DeuflhardNowakWeiser2002, author = {Deuflhard, Peter and Nowak, Ulrich and Weiser, Martin}, title = {Affine Invariant Adaptive Newton Codes for Discretized PDEs}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7005}, number = {02-33}, year = {2002}, abstract = {The paper deals with three different Newton algorithms that have recently been worked out in the general frame of affine invariance. Of particular interest is their performance in the numerical solution of discretized boundary value problems (BVPs) for nonlinear partial differential equations (PDEs). Exact Newton methods, where the arising linear systems are solved by direct elimination, and inexact Newton methods, where an inner iteration is used instead, are synoptically presented, both in affine invariant convergence theory and in numerical experiments. The three types of algorithms are: (a) affine covariant (formerly just called affine invariant) Newton algorithms, oriented toward the iterative errors, (b) affine contravariant Newton algorithms, based on iterative residual norms, and (c) affine conjugate Newton algorithms for convex optimization problems and discrete nonlinear elliptic PDEs.}, language = {en} } @misc{PowellWeiser1999, author = {Powell, Gary and Weiser, Martin}, title = {Container Adaptors}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4308}, number = {SC-99-41}, year = {1999}, abstract = {The C++ standard template library has many useful containers for data. The standard library includes two adpators, queue, and stack. The authors have extended this model along the lines of relational database semantics. Sometimes the analogy is striking, and we will point it out occasionally. An adaptor allows the standard algorithms to be used on a subset or modification of the data without having to copy the data elements into a new container. The authors provide many useful adaptors which can be used together to produce interesting views of data in a container.}, language = {en} } @misc{ErdmannLangSeebass1997, author = {Erdmann, Bodo and Lang, Jens and Seebass, Martin}, title = {Optimization of Temperature Distributions for Regional Hyperthermia Based on a Nonlinear Heat Transfer Model}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3283}, number = {SC-97-59}, year = {1997}, abstract = {We describe an optimization process specially designed for regional hyperthermia of deap seated tumors in order to achieve desired steady--state temperature distributions. A nonlinear three--dimensional heat--transfer model based on temperature--dependent blood perfusion is applied to predict the temperature. Optimal heating is obtained by minimizing an integral object function which measures the distance between desired and model predicted temperatures. Sequential minima are calculated from successively improved constant--rate perfusion models employing a damped Newton method in an inner iteration. Numerical results for a Sigma 60 applicator are presented. This work has been supported by Deutsche Forschungsgemeinschaft (DFG) within the Sonderforschungsbereich 273 \glqq Hyperthermie: Methodik und Klinik \grqq .}, language = {en} } @misc{SchielaWeiser2005, author = {Schiela, Anton and Weiser, Martin}, title = {Superlinear Convergence of the Control Reduced Interior Point Method for PDE Constrained Optimization}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8490}, number = {05-15}, year = {2005}, abstract = {A thorough convergence analysis of the Control Reduced Interior Point Method in function space is performed. This recently proposed method is a primal interior point pathfollowing scheme with the special feature, that the control variable is eliminated from the optimality system. Apart from global linear convergence we show, that this method converges locally almost quadratically, if the optimal solution satisfies a function space analogue to a non-degeneracy condition. In numerical experiments we observe, that a prototype implementation of our method behaves in compliance with our theoretical results.}, language = {en} } @misc{KoberErdmannHellmichetal.2004, author = {Kober, Cornelia and Erdmann, Bodo and Hellmich, Christian and Sader, Robert and Zeilhofer, Hans-Florian}, title = {Anisotropic Simulation of the Human Mandible}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7873}, number = {04-12}, year = {2004}, abstract = {We focus on the role of anisotropic elasticity in the simulation of the load distribution in a human mandible due to a lateral bite on the leftmost premolar. Based on experimental evidence, we adopt ``local''" orthotropy of the elastic properties of the bone tissue. Since the trajectories of anisotropic elasticity are not accessible from Computer Tomographic (CT) data, they will be reconstructed from (i) the organ's geometry and (ii) from coherent structures which can be recognized from the spatial distribution of the CT values. A sensitivity analysis comprising various 3D FE simulations reveals the relevance of elastic anisotropy for the load carrying behavior of a human mandible: Comparison of the load distributions in isotropic and anisotropic simulations indicates that anisotropy seems to ``spare''" the mandible from loading. Moreover, a maximum degree of anisotropy leads to kind of an load minimization of the mandible, expressed by a minimum of different norms of local strain, evaluated throughout the organ. Thus, we may suggest that anisotropy is not only relevant, but also in some sense ``optimal''.}, language = {en} } @misc{KoberErdmannLangetal.2004, author = {Kober, Cornelia and Erdmann, Bodo and Lang, Jens and Sader, Robert and Zeilhofer, Hans-Florian}, title = {Adaptive Finite Element Simulation of the Human Mandible Using a New Physiological Model of the Masticatory Muscles}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7917}, number = {04-16}, year = {2004}, abstract = {Structural mechanics simulation of bony organs is of general medical and biomechanical interest, because of the interdependence of the inner architecture of bone and its functional loading already stated by Wolff in 1892. This work is part of a detailed research project concerning the human mandible. By adaptive finite element techniques, stress/strain profiles occurring in the bony structure under biting were simulated. Estimates of the discretization errors, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method. In general, our simulation requires a representation of the organ's geometry, an appropriate material description, and the load case due to teeth, muscle, or joint forces. In this paper, we want to focus on the influence of the masticatory system. Our goal is to capture the physiological situation as far as possible. By means of visualization techniques developed by the group, we are able to extract individual muscle fibres from computed tomography data. By a special algorithm, the fibres are expanded to fanlike (esp. for the musc. temporalis) coherent vector fields similar to the anatomical reality. The activity of the fibres can be adapted according to compartmentalisation of the muscles as measured by electromyological experiments. A refined sensitivity analysis proved remarkable impact of the presented approach on the simulation results.}, language = {en} } @misc{WeiserSchiela2004, author = {Weiser, Martin and Schiela, Anton}, title = {Function space interior point methods for PDE constrained optimization}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8027}, number = {04-27}, year = {2004}, abstract = {A primal-dual interior point method for optimal control problems with PDE constraints is considered. The algorithm is directly applied to the infinite dimensional problem. Existence and convergence of the central path are analyzed. Numerical results from an inexact continuation method applied to a model problem are shown.}, language = {en} } @article{SchenkWaechterWeiser2008, author = {Schenk, O. and W{\"a}chter, Andreas and Weiser, Martin}, title = {Inertia Revealing Preconditioning For Large-Scale Nonconvex Constrained Optimization}, volume = {31}, journal = {SIAM J. Sci. Comp.}, number = {2}, pages = {939 -- 960}, year = {2008}, language = {en} } @article{WeiserDeuflhard2007, author = {Weiser, Martin and Deuflhard, Peter}, title = {Inexact central path following algorithms for optimal control problems}, volume = {46}, journal = {SIAM J. Control Opt.}, number = {3}, pages = {792 -- 815}, year = {2007}, language = {en} } @article{BergmannGraichenRohlmannetal.2008, author = {Bergmann, G. and Graichen, F. and Rohlmann, A. and Bender, Alwina and Ehrig, Rainald and Westerhoff, P. and Heinlein, B.}, title = {Design and calibration of load sensing orthopaedic implants}, journal = {Journal Biomech. Eng.}, pages = {02109 -- 1-02109-9}, year = {2008}, language = {en} } @article{ColliFranzoneDeuflhardErdmannetal.2006, author = {Colli Franzone, Piero and Deuflhard, Peter and Erdmann, Bodo and Lang, Jens and Pavarino, Luca Franco}, title = {Adaptivity in Space and Time for Reaction-Diffusion Systems in Electrocardiology}, volume = {28}, journal = {SIAM J. Sc. Comp.}, pages = {942 -- 962}, year = {2006}, language = {en} } @article{DahlkeSteidlTeschke2004, author = {Dahlke, S. and Steidl, G. and Teschke, G.}, title = {Coorbit Spaces and Banach Frames on Homogeneous Spaces with Applications to Analyzing Functions on Spheres}, volume = {1-2}, journal = {Adv. Comput.Math.}, number = {21}, pages = {147 -- 180}, year = {2004}, language = {en} } @article{DahlkeSteidlTeschke2004, author = {Dahlke, S. and Steidl, G. and Teschke, G.}, title = {Weighted Coorbit Spaces and Banach Frames on Homogeneous Spaces}, volume = {5}, journal = {Journal of Fourier Analysis and Applications}, number = {10}, pages = {507 -- 539}, year = {2004}, language = {en} } @book{DeuflhardWeiser2011, author = {Deuflhard, Peter and Weiser, Martin}, title = {Numerische Mathematik 3}, publisher = {de Gruyter, Berlin}, year = {2011}, language = {de} } @book{DeuflhardWeiser2012, author = {Deuflhard, Peter and Weiser, Martin}, title = {Adaptive numerical solution of PDEs}, publisher = {de Gruyter}, address = {Berlin}, year = {2012}, language = {en} } @article{DeuflhardWeiserZachow2006, author = {Deuflhard, Peter and Weiser, Martin and Zachow, Stefan}, title = {Mathematics in Facial Surgery}, volume = {53}, journal = {AMS Notices}, number = {9}, pages = {1012 -- 1016}, year = {2006}, language = {en} } @article{DeuflhardHochmuth2004, author = {Deuflhard, Peter and Hochmuth, Reinhard}, title = {Multiscale analysis of thermoregulation in the human microvascular system}, volume = {27}, journal = {Math. Methods Appl. Sci.}, number = {8}, publisher = {Wiley Verlag}, pages = {971 -- 989}, year = {2004}, language = {en} } @incollection{Deuflhard2003, author = {Deuflhard, Peter}, title = {A Comparison of Related Concepts in Computational Chemistry and Mathematics}, volume = {88}, booktitle = {Chemisty and Mathematics}, number = {330}, pages = {51 -- 66}, year = {2003}, language = {en} } @incollection{Deuflhard2003, author = {Deuflhard, Peter}, title = {From Molecular Dynamics to Conformational Dynamics in Drug Design}, booktitle = {Trends in Nonlinear Analysis}, editor = {et al. Kirkilionis, Markus}, publisher = {Springer Verlag Berlin}, pages = {269 -- 288}, year = {2003}, language = {en} } @misc{DeuflhardZachow2012, author = {Deuflhard, Peter and Zachow, Stefan}, title = {Mathematische Therapie- und Operationsplanung}, publisher = {Berliner Wirtschaftsgespr{\"a}che e.V.}, address = {Berlin}, pages = {89 -- 90}, year = {2012}, language = {en} } @incollection{DeuflhardDoesselLouisetal.2010, author = {Deuflhard, Peter and D{\"o}ssel, Olaf and Louis, Alfred and Zachow, Stefan}, title = {More Mathematics into Medicine!}, booktitle = {Production Factor Mathematics}, publisher = {Springer}, pages = {357 -- 378}, year = {2010}, language = {en} } @article{WeiserZachowDeuflhard2010, author = {Weiser, Martin and Zachow, Stefan and Deuflhard, Peter}, title = {Craniofacial Surgery Planning Based on Virtual Patient Models}, volume = {52}, journal = {it - Information Technology}, number = {5}, publisher = {Oldenbourg Verlagsgruppe}, doi = {10.1524/itit.2010.0600}, pages = {258 -- 263}, year = {2010}, language = {en} } @inproceedings{ZachowHierlErdmann2004, author = {Zachow, Stefan and Hierl, Thomas and Erdmann, Bodo}, title = {A quantitative evaluation of 3D soft tissue prediction in maxillofacial surgery planning}, booktitle = {Proc. 3. Jahrestagung der Deutschen Gesellschaft f{\"u}r Computer- und Roboter-assistierte Chirurgie e.V.}, address = {M{\"u}nchen}, year = {2004}, language = {en} } @inproceedings{WeiserScacchi2017, author = {Weiser, Martin and Scacchi, Simone}, title = {Spectral Deferred Correction methods for adaptive electro-mechanical coupling in cardiac simulation}, booktitle = {G. Russo et al.(eds.) Progress in Industrial Mathematics at ECMI 2014}, publisher = {Springer}, doi = {10.1007/978-3-319-23413-7_42}, pages = {321 -- 328}, year = {2017}, abstract = {We investigate spectral deferred correction (SDC) methods for time stepping and their interplay with spatio-temporal adaptivity, applied to the solution of the cardiac electro-mechanical coupling model. This model consists of the Monodomain equations, a reaction-diffusion system modeling the cardiac bioelectrical activity, coupled with a quasi-static mechanical model describing the contraction and relaxation of the cardiac muscle. The numerical approximation of the cardiac electro-mechanical coupling is a challenging multiphysics problem, because it exhibits very different spatial and temporal scales. Therefore, spatio-temporal adaptivity is a promising approach to reduce the computational complexity. SDC methods are simple iterative methods for solving collocation systems. We exploit their flexibility for combining them in various ways with spatio-temporal adaptivity. The accuracy and computational complexity of the resulting methods are studied on some numerical examples.}, language = {en} } @misc{LubkollSchielaWeiser2015, author = {Lubkoll, Lars and Schiela, Anton and Weiser, Martin}, title = {An affine covariant composite step method for optimization with PDEs as equality constraints}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53954}, year = {2015}, abstract = {We propose a composite step method, designed for equality constrained optimization with partial differential equations. Focus is laid on the construction of a globalization scheme, which is based on cubic regularization of the objective and an affine covariant damped Newton method for feasibility. We show finite termination of the inner loop and fast local convergence of the algorithm. We discuss preconditioning strategies for the iterative solution of the arising linear systems with projected conjugate gradient. Numerical results are shown for optimal control problems subject to a nonlinear heat equation and subject to nonlinear elastic equations arising from an implant design problem in craniofacial surgery.}, language = {en} } @misc{WeiserGhosh2016, author = {Weiser, Martin and Ghosh, Sunayana}, title = {Theoretically optimal inexact SDC methods}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53140}, year = {2016}, abstract = {In several inital value problems with particularly expensive right hand side computation, there is a trade-off between accuracy and computational effort in evaluating the right hand sides. We consider inexact spectral deferred correction (SDC) methods for solving such non-stiff initial value problems. SDC methods are interpreted as fixed point iterations and, due to their corrective iterative nature, allow to exploit the accuracy-work-tradeoff for a reduction of the total computational effort. On one hand we derive an error model bounding the total error in terms of the right hand side evaluation errors. On the other hand, we define work models describing the computational effort in terms of the evaluation accuracy. Combining both, a theoretically optimal tolerance selection is worked out by minimizing the total work subject to achieving the requested tolerance.}, language = {en} } @phdthesis{Goetschel2015, author = {G{\"o}tschel, Sebastian}, title = {Adaptive Lossy Trajectory Compression for Optimal Control of Parabolic PDEs}, year = {2015}, abstract = {Optimal control problems governed by nonlinear, time-dependent PDEs on three-dimensional spatial domains are an important tool in many fields, ranging from engineering applications to medicine. For the solution of such optimization problems, methods working on the reduced objective functional are often employed to avoid a full spatio-temporal discretization of the problem. The evaluation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, requiring the storage of a full 4D data set. If Newton-CG methods are used, two additional trajectories have to be stored. To get numerical results that are accurate enough, in many cases very fine discretizations in time and space are necessary, leading to a significant amount of data to be stored and transmitted to mass storage. This thesis deals with the development and analysis of methods for lossy compression of such finite element solutions. The algorithms are based on a change of basis to reduce correlations in the data, combined with quantization. This is achieved by transforming the finite element coefficient vector from the nodal to the hierarchical basis, followed by rounding the coefficients to a prescribed precision. Due to the inexact reconstruction, and thus inexact data for the adjoint equation, the error induced in the reduced gradient, and reduced Hessian, has to be controlled, to not impede convergence of the optimization. Accuracy requirements of different optimization methods are analyzed, and computable error estimates for the influence of lossy trajectory storage are derived. These tools are used to adaptively control the accuracy of the compressed data. The efficiency of the algorithms is demonstrated on several numerical examples, ranging from a simple linear, scalar equation to a semi-linear system of reaction-diffusion equations. In all examples considerable reductions in storage space and bandwidth requirements are achieved, without significantly influencing the convergence behavior of the optimization methods. Finally, to go beyond pointwise error control, the hierarchical basis transform can be replaced by more sophisticated wavelet transforms. Numerical experiments indicate that choosing suitable norms for error control allows higher compression factors.}, language = {en} } @misc{Freytag2017, author = {Freytag, Yvonne}, title = {Optimal Experimental Design to Estimate the Time of Death in a Bayesian Context}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62475}, year = {2017}, abstract = {This thesis is devoted to the interdisciplinary work between mathematicians and forensic experts: the modeling of the human body cooling process after death laying the foundation for the estimation of the time of death. An inverse problem needs to be solved. In this thesis the inverse problem computes the time of death given the measured body temperature and the Forward Model that simulates the body cooling process. The Forward Model is based on the heat equation established by Fourier. This differential equation is numerically solved by the discretization over space by the Finite Element Method and the discretization over time by the Implicit Euler Method. The applications in this thesis demand a fast computation time. A model reduction is achieved by the Proper Orthogonal Decomposition in combination with the Galerkin Method. For reasons of simplification the computations and the measurements are restricted to a cylindrical phantom that is made out of homogeneous polyethylene. The estimate of the time of death is accompanied by an uncertainty. The inverse problem is incorporated by Bayesian inference to interpret the quality of the estimate and the effciency of the experiment. The uncertainty of the estimate of the time of death is minimized by approaching the Optimal Design of the Experiment. An objective function measures the certainty of the data and lays the foundation of the optimization problem. Solving the optimization problem is successfully done by relaxing the complex discrete NP-hard problem and applying a gradient-based method. The results of this thesis clearly show that the design of an experiment has a great in- uence on the outcome of the quality of the estimate. The comparison of the estimate and its properties based on different designs and conditions reveals the effciency of the Design of Experiment in the context of the estimation of the time of death.}, language = {en} } @inproceedings{GoetschelHoehneKolkoorietal.2016, author = {G{\"o}tschel, Sebastian and H{\"o}hne, Christian and Kolkoori, Sanjeevareddy and Mitzscherling, Steffen and Prager, Jens and Weiser, Martin}, title = {Ray Tracing Boundary Value Problems: Simulation and SAFT Reconstruction for Ultrasonic Testing}, booktitle = {Proceedings 19th World Conference on Non-Destructive Testing (WCNDT 2016)}, year = {2016}, language = {en} } @inproceedings{GoetschelMaierhoferMuelleretal.2016, author = {G{\"o}tschel, Sebastian and Maierhofer, Christiane and M{\"u}ller, Jan and Rothbart, Nick and Weiser, Martin}, title = {Quantitative Defect Reconstruction in Active Thermography for Fiber-Reinforced Composites}, booktitle = {Proceedings 19th World Conference on Non-Destructive Testing (WCNDT 2016)}, year = {2016}, language = {en} } @inproceedings{MuellerGoetschelMaierhoferetal.2017, author = {M{\"u}ller, Jan and G{\"o}tschel, Sebastian and Maierhofer, Christiane and Weiser, Martin}, title = {Determining the material parameters for the reconstruction of defects in carbon fiber reinforced polymers from data measured by flash thermography}, volume = {1806}, booktitle = {AIP Conference Proceedings}, doi = {10.1063/1.4974671}, year = {2017}, language = {en} } @misc{Schuh2016, type = {Master Thesis}, author = {Schuh, Hannah}, title = {Identification of Anisotropies in Welding Seams Using Ultrasound Measurements}, year = {2016}, language = {en} } @misc{FischerGoetschelWeiser2017, author = {Fischer, Lisa and G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Lossy data compression reduces communication time in hybrid time-parallel integrators}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63961}, year = {2017}, abstract = {Parallel in time methods for solving initial value problems are a means to increase the parallelism of numerical simulations. Hybrid parareal schemes interleaving the parallel in time iteration with an iterative solution of the individual time steps are among the most efficient methods for general nonlinear problems. Despite the hiding of communication time behind computation, communication has in certain situations a significant impact on the total runtime. Here we present strict, yet no sharp, error bounds for hybrid parareal methods with inexact communication due to lossy data compression, and derive theoretical estimates of the impact of compression on parallel efficiency of the algorithms. These and some computational experiments suggest that compression is a viable method to make hybrid parareal schemes robust with respect to low bandwidth setups.}, language = {en} } @misc{WeiserErdmannSchenkletal.2017, author = {Weiser, Martin and Erdmann, Bodo and Schenkl, Sebastian and Muggenthaler, Holger and Hubig, Michael and Mall, Gita and Zachow, Stefan}, title = {Uncertainty in Temperature-Based Determination of Time of Death}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63818}, year = {2017}, abstract = {Temperature-based estimation of time of death (ToD) can be per- formed either with the help of simple phenomenological models of corpse cooling or with detailed mechanistic (thermodynamic) heat transfer mod- els. The latter are much more complex, but allow a higher accuracy of ToD estimation as in principle all relevant cooling mechanisms can be taken into account. The potentially higher accuracy depends on the accuracy of tissue and environmental parameters as well as on the geometric resolution. We in- vestigate the impact of parameter variations and geometry representation on the estimated ToD based on a highly detailed 3D corpse model, that has been segmented and geometrically reconstructed from a computed to- mography (CT) data set, differentiating various organs and tissue types. From that we identify the most crucial parameters to measure or estimate, and obtain a local uncertainty quantifcation for the ToD.}, language = {en} } @misc{GoetschelMinion2017, author = {G{\"o}tschel, Sebastian and Minion, Michael L.}, title = {Parallel-in-Time for Parabolic Optimal Control Problems Using PFASST}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64989}, year = {2017}, abstract = {In gradient-based methods for parabolic optimal control problems, it is necessary to solve both the state equation and a backward-in-time adjoint equation in each iteration of the optimization method. In order to facilitate fully parallel gradient-type and nonlinear conjugate gradient methods for the solution of such optimal control problems, we discuss the application of the parallel-in-time method PFASST to adjoint gradient computation. In addition to enabling time parallelism, PFASST provides high flexibility for handling nonlinear equations, as well as potential extra computational savings from reusing previous solutions in the optimization loop. The approach is demonstrated here for a model reaction-diffusion optimal control problem.}, language = {en} } @article{FischerGoetschelWeiser2018, author = {Fischer, Lisa and G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Lossy data compression reduces communication time in hybrid time-parallel integrators}, volume = {19}, journal = {Comput. Vis. Sci.}, number = {1}, doi = {10.1007/s00791-018-0293-2}, pages = {19 -- 30}, year = {2018}, abstract = {Parallel in time methods for solving initial value problems are a means to increase the parallelism of numerical simulations. Hybrid parareal schemes interleaving the parallel in time iteration with an iterative solution of the individual time steps are among the most efficient methods for general nonlinear problems. Despite the hiding of communication time behind computation, communication has in certain situations a significant impact on the total runtime. Here we present strict, yet no sharp, error bounds for hybrid parareal methods with inexact communication due to lossy data compression, and derive theoretical estimates of the impact of compression on parallel efficiency of the algorithms. These and some computational experiments suggest that compression is a viable method to make hybrid parareal schemes robust with respect to low bandwidth setups.}, language = {en} } @misc{Freytag2017, type = {Master Thesis}, author = {Freytag, Yvonne}, title = {Optimal Experimental Design to Estimate the Time of Death in a Bayesian Context}, school = {Zuse Institute Berlin (ZIB)}, pages = {81}, year = {2017}, abstract = {This thesis is devoted to the interdisciplinary work between mathematicians and forensic experts: the modeling of the human body cooling process after death laying the foundation for the estimation of the time of death. An inverse problem needs to be solved. In this thesis the inverse problem computes the time of death given the measured body temperature and the Forward Model that simulates the body cooling process. The Forward Model is based on the heat equation established by Fourier. This differential equation is numerically solved by the discretization over space by the Finite Element Method and the discretization over time by the Implicit Euler Method. The applications in this thesis demand a fast computation time. A model reduction is achieved by the Proper Orthogonal Decomposition in combination with the Galerkin Method. For reasons of simplification the computations and the measurements are restricted to a cylindrical phantom that is made out of homogeneous polyethylene. The estimate of the time of death is accompanied by an uncertainty. The inverse problem is incorporated by Bayesian inference to interpret the quality of the estimate and the effciency of the experiment. The uncertainty of the estimate of the time of death is minimized by approaching the Optimal Design of the Experiment. An objective function measures the certainty of the data and lays the foundation of the optimization problem. Solving the optimization problem is successfully done by relaxing the complex discrete NP-hard problem and applying a gradient-based method. The results of this thesis clearly show that the design of an experiment has a great in- uence on the outcome of the quality of the estimate. The comparison of the estimate and its properties based on different designs and conditions reveals the effciency of the Design of Experiment in the context of the estimation of the time of death.}, language = {en} } @misc{Fischer2017, type = {Master Thesis}, author = {Fischer, Lisa}, title = {On the convergence of inexact time parallel time integration}, pages = {61}, year = {2017}, language = {en} }