Changeset 766ae3 in git
- Timestamp:
- Jan 26, 2015, 6:03:39 PM (9 years ago)
- Branches:
- (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
- Children:
- 87a13aa432705563e5a24a51194c2ce8d07ce580
- Parents:
- f999689a0462639ba7b739684fd5d5bf45a2cab3
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/normal.lib
rf999689 r766ae3 2286 2286 if(size(N) < nvars(R0)) 2287 2287 { 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; 2290 2294 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)) 2292 2297 { 2293 2298 attrib(I,"isSB",1); … … 2297 2302 I=std(I); 2298 2303 } 2304 kill @rl0; 2299 2305 } 2300 2306 else … … 2307 2313 if(dim(I)==2) 2308 2314 { 2309 2315 def newR=basering; 2310 2316 } 2311 2317 else 2312 2318 { 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); 2330 2337 } 2331 2338 2332 2339 if((nvars(basering)<=3)&&(size(I)>1)) 2333 2340 { 2334 2341 ERROR("This is not equidimensional"); 2335 2342 } 2336 2343
Note: See TracChangeset
for help on using the changeset viewer.