#1




Question on Jacobi implementation in NR3
From page 573 to 575, some lines of code in implementation part are not so clear to me.
1. tracking the accumulated rotation of eigen vectors matrix v[0..n1][0..n1] should be as in equation 11.1.24, but the actual implementation is rot(v,s,tau,j,ip,j,iq). It is not clear to me how 11.1.24 is translated as rot(...) 2. I cannot see why we need vector b(n),d(n),z(n), since all accumulated rotation will be in the elements of matrix a[0..n1][0..n1] and the diagonal will be the eigenvalues after convergence(the matrix a is modified through rot anyway). Thanks! 
Thread Tools  
Display Modes  

