/usr/include/Cauchy/Eigen3/Norm.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 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 | /*
* 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_NORM_H_
#define _CAUCHY_NORM_H_
#include <Eigen/SVD>
#include "Infinity.h"
namespace Cauchy {
double norm(const Eigen::MatrixXd& m, double p = 2)
{
if (p == 2)
{
Eigen::JacobiSVD<Eigen::MatrixXd> svd(m);
return svd.singularValues()(0,0);
}
else if (p == 1) {
double max = 0.0;
for(int i = 0; i < m.cols(); ++i)
{
double cn = m.col(i).lpNorm<1>();
if(cn > max) max = cn;
}
return max;
} else if ( p == Inf) {
double max = 0.0;
for(int i = 0; i < m.rows(); ++i)
{
double cn = m.row(i).lpNorm<1>();
if(cn > max) max = cn;
}
return max;
} else if (p > 1)
{
// TODO unimplemented (but also unimplemented in matlab)
}
// TODO assert on value p < 1
return -1;
}
}
#endif
|