Changeset 7d24a2 in git for Singular

Timestamp:

Sep 23, 2011, 5:30:15 PM (13 years ago)

Author:

Andreas Steenpaß <steenpas@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

563364d1cd5bd1d1a7baf4123776a2bd37d99f23

Parents:

09723d9937585aae280f53597c2c61b2b3e6304a

Message:

add realclassify.lib

make changes in classify.lib which are needed for realclassify.lib

git-svn-id: file:///usr/local/Singular/svn/trunk@14384 2c84dea3-7e68-4137-9b89-c4e89433aadc

Location:

Singular/LIB

Files:

• classify.lib (modified) (5 diffs)
• realclassify.lib (added)

Singular/LIB/classify.lib

@@ -9 +9 @@
    A library for classifying isolated hypersurface singularities w.r.t. right
    equivalence, based on the determinator of singularities by V.I. Arnold.
+
+SEE ALSO:  realclassify_lib
 
 PROCEDURES:
@@ -25 +27 @@
  debug_log(lev,[]);  print trace and debugging information w.r.t level>@DeBug
  swap(a,b);          swaps the arguments
+ modality(f);        modality of the singularity
+ complexSingType(f); complex type of the singularity as a string
+ prepRealclassify(f);
+                     the modality and the complex type of the singularity at
+                     once
            (parameters in square brackets [] are optional)
 ";
@@ -40 +47 @@
   s = read(l,"VERSION");
   if (s == "" ) {
-    if (printlevel)
+    if (printlevel > 0)
     {
       "(classify.lib): Need to create database...";
@@ -1435 +1442 @@
 static proc printresult (int step, poly f, string typ, list cstn, int m)
 {
+  if(defined(onlyreturninvariants))
+  {
+    if(onlyreturninvariants == 1)
+    {
+      Modality = m;
+      Type = typ;
+    }
+  }
 //---------------------------- initialisation ---------------------------------
   int corank, Mu, K;
@@ -2938 +2953 @@
 
 ///////////////////////////////////////////////////////////////////////////////
+proc prepRealclassify(poly f)
+"
+USAGE:    prepRealclassify(f); f poly
+RETURN:   a list, containing the modality of the singularity and the type of
+          the singularity as a string
+          @* This procedure is needed in realclassify.lib in order to avoid
+          classify() being called more than once.
+EXAMPLE:  example prepRealclassify; shows an example"
+{
+  exportinvariants(f);
+  return(list(Modality, Type));
+}
+example
+{
+  "EXAMPLE"; echo = 2;
+  ring r = 0, (x,y,z), ds;
+  poly f = (x2+3y-2z)^2+xyz-(x-y3+x2z3)^3;
+  prepRealclassify(f);
+}
+
+///////////////////////////////////////////////////////////////////////////////
+proc modality(poly f)
+"
+USAGE:    modality(f); f poly
+RETURN:   the modality of the singularity
+EXAMPLE:  example modality; shows an example"
+{
+  exportinvariants(f);
+  return(Modality);
+}
+example
+{
+  "EXAMPLE"; echo = 2;
+  ring r = 0, (x,y,z), ds;
+  poly f = (x2+3y-2z)^2+xyz-(x-y3+x2z3)^3;
+  modality(f);
+}
+
+///////////////////////////////////////////////////////////////////////////////
+proc complexSingType(poly f)
+"
+USAGE:    complexSingType(f); f poly
+RETURN:   the type of the singularity as a string
+EXAMPLE:  example complexSingType; shows an example"
+{
+  exportinvariants(f);
+  return(Type);
+}
+example
+{
+  "EXAMPLE"; echo = 2;
+  ring r = 0, (x,y,z), ds;
+  poly f = (x2+3y-2z)^2+xyz-(x-y3+x2z3)^3;
+  complexSingType(f);
+}
+
+///////////////////////////////////////////////////////////////////////////////
+/* some of the invariants will be exported during the computation of
+ * classify(f) such that they are accessible at for all procs in the library
+ */
+static proc exportinvariants(poly f)
+{
+  init_debug(-1);
+  if(!defined(onlyreturninvariants))
+  {
+    int onlyreturninvariants;
+    export(onlyreturninvariants);
+  }
+  onlyreturninvariants = 1;
+  if(!defined(Modality))
+  {
+    int Modality;
+    export(Modality);
+  }
+  if(!defined(Type))
+  {
+    string Type;
+    export(Type);
+  }
+  f = classify(f);
+}
+
+///////////////////////////////////////////////////////////////////////////////
 // E n d   O f   F i l e