From 0dafc63e077a52c40d65b8b52a99da3b307aedf8 Mon Sep 17 00:00:00 2001
From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Thu, 21 Feb 2019 12:24:21 +0100
Subject: [PATCH] Fix face orientation in FEFaceEvaluation

---
 include/deal.II/matrix_free/fe_evaluation.h | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index e0143f51f7..ac0e9be18f 100644
--- 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];
           }
-- 
2.39.5