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