/usr/include/root/TMatrixDSymEigen.h is in libroot-math-matrix-dev 5.34.19+dfsg-1.2.
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 63 64 65 66 | // @(#)root/matrix:$Id$
// Authors: Fons Rademakers, Eddy Offermann Dec 2003
/*************************************************************************
* Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
* All rights reserved. *
* *
* For the licensing terms see $ROOTSYS/LICENSE. *
* For the list of contributors see $ROOTSYS/README/CREDITS. *
*************************************************************************/
#ifndef ROOT_TMatrixDSymEigen
#define ROOT_TMatrixDSymEigen
//////////////////////////////////////////////////////////////////////////
// //
// TMatrixDSymEigen //
// //
// Eigenvalues and eigenvectors of a real symmetric matrix. //
// //
//////////////////////////////////////////////////////////////////////////
#ifndef ROOT_TMatrixD
#include "TMatrixD.h"
#endif
#ifndef ROOT_TMatrixDSym
#include "TMatrixDSym.h"
#endif
#ifndef ROOT_TVectorD
#include "TVectorD.h"
#endif
class TMatrixDSymEigen
{
protected :
static void MakeTridiagonal (TMatrixD &v,TVectorD &d,TVectorD &e);
static void MakeEigenVectors(TMatrixD &v,TVectorD &d,TVectorD &e);
TMatrixD fEigenVectors; // Eigen-vectors of matrix
TVectorD fEigenValues; // Eigen-values
public :
enum {kWorkMax = 100}; // size of work array
TMatrixDSymEigen() : fEigenVectors(), fEigenValues() {};
TMatrixDSymEigen(const TMatrixDSym &a);
TMatrixDSymEigen(const TMatrixDSymEigen &another);
virtual ~TMatrixDSymEigen() {}
// If matrix A has shape (rowLwb,rowUpb,rowLwb,rowUpb), then each eigen-vector
// must have an index running between (rowLwb,rowUpb) .
// For convenience, the column index of the eigen-vector matrix
// also runs from rowLwb to rowUpb so that the returned matrix
// has also index/shape (rowLwb,rowUpb,rowLwb,rowUpb) .
// The same is true for the eigen-value vector .
const TMatrixD &GetEigenVectors() const { return fEigenVectors; }
const TVectorD &GetEigenValues () const { return fEigenValues; }
TMatrixDSymEigen &operator= (const TMatrixDSymEigen &source);
ClassDef(TMatrixDSymEigen,1) // Eigen-Vectors/Values of a Matrix
};
#endif
|