/usr/include/Cauchy/Eigen3/Eig.h is in cauchy-dev 0.9.0-0ubuntu1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 | /*
* Copyright (c) 2010 Cyrille Berger <cberger@cberger.net>
*
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public
* License as published by the Free Software Foundation;
* either version 2, or (at your option) any later version of the License.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this library; see the file COPYING.RUNTIME. If not, write to
* the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
* Boston, MA 02110-1301, USA.
*
* The file COPYING.RUNTIME.EXCEPTION contains an exception that allow
* to use this file with other open source license
*/
#ifndef _CAUCHY_EIG_H_
#define _CAUCHY_EIG_H_
#include <Eigen/Eigenvalues>
namespace Cauchy {
Eigen::MatrixXcd eig(const Eigen::MatrixXd& m)
{
Eigen::EigenSolver<Eigen::MatrixXd> es(m);
return es.eigenvalues();
}
Eigen::MatrixXcd eig(const Eigen::MatrixXd& m, Eigen::MatrixXcd* D)
{
Eigen::EigenSolver<Eigen::MatrixXd> es(m);
*D = Eigen::MatrixXcd::Zero(es.eigenvalues().rows(), es.eigenvalues().rows());
for(int i = 0; i < D->rows(); ++i)
{
(*D)(i,i) = es.eigenvalues()[i];
}
return es.eigenvectors();
}
}
#endif
|