@misc{Beck1999, author = {Beck, Rudolf}, title = {Algebraic Multigrid by Component Splitting for Edge Elements on Simplicial Triangulations}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4290}, number = {SC-99-40}, year = {1999}, abstract = {Our focus is on Maxwell's equations in the low frequency range; two specific applications we aim at are time-stepping schemes for eddy current computations and the stationary double-curl equation for time-harmonic fields. We assume that the computational domain is discretized by triangles or tetrahedrons; for the finite element approximation we choose N\'{e}d\'{e}lec's \$H(curl)\$-conforming edge elements of the lowest order. For the solution of the arising linear equation systems we devise an algebraic multigrid preconditioner based on a spatial component splitting of the field. Mesh coarsening takes place in an auxiliary subspace, which is constructed with the aid of a nodal vector basis. Within this subspace coarse grids are created by exploiting the matrix graphs. Additionally, we have to cope with the kernel of the \$curl\$-operator, which comprises a considerable part of the spectral modes on the grid. Fortunately, the kernel modes are accessible via a discrete Helmholtz decomposition of the fields; they are smoothed by additional algebraic multigrid cycles. Numerical experiments are included in order to assess the efficacy of the proposed algorithms.}, language = {en} } @misc{HammerschmidtHerrmannPomplunetal.2015, author = {Hammerschmidt, Martin and Herrmann, Sven and Pomplun, Jan and Zschiedrich, Lin and Burger, Sven and Schmidt, Frank}, title = {Reduced basis method for Maxwell's equations with resonance phenomena}, volume = {9630}, journal = {Proc. SPIE}, issn = {1438-0064}, doi = {10.1117/12.2190425}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55687}, pages = {96300R}, year = {2015}, abstract = {Rigorous optical simulations of 3-dimensional nano-photonic structures are an important tool in the analysis and optimization of scattering properties of nano-photonic devices or parameter reconstruction. To construct geometrically accurate models of complex structured nano-photonic devices the finite element method (FEM) is ideally suited due to its flexibility in the geometrical modeling and superior convergence properties. Reduced order models such as the reduced basis method (RBM) allow to construct self-adaptive, error-controlled, very low dimensional approximations for input-output relationships which can be evaluated orders of magnitude faster than the full model. This is advantageous in applications requiring the solution of Maxwell's equations for multiple parameters or a single parameter but in real time. We present a reduced basis method for 3D Maxwell's equations based on the finite element method which allows variations of geometric as well as material and frequency parameters. We demonstrate accuracy and efficiency of the method for a light scattering problem exhibiting a resonance in the electric field.}, language = {en} } @misc{HammerschmidtPomplunBurgeretal.2014, author = {Hammerschmidt, Martin and Pomplun, Jan and Burger, Sven and Schmidt, Frank}, title = {Adaptive sampling strategies for efficient parameter scans in nano-photonic device simulations}, issn = {1438-0064}, doi = {10.1117/12.2036363}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-50395}, year = {2014}, abstract = {Rigorous optical simulations are an important tool in optimizing scattering properties of nano-photonic devices and are used, for example, in solar cell optimization. The finite element method (FEM) yields rigorous, time-harmonic, high accuracy solutions of the full 3D vectorial Maxwell's equations [1] and furthermore allows for great flexibility and accuracy in the geometrical modeling of these often complex shaped 3D nano-structures. A major drawback of frequency domain methods is the limitation of single frequency evaluations. For example the accurate computation of the short circuit current density of an amorphous silicon / micro-crystalline multi-junction thin film solar cell may require the solution of Maxwell's equations for over a hundred different wavelengths if an equidistant sampling strategy is employed. Also in optical metrology, wavelength scans are frequently used to reconstruct unknown geometrical and material properties of optical systems numerically from measured scatterometric data. In our contribution we present several adaptive numerical integration and sampling routines and study their efficiency in the context of the determination of generation rate profiles of solar cells. We show that these strategies lead to a reduction in the computational effort without loss of accuracy. We discuss the employment of tangential information in a Hermite interpolation scheme to achieve similar accuracy on coarser grids. We explore the usability of these strategies for scatterometry and solar cell simulations.}, language = {en} } @misc{HammerschmidtLockauZschiedrichetal.2014, author = {Hammerschmidt, Martin and Lockau, Daniel and Zschiedrich, Lin and Schmidt, Frank}, title = {Optical modelling of incoherent substrate light-trapping in silicon thin film multi-junction solar cells with finite elements and domain decomposition}, issn = {1438-0064}, doi = {10.1117/12.2036346}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-50410}, year = {2014}, abstract = {In many experimentally realized applications, e.g. photonic crystals, solar cells and light-emitting diodes, nano-photonic systems are coupled to a thick substrate layer, which in certain cases has to be included as a part of the optical system. The finite element method (FEM) yields rigorous, high accuracy solutions of full 3D vectorial Maxwell's equations [1] and allows for great flexibility and accuracy in the geometrical modelling. Time-harmonic FEM solvers have been combined with Fourier methods in domain decomposition algorithms to compute coherent solutions of these coupled system. [2,3] The basic idea of a domain decomposition approach lies in a decomposition of the domain into smaller subdomains, separate calculations of the solutions and coupling of these solutions on adjacent subdomains. In experiments light sources are often not perfectly monochromatic and hence a comparision to simulation results might only be justified if the simulation results, which include interference patterns in the substrate, are spectrally averaged. In this contribution we present a scattering matrix domain decomposition algorithm for Maxwell's equations based on FEM. We study its convergence and advantages in the context of optical simulations of silicon thin film multi-junction solar cells. This allows for substrate light-trapping to be included in optical simulations and leads to a more realistic estimation of light path enhancement factors in thin-film devices near the band edge.}, language = {en} } @misc{HammerschmidtHerrmannBurgeretal.2016, author = {Hammerschmidt, Martin and Herrmann, Sven and Burger, Sven and Pomplun, Jan and Schmidt, Frank}, title = {Reduced basis method for the optimization of nano-photonic devices}, issn = {1438-0064}, doi = {10.1007/s11082-016-0530-1}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57556}, year = {2016}, abstract = {Optical 3D simulations in many-query and real-time contexts require new solution strategies. We study an adaptive, error controlled reduced basis method for solving parametrized time-harmonic optical scattering problems. Application fields are, among others, design and optimization problems of nano-optical devices as well as inverse problems for parameter reconstructions occuring e. g. in optical metrology. The reduced basis method presented here relies on a finite element modeling of the scattering problem with parametrization of materials, geometries and sources.}, language = {en} } @misc{HammerschmidtHerrmannPomplunetal.2016, author = {Hammerschmidt, Martin and Herrmann, Sven and Pomplun, Jan and Burger, Sven and Schmidt, Frank}, title = {Model order reduction for the time-harmonic Maxwell equation applied to complex nanostructures}, volume = {9742}, journal = {Proc. SPIE}, issn = {1438-0064}, doi = {10.1117/12.2212367}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57239}, pages = {97420M}, year = {2016}, abstract = {Fields such as optical metrology and computational lithography require fast and efficient methods for solving the time-harmonic Maxwell's equation. Highly accurate geometrical modeling and numerical accuracy atcomputational costs are a prerequisite for any simulation study of complex nano-structured photonic devices. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem based on the hp-adaptive finite element solver JCMsuite capable of handling geometric and non-geometric parameter dependencies allowing for online evaluations in milliseconds. We apply the RBM to compute light-scatteringoptical wavelengths off periodic arrays of fin field-effect transistors (FinFETs) where geometrical properties such as the width and height of the fin and gate can vary in a large range.}, language = {en} } @misc{HammerschmidtBarthPomplunetal.2016, author = {Hammerschmidt, Martin and Barth, Carlo and Pomplun, Jan and Burger, Sven and Becker, Christiane and Schmidt, Frank}, title = {Reconstruction of photonic crystal geometries using a reduced basis method for nonlinear outputs}, issn = {1438-0064}, doi = {10.1117/12.2212482}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57249}, year = {2016}, abstract = {Maxwell solvers based on the hp-adaptive finite element method allow for accurate geometrical modeling and high numerical accuracy. These features are indispensable for the optimization of optical properties or reconstruction of parameters through inverse processes. High computational complexity prohibits the evaluation of the solution for many parameters. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem allowing to compute solutions for a parameter configuration orders of magnitude faster. The RBM allows to evaluate linear and nonlinear outputs of interest like Fourier transform or the enhancement of the electromagnetic field in milliseconds. We apply the RBM to compute light-scattering off two dimensional photonic crystal structures made of silicon and reconstruct geometrical parameters.}, language = {en} } @misc{KlebanovSikorskiSchuetteetal.2016, author = {Klebanov, Ilja and Sikorski, Alexander and Sch{\"u}tte, Christof and R{\"o}blitz, Susanna}, title = {Prior estimation and Bayesian inference from large cohort data sets}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57475}, year = {2016}, abstract = {One of the main goals of mathematical modelling in systems biology related to medical applications is to obtain patient-specific parameterisations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Using these population data, we propose an iterative algorithm for contructing an informative prior distribution, which then serves as the basis for computing patient-specific posteriors and obtaining individual predictions. We demonsrate the performance of our method by applying it to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling.}, language = {en} } @misc{LieSullivanTeckentrup2018, author = {Lie, Han Cheng and Sullivan, T. J. and Teckentrup, Aretha}, title = {Random forward models and log-likelihoods in Bayesian inverse problems}, volume = {6}, journal = {SIAM/ASA Journal on Uncertainty Quantification}, number = {4}, issn = {1438-0064}, arxiv = {http://arxiv.org/abs/1712.05717}, doi = {10.1137/18M1166523}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66324}, pages = {1600 -- 1629}, year = {2018}, abstract = {We consider the use of randomised forward models and log-likelihoods within the Bayesian approach to inverse problems. Such random approximations to the exact forward model or log-likelihood arise naturally when a computationally expensive model is approximated using a cheaper stochastic surrogate, as in Gaussian process emulation (kriging), or in the field of probabilistic numerical methods. We show that the Hellinger distance between the exact and approximate Bayesian posteriors is bounded by moments of the difference between the true and approximate log-likelihoods. Example applications of these stability results are given for randomised misfit models in large data applications and the probabilistic solution of ordinary differential equations.}, language = {en} } @misc{Beck1999, author = {Beck, Rudolf}, title = {Graph-Based Algebraic Multigrid for Lagrange-Type Finite Elements on Simplicial Meshes}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4100}, number = {SC-99-22}, year = {1999}, abstract = {We present an algebraic multigrid preconditioner which uses only the graphs of system matrices. Some elementary coarsening rules are stated, from which an advancing front algorithm for the selection of coarse grid nodes is derived. This technique can be applied to linear Lagrange-type finite element discretizations; for higher-order elements an extension of the multigrid algorithm is provided. Both two- and three-dimensional second order elliptic problems can be handled. Numerical experiments show that the resulting convergence acceleration is comparable to classical geometric multigrid.}, language = {en} } @misc{NavaYazdaniHegevonTycowicz2021, author = {Nava-Yazdani, Esfandiar and Hege, Hans-Christian and von Tycowicz, Christoph}, title = {A Hierarchical Geodesic Model for Longitudinal Analysis on Manifolds}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-85187}, year = {2021}, abstract = {In many applications, geodesic hierarchical models are adequate for the study of temporal observations. We employ such a model derived for manifold-valued data to Kendall's shape space. In particular, instead of the Sasaki metric, we adapt a functional-based metric, which increases the computational efficiency and does not require the implementation of the curvature tensor. We propose the corresponding variational time discretization of geodesics and employ the approach for longitudinal analysis of 2D rat skulls shapes as well as 3D shapes derived from an imaging study on osteoarthritis. Particularly, we perform hypothesis test and estimate the mean trends.}, language = {en} } @misc{RiberaBorrellQuerRichteretal.2021, author = {Ribera Borrell, Enric and Quer, Jannes and Richter, Lorenz and Sch{\"u}tte, Christof}, title = {Improving control based importance sampling strategies for metastable diffusions via adapted metadynamics}, issn = {1438-0064}, year = {2021}, abstract = {Sampling rare events in metastable dynamical systems is often a computationally expensive task and one needs to resort to enhanced sampling methods such as importance sampling. Since we can formulate the problem of finding optimal importance sampling controls as a stochastic optimization problem, this then brings additional numerical challenges and the convergence of corresponding algorithms might as well suffer from metastabilty. In this article we address this issue by combining systematic control approaches with the heuristic adaptive metadynamics method. Crucially, we approximate the importance sampling control by a neural network, which makes the algorithm in principle feasible for high dimensional applications. We can numerically demonstrate in relevant metastable problems that our algorithm is more effective than previous attempts and that only the combination of the two approaches leads to a satisfying convergence and therefore to an efficient sampling in certain metastable settings.}, language = {en} } @misc{EiflerGleixner2023, author = {Eifler, Leon and Gleixner, Ambros}, title = {Safe and Verified Gomory Mixed Integer Cuts in a Rational MIP Framework}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-90159}, year = {2023}, abstract = {This paper is concerned with the exact solution of mixed-integer programs (MIPs) over the rational numbers, i.e., without any roundoff errors and error tolerances. Here, one computational bottleneck that should be avoided whenever possible is to employ large-scale symbolic computations. Instead it is often possible to use safe directed rounding methods, e.g., to generate provably correct dual bounds. In this work, we continue to leverage this paradigm and extend an exact branch-and-bound framework by separation routines for safe cutting planes, based on the approach first introduced by Cook, Dash, Fukasawa, and Goycoolea in 2009. Constraints are aggregated safely using approximate dual multipliers from an LP solve, followed by mixed-integer rounding to generate provably valid, although slightly weaker inequalities. We generalize this approach to problem data that is not representable in floating-point arithmetic, add routines for controlling the encoding length of the resulting cutting planes, and show how these cutting planes can be verified according to the VIPR certificate standard. Furthermore, we analyze the performance impact of these cutting planes in the context of an exact MIP framework, showing that we can solve 21.5\% more instances and reduce solving times by 26.8\% on the MIPLIB 2017 benchmark test set.}, language = {en} } @misc{BorndoerferDaneckerWeiser2023, author = {Bornd{\"o}rfer, Ralf and Danecker, Fabian and Weiser, Martin}, title = {Convergence Properties of Newton's Method for Globally Optimal Free Flight Trajectory Optimization}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-91309}, year = {2023}, abstract = {The algorithmic efficiency of Newton-based methods for Free Flight Trajectory Optimization is heavily influenced by the size of the domain of convergence. We provide numerical evidence that the convergence radius is much larger in practice than what the theoretical worst case bounds suggest. The algorithm can be further improved by a convergence-enhancing domain decomposition.}, language = {en} } @misc{GriewankHasenfelderRadonsetal.2017, author = {Griewank, Andreas and Hasenfelder, Richard and Radons, Manuel and Lehmann, Lutz and Streubel, Tom}, title = {Integrating Lipschitzian Dynamical Systems using Piecewise Algorithmic Differentiation}, issn = {1438-0064}, doi = {10.1080/10556788.2017.1378653}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64639}, year = {2017}, abstract = {In this article we analyze a generalized trapezoidal rule for initial value problems with piecewise smooth right hand side \(F:R^n \to R^n\) based on a generalization of algorithmic differentiation. When applied to such a problem, the classical trapezoidal rule suffers from a loss of accuracy if the solution trajectory intersects a nondifferentiability of \(F\). The advantage of the proposed generalized trapezoidal rule is threefold: Firstly, we can achieve a higher convergence order than with the classical method. Moreover, the method is energy preserving for piecewise linear Hamiltonian systems. Finally, in analogy to the classical case we derive a third order interpolation polynomial for the numerical trajectory. In the smooth case the generalized rule reduces to the classical one. Hence, it is a proper extension of the classical theory. An error estimator is given and numerical results are presented.}, language = {en} } @misc{GriewankStreubelLehmannetal.2016, author = {Griewank, Andreas and Streubel, Tom and Lehmann, Lutz and Hasenfelder, Richard and Radons, Manuel}, title = {Piecewise linear secant approximation via Algorithmic Piecewise Differentiation}, issn = {1438-0064}, doi = {10.1080/10556788.2017.1387256}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-61642}, year = {2016}, abstract = {It is shown how piecewise differentiable functions \(F: R^n → R^m\) that are defined by evaluation programs can be approximated locally by a piecewise linear model based on a pair of sample points x̌ and x̂. We show that the discrepancy between function and model at any point x is of the bilinear order O(||x - x̌|| ||x - x̂||). This is a little surprising since x ∈ R^n may vary over the whole Euclidean space, and we utilize only two function samples F̌ = F(x̌) and F̂ = F(x̂), as well as the intermediates computed during their evaluation. As an application of the piecewise linearization procedure we devise a generalized Newton's method based on successive piecewise linearization and prove for it sufficient conditions for convergence and convergence rates equaling those of semismooth Newton. We conclude with the derivation of formulas for the numerically stable implementation of the aforedeveloped piecewise linearization methods.}, language = {en} }