diff --git a/tests/test_DiscreteFunctionP0.cpp b/tests/test_DiscreteFunctionP0.cpp
index 90ebb76d8a8d332558dd801e3c38f6b36ca62e90..7eb264566b916944c1800ca27b2e35a88d656561 100644
--- a/tests/test_DiscreteFunctionP0.cpp
+++ b/tests/test_DiscreteFunctionP0.cpp
@@ -2742,7 +2742,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
             parallel_for(
               mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
                 const double x = xj[cell_id][0];
-                Ah[cell_id]    = TinyMatrix<2>{2.3 * x, 1 - x, 2 * x, 1.5 * x};
+                Bh[cell_id]    = TinyMatrix<2>{2.3 * x, 1 - x, 2 * x, 1.5 * x};
               });
 
             CHECK_STD_BINARY_MATH_FUNCTION(Ah, Bh, doubleDot);
@@ -3178,7 +3178,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
             parallel_for(
               mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
                 const double x = xj[cell_id][0];
-                Ah[cell_id]    = TinyMatrix<2>{2.3 * x, 1 - x, 2 * x, 1.5 * x};
+                Bh[cell_id]    = TinyMatrix<2>{2.3 * x, 1 - x, 2 * x, 1.5 * x};
               });
 
             CHECK_STD_BINARY_MATH_FUNCTION(Ah, Bh, doubleDot);
@@ -3567,7 +3567,7 @@ TEST_CASE("DiscreteFunctionP0", "[scheme]")
             parallel_for(
               mesh->numberOfCells(), PUGS_LAMBDA(const CellId cell_id) {
                 const double x = xj[cell_id][0];
-                Ah[cell_id]    = TinyMatrix<2>{2.3 * x, 1 - x, 2 * x, 1.5 * x};
+                Bh[cell_id]    = TinyMatrix<2>{2.3 * x, 1 - x, 2 * x, 1.5 * x};
               });
 
             CHECK_STD_BINARY_MATH_FUNCTION(Ah, Bh, doubleDot);