80 for (
j = 0;
j <
n;
j++)
85 for (i =
n-1; i > 0; i--) {
93 for (k = 0; k < i; k++)
97 for (
j = 0;
j < i;
j++) {
107 for (k = 0; k < i; k++) {
118 for (
j = 0;
j < i;
j++)
123 for (
j = 0;
j < i;
j++) {
128 for (k =
j+1; k <= i-1; k++) {
136 for (
j = 0;
j < i;
j++) {
141 for (
j = 0;
j < i;
j++)
143 for (
j = 0;
j < i;
j++) {
146 for (k =
j; k <= i-1; k++) {
159 for (i = 0; i <
n-1; i++) {
165 for (k = 0; k <= i; k++) {
169 for (
j = 0;
j <= i;
j++) {
171 for (k = 0; k <= i; k++) {
175 for (k = 0; k <= i; k++) {
181 for (k = 0; k <= i; k++) {
186 for (
j = 0;
j <
n;
j++) {
209 for (i = 1; i <
n; i++)
216 for (
l = 0;
l <
n;
l++) {
238 Error(
"MakeEigenVectors",
"too many iterations");
253 for (i =
l+2; i <
n; i++)
266 for (i =
m-1; i >=
l; i--) {
281 for (k = 0; k <
n; k++) {
302 for (i = 0; i <
n-1; i++) {
305 for (
j = i+1;
j <
n;
j++) {
314 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.