LS=cdmath.LinearSolver(divMat,Un,iterGMRESMax, precision, "GMRES","ILU")
else:
LS=cdmath.LinearSolver(divMat,Vn,iterGMRESMax, precision, "GMRES","ILU")
- LS.setComputeConditionNumber()
test_desc["Linear_solver_algorithm"]=LS.getNameOfMethod()
test_desc["Linear_solver_preconditioner"]=LS.getNameOfPc()
test_desc["Linear_system_max_actual_iterations_number"]=max(LS.getNumberOfIter(),test_desc["Linear_system_max_actual_iterations_number"])
test_desc["Linear_system_max_actual_error"]=max(LS.getResidu(),test_desc["Linear_system_max_actual_error"])
- test_desc["Linear_system_max_actual_condition number"]=max(LS.getConditionNumber(),test_desc["Linear_system_max_actual_condition number"])
max_dp=0 ; max_dq=0
for k in range(nbCells):
velocity_field.writeVTK("WaveSystem"+str(dim)+"DStaggered"+meshName+"_velocity",False);
print("Ecart au stationnaire exact : error_p= ",delta_press/p0," error_||u||= ",delta_v.maxVector()[0])
- print
- print"-- Iter: " + str(it) + ", Time: " + str(time) + ", dt: " + str(dt)
+ print( )
+ print("-- Iter: " + str(it) + ", Time: " + str(time) + ", dt: " + str(dt) )
print("Variation temporelle relative : pressure ", max_dp/p0 ,", velocity ", max_dq/rho0 )
- print
+ print( )
if(it>=ntmax):
print("Nombre de pas de temps maximum ntmax= ", ntmax, " atteint")