Thanks Daniel, it worked.
On Wednesday, July 10, 2019 at 9:59:08 PM UTC-4, Daniel Arndt wrote:
>
> Vivek,
>
> just remove the line
> if (i->first >= range.first && i->first < range.second)
> VectorTools::interpolate_boundary_values should only set values for
> locally active degrees of freedom an
Vivek,
just remove the line
if (i->first >= range.first && i->first < range.second)
VectorTools::interpolate_boundary_values should only set values for locally
active degrees of freedom anyway.
As long as the values are consistent between different processes, this
should work just fine.
Best,
Dan
Hi all,
I had a legacy code were PETScWrappers::MPI::Vectors were used to for
parallel computing. For very large problems, I was told to move to
Trilinos. It was fairly easy to do so for most part of the code but I am
stuck in the implementation of boundary values. Currently the boundary
value
