Thank you everyone. I understand your suggestions. From what I learned now in deal.II is, that each cell is owned by a specific processor which means I cannot access the information within a locally owned cell. That's why I have to loop over all cells and work with distributed data without using locally owned entities.
Maybe I describe my aim more properly. I like to find the maximum value within a MPI vector such as our locally_relevant_solution in step-40. Hence, I have to loop over all cells and all nodes to find where the geometrical position and value of the displacement lies. How is this possible within a parallely distributed environment? My approach using local_dof_indices is somehow impossible: for (; cell != endc; ++cell, ++cell_number) { pcout << "CELL ID: " << cell_number << std::endl; // this works surprisingly. cell->get_dof_indices(local_dof_indices); // Output or accessing this won't work I assume. for (unsigned int i = 0, j = 1; i < fe.dofs_per_cell; i += 2, j += 2) { double displacement_x = locally_relevant_solution(local_dof_indices[i]); double displacement_y = locally_relevant_solution(local_dof_indices[j]); double config_forces_norm = sqrt(config_forces_x * config_forces_x + config_forces_y * config_forces_y); } } -- 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.