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