Mercurial > octave
changeset 27845:25f914321e7b
Use different Qhull options for griddataX to get better precision (bug #50494).
* griddata.m: Change Qhull options to {"Qt", "Qbb, "Qc}.
* griddatan.m: When input is 2-D or 3-D, Change Qhull options to
{"Qt", "Qbb, "Qc}.
| author | Rik <rik@octave.org> |
|---|---|
| date | Mon, 16 Dec 2019 14:42:06 -0800 |
| parents | e06ce39f78ad |
| children | 4075474fd9d0 |
| files | scripts/geometry/griddata.m scripts/geometry/griddatan.m |
| diffstat | 2 files changed, 7 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/scripts/geometry/griddata.m Mon Dec 16 14:34:34 2019 -0500 +++ b/scripts/geometry/griddata.m Mon Dec 16 14:42:06 2019 -0800 @@ -81,7 +81,8 @@ z = z(:); ## Triangulate data. - tri = delaunay (x, y); + ## Bug #50494 (loss of precision unless non-default Qhull options used) + tri = delaunay (x, y, {"Qt", "Qbb", "Qc"}); zi = NaN (size (xi)); if (strcmp (method, "cubic"))
--- a/scripts/geometry/griddatan.m Mon Dec 16 14:34:34 2019 -0500 +++ b/scripts/geometry/griddatan.m Mon Dec 16 14:42:06 2019 -0800 @@ -55,6 +55,11 @@ error ("griddatan: dimensional mismatch"); endif + ## Bug #50494 (loss of precision unless non-default Qhull options used) + if (isempty (varargin) && n <= 3) + varargin{1} = {"Qt", "Qbb", "Qc"}; + endif + ## triangulate data tri = delaunayn (x, varargin{:});