@@ 461-470 (lines=10) @@ | ||
458 | $this->ort[$m] -= $g; |
|
459 | // Apply Householder similarity transformation |
|
460 | // H = (I -u * u' / h) * H * (I -u * u') / h) |
|
461 | for ($j = $m; $j < $this->n; ++$j) { |
|
462 | $f = 0.0; |
|
463 | for ($i = $high; $i >= $m; --$i) { |
|
464 | $f += $this->ort[$i] * $this->H[$i][$j]; |
|
465 | } |
|
466 | ||
467 | $f /= $h; |
|
468 | for ($i = $m; $i <= $high; ++$i) { |
|
469 | $this->H[$i][$j] -= $f * $this->ort[$i]; |
|
470 | } |
|
471 | } |
|
472 | ||
473 | for ($i = 0; $i <= $high; ++$i) { |
|
@@ 473-482 (lines=10) @@ | ||
470 | } |
|
471 | } |
|
472 | ||
473 | for ($i = 0; $i <= $high; ++$i) { |
|
474 | $f = 0.0; |
|
475 | for ($j = $high; $j >= $m; --$j) { |
|
476 | $f += $this->ort[$j] * $this->H[$i][$j]; |
|
477 | } |
|
478 | ||
479 | $f = $f / $h; |
|
480 | for ($j = $m; $j <= $high; ++$j) { |
|
481 | $this->H[$i][$j] -= $f * $this->ort[$j]; |
|
482 | } |
|
483 | } |
|
484 | ||
485 | $this->ort[$m] = $scale * $this->ort[$m]; |