OpenXM_contrib2/asir2000/lib/primdec - diff

Return to primdec CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / lib

version 1.6, 2003/12/12 03:08:29

version 1.8, 2006/10/03 08:16:55

Line 45

* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,

* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.

* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec,v 1.5 2003/10/20 00:58:47 takayama Exp $

* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec,v 1.7 2006/02/24 01:15:56 noro Exp $

/* Primary decomposition & Radical decomposition program */

/* written by T.Shimoyama, Fujitsu Lab. Date: 1995.10.12 */

Line 457 def extraction(GF,Pr,VL)

V = minalgdep(Pr,VL,PRIMEORD);

U = listminus(VL,V);

V0 = append(V,U);

#if 0

/* This may be a bug. GF is not a GB w.r.t the elimination order */

if ( V0 != VL ) {

ORD = [[TMPORD1,length(V)],[TMPORD2,length(U)]];

GR(G,GF,V0,ORD);

} else

G = GF;

#else

ORD = [[TMPORD1,length(V)],[TMPORD2,length(U)]];

GR(G,GF,V0,ORD);

#endif

dp_ord(TMPORD1);

for (LL = [],HC = 1,I = 0; I < length(G); I++)

LL = append(LL,cdr(fctr(dp_hc(dp_ptod(G[I],V)))));

Line 1388 def zraddec(F,X)

Line 1394 def zraddec(F,X)

/* contraction from V to X */

def radcont(Q,V,X)

{

dp_ord(irem(PRIMEORD,3));

for (R=[],I=length(Q)-1;I>=0;I--) {

dp_ord(irem(PRIMEORD,3));

G=Q[I];

for (E=1,J=0;J<length(G);J++)

E = lcm(E,dp_hc(dp_ptod(G[J],V)));