Opened 10 years ago

Closed 10 years ago

Last modified 10 years ago

#493 closed bug (invalid)

problems as before

Reported by: gorzel Owned by: somebody
Priority: major Milestone: 3-2-0 and higher
Component: dontKnow Version: 3-1-6
Keywords: factories, absFactorize, readline Cc:

Description

First, the nightly built Singular-3-1-6-x86_64-Linux-static.tar from 15.5.13 gives at start up:

Singular : signal 4 (v: 3160):
current line:>>;return();<<
Segment fault/Bus error occurred at 13eace0 because of 10287 (r:1368634491)
please inform the authors
trying to restart...
fatal flex scanner internal error--end of buffer missed

So I had to start it with option -t

gorzelc@SCHAF17:/tmp/SING316/Singular/3-1-6/x86_64-Linux$ ./Singular -tv
Singular for x86_64-Linux version 3-1-6 (3160)  May 15 2013 07:09:27
with
        factory(@(#) factoryVersion = 3.1.6),libfac(3.1.6,December 2012),
        GMP(5.0),NTL(5.5.2),64bit,static readline,Plural,fan/cone,DBM,

My interest was to check the improved absFactorize with the following example, which seems to give wrong results in older versions, but now it comes to a crash:

>  LIB "absfact.lib";
// ** loaded /data/tmp/SING316/Singular/3-1-6/LIB/absfact.lib (15793,2013-04-30)
>  ring r=0,x,dp;
> poly f = 9*x^6-27*x^5-45*x^4+195*x^3-20*x^2-372*x+276;     
> def S = absFactorize(f);
Singular : signal 4 (v: 3160):
current line:>>  list p_facts = factorize(p);<<
Segment fault/Bus error occurred at 83 because of 10216 (r:1368634610)
please inform the authors
trying to restart...
   ? size(`nothing`) failed
   ? expected size(`bigint`)
   ? expected size(`number`)
   ? expected size(`resolution`)
   ? expected size(`string`)
   ? expected size(`poly`)
   ? expected size(`vector`)
   ? expected size(`ideal`)
   ? expected size(`module`)
   ? expected size(`matrix`)
   ? expected size(`intvec`)
   ? expected size(`intmat`)
   ? expected size(`list`)
   ? expected size(`ring`)
   ? error occurred in or before absfact.lib::smallestProperSimpleFactor line 86: `  int s = size(p_facts[1]);`
   ? expected int-expression. type 'help int;'
   ? leaving absfact.lib::smallestProperSimpleFactor
Auf Wiedersehen.

Change History (2)

comment:1 Changed 10 years ago by hannes

Resolution: invalid
Status: newclosed

not reproducable, hardware not compatible? (signal 4 is SIGILL: Illegal Instruction)

                     SINGULAR                                 /  Development
 A Computer Algebra System for Polynomial Computations       /   version 3-1-6
                                                           0<
 by: W. Decker, G.-M. Greuel, G. Pfister, H. Schoenemann     \   Dec 2012
FB Mathematik der Universitaet, D-67653 Kaiserslautern        \
// ** executing /tmp/hannes/Singular/LIB/.singularrc
>  LIB"absfact.lib";
// ** loaded /tmp/hannes/Singular/LIB/absfact.lib (14191,2011-05-04)
>  ring r=0,x,dp;
>  poly f = 9*x^6-27*x^5-45*x^4+195*x^3-20*x^2-372*x+276;
>  def S = absFactorize(f);

// 'absFactorize' created a ring, in which a list absolute_factors (the
// absolute factors) is stored.
// To access the list of absolute factors, type (if the name S was assigned
// to the return value):
        setring(S); absolute_factors; 
> setring(S); absolute_factors; 
[1]:
   _[1]=9
   _[2]=x+(-a)
[2]:
   1,1
[3]:
   _[1]=(a)
   _[2]=(9a6-27a5-45a4+195a3-20a2-372a+276)
[4]:
   6

comment:2 Changed 10 years ago by anonymous

Hm, I will try to find out where this

(signal 4 is SIGILL: Illegal Instruction)

comes from But I used the same machines with Ubuntu as before.

But note also, whether this realut of absFactorize could not be improved:

> ring r =0,x,dp;
> poly g = x^6 - x^5 + 5*x^4 - 5*x^3 + 10*x^2 - 6*x + 6;
> def Sg = absFactorize(g);

>   setring(Sg); absolute_factors; 
[1]:
   _[1]=1
   _[2]=-x+(a4+3a2-a+3)
[2]:
   1,1
[3]:
   _[1]=(a)
   _[2]=(a6-a5+5a4-5a3+10a2-6a+6)
[4]:
   6

> ring ra4 = (0,a),x,dp; minpoly = a6-a5+5a4-5a3+10a2-6a+6;
> factorize(x6-x5+5x4-5x3+10x2-6x+6);
[1]:
   _[1]=1
   _[2]=x+(-a4-3a2+a-3)
   _[3]=x+(-a)
   _[4]=x2+(-a2+a-2)*x+(a4-a3+4a2-3a+5)
   _[5]=x2+(a4+4a2-a+4)*x+(a5-a4+4a3-5a2+6a-4)
[2]:
   1,1,1,1,1

I would like to see x-a as return value from absfact.

BTW: your lib has not the recent one

Note: See TracTickets for help on using tickets.