@@ 330-339 (lines=10) @@ | ||
327 | $this->ort[$m] -= $g; |
|
328 | // Apply Householder similarity transformation |
|
329 | // H = (I -u * u' / h) * H * (I -u * u') / h) |
|
330 | for ($j = $m; $j < $this->n; ++$j) { |
|
331 | $f = 0.0; |
|
332 | for ($i = $high; $i >= $m; --$i) { |
|
333 | $f += $this->ort[$i] * $this->H[$i][$j]; |
|
334 | } |
|
335 | $f /= $h; |
|
336 | for ($i = $m; $i <= $high; ++$i) { |
|
337 | $this->H[$i][$j] -= $f * $this->ort[$i]; |
|
338 | } |
|
339 | } |
|
340 | for ($i = 0; $i <= $high; ++$i) { |
|
341 | $f = 0.0; |
|
342 | for ($j = $high; $j >= $m; --$j) { |
|
@@ 340-349 (lines=10) @@ | ||
337 | $this->H[$i][$j] -= $f * $this->ort[$i]; |
|
338 | } |
|
339 | } |
|
340 | for ($i = 0; $i <= $high; ++$i) { |
|
341 | $f = 0.0; |
|
342 | for ($j = $high; $j >= $m; --$j) { |
|
343 | $f += $this->ort[$j] * $this->H[$i][$j]; |
|
344 | } |
|
345 | $f = $f / $h; |
|
346 | for ($j = $m; $j <= $high; ++$j) { |
|
347 | $this->H[$i][$j] -= $f * $this->ort[$j]; |
|
348 | } |
|
349 | } |
|
350 | $this->ort[$m] = $scale * $this->ort[$m]; |
|
351 | $this->H[$m][$m-1] = $scale * $g; |
|
352 | } |