Hmm, when doing that I only get the maximum value of the elements located on one processor, but not for all nodes.
Am Donnerstag, 19. April 2018 14:19:00 UTC+2 schrieb Bruno Turcksin: > > On Thursday, April 19, 2018 at 5:56:54 AM UTC-4, Maxi Miller wrote: >> >> Assumed I call that on a LinearAlgebraTrilinos::MPI::Vector, I get the >> error message >> The violated condition was: >> local_index >= 0 >> Additional information: >> You tried to access element 0 of a distributed vector, but this >> element is not stored on the current processor. Note: There are 120 >> elements stored on the current processor from within the range 150 through >> 269 but Trilinos vectors need not store contiguous ranges on each >> processor, and not every element in this range may in fact be stored >> locally. >> when running on multiple cores. How could I access the whole vector >> without running into ghost elements? >> > You can use > http://dealii.org/developer/doxygen/deal.II/classTrilinosWrappers_1_1MPI_1_1Vector.html#a315e9e01462309adce85314a4208c0ff > > to get the IndexSet of the locally owned elements. Then, you can either use > http://dealii.org/developer/doxygen/deal.II/classIndexSet.html#a0db80ed325ca916e10ff7f62ec71e4f5 > > to get the indices or directly loop over the indices of the > locally_owned_elements IndexSet with for (auto index : > locally_owned_elements) > > Best, > > Bruno > -- 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.
