Symplectic Adjoint Method for Exact Gradient of Neural ODE with Minimal Memory
This addresses a key bottleneck in training neural ODEs for continuous-time dynamical systems, offering a more efficient and reliable gradient computation method.
The paper tackles the high memory consumption and computational cost of gradient computation in neural ODEs by proposing the symplectic adjoint method, which achieves exact gradients with memory proportional to the number of uses plus network size, demonstrating reduced memory usage, faster performance, and improved robustness to errors.
A neural network model of a differential equation, namely neural ODE, has enabled the learning of continuous-time dynamical systems and probabilistic distributions with high accuracy. The neural ODE uses the same network repeatedly during a numerical integration. The memory consumption of the backpropagation algorithm is proportional to the number of uses times the network size. This is true even if a checkpointing scheme divides the computation graph into sub-graphs. Otherwise, the adjoint method obtains a gradient by a numerical integration backward in time. Although this method consumes memory only for a single network use, it requires high computational cost to suppress numerical errors. This study proposes the symplectic adjoint method, which is an adjoint method solved by a symplectic integrator. The symplectic adjoint method obtains the exact gradient (up to rounding error) with memory proportional to the number of uses plus the network size. The experimental results demonstrate that the symplectic adjoint method consumes much less memory than the naive backpropagation algorithm and checkpointing schemes, performs faster than the adjoint method, and is more robust to rounding errors.