[GRASS-dev] [GRASS GIS] #1786: r.resamp.rst fails with message "split_and_interpolate() failed

[GRASS GIS]

#1786: r.resamp.rst fails with message "split_and_interpolate() failed
---------------------+------------------------------------------------------
 Reporter:  bko      |       Owner:  grass-dev@…              
     Type:  defect   |      Status:  new                      
 Priority:  normal   |   Milestone:  6.4.3                    
Component:  Default  |     Version:  6.4.1                    
 Keywords:           |    Platform:  Linux                    
      Cpu:  x86-64   |  
---------------------+------------------------------------------------------
 The following command fails on my computer:[[BR]]

 r.resamp.rst input=nearest_hydro_permanent_all_200 at PERMANENT ew_res=10
 ns_res=10 elev=nearest_hydro_permanent_all_200_10m overlap=3 zmult=1.0
 tension=40.

 Here is the output from the command:[[BR]]

 '''65910784 bytes of disk space for temp files.'''
 Temporarily changing the region to desired resolution ...
 Changing back to the original region ...
 Percent complete:
 Cannot fseek elev offset2=-2147383136
 Cannot fseek elev offset2=-2147382860
 Cannot fseek elev offset2=-2147382584
 Cannot fseek elev offset2=-2147382308
 Cannot fseek elev offset2=-2147382032
 Cannot fseek elev offset2=-2147381756
 Cannot fseek elev offset2=-2147381480
 Cannot fseek elev offset2=-2147381204
 Cannot fseek elev offset2=-2147380928
 Cannot fseek elev offset2=-2147380652
 Cannot fseek elev offset2=-2147380376
 Cannot fseek elev offset2=-2147380100
 Cannot fseek elev offset2=-2147379824
 Cannot fseek elev offset2=-2147379548
 Cannot fseek elev offset2=-2147379272
 Cannot fseek elev offset2=-2147378996
 Cannot fseek elev offset2=-2147378720
 Cannot fseek elev offset2=-2147378444
 Cannot fseek elev offset2=-2147378168
 Cannot fseek elev offset2=-2147377892
 Cannot fseek elev offset2=-2147377616
 Cannot fseek elev offset2=-2147377340
 Cannot fseek elev offset2=-2147377064
 Cannot fseek elev offset2=-2147376788
 Cannot fseek elev offset2=-2147376512
 Cannot fseek elev offset2=-2147376236
 Cannot fseek elev offset2=-2147375960
 Cannot fseek elev offset2=-2147375684
 Cannot fseek elev offset2=-2147375408
 Cannot fseek elev offset2=-2147375132
 Cannot fseek elev offset2=-2147374856
 Cannot fseek elev offset2=-2147374576
 interpolate() failed
 dnorm in mainc after grid before out1= 1249.036216
 split_and_interpolate() failed

 The command fails in exactly the same way on two different systems:
 1. Ubuntu 10.04 with grass 6.4 RC with amd64 processor and 64 bit OS
 2. Ubuntu 10.04 with grass 6.4.1 with amd64 processor and 64 bit OS

 This exact command does work on a project that is about 1/10th the size of
 the one on which the command failed.

 I looked at the code and wonder if I had an integer overflow (possibly on
 the variable sdisk). My region is huge and my resolution is high.

 Region:[[BR]]

 projection: 1 (UTM)
 zone:       12
 datum:      nad83
 ellipsoid:  grs80
 north:      6399780
 south:      6058660
 west:       243410
 east:       563010
 nsres:      10
 ewres:      10
 rows:       34112
 cols:       31960
 cells:      1090219520

 I find it strange that the '''65910784 bytes of disk space for temp
 files.''' didn't format to 65.91 MB ....

 Makes me think overflow but I my C knowledge is just about nonexistent.

 Cheers!

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/1786>
GRASS GIS <http://grass.osgeo.org>