diff options
| author | dos-reis <gdr@axiomatics.org> | 2008-04-03 04:23:42 +0000 |
|---|---|---|
| committer | dos-reis <gdr@axiomatics.org> | 2008-04-03 04:23:42 +0000 |
| commit | 001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532 (patch) | |
| tree | da9e2fe5d81ff4cd7709d12e44b8c3e348b8a8e3 /src/algebra/pseudolin.spad.pamphlet | |
| parent | a7bab9a6c2070d05e2dbd256ce455079c8ced385 (diff) | |
| download | open-axiom-001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532.tar.gz | |
Replace `^=' with `~='.
Diffstat (limited to 'src/algebra/pseudolin.spad.pamphlet')
| -rw-r--r-- | src/algebra/pseudolin.spad.pamphlet | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/algebra/pseudolin.spad.pamphlet b/src/algebra/pseudolin.spad.pamphlet index 3eb384b4..fb6d38a5 100644 --- a/src/algebra/pseudolin.spad.pamphlet +++ b/src/algebra/pseudolin.spad.pamphlet @@ -96,19 +96,19 @@ PseudoLinearNormalForm(K:Field): Exports == Implementation where while j <= N and M(i, j) = 0 repeat j := j + 1 if j <= N then -- expand companionblock by lemma 5 - if j ^= i+1 then + if j ~= i+1 then -- perform first a permutation E := permutationMatrix(N, i+1, j) M := changeBase(M, E, sig, der) B := B*E Binv := E*Binv - -- now is M(i, i+1) ^= 0 + -- now is M(i, i+1) ~= 0 E := mulMatrix(N, i+1, siginv inv M(i,i+1)) M := changeBase(M, E, sig, der) B := B*E Binv := inv(E)*Binv for j in 1..N repeat - if j ^= i+1 then + if j ~= i+1 then E := addMatrix(N, i+1, j, siginv(-M(i,j))) M := changeBase(M, E, sig, der) B := B*E |