Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
TMatrixDEigen.h
Go to the documentation of this file.
1// @(#)root/matrix:$Id$
2// Authors: Fons Rademakers, Eddy Offermann Dec 2003
3
4/*************************************************************************
5 * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6 * All rights reserved. *
7 * *
8 * For the licensing terms see $ROOTSYS/LICENSE. *
9 * For the list of contributors see $ROOTSYS/README/CREDITS. *
10 *************************************************************************/
11
12#ifndef ROOT_TMatrixDEigen
13#define ROOT_TMatrixDEigen
14
15//////////////////////////////////////////////////////////////////////////
16// //
17// TMatrixDEigen //
18// //
19// Eigenvalues and eigenvectors of a real matrix. //
20// //
21//////////////////////////////////////////////////////////////////////////
22
23#include "TMatrixD.h"
24#include "TVectorD.h"
25
27{
28protected :
29
30 static void MakeHessenBerg (TMatrixD &v,TVectorD &ortho,TMatrixD &H);
31 static void MakeSchurr (TMatrixD &v,TVectorD &d, TVectorD &e,TMatrixD &H);
32 static void Sort (TMatrixD &v,TVectorD &d, TVectorD &e);
33
34 TMatrixD fEigenVectors; // Eigen-vectors of matrix
35 TVectorD fEigenValuesRe; // Eigen-values
36 TVectorD fEigenValuesIm; // Eigen-values
37
38public :
39
40 enum {kWorkMax = 100}; // size of work array
41
44 TMatrixDEigen(const TMatrixD &a);
45 TMatrixDEigen(const TMatrixDEigen &another);
46 virtual ~TMatrixDEigen() {}
47
48// If matrix A has shape (rowLwb,rowUpb,rowLwb,rowUpb), then each eigen-vector
49// must have an index running between (rowLwb,rowUpb) .
50// For convenience, the column index of the eigen-vector matrix
51// also runs from rowLwb to rowUpb so that the returned matrix
52// has also index/shape (rowLwb,rowUpb,rowLwb,rowUpb) .
53// The same is true for the eigen-value vectors an matrix .
54
55 const TMatrixD &GetEigenVectors () const { return fEigenVectors; }
56 const TVectorD &GetEigenValuesRe() const { return fEigenValuesRe; }
57 const TVectorD &GetEigenValuesIm() const { return fEigenValuesIm; }
58 const TMatrixD GetEigenValues () const;
59
61
62 ClassDef(TMatrixDEigen,1) // Eigen-Vectors/Values of a Matrix
63};
64#endif
#define d(i)
Definition RSha256.hxx:102
#define a(i)
Definition RSha256.hxx:99
#define e(i)
Definition RSha256.hxx:103
#define ClassDef(name, id)
Definition Rtypes.h:342
TMatrixDEigen.
TVectorD fEigenValuesIm
static void MakeHessenBerg(TMatrixD &v, TVectorD &ortho, TMatrixD &H)
Nonsymmetric reduction to Hessenberg form.
TVectorD fEigenValuesRe
const TMatrixD & GetEigenVectors() const
virtual ~TMatrixDEigen()
TMatrixD fEigenVectors
const TMatrixD GetEigenValues() const
Computes the block diagonal eigenvalue matrix.
TMatrixDEigen & operator=(const TMatrixDEigen &source)
Assignment operator.
static void Sort(TMatrixD &v, TVectorD &d, TVectorD &e)
Sort eigenvalues and corresponding vectors in descending order of Re^2+Im^2 of the complex eigenvalue...
const TVectorD & GetEigenValuesRe() const
static void MakeSchurr(TMatrixD &v, TVectorD &d, TVectorD &e, TMatrixD &H)
Nonsymmetric reduction from Hessenberg to real Schur form.
const TVectorD & GetEigenValuesIm() const
#define H(x, y, z)