Changeset 20057b in git
- Timestamp:
- Jun 8, 1999, 5:54:46 PM (24 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- 6047f66a1c9dd5e774d56bb4724616bc964ab76c
- Parents:
- 3375e2707f3e056d897c66030470921db73ca942
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/primdec.lib
r3375e2 r20057b 1 // $Id: primdec.lib,v 1.3 1 1999-03-16 14:18:51 Singular Exp $1 // $Id: primdec.lib,v 1.32 1999-06-08 15:54:46 pfister Exp $ 2 2 //////////////////////////////////////////////////////////////////////////////// 3 3 // primdec.lib // … … 11 11 //////////////////////////////////////////////////////////////////////////////// 12 12 13 version="$Id: primdec.lib,v 1.3 1 1999-03-16 14:18:51 Singular Exp $";13 version="$Id: primdec.lib,v 1.32 1999-06-08 15:54:46 pfister Exp $"; 14 14 info=" 15 15 LIBRARY: primdec.lib: PROCEDURE FOR PRIMARY DECOMPOSITION … … 491 491 "Warning: The characteristic ist too small to use"; 492 492 "the Algorithm of Gianni/Trager/Zacharias."; 493 "The result may be wrong.";493 "The result may be incomplete."; 494 494 "Try primdecSY"; 495 "The result may be also imcomplete."; 495 496 } 496 497 if (reduce(i[m]-t^e,prm,1) !=0) … … 714 715 { 715 716 def @P = basering; 717 int uytrewq; 716 718 int nva = nvars(basering); 717 719 int @k,@s,@n,@k1,zz; … … 731 733 732 734 j=interred(j); 735 733 736 attrib(j,"isSB",1); 734 737 if(vdim(j)==deg(j[1])) … … 950 953 attrib(@qh,"isSB",1); 951 954 ser1=phi1(ser); 955 956 952 957 @lh=zero_decomp (@qh,phi(ser1),@wr); 953 958 // @lh=zero_decomp (@qh,psi(ser),@wr); … … 970 975 else 971 976 { 972 act=factor(@qh[1]); 973 if(2*size(act[1])==size(@lh)) 974 { 975 for(@n=1;@n<=size(act[1]);@n++) 976 { 977 @f=act[1][@n]^act[2][@n]; 978 ser1=psi(@f); 979 lres0[2*@n-1]=interred(primary[2*@k-1]+psi1(ser1)); 980 helpprim=@lh[2*@n]; 981 ser1=psi(helpprim); 982 lres0[2*@n]=psi1(ser1); 983 } 984 } 985 else 986 { 977 //act=factor(@qh[1]); 978 //if(2*size(act[1])==size(@lh)) 979 //{ 980 981 // for(@n=1;@n<=size(act[1]);@n++) 982 // { 983 // @f=act[1][@n]^act[2][@n]; 984 // ser1=psi(@f); 985 // lres0[2*@n-1]=interred(primary[2*@k-1]+psi1(ser1)); 986 // helpprim=@lh[2*@n]; 987 // ser1=psi(helpprim); 988 // lres0[2*@n]=psi1(ser1); 989 // } 990 // } 991 // else 992 // { 987 993 lres1=psi(@lh); 988 994 lres0=psi1(lres1); 989 995 //} 990 996 } 991 997 if(npars(@P)>0)
Note: See TracChangeset
for help on using the changeset viewer.