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