Changeset 766ae3 in git


Ignore:
Timestamp:
Jan 26, 2015, 6:03:39 PM (9 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
Children:
87a13aa432705563e5a24a51194c2ce8d07ce580
Parents:
f999689a0462639ba7b739684fd5d5bf45a2cab3
Message:
ring definitions for Normal::genus
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/normal.lib

    rf999689 r766ae3  
    22862286   if(size(N) < nvars(R0))
    22872287   {
    2288      string newvar=string(N);
    2289      execute("ring R=("+charstr(R0)+"),("+newvar+"),dp;");
     2288     list newvar=N[1..ncols(N)];
     2289     list @rl0=ringlist(R0);
     2290     @rl0[2]=list(N[1..ncols(N)]);
     2291     @rl0[3]=list(list("dp",1:ncols(N)),list("C",0));
     2292     def R=ring(@rl0); kill @rl0;
     2293     setring R;
    22902294     ideal I =imap(R0,I);
    2291      if((ringlist(R0)[3][1][1]=="dp")&&(size(ringlist(R0)[3])==2))
     2295     list @rl0=ringlist(R0);
     2296     if((@rl0[3][1][1]=="dp")&&(size(@rl0[3])==2))
    22922297     {
    22932298       attrib(I,"isSB",1);
     
    22972302       I=std(I);
    22982303     }
     2304     kill @rl0;
    22992305   }
    23002306   else
     
    23072313   if(dim(I)==2)
    23082314   {
    2309       def newR=basering;
     2315     def newR=basering;
    23102316   }
    23112317   else
    23122318   {
    2313       short=0;
    2314       string smp = string(minpoly);
    2315       if(dim(I)==0)
    2316       {
    2317          execute("ring Rhelp=("+charstr(R0)+"),(@s,@t),dp;");
    2318       }
    2319       else
    2320       {
    2321          execute("ring Rhelp=("+charstr(R0)+"),(@s),dp;");
    2322       }
    2323       if (smp!="0")
    2324       { execute("minpoly = "+smp+";");}
    2325       def newR=R+Rhelp;
    2326       setring newR;
    2327       ideal I=imap(R,I);
    2328       I=homog(I,@s);
    2329       attrib(I,"isSB",1);
     2319     list @rl0=ringlist(R0);
     2320     if(dim(I)==0)
     2321     {
     2322       @rl0[2]=list("@s","@t");
     2323       @rl0[3]=list(list("dp",1:2),list("C",0));
     2324     }
     2325     else
     2326     {
     2327       @rl0[2]=list("@s");
     2328       @rl0[3]=list(list("dp",1),list("C",0));
     2329     }
     2330     def Rhelp=ring(@rl0);
     2331     kill @rl0;
     2332     def newR=R+Rhelp;
     2333     setring newR;
     2334     ideal I=imap(R,I);
     2335     I=homog(I,@s);
     2336     attrib(I,"isSB",1);
    23302337   }
    23312338
    23322339   if((nvars(basering)<=3)&&(size(I)>1))
    23332340   {
    2334        ERROR("This is not equidimensional");
     2341     ERROR("This is not equidimensional");
    23352342   }
    23362343
Note: See TracChangeset for help on using the changeset viewer.