Numerical Recipes Forum  

Go Back   Numerical Recipes Forum > Obsolete Editions Forum > Methods: Chapters 2, 11, and 18

Thread Tools Display Modes
Old 05-15-2009, 09:32 AM
Devinzhang Devinzhang is offline
Registered User
Join Date: May 2009
Posts: 5
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..n-1][0..n-1] 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..n-1][0..n-1] and the diagonal will be the eigenvalues after convergence(the matrix a is modified through rot anyway).

Reply With Quote

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is Off
HTML code is Off

Forum Jump

All times are GMT -5. The time now is 01:05 AM.

Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.