@article{GriewankHasenfelderRadonsetal.2018, author = {Griewank, Andreas and Hasenfelder, Richard and Radons, Manuel and Lehmann, Lutz and Streubel, Tom}, title = {Integrating Lipschitzian dynamical systems using piecewise algorithmic differentiation}, series = {Optimization Methods and Software}, volume = {33}, journal = {Optimization Methods and Software}, publisher = {Taylor \& Francis}, doi = {10.1080/10556788.2017.1378653}, pages = {1089 -- 1107}, year = {2018}, abstract = {In this article we analyse a generalized trapezoidal rule for initial value problems with piecewise smooth right-hand side F : IR^n -> IR^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} } @article{GriewankStreubelLehmannetal.2017, author = {Griewank, Andreas and Streubel, Tom and Lehmann, Lutz and Radons, Manuel and Hasenfelder, Richard}, title = {Piecewise linear secant approximation via algorithmic piecewise differentiation}, series = {Optimization Methods and Software}, volume = {33}, journal = {Optimization Methods and Software}, number = {4-6}, publisher = {Taylor \& Francis}, doi = {10.1080/10556788.2017.1387256}, pages = {1108 -- 1126}, year = {2017}, abstract = {It is shown how piecewise differentiable functions F : IR^n -> IR^m that are defined by evaluation programmes can be approximated locally by a piecewise linear model based on a pair of sample points \check x and \hat x. We show that the discrepancy between function and model at any point x is of the bilinear order O(||x - \check x||*||x - \hat x||). 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 equalling those of semismooth Newton. We conclude with the derivation of formulas for the numerically stable implementation of the aforedeveloped piecewise linearization methods.}, language = {en} } @misc{HasenfelderLehmannRadonsetal., author = {Hasenfelder, Richard and Lehmann, Lutz and Radons, Manuel and Streubel, Tom and Strohm, Christian and Griewank, Andreas}, title = {Computational aspects of the Generalized Trapezoidal Rule}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-68615}, abstract = {In this article we analyze a generalized trapezoidal rule for initial value problems with piecewise smooth right hand side F:IR^n -> IR^n. When applied to such a problem, the classical trapezoidal rule suffers from a loss of accuracy if the solution trajectory intersects a non-differentiability of F. In such a situation the investigated generalized trapezoidal rule achieves a higher convergence order than the classical method. While the asymptotic behavior of the generalized method was investigated in a previous work, in the present article we develop the algorithmic structure for efficient implementation strategies and estimate the actual computational cost of the latter. Moreover, energy preservation of the generalized trapezoidal rule is proved for Hamiltonian systems with piecewise linear right hand side.}, language = {en} } @misc{GriewankStreubelLehmannetal., 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}, 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} } @misc{GriewankHasenfelderRadonsetal., 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}, 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} }