complexsolver.cxx

Go to the documentation of this file.

 
#include <Eigen/Eigenvalues>
 
#include "complexsolver.h"
 
//.............................................................................
void complexsolver(const Eigen::Matrix3cd& A, OscProb::vectorD& w)
{
  Eigen::ComplexEigenSolver<Eigen::Matrix3cd> eigensolver;
  eigensolver.compute(A);
 
  for (int t = 0; t < w.size(); t++) {
    w[t] = eigensolver.eigenvalues()(t).real();
  }
}