[CIG-SHORT] Why Pylith sometimes cannot use linear interpolation?

[Matthew Knepley]

On Thu, Jan 17, 2019 at 1:32 PM tu xiang <tuxiang2016 at outlook.com> wrote:

> Dear Charles
> The old problem,When I use linear interpolation Pylith reports that it could not find final slip
> at some point using spatial database. You told me that It should haven’t specified enough
> points to interpolate the slip. I tried sparse mesh grid and dense spatial database, but it still
> reports errors. I read the examples provided by Pylith. I found that, the even the example
> provide sparse spatial database, it can use linear interpolation. So, I do not think "It should haven’t specified enough
> points to interpolate the slip. " is the reason.
> Do you think so?
>
> spatialdata can only do interpolation, not extrapolation (and we would not
want to do that since you do
not have good error bounds). That means if a point lies outside the convex
hull of your specified points,
then it will complain, even if you have "a lot" of points. So, make sure
you have points specified such that
the box they define encloses the entire simulation domain.

  Thanks,

    Matt


> Best regards
> Tu
>
> _________________________________________________
> CIG-SHORT has MOVED to our new forum.  Please
> head to https://community.geodynamics.org.
>
> This mailing list will be disabled at the end of
> January 2019.
> _________________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20190117/28579dd6/attachment.html>