使用 CoDiPack 自动区分

Question

以下代码：

#include <codi.hpp>
...
codi::RealForward Gcodi[l];
for (int p = 0; p < l; p++) 
{
     ...
     double a = Gcodi[p];
}

给我编译错误：

nnBFAD.cpp：在函数“void OptBF()”中：

nnBFAD.cpp:156:25: 错误：无法将 'codi::RealForward {aka codi::ActiveReal >}' 转换为 'double' in initialization double

a = Gcodi[p];

提示？

Answer 1

根据此处的官方文档， RealForward是一种重载了赋值运算符的类型，因此您可以使用 double..

喜欢做

codi::RealForward a = 3.0;

相反的方向当然没有定义，

这就是为什么你不能通过以下方式直接将codi::RealForward转换为 double 的原因：

double a = Gcodi[p];

但你可以调用函数，即

double a = Gcodi[p].getGradient();

---

更新：

然后你可以分配一个带有双精度的 RealForward 对象

double myDouble{3.3};
RealForward a = myDouble;

但直接从 REalForwad 分配双精度值是不合法的：

RealForward a = ...;
double myDouble = a; //not valid!

其他例子

RealForward b = a * a; //this is ok because a * a is a double

Answer 2

我找到了一个函数：getValue()

双 a = Gcodi[p].getValue();

如中所述： https ://www.scicomp.uni-kl.de/codi/d5/dbb/Tutorial1.html