Changeset 5895e8 in git
- Timestamp:
- Sep 23, 2010, 3:51:15 PM (13 years ago)
- Branches:
- (u'spielwiese', '828514cf6e480e4bafc26df99217bf2a1ed1ef45')
- Children:
- 807759680bf2373d9851bbe0a7f175e190f4c18a
- Parents:
- fc2d16d75fa03d73ea81fee4ff230ad4158b821e
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/resjung.lib
rfc2d16 r5895e8 1 // 1 2 version="$Id: jung.lib,v 1.0 2008/13/06 10:32:00 Singular Exp $"; 2 3 category="Commutative Algebra"; … … 7 8 Anne Fruehbis-Krueger, anne@math.uni-hannover.de 8 9 9 MAINPROCEDURES:10 PROCEDURES: 10 11 jungresolve(J,i) computes a resolution of the surface given by the 11 12 ideal J using Jungs Method 12 13 AUXILLARY PROCEDURES:14 13 clocus(J) computes the critical locus of the projection of V(J) 15 14 onto the coordinate plane of the last two coordinates … … 90 89 if(attrib(J,"isPrincipal")==0){ 91 90 setring R; 92 for(j = 1;j<=size(J);j++){ //determines the monic polynomial in var(i) with coefficents in C2 91 for(j = 1;j<=size(J);j++) 92 { //determines the monic polynomial in var(i) with coefficents in C2 93 93 if(defined(w)) {kill w;} 94 94 intvec w = leadexp(J[j]); … … 108 108 setring A; 109 109 } 110 product = product*resultant(J[index],diff(J[index],var(i)),var(i)); //Product of the discriminants, which lies in C2 110 //Product of the discriminants, which lies in C2: 111 product = product*resultant(J[index],diff(J[index],var(i)),var(i)); 111 112 kill index,J,v; 112 113 } … … 115 116 ideal Clocus = imap(A,product); //the critical locus is contained in this 116 117 ideal I = preimage(A,LastTwo,lowdim); 117 Clocus = radical(intersect(Clocus,I)); //radical is necessary since the resultant is not reduced in general 118 //radical is necessary since the resultant is not reduced in general: 119 Clocus = radical(intersect(Clocus,I)); 118 120 export(Clocus); 119 121 return(C2);
Note: See TracChangeset
for help on using the changeset viewer.