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 af9d004f543434ec3c01f99a46a40d5181dc9d0e..569cda39a37133e7d400d1cb70b2581990bb6bcf 100644 (file)
--- 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);
       }
   }