php-ai /
php-ml
| @@ 778-786 (lines=9) @@ | ||
| 775 | $q = $q / $p; |
|
| 776 | $r = $r / $p; |
|
| 777 | // Row modification |
|
| 778 | for ($j = $k; $j < $nn; ++$j) { |
|
| 779 | $p = $this->H[$k][$j] + $q * $this->H[$k + 1][$j]; |
|
| 780 | if ($notlast) { |
|
| 781 | $p = $p + $r * $this->H[$k + 2][$j]; |
|
| 782 | $this->H[$k + 2][$j] = $this->H[$k + 2][$j] - $p * $z; |
|
| 783 | } |
|
| 784 | ||
| 785 | $this->H[$k][$j] = $this->H[$k][$j] - $p * $x; |
|
| 786 | $this->H[$k + 1][$j] = $this->H[$k + 1][$j] - $p * $y; |
|
| 787 | } |
|
| 788 | ||
| 789 | // Column modification |
|
| @@ 790-798 (lines=9) @@ | ||
| 787 | } |
|
| 788 | ||
| 789 | // Column modification |
|
| 790 | for ($i = 0; $i <= min($n, $k + 3); ++$i) { |
|
| 791 | $p = $x * $this->H[$i][$k] + $y * $this->H[$i][$k + 1]; |
|
| 792 | if ($notlast) { |
|
| 793 | $p = $p + $z * $this->H[$i][$k + 2]; |
|
| 794 | $this->H[$i][$k + 2] = $this->H[$i][$k + 2] - $p * $r; |
|
| 795 | } |
|
| 796 | ||
| 797 | $this->H[$i][$k] = $this->H[$i][$k] - $p; |
|
| 798 | $this->H[$i][$k + 1] = $this->H[$i][$k + 1] - $p * $q; |
|
| 799 | } |
|
| 800 | ||
| 801 | // Accumulate transformations |
|
| @@ 802-810 (lines=9) @@ | ||
| 799 | } |
|
| 800 | ||
| 801 | // Accumulate transformations |
|
| 802 | for ($i = $low; $i <= $high; ++$i) { |
|
| 803 | $p = $x * $this->V[$i][$k] + $y * $this->V[$i][$k + 1]; |
|
| 804 | if ($notlast) { |
|
| 805 | $p = $p + $z * $this->V[$i][$k + 2]; |
|
| 806 | $this->V[$i][$k + 2] = $this->V[$i][$k + 2] - $p * $r; |
|
| 807 | } |
|
| 808 | ||
| 809 | $this->V[$i][$k] = $this->V[$i][$k] - $p; |
|
| 810 | $this->V[$i][$k + 1] = $this->V[$i][$k + 1] - $p * $q; |
|
| 811 | } |
|
| 812 | } // ($s != 0) |
|
| 813 | } // k loop |
|
