@@ 280-282 (lines=3) @@ | ||
277 | $this->V[$i][$i] = 1.0; |
|
278 | $h = $this->d[$i + 1]; |
|
279 | if ($h != 0.0) { |
|
280 | for ($k = 0; $k <= $i; ++$k) { |
|
281 | $this->d[$k] = $this->V[$k][$i + 1] / $h; |
|
282 | } |
|
283 | ||
284 | for ($j = 0; $j <= $i; ++$j) { |
|
285 | $g = 0.0; |
|
@@ 286-288 (lines=3) @@ | ||
283 | ||
284 | for ($j = 0; $j <= $i; ++$j) { |
|
285 | $g = 0.0; |
|
286 | for ($k = 0; $k <= $i; ++$k) { |
|
287 | $g += $this->V[$k][$i + 1] * $this->V[$k][$j]; |
|
288 | } |
|
289 | ||
290 | for ($k = 0; $k <= $i; ++$k) { |
|
291 | $this->V[$k][$j] -= $g * $this->d[$k]; |
|
@@ 290-292 (lines=3) @@ | ||
287 | $g += $this->V[$k][$i + 1] * $this->V[$k][$j]; |
|
288 | } |
|
289 | ||
290 | for ($k = 0; $k <= $i; ++$k) { |
|
291 | $this->V[$k][$j] -= $g * $this->d[$k]; |
|
292 | } |
|
293 | } |
|
294 | } |
|
295 |