primdecZ() fails for ideal I=1

[kroeker@…]

minimal failing example:

ring rng = integer,x,dp;
ideal I=1;
minAssZ(I);
primdecZ(I); //error: "identifier A not found in Rhelp"

[kroeker@…]

forgot to load "primdecint.lib" in example. Updated failing example:

LIB("primdecint.lib");
ring rng = integer,x,dp;
ideal I=1;
minAssZ(I);
primdecZ(I); //error: "identifier A not found in Rhelp"

[Oleksandr]

i confirmed this report with the current SW. The output is as follows:

>  minAssZ(I);
empty list
>  primdecZ(I); //error: "identifier A not found in Rhelp"
// ** redefining A **
   ? identifier A not found in Rhelp
   ? error occurred in or before primdecint.lib::primdecZ line 180: `      list A = imap(Rhelp,A);`
   ? leaving primdecint.lib::primdecZ
   skipping text from `;` error at token `)`

[anonymous]

The list A is initialized in primdecZ as a list without assigning a ring-related entry. Hence the list does not belong to the ring, but to the ring-independent variables as long as no ring-related entry is assigned.

Suggested Bug-Fix: list A=ideal(0); in line 97 (formerly list A;) and then make sure to always skip first entry of list A. This glues it to the ring Rhelp.

Best regards, Anne

[kroeker@…]

see pull request ​https://github.com/Singular/Sources/pull/514