Oct, 2023

ODE 系统的时间向量化数值积分

TL;DR将刚性系统的常微分方程通过时间的积分与伴随方法计算参数梯度的高效、隐式、矢量化方法在本文中进行了描述。创新之处在于在独立时间序列的数量和连续时间步骤的批次或 “块” 上对问题进行矢量化,从而有效地进行隐式常微分方程系统的组装。线性化的隐式反向欧拉法的块双对角结构使用并行循环消减进一步进行矢量化。在输入数据的两个轴上进行矢量化,提供了更高的计算带宽给计算设备,即使是相对稀疏数据的问题也能充分利用现代 GPU,实现了超过 100 倍的加速,与标准的顺序时间积分相比。我们通过几个示例问题,包括来自分析刚性和非刚性 ODE 模型以及神经 ODE 模型的问题,展示了隐式的矢量化时间积分的优势,并描述并提供了此处开发的方法的免费开源实现。