search

jakobkroeker / test_singular

manage a personal issue list for Singular CAS
0 stars 0 forks source link

incorrect primdecGTZ() example over prime coefficient field #105

Open jakobkroeker opened 9 years ago

jakobkroeker commented 9 years ago

origin: https://github.com/surface-smoothers/lipmanresolve/issues/382

fails on andromeda: not reproducible on groebner not reproducible at home (fc17 32 bit):

Singular for i686-Linux version 4.0.1 (4010, 32 bit) Oct 16 2014 01:18:32 #89d78ad
with
    GMP(5.0.5),NTL(5.5.2),factory(@(#) factoryVersion = 4.0.1),
    static readline,Plural,DBM,
    dynamic modules,dynamic p_Procs,OM_NDEBUG,SING_NDEBUG,CC:OPTIMIZE,eigenvalues,Gauss-Manin system,random=1413547431
built-in modules: {}
AC_CONFIGURE_ARGS =  '--prefix=/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/release' '--enable-gfanlib' '--enable-omalloc' 'OMALLOC_LIBS=/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/omalloc/libomalloc.la' 'OMALLOC_INCLUDES=-I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/.. -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD' '--with-Singular' 'RESOURCES_LIBS=/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/resources/libresources.la' 'RESOURCES_INCLUDES=-I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/.. ' 'FACTORY_LIBS=/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/factory/libfactory.la' 'FACTORY_INCLUDES=-I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/.. -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/../factory/include -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/factory/include'  OMALLOC_LIBS='/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/omalloc/libomalloc.la' OMALLOC_INCLUDES='-I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/.. -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD' --with-Singular RESOURCES_LIBS='/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/resources/libresources.la' RESOURCES_INCLUDES='-I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/.. ' FACTORY_LIBS='/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/factory/libfactory.la' FACTORY_INCLUDES='-I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/.. -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/../factory/include -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD -I/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/factory/include',
CC = gcc,FLAGS : -g -O2  -pipe -fno-common -O3 -Wno-unused-function -Wno-trigraphs -Wno-unused-parameter -Wno-unused-variable -fomit-frame-pointer -fwrapv -fvisibility=default -finline-functions -fno-exceptions -funroll-loops ,
CXX = g++,FLAGS :  -pipe -fno-common -O3 -Wno-unused-function -Wno-trigraphs -Wno-unused-parameter -Wno-unused-variable -fomit-frame-pointer -fwrapv -fvisibility=default -finline-functions -fno-exceptions -fno-rtti -fno-threadsafe-statics -fno-enforce-eh-specs -fconserve-space -funroll-loops,
DEFS : ,CPPFLAGS : ,
LDFLAGS :  -pipe -fno-common -O3 -Wno-unused-function -Wno-trigraphs -Wno-unused-parameter -Wno-unused-variable -fomit-frame-pointer -fwrapv -fvisibility=default -finline-functions -fno-exceptions -funroll-loops  -rdynamic -Wl,-undefined,dynamic_lookup,LIBS : -lreadline -lncurses -lrt -lpthread  (ver: 4.7.2 20121109 (Red Hat 4.7.2-8))
argv[0]   : ./singular-spielwiese
SearchPath: /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/LIB:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/../factory:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/MOD:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/release/libexec/singular/MOD:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular
Singular  : /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/Singular
BinDir    : /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular
ProcDir   : /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/MOD:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/release/libexec/singular/MOD
RootDir   : /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/..
DataDir   : 
DefaultDir: /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/release
InfoFile  : 
IdxFile   : 
HtmlDir   : 
ManualUrl : http://www.singular.uni-kl.de/Manual/
ExDir     : 
Path      : /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/MOD:/home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/release/libexec/singular/MOD:/usr/lib/qt-3.3/bin:/usr/lib/ccache:/usr/local/bin:/usr/bin:/bin
emacs     : /usr/bin/emacs-24.2
xemacs    : /usr/bin/xemacs-21.5-b33
SingularEmacs:  /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/Singular
EmacsLoad : 
EmacsDir  : 
SingularXterm:  /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/Singular
xterm     : /usr/bin/xterm
EmacsDir  : 
Available HelpBrowsers: firefox-www, konqueror-www, lynx-www, dummy, emacs, 
Current HelpBrowser: firefox-www 
                     SINGULAR                                 /  Development
 A Computer Algebra System for Polynomial Computations       /   version 4.0.1
                                                           0<
 by: W. Decker, G.-M. Greuel, G. Pfister, H. Schoenemann     \   Sep 2014
FB Mathematik der Universitaet, D-67653 Kaiserslautern        \
// ** executing /home/jakob/Projects/lipmanresolve/Singular-spielwiese/BUILD/Singular/LIB/.singularrc

option(redSB);
LIB("wrappers.lib");
LIB("normal.lib"); 

system("random",1467556385);
// system("random",158784614 );

ring rng = (3),(x,y,z),(dp(2),dp(1),C);
ideal J = -x2yz2+z-1,xy2z+z,x3z-xy2+1;

def L0 = primdecGTZ(J);
dbprint(1, "L0",L0);

proc level3 ()
{
   dbprint(1,"voice",voice);

    def L = primdecGTZ(J);
    dbprint(1, "L",L);

    ideal a1 = std(L[2][2]);
    a1;
    ideal a2 = std( radical(L[2][2]) );
    a2;
    idealsEqual(a1,a2); // = 0!

}

proc level2 ()
{
    level3();
}

level2();

output:


L0:
[1]
   [1]
      _[1]=z+1
      _[2]=y3-1
      _[3]=x+yz2-yz-y
   [2]
      _[1]=z+1
      _[2]=y-1
      _[3]=x+yz2-yz-y
[2]
   [1]
      _[1]=z2+1
      _[2]=y3+z+1
      _[3]=x+yz2-yz-y
   [2]
      _[1]=z2+1
      _[2]=y-z+1
      _[3]=x+yz2-yz-y

voice: 3
L:
[1]
   [1]
      _[1]=z+1
      _[2]=y3-1
      _[3]=x+yz2-yz-y
   [2]
      _[1]=z+1
      _[2]=y-1
      _[3]=x+yz2-yz-y
[2]
   [1]
      _[1]=z2+1
      _[2]=y3+z+1
      _[3]=x+yz2-yz-y
   [2]
      _[1]=z2+1
      _[2]=y3+z+1
      _[3]=x+yz2-yz-y

a1[1]=z2+1
a1[2]=x-yz+y
a1[3]=y3+z+1

a2[1]=z2+1
a2[2]=y-z+1
a2[3]=x-z

jakobkroeker jakobkroeker added this to the Singular 4.1.0 milestone on Oct 1 jakobkroeker jakobkroeker added the bug label on Oct 1 jakobkroeker jakobkroeker self-assigned this on Oct 1 jakobkroeker jakobkroeker added the minor label on Oct 1 jakobkroeker jakobkroeker commented on Oct 6

cannot reproduce on my fc14 laptop