import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
def test_validation2DWaveSystemCentered_triangles(scaling):
start = time.time()
#### 2D Delaunay triangles meshes
- meshList=['squareWithTriangles_1','squareWithTriangles_2','squareWithTriangles_3','squareWithTriangles_4']
+ meshList=['squareWithTriangles_1','squareWithTriangles_2','squareWithTriangles_3']#,'squareWithTriangles_4'
meshType="Unstructured_Delaunay_triangles"
testColor="Green"
nbMeshes=len(meshList)
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
def test_validation2DWaveSystemPStagCrossTriangles(scaling):
start = time.time()
#### 2D cross triangle meshes
- meshList=['squareWithCrossTriangles_0','squareWithCrossTriangles_1','squareWithCrossTriangles_2','squareWithCrossTriangles_3']#,'squareWithCrossTriangles_4'
+ meshList=['squareWithCrossTriangles_0','squareWithCrossTriangles_1','squareWithCrossTriangles_2']#,'squareWithCrossTriangles_3','squareWithCrossTriangles_4'
meshType="Unstructured_cross_triangles"
testColor="Green"
nbMeshes=len(meshList)
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
def test_validation2DWaveSystemPStagTriangles(scaling):
start = time.time()
#### 2D Delaunay triangles meshes
- meshList=['squareWithTriangles_1','squareWithTriangles_2','squareWithTriangles_3','squareWithTriangles_4']
+ meshList=['squareWithTriangles_1','squareWithTriangles_2','squareWithTriangles_3']
meshType="Unstructured_Delaunay_triangles"
testColor="Green"
nbMeshes=len(meshList)
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
-import cdmath
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
def test_validation2DWaveSystemPStag_squares(scaling):
start = time.time()
#### 2D square mesh
- #meshList=[7,15,31,51]#,81,101
meshList=['squareWithSquares_1','squareWithSquares_2','squareWithSquares_3','squareWithSquares_4']#,'squareWithSquares_5'
mesh_path='../../../ressources/2DCartesien/'
meshType="Regular squares"
cfl=0.5
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
- #for nx in meshList:
- #my_mesh=cdmath.Mesh(0,1,nx,0,1,nx)
- #error_p_tab[i], error_u_tab[i], mesh_size_tab[i], t_final[i], ndt_final[i], max_vel[i], diag_data_press[i], diag_data_vel[i], time_tab[i], cond_number[i] =WaveSystemPStag.solve(my_mesh,"square"+str(nx)+'x'+str(nx),resolution,scaling,meshType,testColor,cfl,"Periodic",)
error_p_tab[i], error_u_tab[i], mesh_size_tab[i], t_final[i], ndt_final[i], max_vel[i], diag_data_press[i], diag_data_vel[i], time_tab[i], cond_number[i] =WaveSystemPStag.solve_file(mesh_path+filename, mesh_name, resolution,scaling,meshType,testColor,cfl,"Periodic")
assert max_vel[i]>0.5 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
import cdmath
import WaveSystemStaggered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
def test_validation2DWaveSystemUpwindDelaunayTriangles(bctype,scaling):
start = time.time()
#### 2D Delaunay triangles mesh
- meshList=['squareWithTriangles_1','squareWithTriangles_2','squareWithTriangles_3','squareWithTriangles_4','squareWithTriangles_5']
+ meshList=['squareWithTriangles_1','squareWithTriangles_2','squareWithTriangles_3']#,'squareWithTriangles_4','squareWithTriangles_5'
meshType="Unstructured_Delaunay_triangles"
testColor="Green"
nbMeshes=len(meshList)
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemCentered
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemPStag
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt
import cdmath
import WaveSystemUpwind
+import matplotlib
+matplotlib.use("Agg")
import matplotlib.pyplot as plt
import numpy as np
from math import log10, sqrt