# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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] = WaveSystemCentered.solve_file(mesh_path+filename, mesh_name, resolution,scaling,meshType,testColor,cfl,"Periodic")
- assert max_vel[i]>0.94 and max_vel[i]<1.1
+ assert max_vel[i]>0.2 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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] =WaveSystemCentered.solve_file(mesh_path+filename, mesh_name, resolution,scaling,meshType,testColor,cfl,"Periodic")
- assert max_vel[i]>0.87 and max_vel[i]<1
+ assert max_vel[i]>0.07 and max_vel[i]<1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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] = WaveSystemCentered.solve_file(mesh_path+filename, mesh_name, resolution,scaling,meshType,testColor,cfl,"Periodic")
- assert max_vel[i]>0.94 and max_vel[i]<1.1
+ assert max_vel[i]>0.09 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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] = WaveSystemCentered.solve_file(mesh_path+filename, mesh_name, resolution,scaling,meshType,testColor,cfl,"Periodic")
- assert max_vel[i]>0.94 and max_vel[i]<1.1
+ assert max_vel[i]>0.3 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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] = WaveSystemCentered.solve_file(mesh_path+filename, mesh_name, resolution,scaling,meshType,testColor,cfl,"Periodic")
- assert max_vel[i]>0.94 and max_vel[i]<1.1
+ assert max_vel[i]>0.5 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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.94 and max_vel[i]<1.2
+ assert max_vel[i]>0.3 and max_vel[i]<1.2
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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.98 and max_vel[i]<1.3
+ assert max_vel[i]>0.2 and max_vel[i]<1.3
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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.94 and max_vel[i]<1.1
+ assert max_vel[i]>0.18 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
# Storing of numerical errors, mesh sizes and diagonal values
for filename in meshList:
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.94 and max_vel[i]<1.1
+ assert max_vel[i]>0.8 4 and max_vel[i]<1.1
error_p_tab[i]=log10(error_p_tab[i])
error_u_tab[i]=log10(error_u_tab[i])
time_tab[i]=log10(time_tab[i])
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV Centered on 2D DISK with square meshes : scheme order for pressure is ", -ap
+ print( "FV Centered on 2D DISK with square meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV Centered on 2D DISK with square meshes : scheme order for velocity is ", -au
+ print( "FV Centered on 2D DISK with square meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV pseudo staggered on 2D DISK with square meshes : scheme order for pressure is ", -ap
+ print( "FV pseudo staggered on 2D DISK with square meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV pseudo staggered on 2D DISK with square meshes : scheme order for velocity is ", -au
+ print( "FV pseudo staggered on 2D DISK with square meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D DISK with square meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D DISK with square meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D DISK with square meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D DISK with square meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVCentered on 2D triangular meshes : scheme order for velocity without scaling is ", -au
+ print( "FVCentered on 2D triangular meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVCentered on 2D triangular meshes : scheme order for velocity with scaling is ", -au
+ print( "FVCentered on 2D triangular meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
bp=(-a2*b1p+a1*b2p)/det
if(scaling==0):
- print "FV Centered on 2D checkerboard meshes : scheme order for pressure without scaling is ", -ap
+ print( "FV Centered on 2D checkerboard meshes : scheme order for pressure without scaling is ", -ap)
else:
- print "FV Centered on 2D checkerboard meshes : scheme order for pressure with scaling is ", -ap
+ print( "FV Centered on 2D checkerboard meshes : scheme order for pressure with scaling is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FV Centered on 2D checkerboard meshes : scheme order for velocity without scaling is ", -au
+ print( "FV Centered on 2D checkerboard meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FV Centered on 2D checkerboard meshes : scheme order for velocity with scaling is ", -au
+ print( "FV Centered on 2D checkerboard meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVCentered on 2D cross triangle meshes : scheme order for velocity without scaling is ", -au
+ print( "FVCentered on 2D cross triangle meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVCentered on 2D cross triangle meshes : scheme order for velocity with scaling is ", -au
+ print( "FVCentered on 2D cross triangle meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV Centered on 2D deformed quadrangles meshes : scheme order for pressure is ", -ap
+ print( "FV Centered on 2D deformed quadrangles meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV Centered on 2D deformed quadrangles meshes : scheme order for velocity is ", -au
+ print( "FV Centered on 2D deformed quadrangles meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVCentered on 2D Delaunay triangles meshes : scheme order for velocity without scaling is ", -au
+ print( "FVCentered on 2D Delaunay triangles meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVCentered on 2D Delaunay triangles meshes : scheme order for velocity with scaling is ", -au
+ print( "FVCentered on 2D Delaunay triangles meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVCentered on 2D flat cross triangle meshes : scheme order for velocity without scaling is ", -au
+ print( "FVCentered on 2D flat cross triangle meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVCentered on 2D flat cross triangle meshes : scheme order for velocity with scaling is ", -au
+ print( "FVCentered on 2D flat cross triangle meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV Centered on 2D hexagonal meshes : scheme order for pressure is ", -ap
+ print( "FV Centered on 2D hexagonal meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV Centered on 2D hexagonal meshes : scheme order for velocity is ", -au
+ print( "FV Centered on 2D hexagonal meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
bp=(-a2*b1p+a1*b2p)/det
if(scaling==0):
- print "FV Centered on 2D square meshes : scheme order for pressure without scaling is ", -ap
+ print( "FV Centered on 2D square meshes : scheme order for pressure without scaling is ", -ap)
else:
- print "FV Centered on 2D square meshes : scheme order for pressure with scaling is ", -ap
+ print( "FV Centered on 2D square meshes : scheme order for pressure with scaling is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVCentered on 2D square meshes : scheme order for velocity without scaling is ", -au
+ print( "FVCentered on 2D square meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVCentered on 2D square meshes : scheme order for velocity with scaling is ", -au
+ print( "FVCentered on 2D square meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVPStag on 2D triangular meshes : scheme order for velocity without scaling is ", -au
+ print( "FVPStag on 2D triangular meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVPStag on 2D triangular meshes : scheme order for velocity with scaling is ", -au
+ print( "FVPStag on 2D triangular meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
if(scaling==0):
print( "FV pseudo staggered on 2D checkerboard meshes : scheme order for pressure without scaling is ", -ap)
else:
- print( "FV pseudo staggered on 2D checkerboard meshes : scheme order for pressure with scaling is ", -ap)
+ print( "FV pseudo staggered on 2D checkerboard meshes : scheme order for pressure with scaling is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVPStag on 2D cross triangle meshes : scheme order for velocity without scaling is ", -au
+ print( "FVPStag on 2D cross triangle meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVPStag on 2D cross triangle meshes : scheme order for velocity with scaling is ", -au
+ print( "FVPStag on 2D cross triangle meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV PStag on 2D deformed quadrangles meshes : scheme order for pressure is ", -ap
+ print( "FV PStag on 2D deformed quadrangles meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV PStag on 2D deformed quadrangles meshes : scheme order for velocity is ", -au
+ print( "FV PStag on 2D deformed quadrangles meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVPStag on 2D Delaunay triangles meshes : scheme order for velocity without scaling is ", -au
+ print( "FVPStag on 2D Delaunay triangles meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVPStag on 2D Delaunay triangles meshes : scheme order for velocity with scaling is ", -au
+ print( "FVPStag on 2D Delaunay triangles meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVPStag on 2D flat cross triangle meshes : scheme order for velocity without scaling is ", -au
+ print( "FVPStag on 2D flat cross triangle meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVPStag on 2D flat cross triangle meshes : scheme order for velocity with scaling is ", -au
+ print( "FVPStag on 2D flat cross triangle meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV PStag on 2D hexagonal meshes : scheme order for pressure is ", -ap
+ print( "FV PStag on 2D hexagonal meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV PStag on 2D hexagonal meshes : scheme order for velocity is ", -au
+ print( "FV PStag on 2D hexagonal meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
bp=(-a2*b1p+a1*b2p)/det
if(scaling==0):
- print "FV PStag on 2D square meshes : scheme order for pressure without scaling is ", -ap
+ print( "FV PStag on 2D square meshes : scheme order for pressure without scaling is ", -ap)
else:
- print "FV PStag on 2D square meshes : scheme order for pressure with scaling is ", -ap
+ print( "FV PStag on 2D square meshes : scheme order for pressure with scaling is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
bu=(-a2*b1u+a1*b2u)/det
if(scaling==0):
- print "FVPStag on 2D square meshes : scheme order for velocity without scaling is ", -au
+ print( "FVPStag on 2D square meshes : scheme order for velocity without scaling is ", -au)
else:
- print "FVPStag on 2D square meshes : scheme order for velocity with scaling is ", -au
+ print( "FVPStag on 2D square meshes : scheme order for velocity with scaling is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D brick wall meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D brick wall meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D brick wall meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D brick wall meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D checkerboard meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D checkerboard meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D checkerboard meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D checkerboard meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D cross triangles meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D cross triangles meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D cross triangles meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D cross triangles meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D deformed quadrangles meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D deformed quadrangles meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D deformed quadrangles meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D deformed quadrangles meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D triangle meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D triangle meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D triangle meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D triangle meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D flat cross triangles meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D flat cross triangles meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D flat cross triangles meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D flat cross triangles meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D hexagonal meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D hexagonal meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV upwind on 2D hexagonal meshes : scheme order for velocity is ", -au
+ print( "FV upwind on 2D hexagonal meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
ap=( a3*b1p-a2*b2p)/det
bp=(-a2*b1p+a1*b2p)/det
- print "FV upwind on 2D square meshes : scheme order for pressure is ", -ap
+ print( "FV upwind on 2D square meshes : scheme order for pressure is ", -ap)
b1u=np.dot(error_u_tab,mesh_size_tab)
b2u=np.sum(error_u_tab)
au=( a3*b1u-a2*b2u)/det
bu=(-a2*b1u+a1*b2u)/det
- print "FV on 2D square meshes : scheme order for velocity is ", -au
+ print( "FV on 2D square meshes : scheme order for velocity is ", -au)
# Plot of convergence curves
plt.close()
#Sauvegardes
if(it==1 or it%output_freq==0 or it>=ntmax or isStationary or time >=tmax):
- print"-- Iter: " + str(it) + ", Time: " + str(time) + ", dt: " + str(dt)
+ print("-- Iter: " + str(it) + ", Time: " + str(time) + ", dt: " + str(dt))
if(with_source):
print("Variation temporelle relative : pressure ", maxVector[0] ,", velocity x", maxVector[1]/rho0 ,", velocity y", maxVector[2]/rho0)
else: