Robotics, Finite Element Analysis, Programming, Cpp.
matrix = mff.inverse() * kff; Eigen::EigenSolver<Eigen::MatrixXd> es(matrix); Eigen::MatrixXcd eigValeues = es.eigenvalues(); std::vector<double> nfreq; for (int ii = 0; ii < eigValeues.rows(); ii++) { double w = std::sqrt((real(eigValeues(ii)))); nfreq.push_back(w / (2 * M_PI)); }