Josh, The key issue for my application is whether the make_trilinos_map method of > the IndexSet class is called with overlapping equal to true or false. At > the moment, when make_trilinos_map is called from the reinit method for > Vector in the TrilinosWrappers∷MPI namespace, overlapping is always set > to true. > > > I believe that, at lease with my simple grid and running on only one > processor, that overlapping can be set to false. My question is, does the > following code change make any sense? In my case, has_ghosts is evaluating > to false and so when I make this code change, I have linear epetra vectors. > I'll need to think about how, or if, I could make this work with Ifpack and > Hypre when there are Ghost elements. > [...] >
Yes, that change makes sense. It would be great if you could turn it into a pull request. Are you willing to do so? Instructions can be found at https://github.com/dealii/dealii/wiki/Contributing. If you need help, just ask. Best, Daniel -- 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.
