18#include <adolc/adolc.h>
76 { forward(tag,1,1,
i,
i+1,X,Y);
78 cout << Y[0][
i] <<
" - " << y.value() <<
" = " << Y[0][
i]-y.value()
79 <<
" (should be 0)\n";
81 Z[0][
i] = Z[0][
i-1]/
i;
82 cout << Y[0][
i] <<
" - " << Z[0][
i] <<
" = " << Y[0][
i]-Z[0][
i]
83 <<
" (should be 0)\n";
85 reverse(tag,1,1,
i,u,Z);
adouble power(adouble x, int n)
FTensor::Index< 'n', SPACE_DIM > n
FTensor::Index< 'i', SPACE_DIM > i