Seyed, your results look like the ghost values are wrong. My first guess would be that you are still using a completely distributed vector (one initialized with locally_owned_dofs only) instead of a ghosted one (initialized with locally_owned_dofs and locally_relevant_dofs). Can you confirm that after solving you do something like
configurational_forces = local_configurational_forces; and that configurational_forces.has_ghost_values() == true local_configurational_forces.has_ghost_values() == false holds? Can you confirm that you use configurational_forces and not local_configurational_forces for output? What is the output of configurational_forces.print(std::cout)? What kind of error do you get when trying with pcout instead of std::cout? > Maybe this helps: > > How do you compute "reaction forces" within deal.II? > What exactly do you mean by that? We are much better if you can provide us with formulas rather than application-specific names. >From my understanding, you mean a projection, isn't it? So solving for u such that (u,v)=(f,v) for all v and given f? This is implemented in the VectorTools::project functions.[1] Best, Daniel [1] https://www.dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html -- 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.