[GRASS-dev] [GRASS GIS] #2742: r.fill.dir crash with larger dataset

[GRASS GIS]

#2742: r.fill.dir crash with larger dataset
-------------------------+---------------------------------
 Reporter:  neteler      |      Owner:  grass-dev@…
     Type:  defect       |     Status:  new
 Priority:  normal       |  Milestone:  7.0.2
Component:  Raster       |    Version:  svn-releasebranch70
 Keywords:  r.fill.dir   |        CPU:  Unspecified
 Platform:  Unspecified  |
-------------------------+---------------------------------
 With a larger but not huge dataset r.fill.dir crashes:

 {{{
 GRASS 7.1.svn (eu_laea):> g.region -p n=2699750 s=2377750 w=4126750
 e=4536000
 projection: 99 (Lambert Azimuthal Equal Area)
 zone:       0
 datum:      etrs89
 ellipsoid:  grs80
 north:      2699750
 south:      2377750
 west:       4126750
 east:       4536000
 nsres:      25
 ewres:      25
 rows:       12880
 cols:       16370
 cells:      210845600

 GRASS 7.1.svn (eu_laea):>  gdb r.fill.dir
 GNU gdb (GDB) Fedora 7.9.1-17.fc22
 ...
 (gdb) r input=eu_dem_25m_northern_italy output=depressionless
 direction=direction
 Starting program: /home/neteler/software/grass71/dist.x86_64-unknown-
 linux-gnu/bin/r.fill.dir input=eu_dem_25m_northern_italy
 output=depressionless direction=direction
 Missing separate debuginfos, use: dnf debuginfo-install
 glibc-2.21-7.fc22.x86_64
 [Thread debugging using libthread_db enabled]
 Using host libthread_db library "/lib64/libthread_db.so.1".
 Reading input elevation raster map...
  100%
 Filling sinks...
 Determining flow directions for ambiguous cases...

 Program received signal SIGSEGV, Segmentation fault.
 0x0000000000405251 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1995420)
     at dopolys.c:15
 15          for (cnt = 0; cnt < sz; cnt += 3) {
 ...

 (gdb) bt
 #0  0x0000000000405251 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1995420)
     at dopolys.c:15
 #1  0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1995423)
     at dopolys.c:26
 #2  0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1995426)
     at dopolys.c:26
 #3  0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1995429)
     at dopolys.c:26
 ...

 #11759 0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1932495)
     at dopolys.c:26
 ...


 gdb > bt full
 ...
 #1189 0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1977663)
     at dopolys.c:26
         cnt = 1977666
         i = 6227
         j = 5510
         ii = <optimized out>
         jj = <optimized out>
 #1190 0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1977660)
     at dopolys.c:26
         cnt = 1977663
         i = 6227
         j = 5509
         ii = <optimized out>
         jj = <optimized out>
 ...


 #3143 0x00000000004052e1 in recurse_list (flag=flag at entry=53936,
 cells=cells at entry=0x7fffeee93010, sz=sz at entry=6687981,
 start=start at entry=1974657)
     at dopolys.c:26
         cnt = 1974654
         i = 6223
         j = 5435
         ii = <optimized out>
         jj = <optimized out>
 ...

 }}}


 Integer overflow?

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/2742>
GRASS GIS <https://grass.osgeo.org>