(Vec b)
| 86 | } |
| 87 | |
| 88 | public Vec solve(Vec b) |
| 89 | { |
| 90 | //Solve P A x = L U x = P b, for x |
| 91 | |
| 92 | //First solve L y = P b |
| 93 | Vec y = forwardSub(L, P.multiply(b)); |
| 94 | //Sole U x = y |
| 95 | Vec x = backSub(U, y); |
| 96 | |
| 97 | return x; |
| 98 | } |
| 99 | |
| 100 | public Matrix solve(Matrix B) |
| 101 | { |