Abstract
The reliability of numerical results provided by conventional numerical integration methods for ODEs or DAEs depends on the properties of the problem. In particular, since the solution may not be unique at singular points, arbitrary solutions may be obtained. For DAEs, such singularities may occur, if the structure or the dimension of the spaces related to the DAE change. Moreover, even though a numerical singularity is not given in a strict mathematical sense, the numerical behavior may be analogous for sensitive problems. In this contribution, we aim at a characterization of this sensitivity, considering the condition number of a suitable matrix related to the DAE that is constructed using automatic differentiation. We show how this approach, which builds on the projector based analysis, can be applied to properly stated DAEs of index up to 2.
Similar content being viewed by others
References
Deuflhard, P., Bornemann, F.: Numerical mathematics 2. Ordinary differential equations. (Numerische Mathematik 2. Gewöhnliche Differentialgleichungen.) 4th revised and augmented ed. de Gruyter Studium. Berlin (2013)
Estévez Schwarz, D.: Consistent initialization for index-2 differential algebraic equations and its application to circuit simulation. PhD thesis, Berlin: Humboldt-University, Mathematisch-Naturwissenschaftliche Fakultät II, http://edoc.hu-berlin.de/docviews/abstract.php?id=10218 (2000)
Estévez Schwarz, D., Lamour, R.: Diagnosis of singular points of structured DAEs using automatic differentiation. Numerical Algorithms (2014). doi:10.1007/s11075-014-9919-8
Estévez Schwarz, D., Lamour, R.: Progress in Differential-Algebraic Equations, Deskriptor 2013, Differential-Algebraic Equations Forum. Springer. In: Schöps, S., Bartel, A., Günther, M., Ter Maten, E.J.W., Müller, P.C. (eds.) (2014)
Estévez Schwarz, D., Lamour, R.: Projector based integration of DAEs with the Taylor series method using automatic differentiation. J. Comput. Appl. Math. 262, 62–72 (2014)
Golub, G.H., van Loan, C.F.: Matrix Computations. John Hopkins University Press, Baltimore and London (1996)
Jansen, L., Tischendorf, C.: Private communication
Lamour, R., März, R., Tischendorf, C.: Differential-algebraic equations: A projector based analysis. Differential-Algebraic Equations Forum 1. Springer, Berlin (2013)
Sieber, J.: Local error control for general index-1 and index-2 differential-algebraic equations. Technical Report 21, Institut für Mathematik, Humboldt-Universität zu Berlin. http://www2.math.hu-berlin.de/publ/pre/1997/P-97-21.ps (1997)
Silvester, J.R.: Determinants of block matrices. Math. Gaz. 84(501), 460–467 (2000)
Tischendorf, C.: Solution of index-2-DAEs and its application in circuit simulation. PhD thesis, Humboldt-University zu Berlin 1996. appeared in Logos-Verlag Berlin (1996)
Walter, S.F., Lehmann, L.: Algorithmic differentiation in Python with AlgoPy. J. Comput. Sci. 4(5), 334–344 (2013)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Estévez Schwarz, D., Lamour, R. Diagnosis of singular points of properly stated DAEs using automatic differentiation. Numer Algor 70, 777–805 (2015). https://doi.org/10.1007/s11075-015-9973-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-015-9973-x