| @@ 633-637 (lines=5) @@ | ||
| 630 | $p = $p / $r; |
|
| 631 | $q = $q / $r; |
|
| 632 | // Row modification |
|
| 633 | for ($j = $n - 1; $j < $nn; ++$j) { |
|
| 634 | $z = $this->H[$n - 1][$j]; |
|
| 635 | $this->H[$n - 1][$j] = $q * $z + $p * $this->H[$n][$j]; |
|
| 636 | $this->H[$n][$j] = $q * $this->H[$n][$j] - $p * $z; |
|
| 637 | } |
|
| 638 | ||
| 639 | // Column modification |
|
| 640 | for ($i = 0; $i <= $n; ++$i) { |
|
| @@ 640-644 (lines=5) @@ | ||
| 637 | } |
|
| 638 | ||
| 639 | // Column modification |
|
| 640 | for ($i = 0; $i <= $n; ++$i) { |
|
| 641 | $z = $this->H[$i][$n - 1]; |
|
| 642 | $this->H[$i][$n - 1] = $q * $z + $p * $this->H[$i][$n]; |
|
| 643 | $this->H[$i][$n] = $q * $this->H[$i][$n] - $p * $z; |
|
| 644 | } |
|
| 645 | ||
| 646 | // Accumulate transformations |
|
| 647 | for ($i = $low; $i <= $high; ++$i) { |
|
| @@ 647-651 (lines=5) @@ | ||
| 644 | } |
|
| 645 | ||
| 646 | // Accumulate transformations |
|
| 647 | for ($i = $low; $i <= $high; ++$i) { |
|
| 648 | $z = $this->V[$i][$n - 1]; |
|
| 649 | $this->V[$i][$n - 1] = $q * $z + $p * $this->V[$i][$n]; |
|
| 650 | $this->V[$i][$n] = $q * $this->V[$i][$n] - $p * $z; |
|
| 651 | } |
|
| 652 | ||
| 653 | // Complex pair |
|
| 654 | } else { |
|