| @@ 89-91 (lines=3) @@ | ||
| 86 | // Outer loop. |
|
| 87 | for ($j = 0; $j < $this->n; ++$j) { |
|
| 88 | // Make a copy of the j-th column to localize references. |
|
| 89 | for ($i = 0; $i < $this->m; ++$i) { |
|
| 90 | $LUcolj[$i] = &$this->LU[$i][$j]; |
|
| 91 | } |
|
| 92 | // Apply previous transformations. |
|
| 93 | for ($i = 0; $i < $this->m; ++$i) { |
|
| 94 | $LUrowi = $this->LU[$i]; |
|
| @@ 123-125 (lines=3) @@ | ||
| 120 | } |
|
| 121 | // Compute multipliers. |
|
| 122 | if (($j < $this->m) && ($this->LU[$j][$j] != 0.0)) { |
|
| 123 | for ($i = $j+1; $i < $this->m; ++$i) { |
|
| 124 | $this->LU[$i][$j] /= $this->LU[$j][$j]; |
|
| 125 | } |
|
| 126 | } |
|
| 127 | } |
|
| 128 | } // function __construct() |
|