From da1d8f941525124e44f329139688bcdf9d319aff Mon Sep 17 00:00:00 2001 From: ndjinga Date: Wed, 21 Oct 2020 11:09:39 +0200 Subject: [PATCH] Updated documentation (setValue->addValue) --- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../2DDiffusionVF_brickwall/FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../2DDiffusionVF_hexagons/FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../2DDiffusionVF_loc_ref/FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../FiniteVolumes2DDiffusion_SQUARE.py | 2 +- .../2DDiffusionVF_squares/FiniteVolumes2DDiffusion_SQUARE.py | 2 +- CDMATH/tests/doc/2DDiffusionVF/CodeFV5DiffusionProblem.md | 2 +- .../2DPoissonVF_brickwall/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_checkerboard/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../FiniteVolumes2DWithCDMATH.py | 2 +- .../FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_flat_triangles/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_hexagons/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_loc_ref/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_loc_ref/FiniteVolumes2DWithCDMATH.py | 2 +- .../FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_long_triangles/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_squares/FiniteVolumes2DPoisson_SQUARE.py | 2 +- .../2DPoissonVF_triangles/FiniteVolumes2DPoisson_SQUARE.py | 2 +- CDMATH/tests/doc/2DPoissonVF/CodeFV5PoissonProblem.md | 2 +- .../FiniteVolumes2DPoissonStiffBC_DISK.py | 2 +- .../FiniteVolumes2DPoissonStiffBC_DISK.py | 2 +- .../FiniteVolumes2DPoissonStiffBC_DISK.py | 2 +- .../FiniteVolumes2DPoissonStiffBC_DISK.py | 2 +- .../FiniteVolumes2DPoissonStiffBC_SQUARE.py | 2 +- .../FiniteVolumes2DPoissonStiffBC_SQUARE.py | 2 +- .../FiniteVolumes2DPoissonStiffBC_SQUARE.py | 2 +- 35 files changed, 35 insertions(+), 35 deletions(-) diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_Delaunay_triangles/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_Delaunay_triangles/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_Delaunay_triangles/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_Delaunay_triangles/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_brickwall/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_brickwall/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_brickwall/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_brickwall/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_checkerboard/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_checkerboard/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_checkerboard/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_checkerboard/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_deformedQuadrangles/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_deformedQuadrangles/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_deformedQuadrangles/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_deformedQuadrangles/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_flat_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_flat_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_flat_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_flat_cross_triangles/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_hexagons/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_hexagons/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_hexagons/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_hexagons/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_loc_ref/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_loc_ref/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_loc_ref/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_loc_ref/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_long_rectangles/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_long_rectangles/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_long_rectangles/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_long_rectangles/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_skinny_triangles/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_skinny_triangles/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_skinny_triangles/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_skinny_triangles/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_squares/FiniteVolumes2DDiffusion_SQUARE.py b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_squares/FiniteVolumes2DDiffusion_SQUARE.py index 71e6155..f7adba8 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_squares/FiniteVolumes2DDiffusion_SQUARE.py +++ b/CDMATH/tests/doc/2DDiffusionVF/2DDiffusionVF_squares/FiniteVolumes2DDiffusion_SQUARE.py @@ -88,7 +88,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DDiffusionVF/CodeFV5DiffusionProblem.md b/CDMATH/tests/doc/2DDiffusionVF/CodeFV5DiffusionProblem.md index 531d473..2404c5c 100644 --- a/CDMATH/tests/doc/2DDiffusionVF/CodeFV5DiffusionProblem.md +++ b/CDMATH/tests/doc/2DDiffusionVF/CodeFV5DiffusionProblem.md @@ -35,7 +35,7 @@ for i in range(nbCells): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance*(normal[0]*normal[0] + K*normal[1]*normal[1]) - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter())*(normal[0]*normal[0] + K*normal[1]*normal[1]) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_brickwall/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_brickwall/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_brickwall/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_brickwall/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_checkerboard/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_checkerboard/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_checkerboard/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_checkerboard/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_cross_triangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_cross_triangles/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_cross_triangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_cross_triangles/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DWithCDMATH.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DWithCDMATH.py index 21b2c16..7796678 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DWithCDMATH.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_deformedQuadrangles/FiniteVolumes2DWithCDMATH.py @@ -83,7 +83,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) Rigidite.addValue(i,i,coeff) # terme diagonal diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_equilateral_triangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_equilateral_triangles/FiniteVolumes2DPoisson_SQUARE.py index a8aee40..f2b44a9 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_equilateral_triangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_equilateral_triangles/FiniteVolumes2DPoisson_SQUARE.py @@ -83,7 +83,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_flat_triangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_flat_triangles/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_flat_triangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_flat_triangles/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_hexagons/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_hexagons/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_hexagons/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_hexagons/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DWithCDMATH.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DWithCDMATH.py index 21b2c16..7796678 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DWithCDMATH.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_loc_ref/FiniteVolumes2DWithCDMATH.py @@ -83,7 +83,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) Rigidite.addValue(i,i,coeff) # terme diagonal diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_rectangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_rectangles/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_rectangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_rectangles/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_triangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_triangles/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_triangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_long_triangles/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_right_triangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_right_triangles/FiniteVolumes2DPoisson_SQUARE.py index 83c5f47..dbc9a84 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_right_triangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_right_triangles/FiniteVolumes2DPoisson_SQUARE.py @@ -83,7 +83,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_squares/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_squares/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_squares/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_squares/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_triangles/FiniteVolumes2DPoisson_SQUARE.py b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_triangles/FiniteVolumes2DPoisson_SQUARE.py index cb857fe..5203cc7 100644 --- a/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_triangles/FiniteVolumes2DPoisson_SQUARE.py +++ b/CDMATH/tests/doc/2DPoissonVF/2DPoissonVF_triangles/FiniteVolumes2DPoisson_SQUARE.py @@ -80,7 +80,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoissonVF/CodeFV5PoissonProblem.md b/CDMATH/tests/doc/2DPoissonVF/CodeFV5PoissonProblem.md index 155e226..98cf392 100644 --- a/CDMATH/tests/doc/2DPoissonVF/CodeFV5PoissonProblem.md +++ b/CDMATH/tests/doc/2DPoissonVF/CodeFV5PoissonProblem.md @@ -32,7 +32,7 @@ for i in range(nbCells): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_DelaunayTriangles/FiniteVolumes2DPoissonStiffBC_DISK.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_DelaunayTriangles/FiniteVolumes2DPoissonStiffBC_DISK.py index 2554b72..d92d8f9 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_DelaunayTriangles/FiniteVolumes2DPoissonStiffBC_DISK.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_DelaunayTriangles/FiniteVolumes2DPoissonStiffBC_DISK.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_hexagons/FiniteVolumes2DPoissonStiffBC_DISK.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_hexagons/FiniteVolumes2DPoissonStiffBC_DISK.py index ff3c673..8e6e967 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_hexagons/FiniteVolumes2DPoissonStiffBC_DISK.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_hexagons/FiniteVolumes2DPoissonStiffBC_DISK.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshName, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_spiderWeb/FiniteVolumes2DPoissonStiffBC_DISK.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_spiderWeb/FiniteVolumes2DPoissonStiffBC_DISK.py index 2554b72..d92d8f9 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_spiderWeb/FiniteVolumes2DPoissonStiffBC_DISK.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_spiderWeb/FiniteVolumes2DPoissonStiffBC_DISK.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_squares/FiniteVolumes2DPoissonStiffBC_DISK.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_squares/FiniteVolumes2DPoissonStiffBC_DISK.py index 2554b72..d92d8f9 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_squares/FiniteVolumes2DPoissonStiffBC_DISK.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_DISK_squares/FiniteVolumes2DPoissonStiffBC_DISK.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshType, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_equilateralTriangles/FiniteVolumes2DPoissonStiffBC_SQUARE.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_equilateralTriangles/FiniteVolumes2DPoissonStiffBC_SQUARE.py index c2cbccf..7822a02 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_equilateralTriangles/FiniteVolumes2DPoissonStiffBC_SQUARE.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_equilateralTriangles/FiniteVolumes2DPoissonStiffBC_SQUARE.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshName, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_hexagons/FiniteVolumes2DPoissonStiffBC_SQUARE.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_hexagons/FiniteVolumes2DPoissonStiffBC_SQUARE.py index c2cbccf..7822a02 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_hexagons/FiniteVolumes2DPoissonStiffBC_SQUARE.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_hexagons/FiniteVolumes2DPoissonStiffBC_SQUARE.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshName, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary diff --git a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_squares/FiniteVolumes2DPoissonStiffBC_SQUARE.py b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_squares/FiniteVolumes2DPoissonStiffBC_SQUARE.py index c2cbccf..7822a02 100644 --- a/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_squares/FiniteVolumes2DPoissonStiffBC_SQUARE.py +++ b/CDMATH/tests/doc/2DPoisson_StiffBC/2DPoissonVF_StiffBC_DISK/2DPoissonVF_StiffBC_SQUARE_squares/FiniteVolumes2DPoissonStiffBC_SQUARE.py @@ -97,7 +97,7 @@ def solve(my_mesh,filename,resolution, meshName, testColor): Ck=my_mesh.getCell(k) distance=Ci.getBarryCenter().distance(Ck.getBarryCenter()) coeff=Fj.getMeasure()/Ci.getMeasure()/distance - Rigidite.setValue(i,k,-coeff) # terme extradiagonal + Rigidite.addValue(i,k,-coeff) # terme extradiagonal else: coeff=Fj.getMeasure()/Ci.getMeasure()/Ci.getBarryCenter().distance(Fj.getBarryCenter()) #For the particular case where the mesh boundary does not coincide with the domain boundary -- 2.39.2