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