add solver validation + update comments

python/_tools.py

@@ -50,7 +50,7 @@ except :
 
 def get_linear_system_package(choices=None):
     if choices is None:
-        choices = ["mumps","petsc","scipy"]
+        choices = ["mumps","petsc","scipy","petsc4py"]
     if isinstance(choices, str):
         choices = [choices]
     for choice in choices:

python/fluid.py

@@ -139,8 +139,8 @@ class FluidProblem :
         density -- List of fluid phases densities (this should be a vector whom dimension specify if there is one or two fluid)
         sigma -- Surface tension (only when there are two fluids)
         coeff_stab -- Optional argument used as coefficient for extra diffusivity added to stabilise the advection of concentration (only for two fluid problem) 
-        solver -- possible solver are "mumps", "petsc", "scipy"
-        solver_options -- Optional argument to specify solver option (only when petsc is used)
+        solver -- possible solver are "mumps", "petsc", "petsc4py", "scipy"
+        solver_options -- Optional argument to specify solver option (petsc options or mumps parameters)
         drag_in_stab -- States if the drag force is in the stabilisation term
         drag_coefficient_factor -- Factor multiplicating the drag coefficient
         temporal -- Enables d/dt (i.e. False = steady)

validation/cavity/cavity.py

@@ -33,7 +33,7 @@ import unittest
 #Physical parameters for the drops are the ones presented by Metzger et al. (2007) "Falling clouds of particles in viscous fluids"
 
 
-class Poiseuille(unittest.TestCase) :
+class Cavity(unittest.TestCase) :
     def runTest(self) :
         dir_path = os.path.dirname(os.path.realpath(__file__))
         os.chdir(dir_path)

validation/periodic2d/periodic2d.py

@@ -73,7 +73,7 @@ class Periodic2d(unittest.TestCase):
         p.read_vtk(outputdir,0)
 
         #Object fluid creation + Boundary condition of the fluid (field 0 is horizontal velocity; field 1 is vertical velocity; field 2 is pressure)
-        fluid = mbfluid.FluidProblem(2,g,mu,rho)
+        fluid = mbfluid.FluidProblem(2,g,mu,rho,solver="petsc4py")
         fluid.load_msh("mesh.msh")
         fluid.set_wall_boundary("Bottom",velocity=[0,0], pressure = 0)
         fluid.set_wall_boundary("Top",velocity=[0,0])

validation/poiseuille/poiseuille.py

@@ -64,7 +64,7 @@ class Poiseuille(unittest.TestCase) :
         outf = 1                                       # number of iterations between output files
 
         #Object fluid creation + Boundary condition of the fluid (field 0 is horizontal velocity; field 1 is vertical velocity; field 2 is pressure)
-        fluid = mbfluid.FluidProblem(2,g,nu*rho,rho)
+        fluid = mbfluid.FluidProblem(2,g,nu*rho,rho,solver="scipy")
         fluid.load_msh("mesh.msh")
         fluid.set_open_boundary("LeftUp",velocity=[lambda x : 1/(20*mu)*x[:,1]*(1-x[:, 1]),0])
         fluid.set_open_boundary("LeftDown",velocity=[lambda x : 1/(20*mu)*x[:,1]*(1-x[:, 1]),0])