diff options
Diffstat (limited to 'src/algebra/odealg.spad.pamphlet')
| -rw-r--r-- | src/algebra/odealg.spad.pamphlet | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/algebra/odealg.spad.pamphlet b/src/algebra/odealg.spad.pamphlet index 32425c96..03901c7c 100644 --- a/src/algebra/odealg.spad.pamphlet +++ b/src/algebra/odealg.spad.pamphlet @@ -92,7 +92,7 @@ SystemODESolver(F, LO): Exports == Implementation where for sol in sols repeat m := m + count(#1 ~= 0, sol.basis)$List(F) SolMatrix:MF := new(n, m, 0) m := 0 - for sol in reverse_! sols repeat + for sol in reverse! sols repeat i := i+1 er := rec.eqs.i nn := #(er.g) -- size of active companionblock @@ -228,8 +228,8 @@ SystemODESolver(F, LO): Exports == Implementation where if (x(k, j) ~= 0) and ((rown = minr - 1) or degree x(k,j) < degree x(rown,j)) then rown := k rown = minr - 1 => "enuf" - x := swapRows_!(x, i, rown) - swap_!(w, i + offset, rown + offset) + x := swapRows!(x, i, rown) + swap!(w, i + offset, rown + offset) for k in i+1 .. nrows | x(k, j) ~= 0 repeat l := rightLcm(x(i,j), x(k,j)) a := rightQuotient(l, x(i, j)) |