CRA-861: Fixed issue with reading xyz-surfaces in the case when area …

[hgolsen]

…is taken from surface.

[alfbr]

Is this related to #254 or #255 ?

[hgolsen]

This is related to both. This is a general fix for reading xyz-surfaces, but tested on a simpler case. It is not certain that it fixes both these problems, so please re-run both cases when this pull-request is merged.

[alfbr]

Great, already downloaded and compiled. Will try to get it tested over here.

[alfbr]

PR tested. The inversion still fails, at the exact same location. This is the error message:

cravarun: ../libs/nrlib/surface/../grid/grid2d.hpp:147: size_t NRLib::Grid2D::GetIndex(size_t, size_t) const [with A = double]: Assertion `j < nj_' failed.

[hgolsen]

Sorry. This was not merged (this was only a request for merging), so you did not have the updates. (There is usually a small delay from when a request is made to when it is merged into the master branch.)

I have merged the commit into the master branch now. Can you test again tomorrow?

[alfbr]

I am not completely ignorant ;-) I forked your repository, and double checked that this commit was in there before I compiled:

$ git log
commit 403f6c5
Author: Håvard Goodwin Olsen [email protected]
Date: Tue Feb 28 13:06:24 2017 +0100

CRA-861: Fixed issue with reading xyz-surfaces in the case when area is taken from surface.

[hgolsen]

Ah, okey :) I didn't think of that option. Then these issues are still unresolved.

[alfbr]

Yes, I believe so. You have the datasets, so please see if you can reproduce at your end.