From: Marco Tezzele <marcotez@gmail.com>
Date: Wed, 17 Jun 2015 14:51:22 +0000 (+0200)
Subject: fix of get_position_vector
X-Git-Tag: v8.3.0-rc1~105^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=1a8d9ff42bf32880c22b9382bcc8aea5d07cfc59;p=dealii.git

fix of get_position_vector
---

diff --git a/include/deal.II/numerics/vector_tools.templates.h b/include/deal.II/numerics/vector_tools.templates.h
index af9d004f54..569cda39a3 100644
--- a/include/deal.II/numerics/vector_tools.templates.h
+++ b/include/deal.II/numerics/vector_tools.templates.h
@@ -7195,9 +7195,11 @@ namespace VectorTools
               for (unsigned int j=0; j<points.size(); ++j)
                 {
                   if ( fe_to_real[comp_i] == feq.system_to_component_index(j).first)
-                    transfer(i, j) = fe.shape_value(i, points[j]);
+                    transfer(j, i) = fe.shape_value(i, points[j]);
                 }
           }
+
+        transfer.invert(transfer);
         interpolate(dhq, dh, transfer, eulerq, vector);
       }
   }