@misc{WeiserChegini2022, author = {Weiser, Martin and Chegini, Fatemeh}, title = {Higher-order time integration using spectral deferred correction method (SDC) in a cell by cell discretization of cardiac excitation}, doi = {10.35097/716}, year = {2022}, abstract = {This C++ code implements a cell-by-cell model of cardiac excitation using a piecewise-continuous finite element discretization and spectral deferred correction time stepping. The code is based on the Kaskade 7 finite element toolbox and forms a prototype for the µCarp code to be implemented in the Microcard project.}, 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} } @inproceedings{GanderKrauseWeiseretal.2023, author = {Gander, Lia and Krause, Rolf and Weiser, Martin and Costabal, Francisco and Pezzuto, Simone}, title = {On the Accuracy of Eikonal Approximations in Cardiac Electrophysiology in the Presence of Fibrosis}, volume = {13958}, booktitle = {Functional Imaging and Modeling of the Heart. FIMH 2023.}, publisher = {Springer, Cham}, doi = {10.1007/978-3-031-35302-4_14}, year = {2023}, abstract = {Fibrotic tissue is one of the main risk factors for cardiac arrhythmias. It is therefore a key component in computational studies. In this work, we compare the monodomain equation to two eikonal models for cardiac electrophysiology in the presence of fibrosis. We show that discontinuities in the conductivity field, due to the presence of fibrosis, introduce a delay in the activation times. The monodomain equation and eikonal-diffusion model correctly capture these delays, contrarily to the classical eikonal equation. Importantly, a coarse space discretization of the monodomain equation amplifies these delays, even after accounting for numerical error in conduction velocity. The numerical discretization may also introduce artificial conduction blocks and hence increase propagation complexity. Therefore, some care is required when comparing eikonal models to the discretized monodomain equation.}, 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} } @inproceedings{CheginiSteinkeWeiser2022, author = {Chegini, Fatemeh and Steinke, Thomas and Weiser, Martin}, title = {Efficient adaptivity for simulating cardiac electrophysiology with spectral deferred correction methods}, arxiv = {http://arxiv.org/abs/2311.07206}, year = {2022}, abstract = {The locality of solution features in cardiac electrophysiology simulations calls for adaptive methods. Due to the overhead incurred by established mesh refinement and coarsening, however, such approaches failed in accelerating the computations. Here we investigate a different route to spatial adaptivity that is based on nested subset selection for algebraic degrees of freedom in spectral deferred correction methods. This combination of algebraic adaptivity and iterative solvers for higher order collocation time stepping realizes a multirate integration with minimal overhead. This leads to moderate but significant speedups in both monodomain and cell-by-cell models of cardiac excitation, as demonstrated at four numerical examples.}, 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{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} }