But VectorTools::project needs a vector without ghosts, right ? And thats what the error is saying.
In step-40 solve function, we already use a vector without ghosts when we solve the linear system. Similar thing has to be done for project. best praveen > On 04-Oct-2021, at 10:06 AM, Wolfgang Bangerth <[email protected]> wrote: > > On 10/2/21 2:31 AM, Toddy Liu wrote: >> When running my code, the error occurred(the details are as follows) and I >> think the errors might result from the VectorTools::project. So what is the >> feasible way to set the values of my "old_locally_relevant_solution" using >> my "InitialValues" function. > > Toddy: It would be interesting to know whether the error is also in the > latest deal.II release, 9.3. But regardless of this, for parallel > computations it is often simpler to use interpolation (via > VectorTools::interpolate()) rather than projection because the former does > not require solving a linear system. > > Best > W. -- 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/2E9BA17D-B029-4FB6-8C6F-56D17AD77C87%40gmx.net.
