Vachan, No there should be no need to call PETScWrappers::VectorBase::compress() after PETScWrappers::VectorBase::scale() (otherwise that would clearly be a bug in deal.II) since scale() calls VecPointwiseMult which is already collective. Does adding compress() in fact change anything for you? Otherwise, it would be good if you can share a minimal example showing the problem.
Best, Daniel Am So., 24. Nov. 2019 um 08:35 Uhr schrieb vachan potluri < [email protected]>: > Hello, > > I am facing a weird problem. At a point in code, I > have PETScWrappers::VectorBase::scale() called for few distributed vectors. > Subsequently, I have assignment operator on ghosted versions of these > vectors for parallel communication. When I launch the code with 2 or 4 > processes, it works fine. But with 3 processes, the code halts after the > scaling operations and before the first assignment. I am limited to 4 > processes. > > 1. Is a compress() required after scale()? With what operation as > argument? > 2. Why does this behaviour occur only when 3 processes are launched? > Has anyone experienced this before? > > Thanks > > -- > 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 [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/31558d55-94ca-4ab8-b668-9c231788a99d%40googlegroups.com > <https://groups.google.com/d/msgid/dealii/31558d55-94ca-4ab8-b668-9c231788a99d%40googlegroups.com?utm_medium=email&utm_source=footer> > . > -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAOYDWbKg13CH6MnjV3pqjnckf-O0dMaO4-UqnzGvRjodq8268Q%40mail.gmail.com.
