Open landam opened 4 months ago
Describe the bug r.terraflow fails with "grid larger than memory" error on very small regions
r.terraflow
north: 21 south: -1.5 west: -1.5 east: 5 nsres: 0.5 ewres: 0.5 rows: 45 cols: 13 cells: 585
To Reproduce Steps to reproduce the behavior:
grass -c EPSG:5514 /tmp/r_terraflow
r.unpack input=dem.pack out=dem
g.region rast=dem
r.terraflow elevation=dem@PERMANENT filled=dem_filled direction=dem_flowdir accumulation=dem_flowacc --o -s
WARNING: Raster map <dem> is of type DCELL_TYPE -- precision may be lost Reading input data... 100% WARNING: grid larger than memory (ignored) r.terraflow: grid.cpp:33: grid::grid(dimension_type, dimension_type, dimension_type, dimension_type, long int, cclabel_type): Assertion `(size_t)width * height * sizeof(gridElement) < getAvailableMemory()' failed. Aborted
Expected behavior r.terraflow should run successfully.
System description (please complete the following information):
It's due ps->size which is very large (140720409312672), see https://github.com/OSGeo/grass/blob/main/raster/r.terraflow/fill.cpp#L454
ps->size
Describe the bug
r.terraflow
fails with "grid larger than memory" error on very small regionsTo Reproduce Steps to reproduce the behavior:
grass -c EPSG:5514 /tmp/r_terraflow
r.unpack input=dem.pack out=dem
g.region rast=dem
r.terraflow elevation=dem@PERMANENT filled=dem_filled direction=dem_flowdir accumulation=dem_flowacc --o -s
Expected behavior
r.terraflow
should run successfully.System description (please complete the following information):