From: Jack Wu
Date: Sun, 18 Dec 2011 02:06:11 -0800 (PST)
I have a question about the definition of 'INDXQ' in routine SLAED1.f. It it
so hard for me to comprehend what you mean by ''the permutation which
separately sorts the two subproblems in D into ascending order''!
Let's take a simple example: suppose D = diag([1 2 3 4]), u = Z = transpose([0
1 1 0]), i.e. Q(in) = identity matrix and RHO = 1, what 'INDXQ' should be for
matrix X = D + u*u^T = diag([1 2 3 4]) + transpose([0 1 1 0])*[0 1 1 0]?
Thank you very much for your attention and look forward to your reply!
Jie Wu, Jack
Graduate Student of Applied Physics
Physical Measurement and Control Group (http://minty.stanford.edu/)