diff --git a/matrix/src/TDecompBase.cxx b/matrix/src/TDecompBase.cxx index 939f349266c61eb667084e9ea3fb31554b6be02b..6b0c396ca214230859cd2f2cbb1e0f7cbfe52c71 100644 --- a/matrix/src/TDecompBase.cxx +++ b/matrix/src/TDecompBase.cxx @@ -1,4 +1,4 @@ -// @(#)root/matrix:$Name: $:$Id: TDecompBase.cxx,v 1.1 2004/01/25 20:33:32 brun Exp $ +// @(#)root/matrix:$Name: $:$Id: TDecompBase.cxx,v 1.2 2004/01/26 12:08:58 brun Exp $ // Authors: Fons Rademakers, Eddy Offermann Dec 2003 /************************************************************************* @@ -385,12 +385,13 @@ void ApplyHouseHolder(const TVectorD &vc,Double_t up,Double_t beta, Double_t * cp = cr.GetPtr(); Double_t s = cp[lp*inc_c]*up; - for (Int_t i = l; i < nv; i++) + Int_t i; + for (i = l; i < nv; i++) s += cp[i*inc_c]*vp[i]; s = s*beta; cp[lp*inc_c] += s*up; - for (Int_t i = l; i < nv; i++) + for (i = l; i < nv; i++) cp[i*inc_c] += s*vp[i]; } diff --git a/matrix/src/TDecompSVD.cxx b/matrix/src/TDecompSVD.cxx index 4664f2edf76c281d2cd2c859d29dac05f5226747..6e1ea6de4c91a01463ef44eeceab2e9691dfc363 100644 --- a/matrix/src/TDecompSVD.cxx +++ b/matrix/src/TDecompSVD.cxx @@ -1,4 +1,4 @@ -// @(#)root/matrix:$Name: $:$Id: TDecompSVD.cxx,v 1.1 2004/01/25 20:33:32 brun Exp $ +// @(#)root/matrix:$Name: $:$Id: TDecompSVD.cxx,v 1.2 2004/01/26 12:08:58 brun Exp $ // Authors: Fons Rademakers, Eddy Offermann Dec 2003 /************************************************************************* @@ -132,7 +132,7 @@ Bool_t TDecompSVD::Bidiagonalize(TMatrixD &v,TMatrixD &u,TVectorD &sDiag,TVector TArrayD ups(nCol_v); TArrayD betas(nCol_v); - Int_t i; + Int_t i,j; for (i = 0; i < nCol_v; i++) { // Set up Householder Transformation q(i) @@ -144,13 +144,13 @@ Bool_t TDecompSVD::Bidiagonalize(TMatrixD &v,TMatrixD &u,TVectorD &sDiag,TVector DefHouseHolder(vc_i,i,i+1,up,beta); // Apply q(i) to v - for (Int_t j = i; j < nCol_v; j++) { + for (j = i; j < nCol_v; j++) { TMatrixDColumn vc_j = TMatrixDColumn(v,j); ApplyHouseHolder(vc_i,up,beta,i,i+1,vc_j); } // Apply q(i) to u - for (Int_t j = 0; j < nCol_u; j++) + for (j = 0; j < nCol_u; j++) { TMatrixDColumn uc_j = TMatrixDColumn(u,j); ApplyHouseHolder(vc_i,up,beta,i,i+1,uc_j); @@ -169,7 +169,7 @@ Bool_t TDecompSVD::Bidiagonalize(TMatrixD &v,TMatrixD &u,TVectorD &sDiag,TVector betas[i] = beta; // apply h(i) to v - for (Int_t j = i; j < nRow_v; j++) { + for (j = i; j < nRow_v; j++) { TMatrixDRow vr_j = TMatrixDRow(v,j); ApplyHouseHolder(vr_i,up,beta,i+1,i+2,vr_j); diff --git a/matrix/src/TMatrixDSymEigen.cxx b/matrix/src/TMatrixDSymEigen.cxx index 5b335df64045b1b13a9f67c0c74681ff76416cec..ed460189d7e14504c5ac20bafed62f0403bace73 100644 --- a/matrix/src/TMatrixDSymEigen.cxx +++ b/matrix/src/TMatrixDSymEigen.cxx @@ -1,4 +1,4 @@ -// @(#)root/matrix:$Name: $:$Id: TMatrixDSymEigen.cxx,v 1.1 2004/01/25 20:33:32 brun Exp $ +// @(#)root/matrix:$Name: $:$Id: TMatrixDSymEigen.cxx,v 1.2 2004/01/26 12:08:58 brun Exp $ // Authors: Fons Rademakers, Eddy Offermann Dec 2003 /************************************************************************* @@ -65,13 +65,13 @@ void TMatrixDSymEigen::MakeTridiagonal(TMatrixD &v,TVectorD &d,TVectorD &e) const Int_t n = v.GetNrows(); + Int_t i,j,k; Int_t off_n1 = (n-1)*n; - for (Int_t j = 0; j < n; j++) + for (j = 0; j < n; j++) pD[j] = pV[off_n1+j]; // Householder reduction to tridiagonal form. - Int_t i,j,k; for (i = n-1; i > 0; i--) { const Int_t off_i1 = (i-1)*n; const Int_t off_i = i*n;