隐式欧拉法的微分反演:符号分析

Uwe Naumann
{"title":"隐式欧拉法的微分反演:符号分析","authors":"Uwe Naumann","doi":"arxiv-2409.05445","DOIUrl":null,"url":null,"abstract":"The implicit Euler method integrates systems of ordinary differential\nequations $$\\frac{d x}{d t}=G(t,x(t))$$ with differentiable right-hand side $G\n: R \\times R^n \\rightarrow R^n$ from an initial state $x=x(0) \\in R^n$ to a\ntarget time $t \\in R$ as $x(t)=E(t,m,x)$ using an equidistant discretization of\nthe time interval $[0,t]$ yielding $m>0$ time steps. We aim to compute the\nproduct of its inverse Jacobian $$ (E')^{-1} \\equiv \\left (\\frac{d E}{d x}\\right )^{-1} \\in R^{n \\times n} $$ with a given vector efficiently. We show that the differential inverse\n$(E')^{-1} \\cdot v$ can be evaluated for given $v \\in R^n$ with a computational\ncost of $\\mathcal{O}(m \\cdot n^2)$ as opposed to the standard $\\mathcal{O}(m\n\\cdot n^3)$ or, naively, even $\\mathcal{O}(m \\cdot n^4).$ The theoretical\nresults are supported by actual run times. A reference implementation is\nprovided.","PeriodicalId":501162,"journal":{"name":"arXiv - MATH - Numerical Analysis","volume":"32 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Differential Inversion of the Implicit Euler Method: Symbolic Analysis\",\"authors\":\"Uwe Naumann\",\"doi\":\"arxiv-2409.05445\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"The implicit Euler method integrates systems of ordinary differential\\nequations $$\\\\frac{d x}{d t}=G(t,x(t))$$ with differentiable right-hand side $G\\n: R \\\\times R^n \\\\rightarrow R^n$ from an initial state $x=x(0) \\\\in R^n$ to a\\ntarget time $t \\\\in R$ as $x(t)=E(t,m,x)$ using an equidistant discretization of\\nthe time interval $[0,t]$ yielding $m>0$ time steps. We aim to compute the\\nproduct of its inverse Jacobian $$ (E')^{-1} \\\\equiv \\\\left (\\\\frac{d E}{d x}\\\\right )^{-1} \\\\in R^{n \\\\times n} $$ with a given vector efficiently. We show that the differential inverse\\n$(E')^{-1} \\\\cdot v$ can be evaluated for given $v \\\\in R^n$ with a computational\\ncost of $\\\\mathcal{O}(m \\\\cdot n^2)$ as opposed to the standard $\\\\mathcal{O}(m\\n\\\\cdot n^3)$ or, naively, even $\\\\mathcal{O}(m \\\\cdot n^4).$ The theoretical\\nresults are supported by actual run times. A reference implementation is\\nprovided.\",\"PeriodicalId\":501162,\"journal\":{\"name\":\"arXiv - MATH - Numerical Analysis\",\"volume\":\"32 1\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2024-09-09\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - MATH - Numerical Analysis\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2409.05445\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - MATH - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.05445","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

摘要

隐式欧拉法积分常微分方程系统 $$\frac{d x}{d t}=G(t,x(t))$$ 具有可微分右边$G:R \times R^n \rightarrow R^n$ 从初始状态 $x=x(0) \in R^n$ 到目标时间 $t \in R$ 为 $x(t)=E(t,m,x)$,使用时间区间 $[0,t]$ 的等距离散化,产生 $m>0$ 的时间步长。我们的目标是在 R^{n \times n} $$中用给定矢量高效地计算其逆雅各布值 $ (E')^{-1} \equiv \left (\frac{d E}{d x}\right )^{-1} \$ 的乘积。我们证明,对于 R^n$ 中的给定 $v ,可以用 $\mathcal{O}(m\cdot n^2)$ 的计算成本求出微分逆 $(E')^{-1} \cdot v$,而不是标准的 $\mathcal{O}(m\cdot n^3)$ ,甚至不是 $\mathcal{O}(m\cdot n^4) 。我们提供了一个参考实现。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
Differential Inversion of the Implicit Euler Method: Symbolic Analysis
The implicit Euler method integrates systems of ordinary differential equations $$\frac{d x}{d t}=G(t,x(t))$$ with differentiable right-hand side $G : R \times R^n \rightarrow R^n$ from an initial state $x=x(0) \in R^n$ to a target time $t \in R$ as $x(t)=E(t,m,x)$ using an equidistant discretization of the time interval $[0,t]$ yielding $m>0$ time steps. We aim to compute the product of its inverse Jacobian $$ (E')^{-1} \equiv \left (\frac{d E}{d x}\right )^{-1} \in R^{n \times n} $$ with a given vector efficiently. We show that the differential inverse $(E')^{-1} \cdot v$ can be evaluated for given $v \in R^n$ with a computational cost of $\mathcal{O}(m \cdot n^2)$ as opposed to the standard $\mathcal{O}(m \cdot n^3)$ or, naively, even $\mathcal{O}(m \cdot n^4).$ The theoretical results are supported by actual run times. A reference implementation is provided.
求助全文
通过发布文献求助,成功后即可免费获取论文全文。 去求助
来源期刊
自引率
0.00%
发文量
0
期刊最新文献
A Lightweight, Geometrically Flexible Fast Algorithm for the Evaluation of Layer and Volume Potentials Adaptive Time-Step Semi-Implicit One-Step Taylor Scheme for Stiff Ordinary Differential Equations Conditions aux limites fortement non lin{é}aires pour les {é}quations d'Euler de la dynamique des gaz Fully guaranteed and computable error bounds on the energy for periodic Kohn-Sham equations with convex density functionals A novel Mortar Method Integration using Radial Basis Functions
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1