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