Project

General

Profile

Feature #1212

New function: GBasisByHomog

Added by Anna Maria Bigatti over 5 years ago. Updated about 1 month ago.

Status:
Closed
Priority:
Normal
Category:
New Function
Target version:
Start date:
05 Aug 2018
Due date:
% Done:

100%

Estimated time:
8.00 h
Spent time:

Description

Implement GBasis computation via homogenization


Related issues

Related to CoCoALib - Slug #1105: Primary Decompositon (zero-dim) slow casesClosed2017-10-02

Related to CoCoALib - Feature #1094: Revive code for SelfSaturating GBasisIn Progress2017-08-09

Related to CoCoA-5 - Slug #1114: Some other examples for 0-dim radicalClosed2017-10-31

Related to CoCoA-5 - Bug #1230: PrimaryDecomposition with lex orderingClosed2018-10-19

Related to CoCoALib - Bug #1379: Fails to recognize zero-dim idealClosed2019-12-27

Related to CoCoALib - Bug #1423: Radical 0-dim NYI? missing case for GBasisByHomogClosed2020-02-25

Related to CoCoALib - Feature #1434: GBasisByHomog for DegLexClosed2020-03-04

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

#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

Also available in: Atom PDF