Re: [Freepooma-devel] out-of- field domain interpolation

[Richard Guenther]

Roman Krylov wrote:
> Another problem I encountered is: if I do gather or scatter operations 
> with positions to be gathered_to/scatter_from are out of mesh domain,
> SEGFAULT emerges.
> Is this documented somewhere?

It's probably cellContaining that segfaults, it should hit a PInsist if
you compile with --bounds checking enabled.  Maybe the cellContaining
code could be properly guarded agains out-of-bounds queries, but the
question remains what cell to return for such positions :)

So, I think you should avoid positions outside the mesh.

Richard.