rebase/conflicts

PHAREHUB · Jan 22, 2022 · a47bdf3 · a47bdf3
1 parent 861da38
commit a47bdf3
Show file tree

Hide file tree

Showing 7 changed files with 10 additions and 59 deletions.
diff --git a/src/amr/data/particles/refine/split.h b/src/amr/data/particles/refine/split.h
diff --git a/src/amr/data/particles/refine/split.hpp b/src/amr/data/particles/refine/split.hpp
@@ -3,7 +3,7 @@
 
 #include "amr/data/particles/refine/split_1d.hpp"
 #include "amr/data/particles/refine/split_2d.hpp"
-#include "amr/data/particles/refine/split_3d.h"
+#include "amr/data/particles/refine/split_3d.hpp"
 
 
 #endif // endif PHARE_SPLIT_HPP
diff --git a/src/amr/data/particles/refine/split_3d.hpp b/src/amr/data/particles/refine/split_3d.hpp
@@ -8,9 +8,9 @@ Splitting reference material can be found @
 
 #include <array>
 #include <cstddef>
-#include "core/utilities/point/point.h"
-#include "core/utilities/types.h"
-#include "splitter.h"
+#include "core/utilities/point/point.hpp"
+#include "core/utilities/types.hpp"
+#include "splitter.hpp"
 
 namespace PHARE::amr
 {

diff --git a/tests/simulator/advance/test_fields_advance_2d.py b/tests/simulator/advance/test_fields_advance_2d.py
@@ -63,7 +63,7 @@ def test_overlaped_fields_are_equal_with_min_max_patch_size_of_max_ghosts(self,
        *per_interp(({"L0": {"B0": Box2D(6, 23)}})),
        *per_interp(({"L0": {"B0": Box2D( 2, 12), "B1": Box2D(13, 25)}})),
        *per_interp(({"L0": {"B0": Box2D( 5, 20)}, "L1": {"B0": Box2D(15, 19)}})),
-       *per_interp(({"L0": {"B0": Box2D( 5, 20)}, "L1": {"B0": Box2D(12, 38)}, "L2": {"B0": Box2D(30, 52)} })),
+       *per_interp(({"L0": {"B0": Box2D( 5, 20)}, "L1": {"B0": Box2D(12, 37)}, "L2": {"B0": Box2D(30, 51)} })),
     )
     @unpack
     def test_field_coarsening_via_subcycles(self, interp_order, refinement_boxes):

diff --git a/tests/simulator/advance/test_particles_advance_3d.py b/tests/simulator/advance/test_particles_advance_3d.py
@@ -14,7 +14,7 @@
 
 ndim = 3
 interp_orders = [1, 2, 3]
-ppc = 10
+ppc = 5
 
 def per_interp(dic):
     return [(interp, dic) for interp in interp_orders]

diff --git a/tests/simulator/test_advance.py b/tests/simulator/test_advance.py
@@ -36,7 +36,6 @@ def getHierarchy(self, interp_order, refinement_boxes, qty,
                      dl=0.2, extra_diag_options={}, time_step_nbr=1, timestamps=None, ndim=1):
         diag_outputs = f"phare_outputs/advance/{diag_outputs}"
         from pyphare.pharein import global_vars
-        clean_up_diags(global_vars.sim)
         global_vars.sim = None
         if smallest_patch_size is None:
             from pyphare.pharein.simulation import check_patch_size

diff --git a/tests/simulator/test_initialization.py b/tests/simulator/test_initialization.py
@@ -39,7 +39,6 @@ def getHierarchy(self, interp_order, refinement_boxes, qty,
                      cells=120, dl=0.1, ndim=1):
         diag_outputs = f"phare_outputs/init/{diag_outputs}"
         from pyphare.pharein import global_vars
-        clean_up_diags(global_vars.sim)
         global_vars.sim =None
 
         if smallest_patch_size is None:
@@ -328,16 +327,14 @@ def domain(patch_data):
 
 
 
-<<<<<<< HEAD
-
 
     def _test_density_is_as_provided_by_user(self, dim, interp_order, **kwargs):
         print(f"test_density_is_as_provided_by_user : dim {dim} interp_order {interp_order}")
 
         empirical_dim_devs = {
               1: 6e-3,
               2: 3e-2,
-              3: 3e-2,
+              3: 2e-1,
         }
         nbParts = {1 : 10000, 2: 1000, 3: 20}
         hier = self.getHierarchy(interp_order, {"L0": {"B0": nDBox(dim, 5, 14)}},
@@ -364,45 +361,9 @@ def _test_density_is_as_provided_by_user(self, dim, interp_order, **kwargs):
                 nbrGhosts = layout.nbrGhosts(interp_order, centering)
                 select = tuple([slice(nbrGhosts,-nbrGhosts) for i in range(dim)])
 
-                xyz = patch.patch_datas["rho"].meshgrid()
-                xyz0 = [v[select] for v in xyz]
-
-                # if dim == 1:
-                #     x0 = xx[select]
-                #     protons_expected = proton_density_fn(x0)
-                #     beam_expected    = beam_density_fn(x0)
-
-                #     # ion_expected     = protons_expected + beam_expected
-                #     # ion_actual     = ion_density[nbrGhosts:-nbrGhosts]
-                #     # beam_actual    = beam_density[nbrGhosts:-nbrGhosts]
-                #     # protons_actual = proton_density[nbrGhosts:-nbrGhosts]
-
-
-                # if dim == 2:
-                #     xx, yy = patch.patch_datas["rho"].meshgrid()
-
-                #     x0 = xx[select]
-                #     y0 = yy[select]
-
-                #     protons_expected = proton_density_fn(x0, y0)
-                #     beam_expected    = beam_density_fn(x0, y0)
-
-                #     # ion_expected     = protons_expected + beam_expected
-                #     # ion_actual     = ion_density[select]
-                #     # beam_actual    = beam_density[select]
-                #     # protons_actual = proton_density[select]
-
-
-
-                # if dim == 3:
-                #     xx, yy, zz = patch.patch_datas["rho"].meshgrid()
-
-                #     x0 = xx[select]
-                #     y0 = yy[select]
-                #     z0 = zz[select]
-
-                protons_expected = proton_density_fn(*xyz0)
-                beam_expected    = beam_density_fn(*xyz0)
+                mesh = patch.patch_datas["rho"].meshgrid(select)
+                protons_expected = proton_density_fn(*mesh)
+                beam_expected    = beam_density_fn(*mesh)
                 ion_expected   = protons_expected + beam_expected
 
                 protons_actual = proton_density[select]