Added a new test which can crash

Author: mesonepigreco

tests/test_graphene_symmetries/dyn_mono_10x10x1_full1

@@ -0,0 +1,55 @@
+Dynamical matrix file
+                                                                           
+  1    2  4  4.6714030  0.0000000  7.0125405  0.0000000  0.0000000  0.0000000
+           1  'C   '    10947.083358978085     
+    1    1      0.0000000000      0.0000000000      0.0000000000
+    2    1      0.4999999950     -0.2886751317      0.0000000000
+
+     Dynamical  Matrix in cartesian axes bla
+
+     q = (    0.000000000   0.000000000   0.000000000 ) 
+
+    1    1
+    1.08322136     0.00000000       0.00000000     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000       1.08322136     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000       0.00000000     0.00000000       0.34072341     0.00000000
+    1    2
+   -1.08322136     0.00000000       0.00000000     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000      -1.08322136     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000       0.00000000     0.00000000      -0.34072341     0.00000000
+    2    1
+   -1.08322136     0.00000000       0.00000000     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000      -1.08322136     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000       0.00000000     0.00000000      -0.34072341     0.00000000
+    2    2
+    1.08322136     0.00000000       0.00000000     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000       1.08322136     0.00000000       0.00000000     0.00000000
+    0.00000000     0.00000000       0.00000000     0.00000000       0.34072341     0.00000000
+
+     Diagonalizing the dynamical matrix
+
+     q = (    0.000000000   0.000000000   0.000000000 ) 
+
+ **************************************************************************
+     freq (    1) =       0.000000 [THz] =       0.000000 [cm-1]
+ (  0.000000  0.000000  0.000000  0.000000 -0.707107  0.000000 ) 
+ (  0.000000  0.000000  0.000000  0.000000 -0.707107  0.000000 ) 
+     freq (    2) =       0.000000 [THz] =       0.000015 [cm-1]
+ ( -0.371080  0.000000  0.601913  0.000000  0.000000  0.000000 ) 
+ ( -0.371080  0.000000  0.601913  0.000000  0.000000  0.000000 ) 
+     freq (    3) =       0.000001 [THz] =       0.000030 [cm-1]
+ ( -0.601913  0.000000 -0.371080  0.000000  0.000000  0.000000 ) 
+ ( -0.601913  0.000000 -0.371080  0.000000  0.000000  0.000000 ) 
+     freq (    4) =      25.956242 [THz] =     865.807028 [cm-1]
+ (  0.000000  0.000000  0.000000  0.000000 -0.707107  0.000000 ) 
+ (  0.000000  0.000000  0.000000  0.000000  0.707107  0.000000 ) 
+     freq (    5) =      46.280648 [THz] =    1543.756236 [cm-1]
+ (  0.000000  0.000000  0.707107  0.000000  0.000000  0.000000 ) 
+ (  0.000000  0.000000 -0.707107  0.000000  0.000000  0.000000 ) 
+     freq (    6) =      46.280648 [THz] =    1543.756236 [cm-1]
+ ( -0.707107  0.000000 -0.000000  0.000000  0.000000  0.000000 ) 
+ (  0.707107  0.000000 -0.000000  0.000000  0.000000  0.000000 ) 
+ **************************************************************************
+
+File generated with q2qstar.x from 'dyn_mono_10x10x1_1'
+

tests/test_graphene_symmetries/test_graphene.py

@@ -0,0 +1,22 @@
+import cellconstructor as CC
+import cellconstructor.Phonons
+
+import sscha, sscha.Ensemble, sscha.SchaMinimizer
+
+import numpy as np
+import sys, os
+
+
+# A mock test that should not crash when we try to impose symmetries on a valid dynamical matrix
+def test_root_step_identity():
+    total_path = os.path.dirname(os.path.abspath(__file__))
+    os.chdir(total_path)
+    
+    dyn = CC.Phonons.Phonons("dyn_mono_10x10x1_full")
+    dyn.Symmetrize()
+
+     
+if __name__ == "__main__":
+    test_root_step_identity()
+
+