期刊
JOURNAL OF COMPUTATIONAL SCIENCE
卷 4, 期 5, 页码 334-344出版社
ELSEVIER
DOI: 10.1016/j.jocs.2011.10.007
关键词
Automatic differentiation; Cholesky decomposition; Hierarchical approach; Higher-order derivatives; Numerical linear algebra; NumPy; Taylor arithmetic
资金
- Bundesministerium fur Bildung und Forschung (BMBF) [03GRPAL3]
- Humboldt Universitat zu Berlin
Many programs for scientific computing in Python are based on NumPy and therefore make heavy use of numerical linear algebra (NLA) functions, vectorized operations, slicing and broadcasting. AlgoPy provides the means to compute derivatives of arbitrary order and Taylor approximations of such programs. The approach is based on a combination of univariate Taylor polynomial arithmetic and matrix calculus in the (combined) forward/reverse mode of Algorithmic Differentiation (AD). In contrast to existing AD tools, vectorized operations and NLA functions are not considered to be a sequence of scalar elementary functions. Instead, dedicated algorithms for the matrix product, matrix inverse and the Cholesky, QR, and symmetric eigenvalue decomposition are implemented in AlgoPy. We discuss the reasons for this alternative approach and explain the underlying idea. Examples illustrate how AlgoPy can be used from a user's point of view. (C) 2011 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据