@@ 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 { |