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[i1]=e[i]; 
...  ...  
301  301 
g=c*rb; 
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