×

关于自适应时间步长算法的离散伴随。 (英语) Zbl 1177.65098号

摘要:我们研究了自适应时间步长数值算法在自动微分(AD)反向模式下的行为。通过区分原算法的时间步长控制器和误差估计器,反向模式AD生成时间步长的伪伴随导数。由此产生的离散伴随模型与伴随常微分方程不一致,并产生不正确的导数。为了恢复一致性,必须抵消离散伴随模型中非物理导数的贡献。
我们证明了一步显式自适应算法(如Runge-Kutta格式)的离散伴随模型可以通过简单的代码修改使其与连续对应的模型保持一致。此外,我们将分析扩展到涵盖通过离散伴随码的额外前向模式微分导出的二阶伴随模型。几个数值例子支持数学推导。

MSC公司:

65升05 常微分方程初值问题的数值方法
65升06 常微分方程的多步、Runge-Kutta和外推方法
34A34飞机 非线性常微分方程和系统
65磅50 常微分方程的网格生成、细化和自适应方法
65L20英寸 常微分方程数值方法的稳定性和收敛性
65升70 常微分方程数值方法的误差界
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Griewank,A.公司。;Walther,A.,《评估衍生品:算法区分的原理和技术》(2008),SIAM:SIAM Philadelphia,PA,USA·Zbl 1159.65026号
[2] Griewank,A.公司。;朱德斯,D。;Utke,J.,《算法755:ADOL-C:用C/C++编写的算法的自动微分包》,ACM数学软件事务,22,2,131-167(1996)·Zbl 0884.65015号
[3] 科尔曼,T.F。;Verma,A.,ADMIT-1:自动微分和MATLAB接口工具箱,ACM数学软件汇刊,26,1,150-175(2000)·Zbl 1137.65329号
[4] 比肖夫,C.H。;卢武铉。;Mauer-Oats,A.J.,ADIC:ANSI-C的可扩展自动区分工具,软件:实践与经验,27,12,1427-1456(1997)
[5] C.H.Bischof,A.Carle,P.M.Khademi,A.Mauer,Fortran 77程序自动区分的ADIFOR 2.0系统,技术代表MCS-P481-1194,美国伊利诺伊州阿贡,1994年;C.H.Bischof,A.Carle,P.M.Khademi,A.Mauer,Fortran 77程序自动区分的ADIFOR 2.0系统,技术代表MCS-P481-1194,美国伊利诺伊州阿贡,1994年
[6] Giering,R。;Kaminski,T.,应用TAF生成Fortran 77-95程序的高效派生代码,应用数学与力学学报,2,1,54-57(2003)
[7] R.Giering,切线-线性和伴随模型编译器,用户手册1.41999;R.Giering,切线-线性和伴随模型编译器,用户手册1.41999
[8] L.Hascöet,V.Pascual,TAPENADE 2.1用户指南,技术代表0300,INRIA,Sophia Antipolis,法国,2004;L.Hascöet,V.Pascual,TAPENADE 2.1用户指南,技术代表0300,INRIA,Sophia Antipolis,法国
[9] Beck,T.,迭代过程的自动微分,(ICCAM’92:第五届国际计算与应用数学会议论文集(1994),Elsevier Science Publishers B.V.:荷兰阿姆斯特丹Elsevie Science出版商B.V.)·Zbl 0806.65012号
[10] Griewank,A.公司。;比肖夫,C.H。;Corliss,G.F。;A.卡尔。;Williamson,K.,迭代方程求解器的导数收敛性,优化方法和软件,2321-355(1993)
[11] Hager,W.W.,最优控制和变换伴随系统中的Runge-Kutta方法,数值数学,87,2,247-282(2000)·Zbl 0991.49020号
[12] Walther,A.,最优控制显式Runge-Kutta方法的自动微分,计算优化与应用,36,1,83-108(2007)·Zbl 1278.49037号
[13] Sandu,A.,《关于Runge-Kutta离散伴随的性质》,(国际计算科学会议,第4卷(2006年))·Zbl 1157.65421号
[14] Eberhard,P。;Bischof,C.,数值积分算法的自动微分,计算数学,68,226,717-731(1999)·Zbl 1017.65062号
[15] 比肖夫,C.H。;Bücker,H.M。;Rasch,A.,使用自动微分的湍流模型敏感性分析,SIAM科学计算杂志,26,2,510-522(2005)·Zbl 1077.68932号
[16] 托伊瓦宁,J.I。;Mäkinen,R.A.E。;Järvenpää,S.(美国)。;Ylä-Oijala,P。;Rahola,J.,使用矩量法和自动微分法进行电磁敏感性分析和形状优化,IEEE天线和传播学报,57,1,168-175(2009)·Zbl 1369.78951号
[17] 塔朱迪丁,M。;福斯,S.A。;Keane,A.J.,结构动力学求解器的伴随微分,(Bücker,H.M.;Corliss,G.;Hovland,P.;Naumann,U.;Norris,B.,《自动微分:应用、理论和实现》,计算科学与工程讲义(2005),施普林格)·Zbl 1270.74203号
[18] 比肖夫,C.H。;Bücker,H.M。;an Mey,D.,《计算微分应用于中子散射的案例研究》(Corliss,G.;Faure,C.;Griewank,A.;Hascoöt,L.;Naumann,U.,《算法的自动微分:从模拟到优化》,计算机与信息科学(2002),Springer:Springer New York,NY),69-74,(第6章)
[19] 茴香,K。;Losch,M。;施罗德,J。;Wenzel,M.,《测试海洋生态系统模型:敏感性分析和参数优化》,《海洋系统杂志》,28,45-63(2001)
[20] 三都。;Daescu,D。;卡迈克尔,G.R。;Chai,T.,区域空气质量模型的伴随敏感性分析,计算物理杂志,204,1,222-252(2005)·Zbl 1061.92061号
[21] A.Sandu,L.Zhang,大气化学传输建模中的离散二阶伴随,计算物理杂志,印刷版;A.Sandu,L.Zhang,大气化学传输建模中的离散二阶伴随,计算物理杂志,印刷版·Zbl 1144.92043号
[22] 曹毅。;李,S。;佩佐德,L。;Serban,R.,微分代数方程的伴随灵敏度分析:伴随DAE系统及其数值解,SIAM科学计算杂志,24,3,1076-1089(2002)·兹比尔1034.65066
[23] 三都。;Daescu,D。;Carmichael,G.R.,《利用KPP对化学动力学系统进行直接和伴随敏感性分析:第一部分理论和软件工具》,《大气环境》,37,36,5083-5096(2003)
[24] R.Serban,A.C.Hindmarsh,CVODES:SUNDIALS中启用灵敏度的ODE解算器,技术代表UCRL-JP-200037,美国加利福尼亚州利弗莫尔市劳伦斯·利弗莫雷国家实验室,2003年;R.Serban,A.C.Hindmarsh,CVODES:SUNDIALS中启用灵敏度的ODE解算器,技术代表UCRL-JP-200037,美国加利福尼亚州利弗莫尔市劳伦斯·利弗莫雷国家实验室,2003年
[25] LeDimet,F.X。;纳文,I.M。;Daescu,D.,数据同化中的二阶信息,《月度天气评论》,130,3,629-648(2002)
[26] 王,Z。;纳文,I.M。;LeDimet,F.X。;邹欣,二阶伴随分析:理论与应用,气象学与大气物理学,50,1-3,3-20(1992)
[27] Øzyurt,D.B。;巴顿,P.I.,刚性ODE嵌入泛函的廉价二阶方向导数,SIAM科学计算杂志,26,5,1725-1743(2005)·Zbl 1076.65067号
[28] 海尔,E。;诺塞特,S.P。;Wanner,G.,(求解常微分方程:非tiff问题。求解常微分方程:非tiff问题,计算数学,第一卷(1993),Springer Verlag)·Zbl 0789.65048号
[29] A.Sandu,关于离散伴随线性多步方法的一致性性质,技术代表TR-07-40,弗吉尼亚理工学院和州立大学,美国弗吉尼亚州布莱克斯堡,2007年;A.Sandu,关于离散伴随线性多步方法的一致性,技术代表TR-07-40,弗吉尼亚理工学院和州立大学,弗吉尼亚州布莱克斯堡,美国,2007
[30] Giering,R。;Kaminski,T.,伴随码构造的方法,ACM数学软件汇刊,24,4437-474(1998)·Zbl 0934.65027号
[31] Dormand,J.R。;Prince,P.J.,嵌入式Runge-Kutta公式家族,计算与应用数学杂志,6,1,19-26(1980)·Zbl 0448.65045号
[32] Shampine,L.F。;Reichelt,M.W.,《matlab代码套件》,SIAM科学计算杂志,18,1,1-22(1997)·Zbl 0868.65040号
[33] 普罗瑟罗,A。;Robinson,A.,《关于求解刚性常微分方程组的一步法的稳定性和准确性》,《计算数学》,28,125,145-162(1974)·Zbl 0309.65034号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。