]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
try to fix petsc FETools::interpolation_difference
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 1 Mar 2013 16:34:21 +0000 (16:34 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 1 Mar 2013 16:34:21 +0000 (16:34 +0000)
git-svn-id: https://svn.dealii.org/trunk@28693 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/fe/fe_tools.cc

index 26ed3a42fa02dd22feb85a5f149ee4766192e3c9..dd0e3dd74c78f93fa392e9e65be82eb2fec3ef62 100644 (file)
@@ -1671,10 +1671,13 @@ namespace FETools
             DoFTools::extract_locally_relevant_dofs (dof2,
                                                      dof2_locally_relevant_dofs);
 
+            PETScWrappers::MPI::Vector  u2_out (dynamic_cast<const PETScWrappers::MPI::Vector *> (&u1)->get_mpi_communicator(),
+                dof2_locally_owned_dofs);
+            interpolate(dof1, u1, dof2, constraints2, u2_out);
             PETScWrappers::MPI::Vector  u2 (dynamic_cast<const PETScWrappers::MPI::Vector *> (&u1)->get_mpi_communicator(),
                                             dof2_locally_owned_dofs,
                                             dof2_locally_relevant_dofs);
-            interpolate(dof1, u1, dof2, constraints2, u2);
+            u2 = u2_out;
             interpolate(dof2, u2, dof1, constraints1, u1_interpolated);
           }
         else

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.