Work around an issue with the clang CUDA compiler.

diff --git a/include/deal.II/numerics/vector_tools_constraints.templates.h b/include/deal.II/numerics/vector_tools_constraints.templates.h
index b5a2f681869c48a82914c9cfc7707fd4fb096c6e..9785674a98954cf0e2462405cfbf6408ff477c25 100644 (file)
--- a/include/deal.II/numerics/vector_tools_constraints.templates.h
+++ b/include/deal.II/numerics/vector_tools_constraints.templates.h
@@ -242,15 +242,17 @@ namespace VectorTools
                                     constraining_vector[0]) >
                           std::numeric_limits<double>::epsilon())
                         constraint_entries.emplace_back(
-                          dof_indices.dof_indices[1],
-                          -constraining_vector[1] / constraining_vector[0]);
+                          std::make_pair(dof_indices.dof_indices[1],
+                                         -constraining_vector[1] /
+                                           constraining_vector[0]));
 
                       if (std::fabs(constraining_vector[2] /
                                     constraining_vector[0]) >
                           std::numeric_limits<double>::epsilon())
                         constraint_entries.emplace_back(
-                          dof_indices.dof_indices[2],
-                          -constraining_vector[2] / constraining_vector[0]);
+                          std::make_pair(dof_indices.dof_indices[2],
+                                         -constraining_vector[2] /
+                                           constraining_vector[0]));
 
                       const double normalized_inhomogeneity =
                         (std::fabs(inhomogeneity / constraining_vector[0]) >
@@ -275,15 +277,17 @@ namespace VectorTools
                                     constraining_vector[1]) >
                           std::numeric_limits<double>::epsilon())
                         constraint_entries.emplace_back(
-                          dof_indices.dof_indices[0],
-                          -constraining_vector[0] / constraining_vector[1]);
+                          std::make_pair(dof_indices.dof_indices[0],
+                                         -constraining_vector[0] /
+                                           constraining_vector[1]));
 
                       if (std::fabs(constraining_vector[2] /
                                     constraining_vector[1]) >
                           std::numeric_limits<double>::epsilon())
                         constraint_entries.emplace_back(
-                          dof_indices.dof_indices[2],
-                          -constraining_vector[2] / constraining_vector[1]);
+                          std::make_pair(dof_indices.dof_indices[2],
+                                         -constraining_vector[2] /
+                                           constraining_vector[1]));
 
                       const double normalized_inhomogeneity =
                         (std::fabs(inhomogeneity / constraining_vector[1]) >
@@ -305,15 +309,17 @@ namespace VectorTools
                                     constraining_vector[2]) >
                           std::numeric_limits<double>::epsilon())
                         constraint_entries.emplace_back(
-                          dof_indices.dof_indices[0],
-                          -constraining_vector[0] / constraining_vector[2]);
+                          std::make_pair(dof_indices.dof_indices[0],
+                                         -constraining_vector[0] /
+                                           constraining_vector[2]));
 
                       if (std::fabs(constraining_vector[1] /
                                     constraining_vector[2]) >
                           std::numeric_limits<double>::epsilon())
                         constraint_entries.emplace_back(
-                          dof_indices.dof_indices[1],
-                          -constraining_vector[1] / constraining_vector[2]);
+                          std::make_pair(dof_indices.dof_indices[1],
+                                         -constraining_vector[1] /
+                                           constraining_vector[2]));
 
                       const double normalized_inhomogeneity =
                         (std::fabs(inhomogeneity / constraining_vector[2]) >