Numerically Stable Computation
To determine the quantities needed for the update, we must solve the off-diagonal equation for zero (Golub & Van Loan 1996, §8.4). This implies that
Set β to half this quantity,
If akℓ is zero we can stop without performing an update, thus we never divide by zero. Let t be tan ϑ. Then with a few trigonometric identities we reduce the equation to
For stability we choose the solution
From this we may obtain c and s as
Although we now could use the algebraic update equations given previously, it may be preferable to rewrite them. Let
so that ρ = tan(ϑ/2). Then the revised update equations are
As previously remarked, we need never explicitly compute the rotation angle ϑ. In fact, we can reproduce the symmetric update determined by Qkℓ by retaining only the three values k, ℓ, and t, with t set to zero for a null rotation.
Read more about this topic: Jacobi Rotation
Famous quotes containing the words numerically, stable and/or computation:
“... in the happy laughter of a theatre audience one can get the most immediate and numerically impressive guarantee that there is nothing in one’s mind which is not familiar to the mass of persons living at the time.”
—Rebecca West (1892–1983)
“In verity ... we are the poor. This humanity we would claim for ourselves is the legacy, not only of the Enlightenment, but of the thousands and thousands of European peasants and poor townspeople who came here bringing their humanity and their sufferings with them. It is the absence of a stable upper class that is responsible for much of the vulgarity of the American scene. Should we blush before the visitor for this deficiency?”
—Mary McCarthy (1912–1989)
“I suppose that Paderewski can play superbly, if not quite at his best, while his thoughts wander to the other end of the world, or possibly busy themselves with a computation of the receipts as he gazes out across the auditorium. I know a great actor, a master technician, can let his thoughts play truant from the scene ...”
—Minnie Maddern Fiske (1865–1932)