[GRASS-dev] [GRASS GIS] #1262: r.proj failing in GRASS 7

#1262: r.proj failing in GRASS 7
--------------------------+-------------------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: major | Milestone: 7.0.0
Component: Raster | Version: unspecified
Keywords: reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
--------------------------+-------------------------------------------------
I just tried to reproject a raster from a UTM location into a latlon
location. After a silent failure in the GUI, I made a vector area in the
UTM location that matched the raster extents and reprojected it (v.proj)
into the latlon location. This worked. Then I set the region to match the
vector extents and tried to reproject the raster. I set the resolution to
match the source raster that I was trying to reproject. r.proj fails
immediately with the following seg fault error:

{{{
GRASS 7.0.svn (polengrass13):~ > r.proj input=h92_N-33-131-D
location=PolenUTM mapset=PERMANENT output=h92 resolution=0.000039682539683
Input Projection Parameters: +proj=tmerc +lat_0=0 +lon_0=19 +k=0.9993
+x_0=500000 +y_0=-5300000 +no_defs +a=6378137 +rf=298.257222101
+towgs84=0.000,0.000,0.000
Input Unit Factor: 1
Output Projection Parameters: +proj=longlat +a=6378137 +rf=298.257223563
+no_defs
Output Unit Factor: 1
Segmentation fault

}}}

I have not yet tried this in other dev versions, but I guess I'd better do
so.

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj failing in GRASS 7
--------------------------+-------------------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: major | Milestone: 7.0.0
Component: Raster | Version: unspecified
Keywords: reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
--------------------------+-------------------------------------------------

Comment(by cmbarton):

Same thing happens in GRASS 6.4.1. It also doesn't matter if I set the
resolution or not.

source location projection information (not a UTM, but a mercator):

zone: 0
datum: etrs89
ellipsoid: grs80
north: 516344.80163701
south: 497432.30163701
west: 380800.49568881
east: 398208.49568881
nsres: 4.25
ewres: 4.25
rows: 4450
cols: 4096
cells: 1822720

target location projection information:

projection: 3 (Latitude-Longitude)
zone: 0
datum: ** unknown (default: WGS84) **
ellipsoid: wgs84
north: 52:30:12N
south: 52:19:48N
west: 17:14:36E
east: 17:30:21E
nsres: 0:00:03
ewres: 0:00:03
rows: 208
cols: 315
cells: 65520

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:1&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
--------------------------+-------------------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
--------------------------+-------------------------------------------------
Changes (by cmbarton):

  * priority: major => blocker
  * milestone: 7.0.0 => 6.4.1

Comment:

I've tested this and on the Mac at least r.proj is broken in all
development versions: 7.0, 6.5, and 6.4.1. Here is the result from testing
it on the North Carolina demo data set, trying to reproject from State
Plane into Latlon (same thing happens going the other way).

There is a segmentation fault and no map is reprojected.

{{{
r.proj input=aspect location=nc_spm_07 mapset=PERMANENT
Input Projection Parameters: +proj=lcc +lat_1=36.16666666666666
+lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
+no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
Input Unit Factor: 1
Output Projection Parameters: +proj=longlat +a=6378137 +rf=298.257223563
+no_defs
Output Unit Factor: 1
Segmentation fault
}}}

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:2&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
--------------------------+-------------------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
--------------------------+-------------------------------------------------

Comment(by mlennert):

Replying to [comment:2 cmbarton]:
> I've tested this and on the Mac at least r.proj is broken in all
development versions: 7.0, 6.5, and 6.4.1. Here is the result from testing
it on the North Carolina demo data set, trying to reproject from State
Plane into Latlon (same thing happens going the other way).
>
> There is a segmentation fault and no map is reprojected.
>
>
> {{{
> r.proj input=aspect location=nc_spm_07 mapset=PERMANENT
> Input Projection Parameters: +proj=lcc +lat_1=36.16666666666666
+lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
+no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
> Input Unit Factor: 1
> Output Projection Parameters: +proj=longlat +a=6378137
+rf=298.257223563 +no_defs
> Output Unit Factor: 1
> Segmentation fault
> }}}
>

I've just tried to reproduce on Debian GNU/Linux using fairly recent
grass7 and grass65, but I don't get any segmentation fault and
reprojection works well.

Moritz

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:3&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------
Changes (by martinl):

  * keywords: reprojection => r.proj, reprojection

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:4&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by neteler):

Replying to [comment:2 cmbarton]:
> There is a segmentation fault and no map is reprojected.
>
> {{{
> r.proj input=aspect location=nc_spm_07 mapset=PERMANENT
> Input Projection Parameters: +proj=lcc +lat_1=36.16666666666666
+lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
+no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
> Input Unit Factor: 1
> Output Projection Parameters: +proj=longlat +a=6378137
+rf=298.257223563 +no_defs
> Output Unit Factor: 1
> Segmentation fault
> }}}

Michael, please use a debugger and add the backtrace to the ticket.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:5&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by martinl):

Replying to [comment:5 neteler]:

> Michael, please use a debugger and add the backtrace to the ticket.

I can confirm that it's working on Linux. Some debug info is needed.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:6&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

Here is the gdb results for r.proj

{{{
GRASS 6.4.1svn (nc_ll):~ > gdb `which r.proj`
GNU gdb 6.3.50-20050815 (Apple version gdb-1510) (Wed Sep 22 02:45:02 UTC
2010)
Copyright 2004 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you
are
welcome to change it and/or distribute copies of it under certain
conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB. Type "show warranty" for
details.
This GDB was configured as "x86_64-apple-darwin"...Reading symbols for
shared libraries ......
warning: Could not find object file
"/Users/Shared/src/gdal/gdal-1.7.1/frmts/o/.libs/wktrasterdataset.o" - no
debug information available for "wktrasterdataset.cpp".

warning: Could not find object file
"/Users/Shared/src/gdal/gdal-1.7.1/frmts/o/.libs/wktrasterrasterband.o" -
no debug information available for "wktrasterrasterband.cpp".

warning: Could not find object file
"/Users/Shared/src/gdal/gdal-1.7.1/frmts/o/.libs/wktrasterwrapper.o" - no
debug information available for "wktrasterwrapper.cpp".

.. done

(gdb) run input=aspect location=nc_spm_07 mapset=PERMANENT
Starting program:
/Users/Shared/grass_dev/grass64_rb/macosx/dist/GRASS-6.4.app/Contents/MacOS/bin/r.proj
input=aspect location=nc_spm_07 mapset=PERMANENT
Reading symbols for shared libraries
.+++++++.....................................................................................................
done
Input Projection Parameters: +proj=lcc +lat_1=36.16666666666666
+lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
+no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
Input Unit Factor: 1
Output Projection Parameters: +proj=longlat +a=6378137 +rf=298.257223563
+no_defs
Output Unit Factor: 1

Program received signal EXC_BAD_ACCESS, Could not access memory.
Reason: KERN_INVALID_ADDRESS at address: 0x0000000000000000
0x00007fff80f94f2f in __vfprintf ()
(gdb) bt full
#0 0x00007fff80f94f2f in __vfprintf ()
No symbol table info available.
#1 0x00007fff80fc9b17 in vsnprintf ()
No symbol table info available.
#2 0x00007fff8100a4da in __vsprintf_chk ()
No symbol table info available.
#3 0x0000000100020f7d in print_error [inlined] () at
/Users/Shared/grass_dev/grass64_rb/lib/gis/error.c:61
         buffer = "\020\035?\001\001\000\000\000@?
p?\000\000?¿_?\000\000&\022???\000\000\000\000\000\000l??\000\000\000\000\000\000\000?¿_?\000\000@u???\000\000\000Ŀ_?\000\000\000Ŀ_?\000\000\000Ŀ_?\000\000\000Ŀ_?\000\000?Ŀ_?\000\000@?
p?\000\000?Ŀ_?\000\000?U\037p?\000\000?Ŀ_?\000\000h?\003\000\001\000\000\000?$?HZ\f???=\032??ā???\000\000\000\000\000\000\000\000??\bN???,??0j\037p?\000\000??\036O??P?"...
#4 0x0000000100020f7d in vfprint_error (type=0, template=<value
temporarily unavailable, due to optimizations>, ap=<value temporarily
unavailable, due to optimizations>) at error.c:64
         format = 1606392479
         type = Cannot access memory at address 0x0

(gdb) l
67 {p_nearest, "nearest", "nearest neighbor"},
68 {p_bilinear, "bilinear", "bilinear"},
69 {p_cubic, "cubic", "cubic convolution"},
70 {NULL, NULL, NULL}
71 };
72
73 static char *make_ipol_list(void);
74
75 int main(int argc, char **argv)
76 {
(gdb) frame 2
#2 0x00007fff8100a4da in __vsprintf_chk ()
(gdb) l
77 char *mapname, /* ptr to name of output layer */
78 *setname, /* ptr to name of input mapset */
79 *ipolname; /* name of interpolation method */
80
81 int fdi, /* input map file descriptor */
82 fdo, /* output map file descriptor */
83 method, /* position of method in table */
84 permissions, /* mapset permissions */
85 cell_type, /* output celltype */
86 cell_size, /* size of a cell in bytes */
(gdb)

}}}

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:7&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

Note that I am running this on the same system where I compiled it with no
errors, using William Kyngesburye's Mac frameworks.

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:8&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by martinl):

Replying to [comment:7 cmbarton]:
{{{
> Program received signal EXC_BAD_ACCESS, Could not access memory.
> Reason: KERN_INVALID_ADDRESS at address: 0x0000000000000000
> 0x00007fff80f94f2f in __vfprintf ()
> (gdb) bt full
> #0 0x00007fff80f94f2f in __vfprintf ()
> No symbol table info available.
> #1 0x00007fff80fc9b17 in vsnprintf ()
> No symbol table info available.
> #2 0x00007fff8100a4da in __vsprintf_chk ()
> No symbol table info available.
> #3 0x0000000100020f7d in print_error [inlined] () at
/Users/Shared/grass_dev/grass64_rb/lib/gis/error.c:61
> buffer = "\020\035?\001\001\000\000\000@?
p?\000\000?¿_?\000\000&\022???\000\000\000\000\000\000l??\000\000\000\000\000\000\000?¿_?\000\000@u???\000\000\000Ŀ_?\000\000\000Ŀ_?\000\000\000Ŀ_?\000\000\000Ŀ_?\000\000?Ŀ_?\000\000@?
p?\000\000?Ŀ_?\000\000?U\037p?\000\000?Ŀ_?\000\000h?\003\000\001\000\000\000?$?HZ\f???=\032??ā???\000\000\000\000\000\000\000\000??\bN???,??0j\037p?\000\000??\036O??P?"...
> #4 0x0000000100020f7d in vfprint_error (type=0, template=<value
temporarily unavailable, due to optimizations>, ap=<value temporarily
unavailable, due to optimizations>) at error.c:64
> format = 1606392479
> type = Cannot access memory at address 0x0
}}}

Any idea? This bug is marked as a blocker for 6.4.1...

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:9&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

All I know is that this is relatively new, affects all dev versions, and
shows up on other computers beyond mine. It may not affect Linux versions.
I don't know about Windows. r.proj is a key, core module. I don't think we
should release a new version of GRASS where this is non-functional. I'm
afraid that I don't know what the problem is. I'd hoped that the debug
code would help someone with with knowledge of C. I'm happy to test
potential fixes.

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:10&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by marisn):

If backtraces involve suspicious memory addresses (0x0 etc.), it's
necessary to run program under valrgind to look for possible memory
corruption or pointer abuse.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:11&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by mmetz):

Replying to [comment:11 marisn]:
> If backtraces involve suspicious memory addresses (0x0 etc.), it's
necessary to run program under valrgind to look for possible memory
corruption or pointer abuse.

Good hint. The segfault happened when printing a message, also on Windows,
not on Linux.

Should be fixed in r45195, r45196, r45197 (6.4, 6.5, 7.0).

Markus M

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:12&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

Replying to [comment:12 mmetz]:
> Replying to [comment:11 marisn]:
> > If backtraces involve suspicious memory addresses (0x0 etc.), it's
necessary to run program under valrgind to look for possible memory
corruption or pointer abuse.
>
> Good hint. The segfault happened when printing a message, also on
Windows, not on Linux.
>
> Should be fixed in r45195, r45196, r45197 (6.4, 6.5, 7.0).
>
> Markus M

Thanks much. I'll test in the next day or so.

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:13&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

I just did an updated from the svn trunk (after a make distclean). This
has changed the results and error, but has not corrected the problem. Just
to be sure that there is nothing weird with the map, here is what I did. I
opened the NC_spm_07 dataset and set the region to match the
elev_state_500m dem. Then I tried to reproject into the state plane region
the 1 arc second dem from the nc_ll location. I set debug=5 and here is
the result.

{{{
GRASS 7.0.svn (nc_spm_07):~ > g.gisenv set="DEBUG=5"
GRASS 7.0.svn (nc_spm_07):~ > r.proj input=elev_ned_1arcsec location=nc_ll
mapset=PERMANENT output=elev_ned_1arcsec
D1/5: G_find_raster(): name=elev_ned_1arcsec mapset=user1
D2/5: G__read_Cell_head
D2/5: G__read_Cell_head_array
D3/5: region item: proj: 99
D3/5: region item: zone: 0
D3/5: region item: north: 318500
D3/5: region item: south: -16000
D3/5: region item: east: 963000
D3/5: region item: west: 124000
D3/5: region item: cols: 1678
D3/5: region item: rows: 669
D3/5: region item: e-w resol: 500
D3/5: region item: n-s resol: 500
D3/5: region item: top: 1
D3/5: region item: bottom: 0
D3/5: region item: cols3: 1678
D3/5: region item: rows3: 669
D3/5: region item: depths: 1
D3/5: region item: e-w resol3: 500
D3/5: region item: n-s resol3: 500
D3/5: region item: t-b resol: 1
D1/5: G_find_raster(): name=elev_ned_1arcsec mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_ned_1arcsec mapset=PERMANENT
D1/5: G_find_raster(): name=elev_ned_1arcsec mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_ned_1arcsec mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_ned_1arcsec mapset=PERMANENT
D5/5: color buf = [0% yellow]
D5/5: color buf = [20% green]
D5/5: color buf = [40% cyan]
D5/5: color buf = [60% blue]
D5/5: color buf = [80% magenta]
D5/5: color buf = [100% red]
Input Projection Parameters: +proj=longlat +a=6378137 +rf=298.257223563
+no_defs
Input Unit Factor: 1
Output Projection Parameters: +proj=lcc +lat_1=36.16666666666666
+lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
+no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
Output Unit Factor: 1
D2/5: G__read_Cell_head
D2/5: G__read_Cell_head_array
D3/5: region item: proj: 3
D3/5: region item: zone: 0
D3/5: region item: north: 35:54:41.000547N
D3/5: region item: south: 35:35:17.001098N
D3/5: region item: east: 78:27:08.000616W
D3/5: region item: west: 78:49:17.999987W
D3/5: region item: cols: 1330
D3/5: region item: rows: 1164
D3/5: region item: e-w resol: 0:00:01
D3/5: region item: n-s resol: 0:00:01
D3/5: region item: format: -1
D3/5: region item: compressed: 1
D3/5: G_adjust_Cell_head: epsilon_ns: 8.59107e-07, epsilon_ew: 1e-06
D3/5: G_adjust_Cell_head: epsilon_ns: 8.59107e-07, epsilon_ew: 1e-06
D1/5: G_find_raster(): name=MASK mapset=user1
D3/5: Top: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
35.588055
D3/5: Right: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
35.588055
D3/5: Bottom: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
35.588055
D3/5: Left: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
35.588055
D3/5: Extra check: xmin: -78.821528; xmax: -78.452361; ymin: 35.588195;
ymax: 35.911250
D3/5: Final check: xmin: -78.821667; xmax: -78.452222; ymin: 35.588056;
ymax: 35.911389
D3/5: G_adjust_Cell_head: epsilon_ns: 8.59107e-07, epsilon_ew: 1e-06
D1/5: G_find_raster(): name=MASK mapset=user1
D3/5: Top: xmin: 625709.948811; xmax: 659030.331786; ymin: 239768.184363;
ymax: 239890.045778
D3/5: Right: xmin: 625709.948811; xmax: 659228.054367; ymin:
204049.351740; ymax: 239890.045778
D3/5: Bottom: xmin: 625709.948811; xmax: 659228.054367; ymin:
203927.002865; ymax: 239890.045778
D3/5: Left: xmin: 625709.948811; xmax: 659228.054367; ymin: 203927.002865;
ymax: 239890.045778
D3/5: Extra check: xmin: 625709.948811; xmax: 659228.054367; ymin:
203927.002865; ymax: 239890.045778
D3/5: Final check: xmin: 625500.000000; xmax: 659500.000000; ymin:
203500.000000; ymax: 240000.000000
Segmentation fault
}}}

I can run a shell debugger on this again if needed.

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:14&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by marisn):

Replying to [comment:14 cmbarton]:
>
> I can run a shell debugger on this again if needed.
>
tl;dr
valgrind!

Works just fine on my AMD64 Linux.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:15&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by neteler):

For valgrind usage, see
http://grass.osgeo.org/wiki/GRASS_Debugging#Using_Valgrind

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:16&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

Replying to [comment:15 marisn]:
> Replying to [comment:14 cmbarton]:
> >
> > I can run a shell debugger on this again if needed.
> >
> tl;dr
> valgrind!

Are these shell commands? The don't work as typed in my terminal.

>
> Works just fine on my AMD64 Linux.

We know that it works on Linux. It is broken for Mac and Windows.

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:17&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by cmbarton):

Looks like valgrind does not come on my Mac by default. I'll have to find
instructions and time to compile this. If anyone using a Mac already has
it, it might be faster if they tried to run this. I can post a GRASS
binary to try.

Michael

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:18&gt;
GRASS GIS <http://grass.osgeo.org>

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
Reporter: cmbarton | Owner: grass-dev@…
     Type: defect | Status: new
Priority: blocker | Milestone: 6.4.1
Component: Raster | Version: unspecified
Keywords: r.proj, reprojection | Platform: MacOSX
      Cpu: OSX/Intel |
----------------------------------+-----------------------------------------

Comment(by neteler):

Here are instructions to compile it which may be helpful:
http://stackoverflow.com/questions/2285964/how-can-i-compile-valgrind-on-
snow-leopard

--
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:19&gt;
GRASS GIS <http://grass.osgeo.org>