Changeset 7d24a2 in git
- Timestamp:
- Sep 23, 2011, 5:30:15 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '0604212ebb110535022efecad887940825b97c3f')
- Children:
- 563364d1cd5bd1d1a7baf4123776a2bd37d99f23
- Parents:
- 09723d9937585aae280f53597c2c61b2b3e6304a
- Location:
- Singular/LIB
- Files:
-
- 1 added
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/classify.lib
r09723d r7d24a2 9 9 A library for classifying isolated hypersurface singularities w.r.t. right 10 10 equivalence, based on the determinator of singularities by V.I. Arnold. 11 12 SEE ALSO: realclassify_lib 11 13 12 14 PROCEDURES: … … 25 27 debug_log(lev,[]); print trace and debugging information w.r.t level>@DeBug 26 28 swap(a,b); swaps the arguments 29 modality(f); modality of the singularity 30 complexSingType(f); complex type of the singularity as a string 31 prepRealclassify(f); 32 the modality and the complex type of the singularity at 33 once 27 34 (parameters in square brackets [] are optional) 28 35 "; … … 40 47 s = read(l,"VERSION"); 41 48 if (s == "" ) { 42 if (printlevel )49 if (printlevel > 0) 43 50 { 44 51 "(classify.lib): Need to create database..."; … … 1435 1442 static proc printresult (int step, poly f, string typ, list cstn, int m) 1436 1443 { 1444 if(defined(onlyreturninvariants)) 1445 { 1446 if(onlyreturninvariants == 1) 1447 { 1448 Modality = m; 1449 Type = typ; 1450 } 1451 } 1437 1452 //---------------------------- initialisation --------------------------------- 1438 1453 int corank, Mu, K; … … 2938 2953 2939 2954 /////////////////////////////////////////////////////////////////////////////// 2955 proc prepRealclassify(poly f) 2956 " 2957 USAGE: prepRealclassify(f); f poly 2958 RETURN: a list, containing the modality of the singularity and the type of 2959 the singularity as a string 2960 @* This procedure is needed in realclassify.lib in order to avoid 2961 classify() being called more than once. 2962 EXAMPLE: example prepRealclassify; shows an example" 2963 { 2964 exportinvariants(f); 2965 return(list(Modality, Type)); 2966 } 2967 example 2968 { 2969 "EXAMPLE"; echo = 2; 2970 ring r = 0, (x,y,z), ds; 2971 poly f = (x2+3y-2z)^2+xyz-(x-y3+x2z3)^3; 2972 prepRealclassify(f); 2973 } 2974 2975 /////////////////////////////////////////////////////////////////////////////// 2976 proc modality(poly f) 2977 " 2978 USAGE: modality(f); f poly 2979 RETURN: the modality of the singularity 2980 EXAMPLE: example modality; shows an example" 2981 { 2982 exportinvariants(f); 2983 return(Modality); 2984 } 2985 example 2986 { 2987 "EXAMPLE"; echo = 2; 2988 ring r = 0, (x,y,z), ds; 2989 poly f = (x2+3y-2z)^2+xyz-(x-y3+x2z3)^3; 2990 modality(f); 2991 } 2992 2993 /////////////////////////////////////////////////////////////////////////////// 2994 proc complexSingType(poly f) 2995 " 2996 USAGE: complexSingType(f); f poly 2997 RETURN: the type of the singularity as a string 2998 EXAMPLE: example complexSingType; shows an example" 2999 { 3000 exportinvariants(f); 3001 return(Type); 3002 } 3003 example 3004 { 3005 "EXAMPLE"; echo = 2; 3006 ring r = 0, (x,y,z), ds; 3007 poly f = (x2+3y-2z)^2+xyz-(x-y3+x2z3)^3; 3008 complexSingType(f); 3009 } 3010 3011 /////////////////////////////////////////////////////////////////////////////// 3012 /* some of the invariants will be exported during the computation of 3013 * classify(f) such that they are accessible at for all procs in the library 3014 */ 3015 static proc exportinvariants(poly f) 3016 { 3017 init_debug(-1); 3018 if(!defined(onlyreturninvariants)) 3019 { 3020 int onlyreturninvariants; 3021 export(onlyreturninvariants); 3022 } 3023 onlyreturninvariants = 1; 3024 if(!defined(Modality)) 3025 { 3026 int Modality; 3027 export(Modality); 3028 } 3029 if(!defined(Type)) 3030 { 3031 string Type; 3032 export(Type); 3033 } 3034 f = classify(f); 3035 } 3036 3037 /////////////////////////////////////////////////////////////////////////////// 2940 3038 // E n d O f F i l e 2941 3039
Note: See TracChangeset
for help on using the changeset viewer.