This would be really nice to have.  I tried recently to use the unsupported matrix power, and ran into a number of segfaults for very simple matrix powers:

#include <iostream>
#include <unsupported/Eigen/MatrixFunctions>
int main()
{
  Eigen::MatrixXcd A(2,2);
  A << 1, 1, 4, -2;
  // Segfaults:
  //Eigen::MatrixPower<Eigen::MatrixXcd> Apow(A);
  //std::cout << Apow(3.1) << "\n\n";

  // Segfaults:
  //std::cout << A.pow(3.1) << "\n\n";

  // Segfaults:
  //std::cout << A.pow(3.1).eval() << "\n\n";

  // Segfaults:
  //Eigen::MatrixXcd B(2,2);
  //B = A.pow(3.1);
}

I use matrix power a lot, it would be nice if this worked in Eigen.

Actually this entry is about coefficient-wise power, see the function signatures. So let's continue your issue on bug 614.

-- GitLab Migration Automatic Message --

This bug has been migrated to gitlab.com's GitLab instance and has been closed from further activity.

You can subscribe and participate further through the new bug through this link to our GitLab instance: https://gitlab.com/libeigen/eigen/issues/465.