@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{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{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} } @inproceedings{SagnolHegeWeiser2016, author = {Sagnol, Guillaume and Hege, Hans-Christian and Weiser, Martin}, title = {Using sparse kernels to design computer experiments with tunable precision}, booktitle = {22nd Intern. Conf. on Computational Statistics - COMPSTAT 2016, Oviedo, Spain, 23-26 August 2016, Proceedings ISBN 978-90-73592-36-0}, pages = {397 -- 408}, year = {2016}, 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{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} } @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{WeiserScacchi2014, author = {Weiser, Martin and Scacchi, Simone}, title = {Spectral Deferred Correction methods for adaptive electro-mechanical coupling in cardiac simulation}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-50695}, year = {2014}, 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{DeuflhardKornhuberSanderetal.2014, author = {Deuflhard, Peter and Kornhuber, Ralf and Sander, Oliver and Schiela, Anton and Weiser, Martin}, title = {Mathematics cures virtual patients}, volume = {1}, journal = {MATHEON-Mathematics for Key Technologies}, editor = {Deuflhard, Peter and Gr{\"o}tschel, Martin and H{\"o}mberg, Dietmar and Horst, Ulrich and Kramer, J{\"u}rg and Mehrmann, Volker and Polthier, Konrad and Schmidt, Frank and Sch{\"u}tte, Christof and Skutella, Martin and Sprekels, J{\"u}rgen}, publisher = {European Mathematical Society}, pages = {7 -- 25}, year = {2014}, language = {en} } @article{MoualeuNgangueWeiserEhrigetal.2015, author = {Moualeu-Ngangue, Dany Pascal and Weiser, Martin and Ehrig, Rainald and Deuflhard, Peter}, title = {Optimal control for a tuberculosis model with undetected cases in Cameroon}, volume = {20}, journal = {Communications in Nonlinear Science and Numerical Simulation}, number = {3}, doi = {10.1016/j.cnsns.2014.06.037}, pages = {986 -- 1003}, year = {2015}, abstract = {This paper considers the optimal control of tuberculosis through education, diagnosis campaign and chemoprophylaxis of latently infected. A mathematical model which includes important components such as undiagnosed infectious, diagnosed infectious, latently infected and lost-sight infectious is formulated. The model combines a frequency dependent and a density dependent force of infection for TB transmission. Through optimal control theory and numerical simulations, a cost-effective balance of two different intervention methods is obtained. Seeking to minimize the amount of money the government spends when tuberculosis remain endemic in the Cameroonian population, Pontryagin's maximum principle is used to characterize the optimal control. The optimality system is derived and solved numerically using the forward-backward sweep method (FBSM). Results provide a framework for designing cost-effective strategies for diseases with multiple intervention methods. It comes out that combining chemoprophylaxis and education, the burden of TB can be reduced by 80 \% in 10 years.}, language = {en} } @misc{GoetschelWeiserMaierhoferetal.2012, author = {G{\"o}tschel, Sebastian and Weiser, Martin and Maierhofer, Christiane and Richter, Regina}, title = {Data Enhancement for Active Thermography}, journal = {E-book Proceedings, 11th International Conference on Quantitative Infrared Thermography, Naples}, editor = {Cardone, Gennaro}, year = {2012}, abstract = {Pulse thermography is a non-destructive testing method based on infrared imaging of transient thermal patterns. Heating the surface of the structure under test for a short period of time generates a non-stationary temperature distribution and thus a thermal contrast between the defect and the sound material. Due to measurement noise, preprocessing of the experimental data is necessary, before reconstruction algorithms can be applied. We propose a decomposition of the measured temperature into Green's function solutions to eliminate noise.}, language = {en} } @article{Weiser2013, author = {Weiser, Martin}, title = {On goal-oriented adaptivity for elliptic optimal control problems}, volume = {28}, journal = {Opt. Meth. Softw.}, number = {13}, pages = {969 -- 992}, year = {2013}, abstract = {The paper proposes goal-oriented error estimation and mesh refinement for optimal control problems with elliptic PDE constraints using the value of the reduced cost functional as quantity of interest. Error representation, hierarchical error estimators, and greedy-style error indicators are derived and compared to their counterparts when using the all-at-once cost functional as quantity of interest. Finally, the efficiency of the error estimator and generated meshes are demonstrated on numerical examples.}, language = {en} } @article{DeuflhardSchielaWeiser2012, author = {Deuflhard, Peter and Schiela, Anton and Weiser, Martin}, title = {Mathematical Cancer Therapy Planning in Deep Regional Hyperthermia}, volume = {21}, journal = {Acta Numerica}, pages = {307 -- 378}, year = {2012}, abstract = {This paper surveys the required mathematics for a typical challenging problem from computational medicine, the cancer therapy planning in deep regional hyperthermia. In the course of many years of close cooperation with clinics, the medical problem gave rise to quite a number of subtle mathematical problems, part of which had been unsolved when the common project started. Efficiency of numerical algorithms, i.e. computational speed and monitored reliability, play a decisive role for the medical treatment. Off-the-shelf software had turned out to be not sufficient to meet the requirements of medicine. Rather, new mathematical theory as well as new numerical algorithms had to be developed. In order to make our algorithms useful in the clinical environment, new visualization software, a virtual lab, including 3D geometry processing of individual virtual patients had to be designed and implemented. Moreover, before the problems could be attacked by numerical algorithms, careful mathematical modelling had to be done. Finally, parameter identification and constrained optimization for the PDEs had to be newly analyzed and realized over the individual patient's geometry. Our new techniques had an impact on the specificity of the individual patients' treatment and on the construction of an improved hyperthermia applicator.}, language = {en} } @misc{GoetschelWeiserSchiela2012, author = {G{\"o}tschel, Sebastian and Weiser, Martin and Schiela, Anton}, title = {Solving Optimal Control Problems with the Kaskade 7 Finite Element Toolbox}, journal = {Advances in DUNE}, editor = {Dedner, A. and Flemisch, B. and Kl{\"o}fkorn, R.}, publisher = {Springer}, pages = {101 -- 112}, year = {2012}, abstract = {This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study.}, language = {en} } @article{Weiser2015, author = {Weiser, Martin}, title = {Faster SDC convergence on non-equidistant grids by DIRK sweeps}, volume = {55}, journal = {BIT Numerical Mathematics}, number = {4}, doi = {10.1007/s10543-014-0540-y}, pages = {1219 -- 1241}, year = {2015}, abstract = {Spectral deferred correction methods for solving stiff ODEs are known to converge rapidly towards the collocation limit solution on equidistant grids, but show a much less favourable contraction on non-equidistant grids such as Radau-IIa points. We interprete SDC methods as fixed point iterations for the collocation system and propose new DIRK-type sweeps for stiff problems based on purely linear algebraic considerations. Good convergence is recovered also on non-equidistant grids. The properties of different variants are explored on a couple of numerical examples.}, language = {en} } @article{WeiserGoetschel2012, author = {Weiser, Martin and G{\"o}tschel, Sebastian}, title = {State Trajectory Compression for Optimal Control with Parabolic PDEs}, volume = {34}, journal = {SIAM J. Sci. Comput.}, number = {1}, doi = {10.1137/11082172X}, pages = {A161 -- A184}, year = {2012}, abstract = {In optimal control problems with nonlinear time-dependent 3D PDEs, full 4D discretizations are usually prohibitive due to the storage requirement. For this reason gradient and quasi-Newton methods working on the reduced functional are often employed. The computation 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, again requiring the storage of a full 4D data set. We propose a lossy compression algorithm using an inexact but cheap predictor for the state data, with additional entropy coding of prediction errors. As the data is used inside a discretized, iterative algorithm, lossy coding maintaining an error bound is sufficient.}, language = {en} } @misc{NadobnyWeihrauchWeiseretal.2007, author = {Nadobny, Johanna and Weihrauch, Mirko and Weiser, Martin and Gellermann, Johanna and Wlodarczyk, Waldemar and Budach, Volker and Wust, Peter}, title = {Advances in the Planning and Control of the MR-guided Regional Hyperthermia Applications}, journal = {Proc. Int. Conf. Electromagnetics in Advanced Applications, ICEAA 2007, Torino, Italy}, pages = {1010 -- 1013}, year = {2007}, language = {en} } @article{GoetschelWeiser2010, author = {G{\"o}tschel, Sebastian and Weiser, Martin}, title = {State Trajectory Compression in Optimal Control}, volume = {10}, journal = {PAMM}, number = {1}, doi = {10.1002/pamm.201010282}, pages = {579 -- 580}, year = {2010}, abstract = {In optimal control problems with nonlinear time-dependent 3D PDEs, the computation of the reduced gradient by adjoint methods requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. Since the state enters into the adjoint equation, the storage of a 4D discretization is necessary. We propose a lossy compression algorithm using a cheap predictor for the state data, with additional entropy coding of prediction errors. Analytical and numerical results indicate that compression factors around 30 can be obtained without exceeding the FE discretization error.}, language = {en} } @misc{Weiser2010, author = {Weiser, Martin}, title = {Delayed Residual Compensation for Bidomain Equations}, volume = {1281}, journal = {AIP Conference Proceedings}, doi = {10.1063/1.3498495}, pages = {419 -- 422}, year = {2010}, abstract = {The biodomain model of cardioelectric excitation consists of a reaction-diffusion equation, an elliptic algebraic constraint, and a set of pointwise ODEs. Fast reaction enforces small time steps, such that for common mesh sizes the reaction-diffusion equation is easily solved implicitly due to a dominating mass matrix. In contrast, the elliptic constraint does not benefit from small time steps and requires a comparably expensive solution. We propose a delayed residual compensation that improves the solution of the elliptic constraint and thus alleviates the need for long iteration times.}, language = {en} } @misc{WeiserErdmannDeuflhard2010, author = {Weiser, Martin and Erdmann, Bodo and Deuflhard, Peter}, title = {On Efficiency and Accuracy in Cardioelectric Simulation}, journal = {Progress in Industrial Mathematics at ECMI 2008}, editor = {Wilson, E. and Fitt, A. and Ockendon, H. and Norbury, J.}, publisher = {Springer}, pages = {371 -- 376}, year = {2010}, abstract = {Reasons for the failure of adaptive methods to deliver improved efficiency when integrating monodomain models for myocardiac excitation are discussed. Two closely related techniques for reducing the computational complexity of linearly implicit integrators, deliberate sparsing and splitting, are investigated with respect to their impact on computing time and accuracy.}, language = {en} } @article{WilhelmsSeemannWeiseretal.2010, author = {Wilhelms, Mathias and Seemann, Gunnar and Weiser, Martin and D{\"o}ssel, Olaf}, title = {Benchmarking Solvers of the Monodomain Equation in Cardiac Electrophysiological Modeling}, volume = {55}, journal = {Biomed. Engineer.}, doi = {10.1515/BMT.2010.712}, pages = {99 -- 102}, year = {2010}, language = {en} } @misc{WustWeihrauchWeiseretal.2010, author = {Wust, Peter and Weihrauch, Mirko and Weiser, Martin and Gellermann, Johanna and Eisenhardt, Steffen and Chobrok, Thorsten and Budach, Volker}, title = {Optimization of clinical radiofrequency hyperthermia by use of MR-thermography in a hybrid system}, journal = {World Congress on Medical Physics and Biomedical Engineering, September 2009, Munich, Germany}, editor = {D{\"o}ssel, O. and Schlegel, W. and Magjarevic, R.}, publisher = {Springer}, pages = {174 -- 175}, year = {2010}, language = {en} } @article{LubkollSchielaWeiser2014, author = {Lubkoll, Lars and Schiela, Anton and Weiser, Martin}, title = {An optimal control problem in polyconvex hyperelasticity}, volume = {52}, journal = {SIAM J. Control Opt.}, number = {3}, doi = {10.1137/120876629}, pages = {1403 -- 1422}, year = {2014}, abstract = {We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented.}, language = {en} } @article{GuentherLameckerWeiser2013, author = {G{\"u}nther, Andreas and Lamecker, Hans and Weiser, Martin}, title = {Flexible Shape Matching with Finite Element Based LDDMM}, volume = {105}, journal = {International Journal of Computer Vision}, number = {2}, doi = {10.1007/s11263-012-0599-3}, pages = {128 -- 143}, year = {2013}, abstract = {We consider Large Deformation Diffeomorphic Metric Mapping of general \$m\$-currents. After stating an optimization algorithm in the function space of admissable morph generating velocity fields, two innovative aspects in this framework are presented and numerically investigated: First, we spatially discretize the velocity field with conforming adaptive finite elements and discuss advantages of this new approach. Second, we directly compute the temporal evolution of discrete \$m\$-current attributes.}, language = {en} } @misc{GuentherLameckerWeiser2011, author = {G{\"u}nther, Andreas and Lamecker, Hans and Weiser, Martin}, title = {Direct LDDMM of Discrete Currents with Adaptive Finite Elements}, journal = {Proceedings of the Third International Workshop on Mathematical Foundations of Computational Anatomy - Geometrical and Statistical Methods for Modelling Biological Shape Variability}, editor = {Pennec, X. and Joshi, S. and Nielsen, M.}, pages = {1 -- 14}, year = {2011}, abstract = {We consider Large Deformation Diffeomorphic Metric Mapping of general \$m\$-currents. After stating an optimization algorithm in the function space of admissable morph generating velocity fields, two innovative aspects in this framework are presented and numerically investigated: First, we spatially discretize the velocity field with conforming adaptive finite elements and discuss advantages of this new approach. Second, we directly compute the temporal evolution of discrete \$m\$-current attributes.}, language = {en} } @misc{SchielaWeiser2010, author = {Schiela, Anton and Weiser, Martin}, title = {Barrier methods for a control problem from hyperthermia treatment planning}, journal = {Recent Advances in Optimization and its Applications in Engineering (Proceedings of 14th Belgian-French-German Conference on Optimization 2009)}, editor = {Diehl, M. and Glineur, F. and Jarlebring, E. and Michiels, W.}, publisher = {Springer}, pages = {419 -- 428}, year = {2010}, abstract = {We consider an optimal control problem from hyperthermia treatment planning and its barrier regularization. We derive basic results, which lay the groundwork for the computation of optimal solutions via an interior point path-following method. Further, we report on a numerical implementation of such a method and its performance at an example problem.}, language = {en} } @article{WeiserRoelligArndtetal.2010, author = {Weiser, Martin and R{\"o}llig, Mathias and Arndt, Ralf and Erdmann, Bodo}, title = {Development and test of a numerical model for pulse thermography in civil engineering}, volume = {46}, journal = {Heat and Mass Transfer}, number = {11-12}, pages = {1419 -- 1428}, year = {2010}, abstract = {Pulse thermography of concrete structures is used in civil engineering for detecting voids, honeycombing and delamination. The physical situation is readily modeled by Fourier's law. Despite the simplicity of the PDE structure, quantitatively realistic numerical 3D simulation faces two major obstacles. First, the short heating pulse induces a thin boundary layer at the heated surface which encapsulates all information and therefore has to be resolved faithfully. Even with adaptive mesh refinement techniques, obtaining useful accuracies requires an unsatisfactorily fine discretization. Second, bulk material parameters and boundary conditions are barely known exactly. We address both issues by a semi-analytic reformulation of the heat transport problem and by parameter identification. Numerical results are compared with measurements of test specimens.}, language = {en} } @article{RannebergWeiserWeihrauchetal.2010, author = {Ranneberg, Maximilian and Weiser, Martin and Weihrauch, Mirko and Budach, Volker and Gellermann, Johanna and Wust, Peter}, title = {Regularized Antenna Profile Adaptation in Online Hyperthermia Treatment}, volume = {37}, journal = {Medical Physics}, doi = {10.1118/1.3488896}, pages = {5382 -- 5394}, year = {2010}, language = {en} } @article{LubkollSchielaWeiser2017, author = {Lubkoll, Lars and Schiela, Anton and Weiser, Martin}, title = {An affine covariant composite step method for optimization with PDEs as equality constraints}, volume = {32}, journal = {Optimization Methods and Software}, number = {5}, doi = {10.1080/10556788.2016.1241783}, pages = {1132 -- 1161}, year = {2017}, 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{GoetschelMaierhoferMuelleretal.2016, author = {G{\"o}tschel, Sebastian and Maierhofer, Christiane and M{\"u}ller, Jan P. and Rothbart, Nick and Weiser, Martin}, title = {Quantitative Defect Reconstruction in Active Thermography for Fiber-Reinforced Composites}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-58374}, year = {2016}, abstract = {Carbon-fiber reinforced composites are becoming more and more important in the production of light-weight structures, e.g., in the automotive and aerospace industry. Thermography is often used for non-destructive testing of these products, especially to detect delaminations between different layers of the composite. In this presentation, we aim at methods for defect reconstruction from thermographic measurements of such carbon-fiber reinforced composites. The reconstruction results shall not only allow to locate defects, but also give a quantitative characterization of the defect properties. We discuss the simulation of the measurement process using finite element methods, as well as the experimental validation on flat bottom holes. Especially in pulse thermography, thin boundary layers with steep temperature gradients occurring at the heated surface need to be resolved. Here we use the combination of a 1D analytical solution combined with numerical solution of the remaining defect equation. We use the simulations to identify material parameters from the measurements. Finally, fast heuristics for reconstructing defect geometries are applied to the acquired data, and compared for their accuracy and utility in detecting different defects like back surface defects or delaminations.}, language = {en} } @misc{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}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-58386}, year = {2016}, abstract = {The application of advanced imaging techniques for the ultrasonic inspection of inhomogeneous anisotropic materials like austenitic and dissimilar welds requires information about acoustic wave propagation through the material, in particular travel times between two points in the material. Forward ray tracing is a popular approach to determine traveling paths and arrival times but is ill suited for inverse problems since a large number of rays have to be computed in order to arrive at prescribed end points. In this contribution we discuss boundary value problems for acoustic rays, where the ray path between two given points is determined by solving the eikonal equation. The implementation of such a two point boundary value ray tracer for sound field simulations through an austenitic weld is described and its efficiency as well as the obtained results are compared to those of a forward ray tracer. The results are validated by comparison with experimental results and commercially available UT simulation tools. As an application, we discuss an implementation of the method for SAFT (Synthetic Aperture Focusing Technique) reconstruction. The ray tracer calculates the required travel time through the anisotropic columnar grain structure of the austenitic weld. There, the formulation of ray tracing as a boundary value problem allows a straightforward derivation of the ray path from a given transducer position to any pixel in the reconstruction area and reduces the computational cost considerably.}, language = {en} } @book{Weiser2016, author = {Weiser, Martin}, title = {Inside Finite Elements}, publisher = {De Gruyter}, year = {2016}, abstract = {All relevant implementation aspects of finite element methods are discussed in this book. The focus is on algorithms and data structures as well as on their concrete implementation. Theory is covered as far as it gives insight into the construction of algorithms.Throughout the exercises a complete FE-solver for scalar 2D problems will be implemented in Matlab/Octave.}, language = {en} } @misc{VolkweinWeiser2000, author = {Volkwein, Stefan and Weiser, Martin}, title = {Affine Invariant Convergence Analysis for Inexact Augmented Lagrangian-SQP Methods}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6243}, number = {00-56}, year = {2000}, abstract = {An affine invariant convergence analysis for inexact augmented Lagrangian-SQP methods is presented. The theory is used for the construction of an accuracy matching between iteration errors and truncation errors, which arise from the inexact linear system solves. The theoretical investigations are illustrated numerically by an optimal control problem for the Burgers equation.}, language = {en} } @misc{WeiserDeuflhard2001, author = {Weiser, Martin and Deuflhard, Peter}, title = {The Central Path towards the Numerical Solution of Optimal Control Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6380}, number = {01-12}, year = {2001}, abstract = {A new approach to the numerical solution of optimal control problems including control and state constraints is presented. Like hybrid methods, the approach aims at combining the advantages of direct and indirect methods. Unlike hybrid methods, however, our method is directly based on interior-point concepts in function space --- realized via an adaptive multilevel scheme applied to the complementarity formulation and numerical continuation along the central path. Existence of the central path and its continuation towards the solution point is analyzed in some theoretical detail. An adaptive stepsize control with respect to the duality gap parameter is worked out in the framework of affine invariant inexact Newton methods. Finally, the performance of a first version of our new type of algorithm is documented by the successful treatment of the well-known intricate windshear problem.}, language = {en} } @misc{Weiser2002, author = {Weiser, Martin}, title = {Linear convergence of an interior point method for linear control constrained optimal control problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6809}, number = {02-13}, year = {2002}, abstract = {The paper provides a detailed analysis of a short step interior point algorithm applied to linear control constrained optimal control problems. Using an affine invariant local norm and an inexact Newton corrector, the well-known convergence results from finite dimensional linear programming can be extended to the infinite dimensional setting of optimal control. The present work complements a recent paper of Weiser and Deuflhard, where convergence rates have not been derived. The choice of free parameters, i.e. the corrector accuracy and the number of corrector steps, is discussed.}, language = {en} } @misc{Weiser2007, author = {Weiser, Martin}, title = {Pointwise Nonlinear Scaling for Reaction-Diffusion-Equations}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-10493}, number = {07-45}, year = {2007}, abstract = {Parabolic reaction-diffusion systems may develop sharp moving reaction fronts which pose a challenge even for adaptive finite element methods. We propose a method to transform the equation into an equivalent form that usually exhibits solutions which are easier to discretize, giving higher accuracy for a given number of degrees of freedom. The transformation is realized as an efficiently computable pointwise nonlinear scaling that is optimized for prototypical planar travelling wave solutions of the underlying reaction-diffusion equation. The gain in either performance or accuracy is demonstrated on different numerical examples.}, language = {en} } @misc{SchenkWaechterWeiser2007, author = {Schenk, Olaf and W{\"a}chter, Andreas and Weiser, Martin}, title = {Inertia Revealing Preconditioning For Large-Scale Nonconvex Constrained Optimization}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-10314}, number = {07-32}, year = {2007}, abstract = {Fast nonlinear programming methods following the all-at-once approach usually employ Newton's method for solving linearized Karush-Kuhn-Tucker (KKT) systems. In nonconvex problems, the Newton direction is only guaranteed to be a descent direction if the Hessian of the Lagrange function is positive definite on the nullspace of the active constraints, otherwise some modifications to Newton's method are necessary. This condition can be verified using the signs of the KKT's eigenvalues (inertia), which are usually available from direct solvers for the arising linear saddle point problems. Iterative solvers are mandatory for very large-scale problems, but in general do not provide the inertia. Here we present a preconditioner based on a multilevel incomplete \$LBL^T\$ factorization, from which an approximation of the inertia can be obtained. The suitability of the heuristics for application in optimization methods is verified on an interior point method applied to the CUTE and COPS test problems, on large-scale 3D PDE-constrained optimal control problems, as well as 3D PDE-constrained optimization in biomedical cancer hyperthermia treatment planning. The efficiency of the preconditioner is demonstrated on convex and nonconvex problems with \$150^3\$ state variables and \$150^2\$ control variables, both subject to bound constraints.}, language = {en} } @misc{GriesseWeiser2005, author = {Griesse, Roland and Weiser, Martin}, title = {On the Interplay Between Interior Point Approximation and Parametric Sensitivities in Optimal Control}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8773}, number = {05-44}, year = {2005}, abstract = {This paper is concerned with the sensitivities of function space oriented interior point approximations in parameter dependent problems. For an abstract setting that covers control constrained optimal control problems, the convergence of interior point sensitivities to the sensitivities of the optimal solution is shown. Error bounds for \$L_q\$ norms are derived and illustrated with numerical examples.}, language = {en} } @misc{WeiserHubigShanmugamSubramaniam2025, author = {Weiser, Martin and Hubig, Michael and Shanmugam Subramaniam, Jayant}, title = {Reconstructing Ambient Temperature Changes in Death Time Estimation with a Bayesian Double-Exponential Approach}, journal = {Zenodo}, doi = {10.5281/zenodo.17702240}, year = {2025}, abstract = {Code and data for the reconstruction of ambient temperature drop in time of death estimation We provide Octave code and temperature measurement data for - empirircally estimating thermal sensor likelihood - estimating time and amplitude of a single sudden ambient temperature drop from temperature measurement data in two thermally different compartments.}, language = {en} } @misc{VillaniWeiser2024, author = {Villani, Paolo and Weiser, Martin}, title = {Adaptive Gaussian process regression for inverse problems. ALGORITMY 2024}, year = {2024}, abstract = {This submission contains the code used for the proceedings paper for ALGORITMY 2024, concerning surrogate model-based inverse problems. Version 2 modified plot appearances and fixed some minor typos.}, language = {en} } @inproceedings{SemlerWeiser2025, author = {Semler, Phillip and Weiser, Martin}, title = {Adaptive gradient-enhanced Gaussian process surrogates for inverse problems}, booktitle = {Mathematical Optimization for Machine Learning: Proceedings of the MATH+ Thematic Einstein Semester 2023}, editor = {Fackeldey, Konstantin and Kannan, Aswin and Pokutta, Sebastian and Sharma, Kartikey and Walter, Daniel and Walther, Andrea and Weiser, Martin}, publisher = {De Gruyter}, arxiv = {http://arxiv.org/abs/2404.01864}, doi = {10.1515/9783111376776-005}, pages = {59 -- 78}, year = {2025}, abstract = {Generating simulated training data needed for constructing sufficiently accurate surrogate models to be used for efficient optimization or parameter identification can incur a huge computational effort in the offline phase. We consider a fully adaptive greedy approach to the computational design of experiments problem using gradient-enhanced Gaussian process regression as surrogates. Designs are incrementally defined by solving an optimization problem for accuracy given a certain computational budget. We address not only the choice of evaluation points but also of required simulation accuracy, both of values and gradients of the forward model. Numerical results show a significant reduction of the computational effort compared to just position-adaptive and static designs as well as a clear benefit of including gradient information into the surrogate training.}, language = {en} } @inproceedings{VillaniUngerWeiser2024, author = {Villani, Paolo and Unger, J{\"o}rg F. and Weiser, Martin}, title = {Adaptive Gaussian Process Regression for Bayesian inverse problems}, booktitle = {Proceedings of the Conference Algoritmy 2024}, pages = {214 -- 224}, year = {2024}, abstract = {We introduce a novel adaptive Gaussian Process Regression (GPR) methodology for efficient construction of surrogate models for Bayesian inverse problems with expensive forward model evaluations. An adaptive design strategy focuses on optimizing both the positioning and simulation accuracy of training data in order to reduce the computational cost of simulating training data without compromising the fidelity of the posterior distributions of parameters. The method interleaves a goal-oriented active learning algorithm selecting evaluation points and tolerances based on the expected impact on the Kullback-Leibler divergence of surrogated and true posterior with a Markov Chain Monte Carlo sampling of the posterior. The performance benefit of the adaptive approach is demonstrated for two simple test problems.}, language = {en} } @article{AndresArconesWeiserKoutsourelakisetal.2023, author = {Andr{\´e}s Arcones, Daniel and Weiser, Martin and Koutsourelakis, Faidon-Stelios and Unger, J{\"o}rg F.}, title = {A Bayesian Framework for Simulation-based Digital Twins of Bridges}, volume = {6}, journal = {EUROSTRUCT 2023:European Association on Quality Control of Bridges and Structures: Digital Transformation in Sustainability}, number = {5}, doi = {10.1002/cepa.2177}, pages = {734 -- 740}, year = {2023}, abstract = {Simulation-based digital twins have emerged as a powerful tool for evaluating the mechanical response of bridges. As virtual representations of physical systems, digital twins can provide a wealth of information that complements traditional inspection and monitoring data. By incorporating virtual sensors and predictive maintenance strategies, they have the potential to improve our understanding of the behavior and performance of bridges over time. However, as bridges age and undergo regular loading and extreme events, their structural characteristics change, often differing from the predictions of their initial design. Digital twins must be continuously adapted to reflect these changes. In this article, we present a Bayesian framework for updating simulation-based digital twins in the context of bridges. Our approach integrates information from measurements to account for inaccuracies in the simulation model and quantify uncertainties. Through its implementation and assessment, this work demonstrates the potential for digital twins to provide a reliable and up-to-date representation of bridge behavior, helping to inform decision-making for maintenance and management.}, language = {en} } @inproceedings{AndresArconesWeiserKoutsourelakisetal.2023, author = {Andr{\´e}s Arcones, Daniel and Weiser, Martin and Koutsourelakis, Faidon-Stelios and Unger, J{\"o}rg F.}, title = {Evaluation of Model Bias Identification Approaches Based on Bayesian Inference and Applications to Digital Twins}, booktitle = {5th ECCOMAS Thematic Conference on Uncertainty Quantification in Computational Sciences and Engineering}, pages = {1 -- 15}, year = {2023}, abstract = {In recent years, the use of simulation-based digital twins for monitoring and assessment of complex mechanical systems has greatly expanded. Their potential to increase the information obtained from limited data makes them an invaluable tool for a broad range of real-world applications. Nonetheless, there usually exists a discrepancy between the predicted response and the measurements of the system once built. One of the main contributors to this difference in addition to miscalibrated model parameters is the model error. Quantifying this socalled model bias (as well as proper values for the model parameters) is critical for the reliable performance of digital twins. Model bias identification is ultimately an inverse problem where information from measurements is used to update the original model. Bayesian formulations can tackle this task. Including the model bias as a parameter to be inferred enables the use of a Bayesian framework to obtain a probability distribution that represents the uncertainty between the measurements and the model. Simultaneously, this procedure can be combined with a classic parameter updating scheme to account for the trainable parameters in the original model. This study evaluates the effectiveness of different model bias identification approaches based on Bayesian inference methods. This includes more classical approaches such as direct parameter estimation using MCMC in a Bayesian setup, as well as more recent proposals such as stat-FEM or orthogonal Gaussian Processes. Their potential use in digital twins, generalization capabilities, and computational cost is extensively analyzed.}, language = {en} } @article{MaierWeiserConrad2025, author = {Maier, Kristina and Weiser, Martin and Conrad, Tim}, title = {Hybrid PDE-ODE Models for Efficient Simulation of Infection Spread in Epidemiology}, volume = {481}, journal = {Proceedings of the Royal Society A}, number = {2306}, publisher = {Royal Society}, arxiv = {http://arxiv.org/abs/2405.12938}, doi = {10.1098/rspa.2024.0421}, year = {2025}, abstract = {This paper introduces a novel hybrid model combining Partial Differential Equations (PDEs) and Ordinary Differential Equations (ODEs) to simulate infectious disease dynamics across geographic regions. By leveraging the spatial detail of PDEs and the computational efficiency of ODEs, the model enables rapid evaluation of public health interventions. Applied to synthetic environments and real-world scenarios in Lombardy, Italy, and Berlin, Germany, the model highlights how interactions between PDE and ODE regions affect infection dynamics, especially in high-density areas. Key findings reveal that the placement of model boundaries in densely populated regions can lead to inaccuracies in infection spread, suggesting that boundaries should be positioned in areas of lower population density to better reflect transmission dynamics. Additionally, regions with low population density hinder infection flow, indicating a need for incorporating, e.g., jumps in the model to enhance its predictive capabilities. Results indicate that the hybrid model achieves a balance between computational speed and accuracy, making it a valuable tool for policymakers in real-time decision-making and scenario analysis in epidemiology and potentially in other fields requiring similar modeling approaches.}, language = {en} } @article{AndresArconesWeiserKoutsourelakisetal.2024, author = {Andr{\´e}s Arcones, Daniel and Weiser, Martin and Koutsourelakis, Faidon-Stelios and Unger, J{\"o}rg F.}, title = {Model bias identification for Bayesian calibration of stochastic digital twins of bridges}, volume = {41}, journal = {Applied Stochastic Models in Business and Industry}, number = {3}, doi = {10.1002/asmb.2897}, year = {2024}, abstract = {Simulation-based digital twins must provide accurate, robust and reliable digital representations of their physical counterparts. Quantifying the uncertainty in their predictions plays, therefore, a key role in making better-informed decisions that impact the actual system. The update of the simulation model based on data must be then carefully implemented. When applied to complex standing structures such as bridges, discrepancies between the computational model and the real system appear as model bias, which hinders the trustworthiness of the digital twin and increases its uncertainty. Classical Bayesian updating approaches aiming to infer the model parameters often fail at compensating for such model bias, leading to overconfident and unreliable predictions. In this paper, two alternative model bias identification approaches are evaluated in the context of their applicability to digital twins of bridges. A modularized version of Kennedy and O'Hagan's approach and another one based on Orthogonal Gaussian Processes are compared with the classical Bayesian inference framework in a set of representative benchmarks. Additionally, two novel extensions are proposed for such models: the inclusion of noise-aware kernels and the introduction of additional variables not present in the computational model through the bias term. The integration of such approaches in the digital twin corrects the predictions, quantifies their uncertainty, estimates noise from unknown physical sources of error and provides further insight into the system by including additional pre-existing information without modifying the computational model.}, language = {en} } @article{GoebelHuynhCheginietal.2025, author = {G{\"o}bel, Fritz and Huynh, Ngoc Mai Monica and Chegini, Fatemeh and Pavarino, Luca and Weiser, Martin and Scacchi, Simone and Anzt, Hartwig}, title = {A BDDC Preconditioner for the Cardiac EMI Model in three Dimensions}, journal = {SIAM J. Sci. Comput.}, arxiv = {http://arxiv.org/abs/2502.07722}, year = {2025}, abstract = {We analyze a Balancing Domain Decomposition by Constraints (BDDC) preconditioner for the solution of three dimensional composite Discontinuous Galerkin discretizations of reaction-diffusion systems of ordinary and partial differential equations arising in cardiac cell-by-cell models like the Extracellular space, Membrane and Intracellular space (EMI) Model. These microscopic models are essential for the understanding of events in aging and structurally diseased hearts which macroscopic models relying on homogenized descriptions of the cardiac tissue, like Monodomain and Bidomain models, fail to adequately represent. The modeling of each individual cardiac cell results in discontinuous global solutions across cell boundaries, requiring the careful construction of dual and primal spaces for the BDDC preconditioner. We provide a scalable condition number bound for the precondition operator and validate the theoretical results with extensive numerical experiments.}, language = {en} }