@article{AndresArconesWeiserKoutsourelakisetal.2026, author = {Andr{\´e}s Arcones, Daniel and Weiser, Martin and Koutsourelakis, Phaedon-Stelios and Unger, J{\"o}rg F.}, title = {Embedded Model Form Uncertainty Quantification with Measurement Noise for Bayesian Model Calibration}, volume = {7}, journal = {Data-Centric Engineering}, arxiv = {http://arxiv.org/abs/2410.12037}, doi = {10.1017/dce.2025.10035}, year = {2026}, abstract = {A key factor in ensuring the accuracy of computer simulations that model physical systems is the proper calibration of their parameters based on real-world observations or experimental data. Inevitably, uncertainties arise, and Bayesian methods provide a robust framework for quantifying and propagating these uncertainties to model predictions. Nevertheless, Bayesian methods paired with inexact models usually produce predictions unable to represent the observed datapoints. Additionally, the quantified uncertainties of these overconfident models cannot be propagated to other Quantities of Interest (QoIs) reliably. A promising solution involves embedding a model inadequacy term in the inference parameters, allowing the quantified model form uncertainty to influence non-observed QoIs. This paper introduces a more interpretable framework for embedding the model inadequacy compared to existing methods. To overcome the limitations of current approaches, we adapt the existing likelihood models to properly account for noise in the measurements and propose two new formulations designed to address their shortcomings. Moreover, we evaluate the performance of this inadequacy-embedding approach in the presence of discrepancies between measurements and model predictions, including noise and outliers. Particular attention is given to how the uncertainty associated with the model inadequacy term propagates to the QoIs, enabling a more comprehensive statistical analysis of prediction's reliability. Finally, the proposed approach is applied to estimate the uncertainty in the predicted heat flux from a transient thermal simulation using temperature bservations.}, language = {en} } @article{SubramaniamHubigSchenkletal.2026, author = {Subramaniam, Jayant Shanmugam and Hubig, Michael and Schenkl, Sebastian and Muggenthaler, Holger and Springer, Steffen and Weiser, Martin and Sudau, Jakob and Shah, Faisal and Mall, Gita}, title = {Reconstructing sudden ambient temperature changes for forensic death time estimation using temperatures in two closed compartments: Proof of concept}, journal = {International Journal of Legal Medicine}, year = {2026}, abstract = {Ambient temperature TA has a strong impact on temperature-based time since death estimation (TTDE). Frequently TA is lowered instantaneously at some time t0 from a previous value TA0 to TA1 < TA0 by, e.g., opening a window or a door. We aim at reconstructing TA0 and t0. TTDE literature suggests temperature measurements in closed compartments such as cupboards or neighboring rooms, where TA0 could have been 'preserved' after t0. We aim to estimate t0 and TA0 from temperature measurements TZ(t) in closed compartments Z at times t > t0. We obtain promising results assuming Newtonian cooling for boxes filled with air, heaps of clothes, or books in two different experimental scenarios. Two different parameter estimators, (TA0^, t0^) based on four temperature measurements and (TA0*, t0*) for 4N measurements were tested. Our results in a climate chamber were partially appropriate for TTDE input. A decline at time t0 from TA0 = 22.5°C ↓ TA1 = 14°C was reconstructed at t = t0 + 95min with relative deviations ρt0^ = 27\% and ρTA0^ = 19\% relative to t - t0 and TA0 - TA1 respectively, for N = 1 with span Δt = 50min. For N = 200 in a time interval [t0 + 95min, t0 + 295min] we found ρt0^ = 5\% and ρTA0^ = 11\% with the same Δt. Further research is necessary to guarantee applicability in routine casework, in particular with respect to more elaborate cooling models, estimation algorithms, and evaluation localization.}, 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{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} } @article{CardereraPokuttaSchuetteetal.2025, author = {Carderera, Alejandro and Pokutta, Sebastian and Sch{\"u}tte, Christof and Weiser, Martin}, title = {An efficient first-order conditional gradient algorithm in data-driven sparse identification of nonlinear dynamics to solve sparse recovery problems under noise}, volume = {470}, journal = {Journal of Computational and Applied Mathematics}, doi = {10.1016/j.cam.2025.116675}, year = {2025}, abstract = {Governing equations are essential to the study of nonlinear dynamics, often enabling the prediction of previously unseen behaviors as well as the inclusion into control strategies. The discovery of governing equations from data thus has the potential to transform data-rich fields where well-established dynamical models remain unknown. This work contributes to the recent trend in data-driven sparse identification of nonlinear dynamics of finding the best sparse fit to observational data in a large library of potential nonlinear models. We propose an efficient first-order Conditional Gradient algorithm for solving the underlying optimization problem. In comparison to the most prominent alternative framework, the new framework shows significantly improved performance on several essential issues like sparsity-induction, structure-preservation, noise robustness, and sample efficiency. We demonstrate these advantages on several dynamics from the field of synchronization, particle dynamics, and enzyme chemistry.}, 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{BecksLippoldWinkleretal.2024, author = {Becks, Henrik and Lippold, Lukas and Winkler, Paul and Rohrer, Maximilian and Leusmann, Thorsten and Anton, David and Sprenger, Bjarne and K{\"a}hler, Philipp and Rudenko, Iryna and Andr{\´e}s Arcones, Daniel and Koutsourelakis, Phaedon-Stelios and Unger, J{\"o}rg F. and Weiser, Martin and Petryna, Yuri and Schnellenbach-Held, Martina and Lowke, Dirk and Wessels, Henning and Lenzen, Armin and Zabel, Volkmar and K{\"o}nke, Carsten and Claßen, Martin and Hegger, Josef}, title = {Neuartige Konzepte f{\"u}r die Zustands{\"u}berwachung und -analyse von Br{\"u}ckenbauwerken - Einblicke in das Forschungsvorhaben SPP100+}, volume = {99}, journal = {Bauingenieur}, number = {10}, doi = {10.37544/0005-6650-2024-10-63}, pages = {327 -- 338}, year = {2024}, abstract = {Die Br{\"u}ckeninfrastruktur in Deutschland und Europa steht aufgrund steigender Verkehrslasten und alternder Bauwerke vor erheblichen Herausforderungen. Das DFG-Schwerpunktprogramm 2388 „Hundert plus - Verl{\"a}ngerung der Lebensdauer komplexer Baustrukturen durch intelligente Digitalisierung" (SPP100+) strebt an, durch digitale Innovationen und pr{\"a}diktive Instandhaltungsstrategien die Nutzungsdauer bestehender Br{\"u}ckenbauwerke zu verl{\"a}ngern. Der vorliegende Beitrag fokussiert sich auf das SPP100+ zugeh{\"o}rige Cluster „Monitoring und Simulation", das sieben Teilprojekte umfasst. Die Projekte entwickeln fortschrittliche Methoden zur {\"U}berwachung und Zustandsbewertung von Br{\"u}cken mittels Digitaler Zwillinge, hochaufl{\"o}sender Sensortechnik und numerischer Simulationen. Innovative Ans{\"a}tze wie nichtlineare Modellanpassungen, stochastische Methoden und k{\"u}nstliche Intelligenz erm{\"o}glichen eine pr{\"a}zise und fr{\"u}hzeitige Identifizierung potenzieller Sch{\"a}den. Die Kombination aus kontinuierlichem Bauwerksmonitoring und effizienter Datenauswertung ist entscheidend f{\"u}r die langfristige Sicherheit und Langlebigkeit bestehender Br{\"u}cken und tr{\"a}gt dar{\"u}ber hinaus zur Ressourcenschonung bei.}, language = {de} } @article{AndresArconesWeiserKoutsourelakisetal.2024, author = {Andr{\´e}s Arcones, Daniel and Weiser, Martin and Koutsourelakis, Phaedon-Stelios and Unger, J{\"o}rg F.}, title = {Bias Identification Approaches for Model Updating of Simulation-based Digital Twins of Bridges}, volume = {2}, journal = {Structural Health Monitoring in the Light of Climate Impact and Data Science. Research and Review Journal of Nondestructive Testing}, number = {2}, doi = {10.58286/30524}, year = {2024}, abstract = {Simulation-based digital twins of bridges have the potential not only to serve as monitoring devices of the current state of the structure but also to generate new knowledge through physical predictions that allow for better-informed decisionmaking. For an accurate representation of the bridge, the underlying models must be tuned to reproduce the real system. Nevertheless, the necessary assumptions and simplifications in these models irremediably introduce discrepancies between measurements and model response. We will show that quantifying the extent of the uncertainties introduced through the models that lead to such discrepancies provides a better understanding of the real system, enhances the model updating process, and creates more robust and trustworthy digital twins. The inclusion of an explicit bias term will be applied to a representative demonstrator case based on the thermal response of the Nibelungenbr{\"u}cke of Worms. The findings from this work are englobed in the initiative SPP 100+, whose main aim is the extension of the service life of structures, especially through the implementation of digital twins.}, language = {en} } @article{NiemannUramWolfetal.2024, author = {Niemann, Jan-Hendrik and Uram, Samuel and Wolf, Sarah and Conrad, Natasa Djurdjevac and Weiser, Martin}, title = {Multilevel Optimization for Policy Design with Agent-Based Epidemic Models}, volume = {77}, journal = {Computational Science}, arxiv = {http://arxiv.org/abs/2304.02281}, doi = {10.1016/j.jocs.2024.102242}, pages = {102242}, year = {2024}, abstract = {Epidemiological models can not only be used to forecast the course of a pandemic like COVID-19, but also to propose and design non-pharmaceutical interventions such as school and work closing. In general, the design of optimal policies leads to nonlinear optimization problems that can be solved by numerical algorithms. Epidemiological models come in different complexities, ranging from systems of simple ordinary differential equations (ODEs) to complex agent-based models (ABMs). The former allow a fast and straightforward optimization, but are limited in accuracy, detail, and parameterization, while the latter can resolve spreading processes in detail, but are extremely expensive to optimize. We consider policy optimization in a prototypical situation modeled as both ODE and ABM, review numerical optimization approaches, and propose a heterogeneous multilevel approach based on combining a fine-resolution ABM and a coarse ODE model. Numerical experiments, in particular with respect to convergence speed, are given for illustrative examples.}, 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} } @article{SubramaniamHubigMuggenthaleretal.2023, author = {Subramaniam, Jayant S. and Hubig, Michael and Muggenthaler, Holger and Schenkl, Sebastian and Ullrich, Julia and Pourtier, Gr{\´e}groire and Weiser, Martin and Mall, Gita}, title = {Sensitivity of temperature-based time since death estimation on measurement location}, volume = {137}, journal = {International Journal of Legal Medicine}, doi = {10.1007/s00414-023-03040-y}, pages = {1815 -- 1837}, year = {2023}, abstract = {Rectal temperature measurement (RTM) from crime scenes is an important parameter for temperature-based time of death estimation (TDE). Various influential variables exist in TDE methods like the uncertainty in thermal and environmental parameters. Although RTM depends in particular on the location of measurement position, this relationship has never been investigated separately. The presented study fills this gap using Finite Element (FE) simulations of body cooling. A manually meshed coarse human FE model and an FE geometry model developed from the CT scan of a male corpse are used for TDE sensitivity analysis. The coarse model is considered with and without a support structure of moist soil. As there is no clear definition of ideal rectal temperature measurement location for TDE, possible variations in RTM location (RTML) are considered based on anatomy and forensic practice. The maximum variation of TDE caused by RTML changes is investigated via FE simulation. Moreover, the influence of ambient temperature, of FE model change and of the models positioning on a wet soil underground are also discussed. As a general outcome, we notice that maximum TDE deviations of up to ca. 2-3 h due to RTML deviations have to be expected. The direction of maximum influence of RTML change on TDE generally was on the line caudal to cranial.}, language = {en} } @article{SemlerWeiser2023, author = {Semler, Phillip and Weiser, Martin}, title = {Adaptive Gaussian Process Regression for Efficient Building of Surrogate Models in Inverse Problems}, volume = {39}, journal = {Inverse Problems}, number = {12}, arxiv = {http://arxiv.org/abs/2303.05824}, doi = {10.1088/1361-6420/ad0028}, pages = {125003}, year = {2023}, abstract = {In a task where many similar inverse problems must be solved, evaluating costly simulations is impractical. Therefore, replacing the model y with a surrogate model y(s) that can be evaluated quickly leads to a significant speedup. The approximation quality of the surrogate model depends strongly on the number, position, and accuracy of the sample points. With an additional finite computational budget, this leads to a problem of (computer) experimental design. In contrast to the selection of sample points, the trade-off between accuracy and effort has hardly been studied systematically. We therefore propose an adaptive algorithm to find an optimal design in terms of position and accuracy. Pursuing a sequential design by incrementally appending the computational budget leads to a convex and constrained optimization problem. As a surrogate, we construct a Gaussian process regression model. We measure the global approximation error in terms of its impact on the accuracy of the identified parameter and aim for a uniform absolute tolerance, assuming that y(s) is computed by finite element calculations. A priori error estimates and a coarse estimate of computational effort relate the expected improvement of the surrogate model error to computational effort, resulting in the most efficient combination of sample point and evaluation tolerance. We also allow for improving the accuracy of already existing sample points by continuing previously truncated finite element solution procedures.}, language = {en} } @article{UllrichWeiserSubramaniametal.2023, author = {Ullrich, Julia and Weiser, Martin and Subramaniam, Jayant and Schenkl, Sebastian and Muggenthaler, Holger and Hubig, Michael and Mall, Gita}, title = {The impact of anatomy variation on temperature based time of death estimation}, volume = {137}, journal = {International Journal of Legal Medicine}, doi = {10.1007/s00414-023-03026-w}, pages = {1615 -- 1627}, year = {2023}, abstract = {Temperature-based time of death estimation (TTDE) using simulation methods such as the finite element (FE) method promises higher accuracy and broader applicability in nonstandard cooling scenarios than established phenomenological methods. Their accuracy depends crucially on the simulation model to capture the actual situation. The model fidelity in turn hinges on the representation of the corpse's anatomy in form of computational meshes as well as on the thermodynamic parameters. While inaccuracies in anatomy representation due to coarse mesh resolution are known to have a minor impact on the estimated time of death, the sensitivity with respect to larger differences in the anatomy has so far not been studied. We assess this sensitivity by comparing four independently generated and vastly different anatomical models in terms of the estimated time of death in an identical cooling scenario. In order to isolate the impact of shape variation, the models are scaled to a reference size, and the possible impact of measurement location variation is excluded explicitly, which gives a lower bound on the impact of anatomy on the estimated time of death.}, language = {en} } @article{HuynhCheginiPavarinoetal.2023, author = {Huynh, Ngoc and Chegini, Fatemeh and Pavarino, Luca and Weiser, Martin and Scacchi, Simone}, title = {Convergence analysis of BDDC preconditioners for hybrid DG discretizations of the cardiac cell-by-cell model}, volume = {45}, journal = {SIAM Journal on Scientific Computing}, number = {6}, arxiv = {http://arxiv.org/abs/2212.12295}, pages = {A2836 -- A2857}, year = {2023}, abstract = {A Balancing Domain Decomposition by Constraints (BDDC) preconditioner is constructed and analyzed for the solution of hybrid Discontinuous Galerkin discretizations of reaction-diffusion systems of ordinary and partial differential equations arising in cardiac cell-by-cell models. The latter are different from the classical Bidomain and Monodomain cardiac models based on homogenized descriptions of the cardiac tissue at the macroscopic level, and therefore they allow the representation of individual cardiac cells, cell aggregates, damaged tissues and nonuniform distributions of ion channels on the cell membrane. The resulting discrete cell-by-cell models have discontinuous global solutions across the cell boundaries, hence the proposed BDDC preconditioner is based on appropriate dual and primal spaces with additional constraints which transfer information between cells (subdomains) without influencing the overall discontinuity of the global solution. A scalable convergence rate bound is proved for the resulting BDDC cell-by-cell preconditioned operator, while numerical tests validate this bound and investigate its dependence on the discretization parameters.}, language = {en} } @article{BartelsFisikopoulosWeiser2023, author = {Bartels, Tinko and Fisikopoulos, Vissarion and Weiser, Martin}, title = {Fast Floating-Point Filters for Robust Predicates}, volume = {63}, journal = {BIT Numerical Mathematics}, arxiv = {http://arxiv.org/abs/2208.00497}, doi = {10.1007/s10543-023-00975-x}, year = {2023}, abstract = {Geometric predicates are at the core of many algorithms, such as the construction of Delaunay triangulations, mesh processing and spatial relation tests. These algorithms have applications in scientific computing, geographic information systems and computer-aided design. With floating-point arithmetic, these geometric predicates can incur round-off errors that may lead to incorrect results and inconsistencies, causing computations to fail. This issue has been addressed using a combination of exact arithmetic for robustness and floating-point filters to mitigate the computational cost of exact computations. The implementation of exact computations and floating-point filters can be a difficult task, and code generation tools have been proposed to address this. We present a new C++ meta-programming framework for the generation of fast, robust predicates for arbitrary geometric predicates based on polynomial expressions. We combine and extend different approaches to filtering, branch reduction, and overflow avoidance that have previously been proposed. We show examples of how this approach produces correct results for data sets that could lead to incorrect predicate results with naive implementations. Our benchmark results demonstrate that our implementation surpasses state-of-the-art implementations.}, language = {en} } @article{BorndoerferDaneckerWeiser2023, author = {Bornd{\"o}rfer, Ralf and Danecker, Fabian and Weiser, Martin}, title = {Newton's Method for Global Free Flight Trajectory Optimization}, volume = {4}, journal = {Operations Research Forum}, doi = {10.1007/s43069-023-00238-z}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-91846}, year = {2023}, abstract = {Globally optimal free flight trajectory optimization can be achieved with a combination of discrete and continuous optimization. A key requirement is that Newton's method for continuous optimization converges in a sufficiently large neighborhood around a minimizer. We show in this paper that, under certain assumptions, this is the case.}, language = {en} } @article{BorndoerferDaneckerWeiser2023, author = {Bornd{\"o}rfer, Ralf and Danecker, Fabian and Weiser, Martin}, title = {Error Bounds for Discrete-Continuous Free Flight Trajectory Optimization}, volume = {198}, journal = {Journal of Optimization Theory and Applications}, doi = {10.1007/s10957-023-02264-7}, pages = {830 -- 856}, year = {2023}, abstract = {Flight planning, the computation of optimal routes in view of flight time and fuel consumption under given weather conditions, is traditionally done by finding globally shortest paths in a predefined airway network. Free flight trajectories, not restricted to a network, have the potential to reduce the costs significantly, and can be computed using locally convergent continuous optimal control methods. Hybrid methods that start with a discrete global search and refine with a fast continuous local optimization combine the best properties of both approaches, but rely on a good switchover, which requires error estimates for discrete paths relative to continuous trajectories. Based on vertex density and local complete connectivity, we derive localized and a priori bounds for the flight time of discrete paths relative to the optimal continuous trajectory, and illustrate their properties on a set of benchmark problems. It turns out that localization improves the error bound by four orders of magnitude, but still leaves ample opportunities for tighter bounds using a posteriori error estimators.}, language = {en} } @article{AlhaddadFoerstnerGrothetal.2022, author = {Alhaddad, Samer and F{\"o}rstner, Jens and Groth, Stefan and Gr{\"u}newald, Daniel and Grynko, Yevgen and Hannig, Frank and Kenter, Tobias and Pfreundt, F.J. and Plessl, Christian and Schotte, Merlind and Steinke, Thomas and Teich, J. and Weiser, Martin and Wende, Florian}, title = {The HighPerMeshes Framework for Numerical Algorithms on Unstructured Grids}, volume = {34}, journal = {Concurrency and Computation: Practice and Experience}, number = {14}, doi = {10.1002/cpe.6616}, year = {2022}, abstract = {Solving PDEs on unstructured grids is a cornerstone of engineering and scientific computing. Heterogeneous parallel platforms, including CPUs, GPUs, and FPGAs, enable energy-efficient and computationally demanding simulations. In this article, we introduce the HPM C++-embedded DSL that bridges the abstraction gap between the mathematical formulation of mesh-based algorithms for PDE problems on the one hand and an increasing number of heterogeneous platforms with their different programming models on the other hand. Thus, the HPM DSL aims at higher productivity in the code development process for multiple target platforms. We introduce the concepts as well as the basic structure of the HPM DSL, and demonstrate its usage with three examples. The mapping of the abstract algorithmic description onto parallel hardware, including distributed memory compute clusters, is presented. A code generator and a matching back end allow the acceleration of HPM code with GPUs. Finally, the achievable performance and scalability are demonstrated for different example problems.}, language = {en} } @article{GoetschelSchielaWeiser2021, author = {G{\"o}tschel, Sebastian and Schiela, Anton and Weiser, Martin}, title = {Kaskade 7 - a Flexible Finite Element Toolbox}, volume = {81}, journal = {Computers and Mathematics with Applications}, doi = {10.1016/j.camwa.2020.02.011}, pages = {444 -- 458}, year = {2021}, 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 \textsc{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} } @article{AlhaddadFoerstnerGrothetal.2021, author = {Alhaddad, Samer and F{\"o}rstner, Jens and Groth, Stefan and Gr{\"u}newald, Daniel and Grynko, Yevgen and Hannig, Frank and Kenter, Tobias and Pfreundt, Franz-Josef and Plessl, Christian and Schotte, Merlind and Steinke, Thomas and Teich, J{\"u}rgen and Weiser, Martin and Wende, Florian}, title = {HighPerMeshes - A Domain-Specific Language for Numerical Algorithms on Unstructured Grids}, journal = {Euro-Par 2020: Parallel Processing Workshops.}, publisher = {Springer}, doi = {10.1007/978-3-030-71593-9_15}, pages = {185 -- 196}, year = {2021}, abstract = {Solving partial differential equations on unstructured grids is a cornerstone of engineering and scientific computing. Nowadays, heterogeneous parallel platforms with CPUs, GPUs, and FPGAs enable energy-efficient and computationally demanding simulations. We developed the HighPerMeshes C++-embedded Domain-Specific Language (DSL) for bridging the abstraction gap between the mathematical and algorithmic formulation of mesh-based algorithms for PDE problems on the one hand and an increasing number of heterogeneous platforms with their different parallel programming and runtime models on the other hand. Thus, the HighPerMeshes DSL aims at higher productivity in the code development process for multiple target platforms. We introduce the concepts as well as the basic structure of the HighPer-Meshes DSL, and demonstrate its usage with three examples, a Poisson and monodomain problem, respectively, solved by the continuous finite element method, and the discontinuous Galerkin method for Maxwell's equation. The mapping of the abstract algorithmic description onto parallel hardware, including distributed memory compute clusters is presented. Finally, the achievable performance and scalability are demonstrated for a typical example problem on a multi-core CPU cluster.}, language = {en} } @article{SchneckWeiserWende2021, author = {Schneck, Jakob and Weiser, Martin and Wende, Florian}, title = {Impact of mixed precision and storage layout on additive Schwarz smoothers}, volume = {28}, journal = {Numerical Linear Algebra with Applications}, number = {4}, doi = {10.1002/nla.2366}, year = {2021}, abstract = {The growing discrepancy between CPU computing power and memory bandwidth drives more and more numerical algorithms into a bandwidth-bound regime. One example is the overlapping Schwarz smoother, a highly effective building block for iterative multigrid solution of elliptic equations with higher order finite elements. Two options of reducing the required memory bandwidth are sparsity exploiting storage layouts and representing matrix entries with reduced precision in floating point or fixed point format. We investigate the impact of several options on storage demand and contraction rate, both analytically in the context of subspace correction methods and numerically at an example of solid mechanics. Both perspectives agree on the favourite scheme: fixed point representation of Cholesky factors in nested dissection storage.}, language = {en} } @article{CardereraPokuttaSchuetteetal.2021, author = {Carderera, Alejandro and Pokutta, Sebastian and Sch{\"u}tte, Christof and Weiser, Martin}, title = {An efficient first-order conditional gradient algorithm in data-driven sparse identification of nonlinear dynamics to solve sparse recovery problems under noise}, journal = {Journal of Computational and Applied Mathematics}, arxiv = {http://arxiv.org/abs/2101.02630}, year = {2021}, abstract = {Governing equations are essential to the study of nonlinear dynamics, often enabling the prediction of previously unseen behaviors as well as the inclusion into control strategies. The discovery of governing equations from data thus has the potential to transform data-rich fields where well-established dynamical models remain unknown. This work contributes to the recent trend in data-driven sparse identification of nonlinear dynamics of finding the best sparse fit to observational data in a large library of potential nonlinear models. We propose an efficient first-order Conditional Gradient algorithm for solving the underlying optimization problem. In comparison to the most prominent alternative algorithms, the new algorithm shows significantly improved performance on several essential issues like sparsity-induction, structure-preservation, noise robustness, and sample efficiency. We demonstrate these advantages on several dynamics from the field of synchronization, particle dynamics, and enzyme chemistry.}, language = {en} } @article{CheginiKopanicakovaKrauseetal.2021, author = {Chegini, Fatemeh and Kopanicakova, Alena and Krause, Rolf and Weiser, Martin}, title = {Efficient Identification of Scars using Heterogeneous Model Hierarchies}, volume = {23}, journal = {EP Europace}, doi = {10.1093/europace/euaa402}, pages = {i113 -- i122}, year = {2021}, abstract = {Aims. Detection and quantification of myocardial scars are helpful both for diagnosis of heart diseases and for building personalized simulation models. Scar tissue is generally charac­terized by a different conduction of electrical excitation. We aim at estimating conductivity-related parameters from endocardial mapping data, in particular the conductivity tensor. Solving this inverse problem requires computationally expensive monodomain simulations on fine discretizations. Therefore, we aim at accelerating the estimation using a multilevel method combining electrophysiology models of different complexity, namely the mono­domain and the eikonal model. Methods. Distributed parameter estimation is performed by minimizing the misfit between simulated and measured electrical activity on the endocardial surface, subject to the mono­domain model and regularization, leading to a constrained optimization problem. We formulate this optimization problem, including the modeling of scar tissue and different regularizations, and design an efficient iterative solver. We consider monodomain grid hierarchies and monodomain-eikonal model hierarchies in a recursive multilevel trust-region method. Results. From several numerical examples, both the efficiency of the method and the estimation quality, depending on the data, are investigated. The multilevel solver is significantly faster than a comparable single level solver. Endocardial mapping data of realistic density appears to be just sufficient to provide quantitatively reasonable estimates of location, size, and shape of scars close to the endocardial surface. Conclusion. In several situations, scar reconstruction based on eikonal and monodomain models differ significantly, suggesting the use of the more accurate but more expensive monodomain model for this purpose. Still, eikonal models can be utilized to accelerate the computations considerably, enabling the use of complex electrophysiology models for estimating myocardial scars from endocardial mapping data.}, language = {en} } @article{SchielaStoeckleinWeiser2021, author = {Schiela, Anton and St{\"o}cklein, Matthias and Weiser, Martin}, title = {A primal dual projection algorithm for efficient constraint preconditioning}, volume = {43}, journal = {SIAM Journal on Scientific Computing}, number = {6}, doi = {10.1137/20M1380739}, pages = {A4095 -- A4120}, year = {2021}, abstract = {We consider a linear iterative solver for large scale linearly constrained quadratic minimization problems that arise, for example, in optimization with PDEs. By a primal-dual projection (PDP) iteration, which can be interpreted and analysed as a gradient method on a quotient space, the given problem can be solved by computing sulutions for a sequence of constrained surrogate problems, projections onto the feasible subspaces, and Lagrange multiplier updates. As a major application we consider a class of optimization problems with PDEs, where PDP can be applied together with a projected cg method using a block triangular constraint preconditioner. Numerical experiments show reliable and competitive performance for an optimal control problem in elasticity.}, language = {en} } @article{BorndoerferDaneckerWeiser2021, author = {Bornd{\"o}rfer, Ralf and Danecker, Fabian and Weiser, Martin}, title = {A Discrete-Continuous Algorithm for Free Flight Planning}, volume = {14}, journal = {Algorithms}, number = {1}, publisher = {MDPI}, issn = {1438-0064}, doi = {10.3390/a14010004}, pages = {4}, year = {2021}, abstract = {We propose a hybrid discrete-continuous algorithm for flight planning in free flight airspaces. In a first step, our DisCOptER method discrete-continuous optimization for enhanced resolution) computes a globally optimal approximate flight path on a discretization of the problem using the A* method. This route initializes a Newton method that converges rapidly to the smooth optimum in a second step. The correctness, accuracy, and complexity of the method are goverened by the choice of the crossover point that determines the coarseness of the discretization. We analyze the optimal choice of the crossover point and demonstrate the asymtotic superority of DisCOptER over a purely discrete approach.}, language = {en} } @article{OzelKulkarniHasanetal.2019, author = {Ozel, Mehmet Neset and Kulkarni, Abhishek and Hasan, Amr and Brummer, Josephine and Moldenhauer, Marian and Daumann, Ilsa-Maria and Wolfenberg, Heike and Dercksen, Vincent J. and Kiral, Ferdi Ridvan and Weiser, Martin and Prohaska, Steffen and von Kleist, Max and Hiesinger, Peter Robin}, title = {Serial synapse formation through filopodial competition for synaptic seeding factors}, volume = {50}, journal = {Developmental Cell}, number = {4}, doi = {10.1016/j.devcel.2019.06.014}, pages = {447 -- 461}, year = {2019}, abstract = {Following axon pathfinding, growth cones transition from stochastic filopodial exploration to the formation of a limited number of synapses. How the interplay of filopodia and synapse assembly ensures robust connectivity in the brain has remained a challenging problem. Here, we developed a new 4D analysis method for filopodial dynamics and a data-driven computational model of synapse formation for R7 photoreceptor axons in developing Drosophila brains. Our live data support a 'serial synapse formation' model, where at any time point only a single 'synaptogenic' filopodium suppresses the synaptic competence of other filopodia through competition for synaptic seeding factors. Loss of the synaptic seeding factors Syd-1 and Liprin-α leads to a loss of this suppression, filopodial destabilization and reduced synapse formation, which is sufficient to cause the destabilization of entire axon terminals. Our model provides a filopodial 'winner-takes-all' mechanism that ensures the formation of an appropriate number of synapses.}, language = {en} } @article{GoetschelWeiser2019, author = {G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Compression Challenges in Large Scale Partial Differential Equation Solvers}, volume = {12}, journal = {Algorithms}, number = {9}, doi = {10.3390/a12090197}, pages = {197}, year = {2019}, abstract = {Solvers for partial differential equations (PDEs) are one of the cornerstones of computational science. For large problems, they involve huge amounts of data that need to be stored and transmitted on all levels of the memory hierarchy. Often, bandwidth is the limiting factor due to the relatively small arithmetic intensity, and increasingly 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 over the recent decades. This paper surveys data compression challenges and discusses examples of corresponding solution approaches for PDE problems, covering all levels of the memory hierarchy from mass storage up to the main memory. We illustrate concepts for particular methods, with examples, and give references to alternatives.}, language = {en} } @article{WeiserGhosh2018, author = {Weiser, Martin and Ghosh, Sunayana}, title = {Theoretically optimal inexact SDC methods}, journal = {Communications in Applied Mathematics and Computational Science}, number = {13-1}, doi = {10.2140/camcos.2018.13.53}, pages = {53 -- 86}, year = {2018}, abstract = {In several inital value problems with particularly expensive right hand side evaluation or implicit step computation, there is a trade-off between accuracy and computational effort. We consider inexact spectral deferred correction (SDC) methods for solving such 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 error models bounding the total error in terms of the 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 local tolerance selection is worked out by minimizing the total work subject to achieving the requested tolerance. The properties of optimal local tolerances and the predicted efficiency gain compared to simpler heuristics, and a reasonable practical performance, are illustrated on simple numerical examples.}, language = {en} } @article{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}, volume = {34}, journal = {Inverse Problems}, number = {12}, doi = {10.1088/1361-6420/aae7a5}, pages = {125005}, 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} } @article{HammerschmidtSchneiderSantiagoetal.2018, author = {Hammerschmidt, Martin and Schneider, Philipp-Immanuel and Santiago, Xavier Garcia and Zschiedrich, Lin and Weiser, Martin and Burger, Sven}, title = {Solving inverse problems appearing in design and metrology of diffractive optical elements by using Bayesian optimization}, volume = {10694}, journal = {Proc. SPIE}, doi = {10.1117/12.2315468}, pages = {1069407}, year = {2018}, language = {en} } @article{TierneyMiraReinholdetal.2018, author = {Tierney, Nicholas J. and Mira, Antonietta and Reinhold, J. Jost and Weiser, Martin and Burkart, Roman and Benvenuti, Claudio and Auricchio, Angelo}, title = {Novel relocation methods for automatic external defibrillator improve out-of-hospital cardiac arrest coverage under limited resources}, journal = {Resuscitation}, number = {125}, doi = {10.1016/j.resuscitation.2018.01.055}, pages = {83 -- 89}, year = {2018}, abstract = {Background Mathematical optimisation models have recently been applied to identify ideal Automatic External Defibrillator (AED) locations that maximise coverage of Out of Hospital Cardiac Arrest (OHCA). However, these fixed location models cannot relocate existing AEDs in a flexible way, and have nearly exclusively been applied to urban regions. We developed a flexible location model for AEDs, compared its performance to existing fixed location and population models, and explored how these perform across urban and rural regions. Methods Optimisation techniques were applied to AED deployment and OHCA coverage was assessed. A total of 2802 geolocated OHCAs occurred in Canton Ticino, Switzerland, from January 1st 2005 to December 31st 2015. Results There were 719 AEDs in Canton Ticino. 635 (23\%) OHCA events occurred within 100m of an AED, with 306 (31\%) in urban, and 329 (18\%) in rural areas. Median distance from OHCA events to the nearest AED was 224m (168m urban vs. 269m rural). Flexible location models performed better than fixed location and population models, with the cost to deploy 20 new AEDs instead relocating 171 existing AEDs to new locations, improving OHCA coverage to 38\%, compared to 26\% using fixed models, and 24\% with the population based model. Conclusions Optimisation models for AEDs placement are superior to population models and should be strongly considered by communities when selecting areas for AED deployment. Compared to other models, flexible location models increase overall OHCA coverage, and decreases the distance to nearby AEDs, even in rural areas, while saving significant financial resources.}, 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} } @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} } @article{SchenklMuggenthalerHubigetal.2017, author = {Schenkl, Sebastian and Muggenthaler, Holger and Hubig, Michael and Erdmann, Bodo and Weiser, Martin and Zachow, Stefan and Heinrich, Andreas and G{\"u}ttler, Felix Victor and Teichgr{\"a}ber, Ulf and Mall, Gita}, title = {Automatic CT-based finite element model generation for temperature-based death time estimation: feasibility study and sensitivity analysis}, volume = {131}, journal = {International Journal of Legal Medicine}, number = {3}, doi = {doi:10.1007/s00414-016-1523-0}, pages = {699 -- 712}, year = {2017}, abstract = {Temperature based death time estimation is based either on simple phenomenological models of corpse cooling or on detailed physical heat transfer models. The latter are much more complex, but allow a higher accuracy of death time estimation as in principle all relevant cooling mechanisms can be taken into account. Here, a complete work flow for finite element based cooling simulation models is presented. The following steps are demonstrated on CT-phantoms: • CT-scan • Segmentation of the CT images for thermodynamically relevant features of individual geometries • Conversion of the segmentation result into a Finite Element (FE) simulation model • Computation of the model cooling curve • Calculation of the cooling time For the first time in FE-based cooling time estimation the steps from the CT image over segmentation to FE model generation are semi-automatically performed. The cooling time calculation results are compared to cooling measurements performed on the phantoms under controlled conditions. In this context, the method is validated using different CTphantoms. Some of the CT phantoms thermodynamic material parameters had to be experimentally determined via independent experiments. Moreover the impact of geometry and material parameter uncertainties on the estimated cooling time is investigated by a sensitivity analysis.}, 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} } @article{HammerschmidtWeiserSantiagoetal.2017, author = {Hammerschmidt, Martin and Weiser, Martin and Santiago, Xavier Garcia and Zschiedrich, Lin and Bodermann, Bernd and Burger, Sven}, title = {Quantifying parameter uncertainties in optical scatterometry using Bayesian inversion}, volume = {10330}, journal = {Proc. SPIE}, arxiv = {http://arxiv.org/abs/1707.08467}, doi = {10.1117/12.2270596}, pages = {1033004}, year = {2017}, 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} } @article{GoetschelWeiser2015, author = {G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Lossy Compression for PDE-constrained Optimization: Adaptive Error Control}, volume = {62}, journal = {Comput. Optim. Appl.}, number = {1}, publisher = {Springer}, pages = {131 -- 155}, year = {2015}, abstract = {For the solution of optimal control problems governed by nonlinear parabolic PDEs, 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 ad-joint 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 which are accurate enough, in many case very fine discretizations in time and space are necessary, which leads to a significant amount of data to be stored and transmitted to mass storage. Lossy compression methods were developed to overcome the storage problem by reducing the accuracy of the stored trajectories. The inexact data induces errors in the reduced gradient and reduced Hessian. In this paper, we analyze the influence of such a lossy trajectory compression method on Newton-CG methods for optimal control of parabolic PDEs and design an adaptive strategy for choosing appropriate quantization tolerances.}, 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} } @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{GoetschelNagaiahKunischetal.2014, author = {G{\"o}tschel, Sebastian and Nagaiah, Chamakuri and Kunisch, Karl and Weiser, Martin}, title = {Lossy Compression in Optimal Control of Cardiac Defibrillation}, volume = {60}, journal = {J. Sci. Comput.}, number = {1}, doi = {10.1007/s10915-013-9785-x}, pages = {35 -- 59}, year = {2014}, abstract = {This paper presents efficient computational techniques for solving an optimization problem in cardiac defibrillation governed by the monodomain equations. Time-dependent electrical currents injected at different spatial positions act as the control. Inexact Newton-CG methods are used, with reduced gradient computation by adjoint solves. In order to reduce the computational complexity, adaptive mesh refinement for state and adjoint equations is performed. To reduce the high storage and bandwidth demand imposed by adjoint gradient and Hessian-vector evaluations, a lossy compression technique for storing trajectory data is applied. An adaptive choice of quantization tolerance based on error estimates is developed in order to ensure convergence. The efficiency of the proposed approach is demonstrated on numerical examples.}, 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{GoetschelWeiserMaierhoferetal.2013, author = {G{\"o}tschel, Sebastian and Weiser, Martin and Maierhofer, Christiane and Richter, Regina and R{\"o}llig, Mathias}, title = {Fast Defect Shape Reconstruction Based on the Travel Time in Pulse Thermography}, volume = {6}, journal = {Nondestructive Testing of Materials and Structures}, pages = {83 -- 89}, year = {2013}, 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{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} } @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} } @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} } @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} } @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} } @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} }