81 for (
j = 0;
j <
n;
j++)
86 for (i =
n-1; i > 0; i--) {
94 for (k = 0; k < i; k++)
98 for (
j = 0;
j < i;
j++) {
108 for (k = 0; k < i; k++) {
119 for (
j = 0;
j < i;
j++)
124 for (
j = 0;
j < i;
j++) {
129 for (k =
j+1; k <= i-1; k++) {
137 for (
j = 0;
j < i;
j++) {
142 for (
j = 0;
j < i;
j++)
144 for (
j = 0;
j < i;
j++) {
147 for (k =
j; k <= i-1; k++) {
160 for (i = 0; i <
n-1; i++) {
166 for (k = 0; k <= i; k++) {
170 for (
j = 0;
j <= i;
j++) {
172 for (k = 0; k <= i; k++) {
176 for (k = 0; k <= i; k++) {
182 for (k = 0; k <= i; k++) {
187 for (
j = 0;
j <
n;
j++) {
210 for (i = 1; i <
n; i++)
217 for (
l = 0;
l <
n;
l++) {
239 Error(
"MakeEigenVectors",
"too many iterations");
254 for (i =
l+2; i <
n; i++)
267 for (i =
m-1; i >=
l; i--) {
282 for (k = 0; k <
n; k++) {
303 for (i = 0; i <
n-1; i++) {
306 for (
j = i+1;
j <
n;
j++) {
315 for (
j = 0;
j <
n;
j++) {
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
#define R__ASSERT(e)
Checks condition e and reports a fatal error if it's false.
void Error(const char *location, const char *msgfmt,...)
Use this function in case an error occurred.
winID h TVirtualViewer3D TVirtualGLPainter p
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t r
static void MakeEigenVectors(TMatrixD &v, TVectorD &d, TVectorD &e)
Symmetric tridiagonal QL algorithm.
static void MakeTridiagonal(TMatrixD &v, TVectorD &d, TVectorD &e)
This is derived from the Algol procedures tred2 by Bowdler, Martin, Reinsch, and Wilkinson,...
TMatrixDSymEigen & operator=(const TMatrixDSymEigen &source)
Assignment operator.
TMatrixTBase< Element > & ResizeTo(Int_t nrows, Int_t ncols, Int_t=-1) override
Set size of the matrix to nrows x ncols New dynamic elements are created, the overlapping part of the...
TVectorT< Element > & ResizeTo(Int_t lwb, Int_t upb)
Resize the vector to [lwb:upb] .
Short_t Max(Short_t a, Short_t b)
Returns the largest of a and b.
Double_t Sqrt(Double_t x)
Returns the square root of x.
LongDouble_t Power(LongDouble_t x, LongDouble_t y)
Returns x raised to the power y.
Double_t Hypot(Double_t x, Double_t y)
Returns sqrt(x*x + y*y)
Short_t Abs(Short_t d)
Returns the absolute value of parameter Short_t d.