Как найти вторую производную Eigen AutoDiffScalar? - PullRequest
0 голосов
/ 22 июня 2019

Я погружаюсь в модуль AutoDiff Eigen-3.3.7.Рассмотрим этот простой пример нахождения производных x 2 .

#include <eigen3/Eigen/Core>
#include <eigen3/unsupported/Eigen/AutoDiff>
#include <iostream>

int main()
{
  Eigen::AutoDiffScalar<Eigen::Vector2d> x(8.0, Eigen::Vector2d(1,0)), y;

  y = x*x;

  std::cout << "x = " << x << "\n"
            << "y = " << y << "\n"
            << "y' = " << y.derivatives()[0] << "\n"
            << "y'' = " << y.derivatives()[1] << "\n";

  return 0;
}

Выход для второй производной - 0 вместо 2. Как правильно ее найти?

...