Changeset 71a8ba in git


Ignore:
Timestamp:
Apr 11, 2011, 9:58:27 AM (12 years ago)
Author:
Stefan Steidel <steidel@…>
Branches:
(u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')
Children:
0f93f5f4a65af4c121349d23f99a0869a0a9be2d
Parents:
2272157187c82cb5bd15a8c9b7faff7e3b064132
Message:
intersectZ --> intersect

git-svn-id: file:///usr/local/Singular/svn/trunk@14132 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed

  • Singular/LIB/primdecint.lib

    r2272157 r71a8ba  
    2222 heightZ(I);        compute the height of I
    2323 equidimZ(I);       compute the equidimensional part of I
    24  intersectZ(I,J)    compute the intersection of I and J
    2524";
    2625
    … …  
    233232               if(size(reduce(TES,Q))>0)
    234233               {
    235                   TQ=intersectZ(TQ,Q);
    236                   //TQ=intersect(TQ,Q);
     234                  TQ=intersect(TQ,Q);
    237235                  P[size(P)+1]=list(Q,K);
    238236               }
    … …  
    253251               if(size(reduce(TES,Q))>0)
    254252               {
    255                   //TQ=intersectZ(TQ,Q);
     253                  //TQ=intersect(TQ,Q);
    256254                  P[size(P)+1]=list(Q,K);
    257255               }
    … …  
    262260               JJ=JJ,p;
    263261               JJ=stdZ(JJ);
    264                TQ=intersectZ(TQ,JJ);
    265                //TQ=intersect(TQ,JJ);
     262               TQ=intersect(TQ,JJ);
    266263            }
    267264         }
    … …  
    276273      if(!ex){return(P);}
    277274      J=stdZ(J);
    278       TQ=intersectZ(TQ,TES);
    279       //TQ=intersect(TQ,TES);
     275      TQ=intersect(TQ,TES);
    280276      if(size(reduce(TQ,J))!=0)
    281277      {
    … …  
    284280         ideal W=K;
    285281         m++;
    286          while(size(reduce(intersectZ(W,TQ),J))!=0)
    287          //while(size(reduce(intersect(W,TQ),J))!=0)
     282         while(size(reduce(intersect(W,TQ),J))!=0)
    288283         {
    289284            //W=stdZ(addIdealZ(I,K^m));
    … …  
    363358   ideal I5=9c5,6d5;
    364359   ideal I6=17,a15,b15,c15,d15;
    365    ideal I=intersectZ(I1,I2);
    366    I=intersectZ(I,I3);
    367    I=intersectZ(I,I4);
    368    I=intersectZ(I,I5);
    369    I=intersectZ(I,I6);
     360   ideal I=intersect(I1,I2);
     361   I=intersect(I,I3);
     362   I=intersect(I,I4);
     363   I=intersect(I,I5);
     364   I=intersect(I,I6);
    370365   primdecZ(I);
    371    ideal J=intersectZ(ideal(17,a),ideal(17,a2,b));
     366   ideal J=intersect(ideal(17,a),ideal(17,a2,b));
    372367   primdecZ(J);
    373    ideal K=intersectZ(ideal(9,a+3),ideal(9,b+3));
     368   ideal K=intersect(ideal(9,a+3),ideal(9,b+3));
    374369   primdecZ(K);
    375370}
    … …  
    499494   ideal I5=9c5,6d5;
    500495   ideal I6=17,a15,b15,c15,d15;
    501    ideal I=intersectZ(I1,I2);
    502    I=intersectZ(I,I3);
    503    I=intersectZ(I,I4);
    504    I=intersectZ(I,I5);
    505    I=intersectZ(I,I6);
     496   ideal I=intersect(I1,I2);
     497   I=intersect(I,I3);
     498   I=intersect(I,I4);
     499   I=intersect(I,I5);
     500   I=intersect(I,I6);
    506501   minAssZ(I);
    507    ideal J=intersectZ(ideal(17,a),ideal(17,a2,b));
     502   ideal J=intersect(ideal(17,a),ideal(17,a2,b));
    508503   minAssZ(J);
    509    ideal K=intersectZ(ideal(9,a+3),ideal(9,b+3));
     504   ideal K=intersect(ideal(9,a+3),ideal(9,b+3));
    510505   minAssZ(K);
    511506}
    … …  
    595590   ideal I5=9c5,6d5;
    596591   ideal I6=17,a15,b15,c15,d15;
    597    ideal I=intersectZ(I1,I2);
    598    I=intersectZ(I,I3);
    599    I=intersectZ(I,I4);
    600    I=intersectZ(I,I5);
    601    I=intersectZ(I,I6);
     592   ideal I=intersect(I1,I2);
     593   I=intersect(I,I3);
     594   I=intersect(I,I4);
     595   I=intersect(I,I5);
     596   I=intersect(I,I6);
    602597   heightZ(I);
    603598}
    … …  
    652647         B=B,p;
    653648         B=stdZ(B);
    654          K=stdZ(intersectZ(K,B));
    655          //K=stdZ(intersect(K,B));
     649         K=stdZ(intersect(K,B));
    656650         setring Rhelp;
    657651      }
    … …  
    681675   {
    682676      ideal M=radicalZ(J);
    683       K=intersectZ(K,M);
    684       //K=intersect(K,M);
     677      K=intersect(K,M);
    685678   }
    686679   return(K);
    … …  
    695688   ideal I5=9c5,6d5;
    696689   ideal I6=17,a15,b15,c15,d15;
    697    ideal I=intersectZ(I1,I2);
    698    I=intersectZ(I,I3);
    699    I=intersectZ(I,I4);
    700    I=intersectZ(I,I5);
    701    I=intersectZ(I,I6);
     690   ideal I=intersect(I1,I2);
     691   I=intersect(I,I3);
     692   I=intersect(I,I4);
     693   I=intersect(I,I5);
     694   I=intersect(I,I6);
    702695   radicalZ(I);
    703    ideal J=intersectZ(ideal(17,a),ideal(17,a2,b));
     696   ideal J=intersect(ideal(17,a),ideal(17,a2,b));
    704697   radicalZ(J);
    705698}
    … …  
    811804                  N=stdZ(N);
    812805                  Q=extractZ(N,j,IS,B);
    813                   TQ=intersectZ(TQ,Q);
    814                   //TQ=intersect(TQ,Q);
     806                  TQ=intersect(TQ,Q);
    815807               }
    816808               setring Rhelp;
    … …  
    825817               E=E,p;
    826818               E=stdZ(E);
    827                TQ=intersectZ(TQ,E);
    828                //TQ=intersect(TQ,E);
     819               TQ=intersect(TQ,E);
    829820            }
    830821         }
    … …  
    863854       if(he==heightZ(M))
    864855       {
    865           E=intersectZ(M,E);
    866           //E=intersect(M,E);
     856          E=intersect(M,E);
    867857       }
    868858   }
    … …  
    878868   ideal I5=9c5,6d5;
    879869   ideal I6=17,a15,b15,c15,d15;
    880    ideal I=intersectZ(I1,I2);
    881    I=intersectZ(I,I3);
    882    I=intersectZ(I,I4);
    883    I=intersectZ(I,I5);
    884    I=intersectZ(I,I6);
     870   ideal I=intersect(I1,I2);
     871   I=intersect(I,I3);
     872   I=intersect(I,I4);
     873   I=intersect(I,I5);
     874   I=intersect(I,I6);
    885875   equidimZ(I);
    886 }
    887 
    888 ////////////////////////////////////////////////////////////////////////////////
    889 
    890 proc intersectZ(ideal I, ideal J)
    891 "USAGE:  intersectZ(I,J); I,J ideals
    892 RETURN:  the intersection of the input ideals
    893 NOTE:    this is needed because intersect(I,J) does not work, should be replaced
    894          by intersect later
    895 EXAMPLE: example intersectZ; shows an example
    896 {
    897    def R = basering;
    898    execute("ring S=integer,(t,"+varstr(R)+"),(dp(1),dp(nvars(R)));");
    899    ideal I=imap(R,I);
    900    ideal J=imap(R,J);
    901    ideal K=addIdealZ(t*I,(1-t)*J);
    902    K=stdZ(K);
    903    int i;
    904    ideal L;
    905    for(i=1;i<=size(K);i++)
    906    {
    907       if(lead(K[i])/t==0){L[size(L)+1]=K[i];}
    908    }
    909    setring R;
    910    ideal L=imap(S,L);
    911    return(L);
    912 }
    913 example
    914 { "EXAMPLE:";  echo = 2;
    915    ring R=integer,(a,b,c,d),dp;
    916    ideal I1=9,a,b;
    917    ideal I2=3,c;
    918    ideal I3=11,2a,7b;
    919    ideal I4=13a2,17b4;
    920    ideal I5=9c5,6d5;
    921    ideal I6=17,a15,b15,c15,d15;
    922    ideal I=intersectZ(I1,I2); I;
    923    I=intersectZ(I,I3); I;
    924    I=intersectZ(I,I4); I;
    925    I=intersectZ(I,I5); I;
    926    I=intersectZ(I,I6); I;
    927876}
    928877
    … …  
    12701219   def R=basering;
    12711220   int i;
    1272    ideal K=intersectZ(I,ideal(f));
    1273    //ideal K=intersect(I,ideal(f));
     1221   ideal K=intersect(I,ideal(f));
    12741222   //=== K[i]/f; does not work in rings with integer! This should be replaced
    12751223   //=== later
    … …  
    12961244   for(i=2;i<=size(J);i++)
    12971245   {
    1298       K=intersectZ(K,quotientOneZ(I,J[i]));
    1299       //K=intersect(K,quotientOneZ(I,J[i]));
     1246      K=intersect(K,quotientOneZ(I,J[i]));
    13001247   }
    13011248   return(K);
    … …  
    13891336   for(i=2;i<=size(L);i++)
    13901337   {
    1391       K=intersectZ(K,L[i][1]);
    1392       //K=intersect(K,L[i][1]);
     1338      K=intersect(K,L[i][1]);
    13931339   }
    13941340   i=size(reduce(K,stdZ(I)))+size(reduce(I,stdZ(K)));
    … …  
    15411487        x(1)*x(4)^3*x(5),
    15421488        x(4)^4*x(5);
    1543       I=intersectZ(I,ideal(64*181,x(1)^2));
     1489      I=intersect(I,ideal(64*181,x(1)^2));
    15441490
    15451491ring R17=integer,(x,y,z),dp;
Note: See TracChangeset for help on using the changeset viewer.