rge remote-tracking branch 'origin/Package' into Coverage

.gitignore

@@ -4,6 +4,11 @@
 *.ipynb_checkpoints
 # ignore python cache
 *__pycache__*
+*.__pycache__
+*test/__pycache__
+*src/__pycache__
+*examples/__pycache__
+*examples/ReproduableCode/__pycache__
 *.pytest_cache*
 # ignore Sphinx documentation output
 *docs/build

src/Eq02.py

@@ -13,7 +13,7 @@ from dolfinx.nls.petsc import NewtonSolver
 from ufl import TestFunctions, split, dot, grad, dx, inner, Mesh, exp
 from basix.ufl import element, mixed_element
 import matplotlib.pyplot as plt
-from RefinedMesh1D import create_refined_mesh
+from src.RefinedMesh1D import create_refined_mesh
 
 
 def solve_System_2eq(phi_left:float, phi_right:float, p_right:float, z_A:float, z_C:float, y_A_R:float, y_C_R:float, K:float|str, Lambda2:float, a2:float, number_cells:int, solvation:float = 0, PoissonBoltzmann:bool=False, relax_param:float=None, x0:float=0, x1:float=1, refinement_style:str='uniform', return_type:str='Scalar', rtol:float=1e-8, max_iter:float=500):

src/Eq04.py

@@ -13,7 +13,7 @@ from dolfinx.nls.petsc import NewtonSolver
 from ufl import TestFunctions, split, dot, grad, dx, inner, ln, Mesh
 from basix.ufl import element, mixed_element
 import matplotlib.pyplot as plt
-from RefinedMesh1D import create_refined_mesh
+from src.RefinedMesh1D import create_refined_mesh
 
 
 def solve_System_4eq(phi_left:float, phi_right:float, p_right:float, z_A:float, z_C:float, y_A_R:float, y_C_R:float, K:float|str, Lambda2:float, a2:float, number_cells:int, solvation:float = 0, PoissonBoltzmann:bool=False, relax_param:float=None, x0:float=0, x1:float=1, refinement_style:str='uniform', return_type:str='Scalar', rtol:float=1e-8, max_iter:float=500):

src/EqN.py

@@ -13,7 +13,7 @@ from dolfinx.nls.petsc import NewtonSolver
 from ufl import TestFunctions, split, dot, grad, dx, inner, ln, Mesh
 from basix.ufl import element, mixed_element
 import matplotlib.pyplot as plt
-from RefinedMesh1D import create_refined_mesh
+from src.RefinedMesh1D import create_refined_mesh
 
 def solve_System_Neq(phi_left:float, phi_right:float, p_right:float, z_alpha:list, y_R:list, K:float|str, Lambda2:float, a2:float, number_cells:int, solvation:float = 0, PoissonBoltzmann:bool=False, relax_param:float=None, x0:float=0, x1:float=1, refinement_style:str='uniform', return_type:str='Vector', rtol:float=1e-8, max_iter:float=500):
     '''

tests/test_ElectrolyticDiode.py

@@ -28,7 +28,7 @@ refinement_style = 'uniform'
 solvation = 5
 PoissonBoltzmann = False
 rtol = 1e-3
-relax_param = 0.1
+relax_param = 0.08
 max_iter = 1_000