Added visualisation of spectrum in examples

diff --git a/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DEF/SpectrumLaplace2DEF_SQUARE.py b/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DEF/SpectrumLaplace2DEF_SQUARE.py
index 3301a1da1a8331f2f86db4e6ef7989fbf7a9d3a2..984733a7ca60129d11b68391e310e7321501170f 100755 (executable)
--- a/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DEF/SpectrumLaplace2DEF_SQUARE.py
+++ b/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DEF/SpectrumLaplace2DEF_SQUARE.py
@@ -117,7 +117,9 @@ for i in range(nbCells):
                                #else: si condition limite non nulle au bord, ajouter la contribution du bord au second membre de la cellule j
 
 print("Stiffness matrix construction done")
-quit()
+Rigidite.viewMatrix(True, 0, "RigidityMatrix_FiniteElementsOn"+mesh_name+"Laplace")
+Rigidite.plotEigenvalues("FiniteElementsOn"+mesh_name+"Laplace")
+
 # Conditionnement de la matrice de rigidité
 #=================================
 cond = Rigidite.getConditionNumber()

diff --git a/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DVF/SpectrumLaplace2DVF_SQUARE.py b/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DVF/SpectrumLaplace2DVF_SQUARE.py
index a60054c80023255bbadcf1d95c09e7a4e0c1c02f..40eb0ba52a1c519e5372be8bcd896b92c86ea35c 100755 (executable)
--- a/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DVF/SpectrumLaplace2DVF_SQUARE.py
+++ b/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplace2DVF/SpectrumLaplace2DVF_SQUARE.py
@@ -74,6 +74,8 @@ for i in range(nbCells):
         Rigidite.addValue(i,i,coeff) # terme diagonal
 
 print("Stiffness matrix construction done")
+Rigidite.viewMatrix(True, 0, "RigidityMatrix_FiniteVolumesOn"+mesh_name+"Laplace")
+Rigidite.plotEigenvalues("FiniteVolumesOn"+mesh_name+"Laplace")
 
 # Conditionnement de la matrice de rigidité
 #=================================

diff --git a/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplaceBeltrami3DEF/SpectrumFiniteElements3DLaplace-Beltrami.py b/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplaceBeltrami3DEF/SpectrumFiniteElements3DLaplace-Beltrami.py
index a6b3f3efcc6e003210fd3b445917b278ee968970..0b99481a11b5d2b614125ee0c5ce459107c2b082 100755 (executable)
--- a/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplaceBeltrami3DEF/SpectrumFiniteElements3DLaplace-Beltrami.py
+++ b/CDMATH/tests/examples/SpectrumLaplace/SpectrumLaplaceBeltrami3DEF/SpectrumFiniteElements3DLaplace-Beltrami.py
@@ -110,6 +110,8 @@ for i in range(nbCells):
                        Rigidite.addValue(j,k,GradShapeFuncs[j]*GradShapeFuncs[k]/Ci.getMeasure())
 
 print("Linear system matrix building done")
+Rigidite.viewMatrix(True, 0, "RigidityMatrix_FiniteElementsOn"+mesh_name+"LaplaceBeltrami")
+Rigidite.plotEigenvalues("FiniteElementsOn"+mesh_name+"LaplaceBeltrami")
 
 # Conditionnement de la matrice de rigidité
 #==========================================