Feature #1212
New function: GBasisByHomog
Description
Implement GBasis computation via homogenization
Related issues
History
#1 Updated by Anna Maria Bigatti over 5 years ago
- Related to Slug #1105: Primary Decompositon (zero-dim) slow cases added
#2 Updated by Anna Maria Bigatti over 5 years ago
- Related to Feature #1094: Revive code for SelfSaturating GBasis added
#3 Updated by Anna Maria Bigatti over 5 years ago
- Target version changed from CoCoALib-0.99600 to CoCoALib-0.99650 November 2019
Implemented for the case ordering is StdDegRevLex (SparsePolyOps-ideal.C
).
This is enough for PrimaryDecomposition slug #1105.
General case postponed.
#4 Updated by Anna Maria Bigatti over 5 years ago
- % Done changed from 10 to 60
#5 Updated by Anna Maria Bigatti over 5 years ago
- Subject changed from New function: GBasisViaHomog to New function: GBasisByHomog
#6 Updated by Anna Maria Bigatti over 5 years ago
- Related to Slug #1114: Some other examples for 0-dim radical added
#7 Updated by Anna Maria Bigatti over 5 years ago
- Related to Bug #1230: PrimaryDecomposition with lex ordering added
#8 Updated by John Abbott over 4 years ago
- Target version changed from CoCoALib-0.99650 November 2019 to CoCoALib-0.99700
#9 Updated by Anna Maria Bigatti over 4 years ago
- Related to Bug #1379: Fails to recognize zero-dim ideal added
#10 Updated by Anna Maria Bigatti over 4 years ago
I had implemented a trick which I believed would be enough for interreducing (compute the GB twice).
But it is not (see #1379): indeed it is not even minimal (with lex).
Use QQ[x,y,z], Lex; I := ideal(y*z -2*z^5 -2*z^4 +2*z^3, y^2 -16*z^5 -12*z^4 +12*z^3 +12*z^2, x*y -2*x*z -5*y +6*z^5 +2*z^3 -8*z^2 +10*z, x^2 +10*x*z -10*x -51*z^5 +74*z^4 +16*z^3 -14*z^2 -50*z +25, x*z^2 -x*z +5*z^5 -6*z^4 -4*z^3 +5*z, z^6 -z^5 -z^4 +z^2, z^3 -z -1); GBasisByHomog(I); --> [y -4*z -2, x -z^2, x*z -z -1, x^2 -z^2 -z, z^3 -z -1]
#11 Updated by Anna Maria Bigatti over 4 years ago
The trick didn't work ... they usually don't.
Fixed interreduction
#12 Updated by Anna Maria Bigatti over 4 years ago
checked-in. Made new pre-release for Mac.
#13 Updated by Anna Maria Bigatti over 4 years ago
- Estimated time changed from 5.00 h to 8.00 h
#14 Updated by Anna Maria Bigatti over 4 years ago
This works fine now, but only for few orderings.
The general case will not be done for this release. Postponing.
#15 Updated by Anna Maria Bigatti over 4 years ago
- Target version changed from CoCoALib-0.99700 to CoCoALib-0.99800
#16 Updated by Anna Maria Bigatti about 4 years ago
- Related to Bug #1423: Radical 0-dim NYI? missing case for GBasisByHomog added
#17 Updated by Anna Maria Bigatti about 4 years ago
Urgent: DegLex (should be easy)
#18 Updated by Anna Maria Bigatti about 4 years ago
- Related to Feature #1434: GBasisByHomog for DegLex added
#19 Updated by Anna Maria Bigatti about 4 years ago
- % Done changed from 60 to 70
Notes for me. In general: M matrix ordering.
If first row of M all positive
m m m 1
M 0
otherwise
1 1 1 1
1 1 1 0
M 0
and then call MakeTermOrd so make it square
#20 Updated by John Abbott about 2 years ago
- Target version changed from CoCoALib-0.99800 to CoCoALib-0.99850
#21 Updated by Anna Maria Bigatti about 1 month ago
- Status changed from In Progress to Closed
- % Done changed from 70 to 100
Good thing I wrote #1212-19! My first implementation was wrong (embarrassing).
Now all implemented and I added a CoCoA-5 test: test-GB.cocoa5