基于哈密尔顿力学的物理引擎

Question

我正在尝试实现一个基于哈密尔顿力学的物理引擎。 面对几个问题，比如

1. 区分H功能
2. 部分评估H函数的偏导数以有效使用Runge-Kutta方法

我决定将函数表示为struct（Partial evaluation of a function represented as a struct）。

然而，将C ++模板的代码与标准C ++代码相结合证明是困难的。此外，评估结构函数似乎比相应的标准函数至少花费2倍的时间。

所以我问Stackoverflow社区，我应该保留我的方法还是改变它。

修改

我使用g ++并从命令行编译：

g ++ -std = c ++ 14 -IC：\ Hudba \ include Stin.cpp Tma.cpp -LC：\ Hudba \ lib -lglfw3 -lglew32 -lglu32 -lopengl32 -lkernel32 -luser32 -lgdi32 -lws2_32。

为了部分评估struct-function，我需要传递一个在编译时不知道的值数组，使用模板是不可能的。