Changeset e2af3af in git for Singular/LIB/primdec.lib
- Timestamp:
- Jul 1, 2019, 4:07:30 PM (5 years ago)
- Branches:
- (u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')
- Children:
- 2050058328d464a197ff417b127e85fe08b3a3cd
- Parents:
- 77006ca97061dca7e90bdacf1fbaff76644280a24142c906ad243d17ff54542115db04d9e8264de7
- git-author:
- Hans Schoenemann <hannes@mathematik.uni-kl.de>2019-07-01 16:07:30+02:00
- git-committer:
- GitHub <noreply@github.com>2019-07-01 16:07:30+02:00
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/primdec.lib
r77006c re2af3af 2564 2564 2565 2565 indep=maxIndependSet(j); 2566 2567 execute("ring gnir1 = ("+charstr(basering)+"),("+indep[1][1]+"),(" 2568 +indep[1][2]+");"); 2566 ring gnir1 = create_ring(ringlist(basering)[1],"("+indep[1][1]+")","("+indep[1][2]+")","no_minpoly"); 2569 2567 if(homo==1) 2570 2568 { … … 3112 3110 newva[size(newva)]=")"; 3113 3111 newma[size(newma)]=";"; 3114 execute("ring @deirf=("+charstr(gnir)+"),("+newva+",lp;");3112 ring @deirf = create_ring(ringlist(gnir)[1],"("+newva,"lp","no_minpoly"); 3115 3113 execute("map @kappa=gnir,"+newma); 3116 3114 ideal @j= @kappa(@j); … … 3337 3335 @j=fetch(@P,i); 3338 3336 } 3339 execute("ring gnir1 = ("+charstr(basering)+"),("+indep[@m][1]+"),(" 3340 +indep[@m][2]+");"); 3337 ring gnir1 = create_ring(ringlist(basering)[1],"("+indep[@m][1]+")","("+indep[@m][2]+")","no_minpoly"); 3341 3338 execute("map phi=gnir,"+@va+";"); 3342 3339 op=option(get); … … 3646 3643 { 3647 3644 @va=string(maxideal(1)); 3648 execute("ring gnir1 = ("+charstr(basering)+"),("+ 3649 restindep[@m][1]+"),(" +restindep[@m][2]+");"); 3645 ring gnir1 = create_ring(ringlist(basering)[1],"("+restindep[@m][1]+")","("+restindep[@m][2]+")","no_minpoly"); 3650 3646 execute("map phi=gnir,"+@va+";"); 3651 3647 op=option(get); … … 4614 4610 def oldring=basering; 4615 4611 string n=system("neworder",PS); 4616 execute("ring r=("+charstr(oldring)+"),("+n+"),dp;");4612 ring r = create_ring(ringlist(oldring)[1],"("+n+")","dp","no_minpoly"); 4617 4613 ideal PS=imap(oldring,PS); 4618 4614 matrix m=char_series(PS); // We compute an irreducible … … 4629 4625 for(i=1;i<nvars(oldring);i++) { n=n+","+string(var(i)); } 4630 4626 kill r; 4631 execute("ring r=("+charstr(oldring)+"),("+n+"),dp;");4627 ring r = create_ring(ringlist(oldring)[1],"("+n+")","dp","no_minpoly"); 4632 4628 ideal PS=imap(oldring,PS); 4633 4629 matrix m=char_series(PS); … … 5272 5268 5273 5269 U[size(U)]=")"; // we compute the extractor of I (w.r.t. U) 5274 execute("ring RAU=("+charstr(basering)+"),("+A+U+",(dp("+string(a)+"),dp);");5270 ring RAU = create_ring(ringlist(basering)[1],"("+A+U,"(dp("+string(a)+"),dp)","no_minpoly"); 5275 5271 ideal I=imap(R,SI); 5276 5272 //I=std(I,hv); // the standard basis in (R[U])[A] … … 6763 6759 @va = string(maxideal(1)); 6764 6760 6765 execute("ring gnir1 = (" + charstr(basering) + "), (" + indep[@m][1] + "),(" 6766 + indep[@m][2] + ");"); 6761 ring gnir1 = create_ring(ringlist(basering)[1],"(" + indep[@m][1] + ")","(" + indep[@m][2] + ")","no_minpoly"); 6767 6762 execute("map phi = @P," + @va + ";"); 6768 6763 ideal @j = groebner(phi(I)); … … 7325 7320 string(timer-time)); 7326 7321 //---------------- compute characteristic polynomial of matrix -------------- 7327 execute("ring P1=("+charstr(R)+"),T,dp;");7322 ring P1 = create_ring(ringlist(R)[1],"T","dp","no_minpoly"); 7328 7323 matrix n=imap(R,n); 7329 7324 time = timer; … … 7355 7350 //Hence it is better to simultaneously reduce with I. For this we need a new 7356 7351 //ring. 7357 execute("ring P=("+charstr(R)+"),(T,"+varstr(R)+"),(dp(1),dp);");7352 ring P = create_ring(ringlist(R)[1],"(T,"+varstr(R)+")","(dp(1),dp)","no_minpoly"); 7358 7353 list rfac=imap(P1,fac); 7359 7354 intvec ov=option(get);; … … 7652 7647 newva[size(newva)]=")"; 7653 7648 newma[size(newma)]=";"; 7654 execute("ring @deirf=("+charstr(gnir)+"),("+newva+",lp;");7649 ring @deirf = create_ring(ringlist(gnir)[1],"("+newva,"lp","no_minpoly"); 7655 7650 execute("map @kappa=gnir,"+newma); 7656 7651 ideal @j= @kappa(@j); … … 8043 8038 { 8044 8039 @va=string(maxideal(1)); 8045 execute("ring gnir1 = ("+charstr(basering)+"),("+ 8046 restindep[@m][1]+"),(" +restindep[@m][2]+");"); 8040 ring gnir1 = create_ring(ringlist(basering)[1],"(" +restindep[@m][1]+")","(" +restindep[@m][2]+")","no_minpoly"); 8047 8041 execute("map phi=gnir,"+@va+";"); 8048 8042 op=option(get); … … 8387 8381 //} 8388 8382 //======================= 8389 execute("ring gnir1 = ("+charstr(basering)+"),("+indepInfo[1]+"),(" 8390 +indepInfo[2]+");"); 8383 ring gnir1 = create_ring(ringlist(basering)[1],"("+indepInfo[1]+")","(" +indepInfo[2]+")","no_minpoly"); 8391 8384 execute("map phi=gnir,"+@va+";"); 8392 8385 op=option(get);
Note: See TracChangeset
for help on using the changeset viewer.