normalFace1[2]=Ci.getNormalVector(1,2)
normalCell = normalFace0.crossProduct(normalFace1)
- normalCell = normalCell/normalCell.norm()
+ normalCell = normalCell*(1/normalCell.norm())
cellMat=cdmath.Matrix(4)
cellMat[0,0]=N0.x()
#=================================
LS=cdmath.LinearSolver(Rigidite,RHS,100,1.E-6,"GMRES","ILU")
LS.setMatrixIsSingular()#En raison de l'absence de bord
- LS.setComputeConditionNumber()
SolSyst=LS.solve()
print("Preconditioner used : ", LS.getNameOfPc() )
test_desc["Linear_solver_maximum_iterations"]=LS.getNumberMaxOfIter()
test_desc["Linear_system_max_actual_iterations_number"]=LS.getNumberOfIter()
test_desc["Linear_system_max_actual_error"]=LS.getResidu()
- test_desc["Linear_system_max_actual_condition number"]=LS.getConditionNumber()
# Création du champ résultat
#===========================
normalFace1[2]=Ci.getNormalVector(1,2)
normalCell = normalFace0.crossProduct(normalFace1)
- normalCell = normalCell/normalCell.norm()
+ normalCell = normalCell*(1/normalCell.norm())
cellMat=cdmath.Matrix(4)
cellMat[0,0]=N0.x()
#=================================
LS=cdmath.LinearSolver(Rigidite,RHS,100,1.E-6,"GMRES","ILU")
LS.setMatrixIsSingular()#En raison de l'absence de bord
- LS.setComputeConditionNumber()
SolSyst=LS.solve()
print( "Preconditioner used : ", LS.getNameOfPc() )
test_desc["Linear_solver_maximum_iterations"]=LS.getNumberMaxOfIter()
test_desc["Linear_system_max_actual_iterations_number"]=LS.getNumberOfIter()
test_desc["Linear_system_max_actual_error"]=LS.getResidu()
- test_desc["Linear_system_max_actual_condition number"]=LS.getConditionNumber()
# Création du champ résultat
#===========================
normalFace1[2]=Ci.getNormalVector(1,2)
normalCell = normalFace0.crossProduct(normalFace1)
- normalCell = normalCell/normalCell.norm()
+ normalCell = normalCell*(1/normalCell.norm())
cellMat=cdmath.Matrix(4)
cellMat[0,0]=N0.x()
#=================================
LS=cdmath.LinearSolver(Rigidite,RHS,100,1.E-6,"CG","CHOLESKY")
LS.setMatrixIsSingular()#En raison de l'absence de bord
- LS.setComputeConditionNumber()
SolSyst=LS.solve()
print( "Preconditioner used : ", LS.getNameOfPc() )
test_desc["Linear_solver_maximum_iterations"]=LS.getNumberMaxOfIter()
test_desc["Linear_system_max_actual_iterations_number"]=LS.getNumberOfIter()
test_desc["Linear_system_max_actual_error"]=LS.getResidu()
- test_desc["Linear_system_max_actual_condition number"]=LS.getConditionNumber()
# Création du champ résultat
#===========================