Further refinement of automorphisms output

Author: Winfried Bruns

NormalizModule.cpp

@@ -172,11 +172,11 @@ static PyObject* CallPythonFuncOnOneArg(PyObject* function, PyObject* single_arg
 #ifndef NMZ_RELEASE
 static_assert(
     false,
-    "Your Normaliz version (unknown) is to old! Update to 3.8.10 or newer.");
+    "Your Normaliz version (unknown) is to old! Update to 3.9.0 or newer.");
 #endif
-#if NMZ_RELEASE < 30810
+#if NMZ_RELEASE < 30900
 static_assert(false,
-              "Your Normaliz version is to old! Update to 3.8.10 or newer.");
+              "Your Normaliz version is to old! Update to 3.9.0 or newer.");
 #endif
 
 /***************************************************************************
@@ -709,14 +709,18 @@ template < typename Integer >
 static PyObject*
 NmzAutomorphismsToPython(const AutomorphismGroup< Integer >& grp)
 {
+    int list_size = 6;
+    if(grp.IsInput() || grp.IsAmbient())
+        list_size =7;    
 
-    PyObject* list = PyList_New(6);
+    PyObject* list = PyList_New(list_size);
 
     PyList_SetItem(list, 0, NmzToPyNumber(grp.getOrder()));
     PyList_SetItem(list, 1, BoolToPyBool(grp.IsIntegralityChecked()));
-    PyList_SetItem(list, 2, BoolToPyBool(grp.IsIntegral())); 
+    PyList_SetItem(list, 2, BoolToPyBool(grp.IsIntegral()));
 
     if(grp.IsInput() || grp.IsAmbient()){
+        PyList_SetItem(list, 6, NmzMatrixToPyList(grp.getGens().get_elements()));
         PyObject* current = PyList_New(2);        
         PyList_SetItem(current, 0, NmzMatrixToPyList(grp.getGensPerms()));
         PyList_SetItem(current, 1, NmzMatrixToPyList(grp.getGensOrbits()));

PyNormaliz.py

@@ -27,12 +27,31 @@ def print_perms_and_orbits(data, name):
 
 def print_automs(Automs):
     print("order ", Automs[0])
-    if len(Automs[1][0]) >0:
-        print_perms_and_orbits(Automs[1],"extreme rays of (recession) cone")
-    if len(Automs[2][0]) >0:
-        print_perms_and_orbits(Automs[2],"vertices of polyhedron")
+    if Automs[1]:
+        if Automs[2]:
+            print("automorphisms are integral")
+        else:
+            print("automorphisms are not integral")
+    else:
+        print("integrality of automorphisms unknown")
+        
+    gen_name ="extreme rays of (recession) cone"
+    if len(Automs) == 7:
+        gen_name = "input vectors";        
+    lf_name = "support hyperplanes"
+    if len(Automs) == 7:
+        lf_name = "coordinates"
+
     if len(Automs[3][0]) >0:
-        print_perms_and_orbits(Automs[3],"support hyperplanes")
+        print_perms_and_orbits(Automs[3],gen_name)
+    if len(Automs[4][0]) >0:
+        print_perms_and_orbits(Automs[4],"vertices of polyhedron")
+    if len(Automs[5][0]) >0:
+        print_perms_and_orbits(Automs[5],lf_name)
+    
+    if len(Automs) == 7:
+        print("input vectors")
+        print_matrix(Automs[6])
     return
 
 def print_Stanley_dec(dec):

README.md

@@ -13,9 +13,9 @@ A full documentation is conatined in [Appendix E](doc/PyNormaliz.pdf) of the Nor
 ## Requirements
 
 * python 2.7 or higher or python 3.4 or higher
-* Normaliz 3.8.10 or higher <https://github.com/Normaliz/Normaliz/releases>
+* Normaliz 3.9.0 or higher <https://github.com/Normaliz/Normaliz/releases>
 
-(The current version of PyNormaliz is under construction. Normaliz 3.8.10 not yet released. Use current master of Normaliz for current master of PyNormaliz.) The source packages of the Normaliz realeases contain PyNormaliz.
+(The current version of PyNormaliz is under construction. Normaliz 3.9.0 not yet released. Use current master of Normaliz for current master of PyNormaliz.) The source packages of the Normaliz realeases contain PyNormaliz.
 
 ## Installation