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/smith.spad.pamphlet | |
| parent | a7bab9a6c2070d05e2dbd256ce455079c8ced385 (diff) | |
| download | open-axiom-001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532.tar.gz | |
Replace `^=' with `~='.
Diffstat (limited to 'src/algebra/smith.spad.pamphlet')
| -rw-r--r-- | src/algebra/smith.spad.pamphlet | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/algebra/smith.spad.pamphlet b/src/algebra/smith.spad.pamphlet index 8c89d9ef..dd3ff8ca 100644 --- a/src/algebra/smith.spad.pamphlet +++ b/src/algebra/smith.spad.pamphlet @@ -99,7 +99,7 @@ SmithNormalForm(R,Row,Col,M) : Exports == Implementation where m1:= nrows m n1:= ncols m for i in 1..m1 repeat - for j in 1..n1 | (j ^= i) repeat + for j in 1..n1 | (j ~= i) repeat if not zero?(m(i,j)) then return false true @@ -111,14 +111,14 @@ SmithNormalForm(R,Row,Col,M) : Exports == Implementation where m -- elementary operation of second kind: add to row i-- - -- a*row j (i^=j) -- + -- a*row j (i~=j) -- elRow2(m : M,a:R,i:I,j:I) : M == vec:= map(a*#1,row(m,j)) vec:=map("+",row(m,i),vec) setRow!(m,i,vec) m -- elementary operation of second kind: add to column i -- - -- a*column j (i^=j) -- + -- a*column j (i~=j) -- elColumn2(m : M,a:R,i:I,j:I) : M == vec:= map(a*#1,column(m,j)) vec:=map("+",column(m,i),vec) @@ -157,7 +157,7 @@ SmithNormalForm(R,Row,Col,M) : Exports == Implementation where lastStep(sf : SmithForm) : SmithForm == m:=sf.Smith m1:=min(nrows m,ncols m) - for i in 1..m1 while (mii:=m(i,i)) ^=0 repeat + for i in 1..m1 while (mii:=m(i,i)) ~=0 repeat for j in i+1..m1 repeat if (m(j,j) exquo mii) case "failed" then return lastStep(ijDivide(sf,i,j)) |