Context Navigation

Changeset cc93648 in git

Timestamp:

Apr 24, 2014, 7:03:24 PM (10 years ago)

Author:

Jakob Kroeker <kroeker@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

d21b700adb5e8106ab69ce599d924b7818ffb9ad

Parents:

dcd3cea67ee020fad5765aab2ede6d9e49aa2e5c

git-author:

Jakob Kroeker <kroeker@math.uni-hannover.de>2014-04-24 19:03:24+02:00

git-committer:

Jakob Kröker <kroeker@math.uni-hannover.de>2014-04-25 10:51:36+02:00

Message:

fix variable conflict issue in primdecint.lib

Files:

-                      rdcd3ce
+                      rcc93648
 "USAGE:  intersectZ(I,J); I,J ideals
 RETURN:  the intersection of the input ideals
 NOTE:    this is needed because intersect(I,J) does not work, should be replaced
          by intersect later
+NOTE:    this is an alternative to intersect(I,J) over integers,
+         is faster for some examples and should be kept for debug purposes.
 EXAMPLE: example intersectZ; shows an example
+{
    def R = basering;
+   execute("ring S=integer,(t,"+varstr(R)+"),(dp(1),dp(nvars(R)));");
+   ideal I=imap(R,I);
+   ideal J=imap(R,J);
+   ideal K=addIdealZ(t*I,(1-t)*J);
+   K=stdZ(K);
+   execute("ring S=integer,( X(1..nvars(R)) ), ( dp(nvars(R)) );");
+   ideal I = fetch(R,I);
+   ideal J = fetch(R,J);
+   execute("ring St=integer, ( t, X(1..nvars(R))  ), ( dp(1), dp(nvars(R)) );");
+   ideal I = imap(S,I);
+   ideal J = imap(S,J);
+   ideal K = addIdealZ( var(1)*I,(1-var(1))*J);
+   K = stdZ(K);
    int i;
    ideal L;
+   for(i=1;i<=size(K);i++)
+   {
+      if(lead(K[i])/t==0){L[size(L)+1]=K[i];}
+   }
+   for(i=1; i<=size(K); i++)
+   {
+      if ( lead(K[i])/var(1)==0 ) { L[size(L)+1] = K[i]; }
+   }
+   setring S;
+   ideal L = imap(St, L);
    setring R;
    ideal L=imap(S,L);
+   ideal L = fetch(S, L);
    return(L);
+}

Note: See TracChangeset for help on using the changeset viewer.