]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix of get_position_vector
authorMarco Tezzele <marcotez@gmail.com>
Wed, 17 Jun 2015 14:51:22 +0000 (16:51 +0200)
committerMarco Tezzele <marcotez@gmail.com>
Wed, 17 Jun 2015 14:52:18 +0000 (16:52 +0200)
include/deal.II/numerics/vector_tools.templates.h

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.