WIP tests for polynomial reconstruction of degree 3 [ci-skip]

Seems buggy in 3d for degree 3+

src/scheme/PolynomialReconstruction.cpp

@@ -1384,6 +1384,22 @@ PolynomialReconstruction::_build(
           }
         } else {
           Givens::solveCollectionInPlace(A, X, B);
+#warning REMOVE
+          if (cell_j_id == 0) {
+            for (size_t l = 0; l < B.numberOfRows(); ++l) {
+              for (size_t i = 0; i < B.numberOfColumns(); ++i) {
+                std::clog << "B(" << l << ", " << i << ") =" << B(l, i) << '\n';
+                ;
+              }
+            }
+
+            for (size_t l = 0; l < X.numberOfRows(); ++l) {
+              for (size_t i = 0; i < X.numberOfColumns(); ++i) {
+                std::clog << "X(" << l << ", " << i << ") =" << X(l, i) << '\n';
+                ;
+              }
+            }
+          }
         }
 
         column_begin = 0;

tests/CMakeLists.txt

@@ -281,6 +281,7 @@ add_executable (mpi_unit_tests
   test_Partitioner.cpp
   test_PolynomialReconstruction_degree_1.cpp
   test_PolynomialReconstruction_degree_2.cpp
+  test_PolynomialReconstruction_degree_3.cpp
   test_PolynomialReconstructionDescriptor.cpp
   test_RandomEngine.cpp
   test_StencilBuilder_cell2cell.cpp