Fix face orientation in FEFaceEvaluation

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index e0143f51f754a6ec69605edaf4873bca51f1ff46..ac0e9be18f9df86ab3490b91169322f46b8e9be0 100644 (file)
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@ -8562,10 +8562,10 @@ FEFaceEvaluation<dim, fe_degree, n_q_points_1d, n_components, Number>::
         {
           if (integrate)
             for (unsigned int q = 0; q < n_q_points; ++q)
-              tmp_values[orientations[q]] = this->values_quad[c][q];
+              tmp_values[q] = this->values_quad[c][orientations[q]];
           else
             for (unsigned int q = 0; q < n_q_points; ++q)
-              tmp_values[q] = this->values_quad[c][orientations[q]];
+              tmp_values[orientations[q]] = this->values_quad[c][q];
           for (unsigned int q = 0; q < n_q_points; ++q)
             this->values_quad[c][q] = tmp_values[q];
         }
@@ -8574,10 +8574,10 @@ FEFaceEvaluation<dim, fe_degree, n_q_points_1d, n_components, Number>::
           {
             if (integrate)
               for (unsigned int q = 0; q < n_q_points; ++q)
-                tmp_values[orientations[q]] = this->gradients_quad[c][d][q];
+                tmp_values[q] = this->gradients_quad[c][d][orientations[q]];
             else
               for (unsigned int q = 0; q < n_q_points; ++q)
-                tmp_values[q] = this->gradients_quad[c][d][orientations[q]];
+                tmp_values[orientations[q]] = this->gradients_quad[c][d][q];
             for (unsigned int q = 0; q < n_q_points; ++q)
               this->gradients_quad[c][d][q] = tmp_values[q];
           }