| @@ 163-165 (lines=3) @@ | ||
| 160 | $this->V[$i][$i] = 1.0; |
|
| 161 | $h = $this->d[$i+1]; |
|
| 162 | if ($h != 0.0) { |
|
| 163 | for ($k = 0; $k <= $i; ++$k) { |
|
| 164 | $this->d[$k] = $this->V[$k][$i+1] / $h; |
|
| 165 | } |
|
| 166 | for ($j = 0; $j <= $i; ++$j) { |
|
| 167 | $g = 0.0; |
|
| 168 | for ($k = 0; $k <= $i; ++$k) { |
|
| @@ 168-170 (lines=3) @@ | ||
| 165 | } |
|
| 166 | for ($j = 0; $j <= $i; ++$j) { |
|
| 167 | $g = 0.0; |
|
| 168 | for ($k = 0; $k <= $i; ++$k) { |
|
| 169 | $g += $this->V[$k][$i+1] * $this->V[$k][$j]; |
|
| 170 | } |
|
| 171 | for ($k = 0; $k <= $i; ++$k) { |
|
| 172 | $this->V[$k][$j] -= $g * $this->d[$k]; |
|
| 173 | } |
|
| @@ 171-173 (lines=3) @@ | ||
| 168 | for ($k = 0; $k <= $i; ++$k) { |
|
| 169 | $g += $this->V[$k][$i+1] * $this->V[$k][$j]; |
|
| 170 | } |
|
| 171 | for ($k = 0; $k <= $i; ++$k) { |
|
| 172 | $this->V[$k][$j] -= $g * $this->d[$k]; |
|
| 173 | } |
|
| 174 | } |
|
| 175 | } |
|
| 176 | for ($k = 0; $k <= $i; ++$k) { |
|