Refine
Year of publication
- 2018 (2) (remove)
Document Type
- ZIB-Report (2) (remove)
Language
- English (2)
Has Fulltext
- yes (2)
Is part of the Bibliography
- no (2)
Keywords
- Algorithmic Differentiation (1)
- Automatic Differentiation (1)
- Computational Cost (1)
- Implementation (1)
- Lipschitz Continuity (1)
- Nonsmooth (1)
- Piecewise Linearization (1)
- Trapezoidal Rule (1)
Institute
- Mathematical Optimization (2) (remove)
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.
Recent research has shown that piecewise smooth (PS) functions can be approximated by piecewise linear functions with second order error in the distance to
a given reference point. A semismooth Newton type algorithm based on successive application of these piecewise linearizations was subsequently developed
for the solution of PS equation systems. For local bijectivity of the linearization
at a root, a radius of quadratic convergence was explicitly calculated in terms
of local Lipschitz constants of the underlying PS function. In the present work
we relax the criterium of local bijectivity of the linearization to local openness.
For this purpose a weak implicit function theorem is proved via local mapping
degree theory. It is shown that there exist PS functions f:IR^2 --> IR^2 satisfying the weaker
criterium where every neighborhood of the root of f contains a point x such that
all elements of the Clarke Jacobian at x are singular. In such neighborhoods
the steps of classical semismooth Newton are not defined, which establishes
the new method as an independent algorithm. To further clarify the relation between a PS function and its piecewise linearization,
several statements about structure correspondences between the two are proved.
Moreover, the influence of the specific representation of the local piecewise linear models
on the robustness of our method is studied.
An example application from cardiovascular mathematics is given.