#4 Program stops during pass 1 of scale domain 1

[Jimm Domingo]

Several people have reported that the program just silently quits while interpolating the raster surface on the first pass of scale domain 1. Example console output provided by J. Muckenhoupt:

C:\Data\MCC-test\mcc-lidar -s 1.25 -t 0.2 45121E8206.las 45121E8206_mcc.las
Reading 19324586 points from '45121E8206.las' ...
  min X Y : 7.83419e+006 708065
  max X Y : 7.83747e+006 712678
Scale domain: 1 (cell resolution: 0.625)
Pass: 1
  Interpolating 19324586 points...

C:\Data\MCC-test>

[Jimm Domingo]

• owner set to jimm_domingo
• status changed from new to accepted

Determined that the error is an invalid index to a Boost.MultiArray:

/opt/local/include/boost/multi_array/base.hpp:178: failed assertion `size_type(idx - index_bases[0]) < extents[0]'

The invalid index is due to the grid of interpolation regions having an X-Y extent that doesn't match the extent of raster surface. Suspected the cause was cumulative rounding error in calculating X and Y coordinates for intervals and cell centers. This was confirmed by modifying the unit tests for the grid base class to use values from one of the problematic input files (revision r80).

[Jimm Domingo]

• status changed from accepted to closed
• resolution set to fixed

Fixed in revision r82.

[Jimm Domingo]

• milestone set to 1.0 (b2)