I believe operator= will do the right thing. In particular, when you call
close() after setting the values, any communication that needs to be done
will get done. Be sure to run with METHOD=dbg or METHOD=devel to trip any
asserts that you might hit when you're debugging.
Best,
Paul
On Wed, Jan
Hi Folks,
I have a nonlinear system and I want to set its current_local_solution
from a known vector, as follows:
NonlinearImplicitSystem system;
NumericVector v;
Set system.current_local_solution to be v;
How can I do this when the program is to be run in parallel? There are
ghost-node issues