B to primal

Author: PhilipDeegan

src/core/data/grid/gridlayoutdefs.hpp

@@ -1,11 +1,11 @@
 #ifndef PHARE_CORE_GRID_GRIDLAYOUTDEFS_HPP
 #define PHARE_CORE_GRID_GRIDLAYOUTDEFS_HPP
 
-#include <cstddef>
 
-#include "core/hybrid/hybrid_quantities.hpp"
-#include "core/utilities/types.hpp"
 #include "core/utilities/point/point.hpp"
+#include "core/hybrid/hybrid_quantities.hpp"
+
+#include <cstddef>
 
 namespace PHARE
 {

src/core/data/grid/gridlayoutimplyee.hpp

@@ -528,15 +528,16 @@ namespace core
             }
             else if constexpr (dimension == 3)
             {
-                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0, 0}, 0.5};
-                constexpr WeightPoint<dimension> P2{Point<int, dimension>{0, iShift, iShift}, 0.5};
-                return std::array<WeightPoint<dimension>, 2>{P1, P2};
+                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0, 0}, 0.25};
+                constexpr WeightPoint<dimension> P2{Point<int, dimension>{0, iShift, 0}, 0.25};
+                constexpr WeightPoint<dimension> P3{Point<int, dimension>{0, 0, iShift}, 0.25};
+                constexpr WeightPoint<dimension> P4{Point<int, dimension>{iShift, iShift, 0}, 0.25};
+                return std::array<WeightPoint<dimension>, 2>{P1, P2, P3, P4};
             }
         }
 
 
 
-
         NO_DISCARD auto static constexpr ByToMoments()
         {
             // By is dual primal dual
@@ -559,9 +560,11 @@ namespace core
             }
             else if constexpr (dimension == 3)
             {
-                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0, 0}, 0.5};
-                constexpr WeightPoint<dimension> P2{Point<int, dimension>{iShift, 0, iShift}, 0.5};
-                return std::array<WeightPoint<dimension>, 2>{P1, P2};
+                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0, 0}, 0.25};
+                constexpr WeightPoint<dimension> P2{Point<int, dimension>{iShift, 0, 0}, 0.25};
+                constexpr WeightPoint<dimension> P3{Point<int, dimension>{0, 0, iShift}, 0.25};
+                constexpr WeightPoint<dimension> P4{Point<int, dimension>{iShift, 0, iShift}, 0.25};
+                return std::array<WeightPoint<dimension>, 2>{P1, P2, P3, P4};
             }
         }
 
@@ -584,15 +587,19 @@ namespace core
             }
             if constexpr (dimension == 2)
             {
-                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0}, 0.5};
-                constexpr WeightPoint<dimension> P2{Point<int, dimension>{iShift, iShift}, 0.5};
-                return std::array<WeightPoint<dimension>, 2>{P1, P2};
+                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0}, 0.25};
+                constexpr WeightPoint<dimension> P2{Point<int, dimension>{iShift, 0}, 0.25};
+                constexpr WeightPoint<dimension> P3{Point<int, dimension>{0, iShift}, 0.25};
+                constexpr WeightPoint<dimension> P4{Point<int, dimension>{iShift, iShift}, 0.25};
+                return std::array<WeightPoint<dimension>, 2>{P1, P2, P3, P4};
             }
             else if constexpr (dimension == 3)
             {
-                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0, 0}, 0.5};
-                constexpr WeightPoint<dimension> P2{Point<int, dimension>{iShift, iShift, 0}, 0.5};
-                return std::array<WeightPoint<dimension>, 2>{P1, P2};
+                constexpr WeightPoint<dimension> P1{Point<int, dimension>{0, 0, 0}, 0.25};
+                constexpr WeightPoint<dimension> P2{Point<int, dimension>{iShift, 0, 0}, 0.25};
+                constexpr WeightPoint<dimension> P3{Point<int, dimension>{0, iShift, 0}, 0.25};
+                constexpr WeightPoint<dimension> P4{Point<int, dimension>{iShift, iShift, 0}, 0.25};
+                return std::array<WeightPoint<dimension>, 2>{P1, P2, P3, P4};
             }
         }
 
@@ -818,7 +825,8 @@ namespace core
 
 
         NO_DISCARD auto static constexpr ByToEx()
-        { // By is dual primal dual
+        {
+            // By is dual primal dual
             // Ex is dual primal primal
             // operation is thus dpD to dpP
             // shift only in the Z direction