On 12/12/18 8:32 AM, 'Maxi Miller' via deal.II User Group wrote: > I would like to use the function vector_value_list for the interpolation > between two different grids (one externally, one internally) instead of > having > to loop over all points and calling value_list for each point. Now I already > know that not all points are within the range of the current MPI-thread, thus > both functions will throw ExcPointNotAvailableHere sooner or later. If that > happens, I can deal with it when looping over every single point (by simply > hopping over it). If that happens when calling vector_value_list, how does > the > function behave then? Does it throw the error, but all points which can be > set > have been set, or does it throw the error before all possible points have > been > set? If the latter, how can I make sure that I got all possible points in > that > thread?
I had to look at the function's implementation and it turns out that it throws the exception the first time it encounters a cell that is not locally owned. It would be a reasonable idea to restructure the function in such a way that it first fills the data for all points for which this is possible and only then throws the exception. Would that be useful to you? Would you be interested in writing a patch to this end? Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- 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]. For more options, visit https://groups.google.com/d/optout.
