diff --git a/src/mesh/FEM_quadrature.f90 b/src/mesh/FEM_quadrature.f90
index 4cc5e5468..79f02c84d 100644
--- a/src/mesh/FEM_quadrature.f90
+++ b/src/mesh/FEM_quadrature.f90
@@ -224,6 +224,7 @@ pure function permutationStar111(point) result(qPt)
 
   real(pReal), dimension(3,6) :: temp
 
+  temp = 0.0_pReal
   temp(:,1) = [point(1), point(2), 1.0_pReal - point(1) - point(2)]
   temp(:,2) = [point(1), 1.0_pReal - point(1) - point(2), point(2)]
   temp(:,4) = [point(2), 1.0_pReal - point(1) - point(2), point(1)]