@article{GoetschelWeiser2010, author = {G{\"o}tschel, Sebastian and Weiser, Martin}, title = {State Trajectory Compression in Optimal Control}, volume = {10}, journal = {PAMM}, number = {1}, doi = {10.1002/pamm.201010282}, pages = {579 -- 580}, year = {2010}, abstract = {In optimal control problems with nonlinear time-dependent 3D PDEs, the computation of the reduced gradient by adjoint methods requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. Since the state enters into the adjoint equation, the storage of a 4D discretization is necessary. We propose a lossy compression algorithm using a cheap predictor for the state data, with additional entropy coding of prediction errors. Analytical and numerical results indicate that compression factors around 30 can be obtained without exceeding the FE discretization error.}, language = {en} } @misc{Weiser2010, author = {Weiser, Martin}, title = {Delayed Residual Compensation for Bidomain Equations}, volume = {1281}, journal = {AIP Conference Proceedings}, doi = {10.1063/1.3498495}, pages = {419 -- 422}, year = {2010}, abstract = {The biodomain model of cardioelectric excitation consists of a reaction-diffusion equation, an elliptic algebraic constraint, and a set of pointwise ODEs. Fast reaction enforces small time steps, such that for common mesh sizes the reaction-diffusion equation is easily solved implicitly due to a dominating mass matrix. In contrast, the elliptic constraint does not benefit from small time steps and requires a comparably expensive solution. We propose a delayed residual compensation that improves the solution of the elliptic constraint and thus alleviates the need for long iteration times.}, language = {en} } @misc{WeiserErdmannDeuflhard2010, author = {Weiser, Martin and Erdmann, Bodo and Deuflhard, Peter}, title = {On Efficiency and Accuracy in Cardioelectric Simulation}, journal = {Progress in Industrial Mathematics at ECMI 2008}, editor = {Wilson, E. and Fitt, A. and Ockendon, H. and Norbury, J.}, publisher = {Springer}, pages = {371 -- 376}, year = {2010}, abstract = {Reasons for the failure of adaptive methods to deliver improved efficiency when integrating monodomain models for myocardiac excitation are discussed. Two closely related techniques for reducing the computational complexity of linearly implicit integrators, deliberate sparsing and splitting, are investigated with respect to their impact on computing time and accuracy.}, language = {en} } @article{WilhelmsSeemannWeiseretal.2010, author = {Wilhelms, Mathias and Seemann, Gunnar and Weiser, Martin and D{\"o}ssel, Olaf}, title = {Benchmarking Solvers of the Monodomain Equation in Cardiac Electrophysiological Modeling}, volume = {55}, journal = {Biomed. Engineer.}, doi = {10.1515/BMT.2010.712}, pages = {99 -- 102}, year = {2010}, language = {en} } @misc{WustWeihrauchWeiseretal.2010, author = {Wust, Peter and Weihrauch, Mirko and Weiser, Martin and Gellermann, Johanna and Eisenhardt, Steffen and Chobrok, Thorsten and Budach, Volker}, title = {Optimization of clinical radiofrequency hyperthermia by use of MR-thermography in a hybrid system}, journal = {World Congress on Medical Physics and Biomedical Engineering, September 2009, Munich, Germany}, editor = {D{\"o}ssel, O. and Schlegel, W. and Magjarevic, R.}, publisher = {Springer}, pages = {174 -- 175}, year = {2010}, language = {en} } @article{LubkollSchielaWeiser2014, author = {Lubkoll, Lars and Schiela, Anton and Weiser, Martin}, title = {An optimal control problem in polyconvex hyperelasticity}, volume = {52}, journal = {SIAM J. Control Opt.}, number = {3}, doi = {10.1137/120876629}, pages = {1403 -- 1422}, year = {2014}, abstract = {We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented.}, language = {en} } @article{GuentherLameckerWeiser2013, author = {G{\"u}nther, Andreas and Lamecker, Hans and Weiser, Martin}, title = {Flexible Shape Matching with Finite Element Based LDDMM}, volume = {105}, journal = {International Journal of Computer Vision}, number = {2}, doi = {10.1007/s11263-012-0599-3}, pages = {128 -- 143}, year = {2013}, abstract = {We consider Large Deformation Diffeomorphic Metric Mapping of general \$m\$-currents. After stating an optimization algorithm in the function space of admissable morph generating velocity fields, two innovative aspects in this framework are presented and numerically investigated: First, we spatially discretize the velocity field with conforming adaptive finite elements and discuss advantages of this new approach. Second, we directly compute the temporal evolution of discrete \$m\$-current attributes.}, language = {en} } @misc{GuentherLameckerWeiser2011, author = {G{\"u}nther, Andreas and Lamecker, Hans and Weiser, Martin}, title = {Direct LDDMM of Discrete Currents with Adaptive Finite Elements}, journal = {Proceedings of the Third International Workshop on Mathematical Foundations of Computational Anatomy - Geometrical and Statistical Methods for Modelling Biological Shape Variability}, editor = {Pennec, X. and Joshi, S. and Nielsen, M.}, pages = {1 -- 14}, year = {2011}, abstract = {We consider Large Deformation Diffeomorphic Metric Mapping of general \$m\$-currents. After stating an optimization algorithm in the function space of admissable morph generating velocity fields, two innovative aspects in this framework are presented and numerically investigated: First, we spatially discretize the velocity field with conforming adaptive finite elements and discuss advantages of this new approach. Second, we directly compute the temporal evolution of discrete \$m\$-current attributes.}, language = {en} } @misc{SchielaWeiser2010, author = {Schiela, Anton and Weiser, Martin}, title = {Barrier methods for a control problem from hyperthermia treatment planning}, journal = {Recent Advances in Optimization and its Applications in Engineering (Proceedings of 14th Belgian-French-German Conference on Optimization 2009)}, editor = {Diehl, M. and Glineur, F. and Jarlebring, E. and Michiels, W.}, publisher = {Springer}, pages = {419 -- 428}, year = {2010}, abstract = {We consider an optimal control problem from hyperthermia treatment planning and its barrier regularization. We derive basic results, which lay the groundwork for the computation of optimal solutions via an interior point path-following method. Further, we report on a numerical implementation of such a method and its performance at an example problem.}, language = {en} } @article{WeiserRoelligArndtetal.2010, author = {Weiser, Martin and R{\"o}llig, Mathias and Arndt, Ralf and Erdmann, Bodo}, title = {Development and test of a numerical model for pulse thermography in civil engineering}, volume = {46}, journal = {Heat and Mass Transfer}, number = {11-12}, pages = {1419 -- 1428}, year = {2010}, abstract = {Pulse thermography of concrete structures is used in civil engineering for detecting voids, honeycombing and delamination. The physical situation is readily modeled by Fourier's law. Despite the simplicity of the PDE structure, quantitatively realistic numerical 3D simulation faces two major obstacles. First, the short heating pulse induces a thin boundary layer at the heated surface which encapsulates all information and therefore has to be resolved faithfully. Even with adaptive mesh refinement techniques, obtaining useful accuracies requires an unsatisfactorily fine discretization. Second, bulk material parameters and boundary conditions are barely known exactly. We address both issues by a semi-analytic reformulation of the heat transport problem and by parameter identification. Numerical results are compared with measurements of test specimens.}, language = {en} } @article{RannebergWeiserWeihrauchetal.2010, author = {Ranneberg, Maximilian and Weiser, Martin and Weihrauch, Mirko and Budach, Volker and Gellermann, Johanna and Wust, Peter}, title = {Regularized Antenna Profile Adaptation in Online Hyperthermia Treatment}, volume = {37}, journal = {Medical Physics}, doi = {10.1118/1.3488896}, pages = {5382 -- 5394}, year = {2010}, language = {en} } @article{LubkollSchielaWeiser2017, author = {Lubkoll, Lars and Schiela, Anton and Weiser, Martin}, title = {An affine covariant composite step method for optimization with PDEs as equality constraints}, volume = {32}, journal = {Optimization Methods and Software}, number = {5}, doi = {10.1080/10556788.2016.1241783}, pages = {1132 -- 1161}, year = {2017}, abstract = {We propose a composite step method, designed for equality constrained optimization with partial differential equations. Focus is laid on the construction of a globalization scheme, which is based on cubic regularization of the objective and an affine covariant damped Newton method for feasibility. We show finite termination of the inner loop and fast local convergence of the algorithm. We discuss preconditioning strategies for the iterative solution of the arising linear systems with projected conjugate gradient. Numerical results are shown for optimal control problems subject to a nonlinear heat equation and subject to nonlinear elastic equations arising from an implant design problem in craniofacial surgery.}, language = {en} } @misc{GoetschelMaierhoferMuelleretal.2016, author = {G{\"o}tschel, Sebastian and Maierhofer, Christiane and M{\"u}ller, Jan P. and Rothbart, Nick and Weiser, Martin}, title = {Quantitative Defect Reconstruction in Active Thermography for Fiber-Reinforced Composites}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-58374}, year = {2016}, abstract = {Carbon-fiber reinforced composites are becoming more and more important in the production of light-weight structures, e.g., in the automotive and aerospace industry. Thermography is often used for non-destructive testing of these products, especially to detect delaminations between different layers of the composite. In this presentation, we aim at methods for defect reconstruction from thermographic measurements of such carbon-fiber reinforced composites. The reconstruction results shall not only allow to locate defects, but also give a quantitative characterization of the defect properties. We discuss the simulation of the measurement process using finite element methods, as well as the experimental validation on flat bottom holes. Especially in pulse thermography, thin boundary layers with steep temperature gradients occurring at the heated surface need to be resolved. Here we use the combination of a 1D analytical solution combined with numerical solution of the remaining defect equation. We use the simulations to identify material parameters from the measurements. Finally, fast heuristics for reconstructing defect geometries are applied to the acquired data, and compared for their accuracy and utility in detecting different defects like back surface defects or delaminations.}, language = {en} } @misc{GoetschelHoehneKolkoorietal.2016, author = {G{\"o}tschel, Sebastian and H{\"o}hne, Christian and Kolkoori, Sanjeevareddy and Mitzscherling, Steffen and Prager, Jens and Weiser, Martin}, title = {Ray Tracing Boundary Value Problems: Simulation and SAFT Reconstruction for Ultrasonic Testing}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-58386}, year = {2016}, abstract = {The application of advanced imaging techniques for the ultrasonic inspection of inhomogeneous anisotropic materials like austenitic and dissimilar welds requires information about acoustic wave propagation through the material, in particular travel times between two points in the material. Forward ray tracing is a popular approach to determine traveling paths and arrival times but is ill suited for inverse problems since a large number of rays have to be computed in order to arrive at prescribed end points. In this contribution we discuss boundary value problems for acoustic rays, where the ray path between two given points is determined by solving the eikonal equation. The implementation of such a two point boundary value ray tracer for sound field simulations through an austenitic weld is described and its efficiency as well as the obtained results are compared to those of a forward ray tracer. The results are validated by comparison with experimental results and commercially available UT simulation tools. As an application, we discuss an implementation of the method for SAFT (Synthetic Aperture Focusing Technique) reconstruction. The ray tracer calculates the required travel time through the anisotropic columnar grain structure of the austenitic weld. There, the formulation of ray tracing as a boundary value problem allows a straightforward derivation of the ray path from a given transducer position to any pixel in the reconstruction area and reduces the computational cost considerably.}, language = {en} } @book{Weiser2016, author = {Weiser, Martin}, title = {Inside Finite Elements}, publisher = {De Gruyter}, year = {2016}, abstract = {All relevant implementation aspects of finite element methods are discussed in this book. The focus is on algorithms and data structures as well as on their concrete implementation. Theory is covered as far as it gives insight into the construction of algorithms.Throughout the exercises a complete FE-solver for scalar 2D problems will be implemented in Matlab/Octave.}, language = {en} } @misc{Makarenko2017, type = {Master Thesis}, author = {Makarenko, Oleksandra}, title = {Square Root Approximation und statistische Umgewichtungstechnik f{\"u}r das Faltungsverhalten von Polymerb{\"u}rsten}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-65283}, year = {2017}, language = {de} } @misc{VolkweinWeiser2000, author = {Volkwein, Stefan and Weiser, Martin}, title = {Affine Invariant Convergence Analysis for Inexact Augmented Lagrangian-SQP Methods}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6243}, number = {00-56}, year = {2000}, abstract = {An affine invariant convergence analysis for inexact augmented Lagrangian-SQP methods is presented. The theory is used for the construction of an accuracy matching between iteration errors and truncation errors, which arise from the inexact linear system solves. The theoretical investigations are illustrated numerically by an optimal control problem for the Burgers equation.}, language = {en} } @misc{WeiserDeuflhard2001, author = {Weiser, Martin and Deuflhard, Peter}, title = {The Central Path towards the Numerical Solution of Optimal Control Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6380}, number = {01-12}, year = {2001}, abstract = {A new approach to the numerical solution of optimal control problems including control and state constraints is presented. Like hybrid methods, the approach aims at combining the advantages of direct and indirect methods. Unlike hybrid methods, however, our method is directly based on interior-point concepts in function space --- realized via an adaptive multilevel scheme applied to the complementarity formulation and numerical continuation along the central path. Existence of the central path and its continuation towards the solution point is analyzed in some theoretical detail. An adaptive stepsize control with respect to the duality gap parameter is worked out in the framework of affine invariant inexact Newton methods. Finally, the performance of a first version of our new type of algorithm is documented by the successful treatment of the well-known intricate windshear problem.}, language = {en} } @misc{ErdmannKoberLangetal.2001, author = {Erdmann, Bodo and Kober, Cornelia and Lang, Jens and Sader, Robert and Zeilhofer, Hans-Florian and Deuflhard, Peter}, title = {Efficient and Reliable Finite Element Methods for Simulation of the Human Mandible}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6403}, number = {01-14}, year = {2001}, abstract = {By computed tomography data (CT), the individual geometry of the mandible is quite well reproduced, also the separation between cortical and trabecular bone. Using anatomical knowledge about the architecture and the functional potential of the masticatory muscles, realistic situations were approximated. The solution of the underlying partial differential equations describing linear elastic material behaviour is provided by an adaptive finite element method. Estimations of the discretization error, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method.}, language = {en} } @misc{HinzeSchiela2007, author = {Hinze, Michael and Schiela, Anton}, title = {Discretization of Interior Point Methods for State Constrained Elliptic Optimal Control Problems: Optimal Error Estimates and Parameter Adjustment}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-10414}, number = {07-40}, year = {2007}, abstract = {An adjustment scheme for the relaxation parameter of interior point approaches to the numerical solution of pointwise state constrained elliptic optimal control problems is introduced. The method is based on error estimates of an associated finite element discretization of the relaxed problems and optimally selects the relaxation parameter in dependence on the mesh size of discretization. The finite element analysis for the relaxed problems is carried out and a numerical example is presented which confirms our analytical findings.}, language = {en} }