Enhancement to r.cross - add more layers to cross

[dnewcomb]

I'd like r.cross to be able to cross more than 10 layers. I did try changing the NFILES setting to 30 in raster/r.cross/glob.h. This allows me to cross 20 layers, but segfaults at 21 layers. Obviously, there is something else limiting the number of layers besides NFILES which my lack of skill at C programming does not permit me to see.

The rasters being crossed are in a region with 11451 rows and 25977 columns (297462627 cells).

The version of grass that I altered ( NFILES=10 --> 30)was the 6.4 svn weekly snapshot from 10/04/2008

[hamish]

You can follow the instructions on the wiki to run gdb to find out where the segfault comes from:

http://grass.osgeo.org/wiki/GRASS_Debugging

Hamish

[dnewcomb]

Sorry to take so long. I upgraded to 6.4 snapshot of 10/18/08 and segfault is gone. Tested 23 layers and worked fine with NFILES set to 30. Would 30 be a good number to set to? I can test to 30 or more.

Doug

[neteler]

I have increased it to 30 files (7: r34000, 6.4: 34001).

Markus