diff options
Diffstat (limited to 'src/algebra/intclos.spad.pamphlet')
| -rw-r--r-- | src/algebra/intclos.spad.pamphlet | 12 |
1 files changed, 4 insertions, 8 deletions
diff --git a/src/algebra/intclos.spad.pamphlet b/src/algebra/intclos.spad.pamphlet index 802f079d..5e5a13a5 100644 --- a/src/algebra/intclos.spad.pamphlet +++ b/src/algebra/intclos.spad.pamphlet @@ -158,8 +158,7 @@ IntegralBasisTools(R,UP,F): Exports == Implementation where for i in 1..n repeat for j in i..n repeat if not zero?(mij := qelt(mat,i,j)) then d := gcd(d,mij) --- one? d => return d - (d = 1) => return d + one? d => return d d divideIfCan_!(matrix,matrixOut,prime,n) == @@ -654,8 +653,7 @@ NumberFieldIntegralBasis(UP,F): Exports == Implementation where -- these are wrt the original basis for F runningRbden : I := 1 -- runningRbden = denominator for current basis matrix --- one? sing and empty? wilds => [runningRb, runningRbden, runningRbinv] - (sing = 1) and empty? wilds => [runningRb, runningRbden, runningRbinv] + one? sing and empty? wilds => [runningRb, runningRbden, runningRbinv] -- id = basis matrix of the ideal (p-radical) wrt current basis matrixOut : Mat := scalarMatrix(n,0) for p in wilds repeat @@ -721,8 +719,7 @@ NumberFieldIntegralBasis(UP,F): Exports == Implementation where rbinv := UpTriBddDenomInv(rb, rbden) disc := disc0 quo (index * index) indexChange := index quo oldIndex; oldIndex := index --- one? indexChange => return [rb, rbden, rbinv, disc] - (indexChange = 1) => return [rb, rbden, rbinv, disc] + one? indexChange => return [rb, rbden, rbinv, disc] tfm := ((rb * traceMat * transpose rb) exquo (rbden * rbden)) :: Mat iWildLocalIntegralBasis(matrixOut,disc,p) == @@ -755,8 +752,7 @@ NumberFieldIntegralBasis(UP,F): Exports == Implementation where rbinv := UpTriBddDenomInv(rb, rbden) indexChange := index quo oldIndex; oldIndex := index disc := disc quo (indexChange * indexChange) --- one? indexChange or gcd(p2,disc) ~= p2 => - (indexChange = 1) or gcd(p2,disc) ~= p2 => + one? indexChange or gcd(p2,disc) ~= p2 => return [rb, rbden, rbinv, disc] discriminant() == |