diff --git a/src/mesh/FEM_quadrature.f90 b/src/mesh/FEM_quadrature.f90
index 79f02c84d..5c5c3f260 100644
--- a/src/mesh/FEM_quadrature.f90
+++ b/src/mesh/FEM_quadrature.f90
@@ -224,9 +224,9 @@ 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(:,3) = 0.0_pReal
   temp(:,4) = [point(2), 1.0_pReal - point(1) - point(2), point(1)]
   temp(:,5) = [1.0_pReal - point(1) - point(2), point(2), point(1)]
   temp(:,6) = [1.0_pReal - point(1) - point(2), point(1), point(2)]