By the way, it was this collection of lines from the stack trace that made me suspicious that there's a problem thats not actually directly related to the use of the MappingQEulerian class:
[pasha:03624] [18] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) > [0x7fd45a7ebcb0] > [pasha:03624] [19] > /home/pasha/deal/install_8.4.1/lib/libdeal_II.g.so.8.4.1(_ZNK6dealii16TrilinosWrappers10VectorBase5printERSojbb+0x27d) > > [0x7fd46293216f] > [pasha:03624] [20] > /home/pasha/deal/install_8.4.1/lib/libdeal_II.g.so.8.4.1(_ZNK6dealii16TrilinosWrappers3MPI11BlockVector5printERSojbb+0x126) > > [0x7fd462513f60] > [pasha:03624] [21] > ./FEM(_ZN3FEMILi2EE26assemble_large_deformationEb+0x2327) [0x560017] However, maybe this is not the case and some else has some other ideas. On Wednesday, January 18, 2017 at 10:37:48 AM UTC+1, Jean-Paul Pelteret wrote: > > Dear Pasha > > I think that there are two problems with how you've created your Eulerian > mapping: > > dealii::TrilinosWrappers::MPI::BlockVector >> rel_solution2(partition_relevant); >> > > Here you've provided the partitioning, but not an MPI communicator. You > probably only need to use the locally owned partitioning as well. > > >> MappingQEulerian<dim, dealii::TrilinosWrappers::MPI::BlockVector> >> q_mapping(1, dof_handler, rel_solution2); > > > I get the feeling that MappingQEulerian doesn't work with BlockVector's, > but I'm not 100% sure. When the post-processor for step-44 > <https://www.dealii.org/developer/doxygen/deal.II/step_44.html#Solidoutput_results> > > was developed, this certainly was the case. You could try to copy the block > vector into a normal MPI vector and pass that to the mapping. In the mean > time I've made a note <https://github.com/dealii/dealii/issues/3821> to > check whether BlockVector's work with this mapping. > > Regards, > Jean-Paul > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
