Revision 1264 trunk/src/linad99/eigenv.cpp

eigenv.cpp (revision 1264)
255 255
  dvector& e = (dvector&) _e;
256 256
  dmatrix& z = (dmatrix&) _z;
257 257
  int n=d.size();
258
  int m,l,iter,i,k;
258
  int m,l,iter,i;
259 259
  double s,r,p,g,f,dd,c,b;
260 260

  
261 261
  for (i=2;i<=n;i++) e[i-1]=e[i];
......
301 301
          g=c*r-b;
302 302
          /* Next loop can be omitted if eigenvectors not wanted */
303 303
          #ifdef EIGEN_VECTORS
304
            for (k=1;k<=n;k++)
304
            for (int k=1;k<=n;k++)
305 305
            {
306 306
              f=z[k][i+1];
307 307
              z[k][i+1]=s*z[k][i]+c*f;

Also available in: Unified diff