[GRASS-dev] G7: v.loutlier Decomposition failed

2013-12-06 Thread Blumentrath, Stefan
Hi

Maybe it is not of importance, because no one else would use such extreme 
settings, but in case someone considers it as a bug:

I was testing a bit with v.outlier in GRASS 7 on dense LiDAR data (~4 points 
per m2).

When I use a very low Tykhonov regularization parameter (lambda_i) of 
0.001 and rather short spline steps (1.5) at 0.5m resolution I 
get an error-message:
Decomposition failed at row 16499 and col 0
Stefan

___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

[GRASS-dev] segfault on 'r.stream.extract' - debian armh

2013-12-06 Thread epi
Hi,

i’m trying to run  “r.stream.extract” on a little linux machine, i got got 
grass up and running on a small quad-core Arm 1gb ram
OS :  Debian SID ArmHF.

the command i’m using is :

r.stream.extract elevation=elevation@PERMANENT accumulation=accum threshold=40 
stream_rast=stream_network stream_vect=streams --o --q

location :  nc_spm_08_grass7/PERMANENT/


i set the debug level to 5, this the segfault log :

https://gist.github.com/epifanio/7829206

if helpful, this is the log of make clean and make :

https://gist.github.com/epifanio/7829256

On other platform (same grass and r.stream.extract version it wks just fine)
Have you any idea on what’s wrong ?

Thanks,

Massimo.___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Re: [GRASS-dev] segfault on 'r.stream.extract' - debian armh

2013-12-06 Thread epi
i tried to add :
handle SIGILL nostop
to have more debug ..

i can’t understand the log :/


(gdb) handle SIGILL nostop
SignalStop  Print   Pass to program Description
SIGILLNoYes Yes Illegal instruction
(gdb) 




(gdb) exec-file r.stream.extract elevation=elevation@PERMANENT 
accumulation=accum threshold=20 stream_rast=stream_network stream_vect=streams 
--o
(gdb) r
Starting program: /home/epi/Envs/env1/grass-7.0.svn/bin/r.stream.extract 
.stream.extract elevation=elevation50m@PERMANENT accumulation=accum 
threshold=40 stream_rast=stream_network stream_vect=streams --o
[Thread debugging using libthread_db enabled]
Using host libthread_db library /lib/arm-linux-gnueabihf/libthread_db.so.1.

Program received signal SIGILL, Illegal instruction.
D2/5: G_file_name(): path = /home/epi/Envs/grass7data/nc_spm_08_grass7/PERMANENT

Description:
 Performs stream network extraction.

Keywords:
 raster, hydrology, stream network

Usage:
 r.stream.extract elevation=name [accumulation=name] [depression=name]
   threshold=value [d8cut=value] [mexp=value] [stream_length=value]
   [memory=value] [stream_rast=name] [stream_vect=name] [direction=name]
   [--overwrite] [--verbose] [--quiet]

Flags:
 --o   Allow output files to overwrite existing files
 --v   Verbose module output
 --q   Quiet module output

Parameters:
  elevation   Name of input elevation raster map
   accumulation   Name of input accumulation raster map
   Stream extraction will use provided accumulation instead of 
calculating it a new
 depression   Name of raster map with real depressions
   Streams will not be routed out of real depressions
  threshold   Minimum flow accumulation for streams
   Must be  0
  d8cut   Use SFD above this threshold
   If accumulation is larger than d8cut, SFD is used instead of 
MFD. Applies only if no accumulation map is given.
  default: infinity
   mexp   Montgomery exponent for slope, disabled with 0
   Montgomery: accumulation is multiplied with pow(slope,mexp) 
and then compared with threshold.
  default: 0
  stream_length   Delete stream segments shorter than stream_length cells.
   Applies only to first-order stream segments (springs/stream 
heads).
  default: 0
 memory   Maximum memory to be used in MB
  default: 300
stream_rast   Name for output raster map with unique stream ids
stream_vect   Name for output vector map with unique stream ids
  direction   Name for output raster map with flow direction

ERROR: Option elevation does not accept multiple answers
[Inferior 1 (process 7565) exited with code 01]
(gdb) 



On Dec 6, 2013, at 4:32 PM, epi massimodisa...@gmail.com wrote:

 Hi !
 
 this the output of gdb :
 
 Starting program: /home/epi/Envs/env1/grass-7.0.svn/bin/r.stream.extract 
 .stream.extract elevation=elevation50m@PERMANENT accumulation=accum 
 threshold=40 stream_rast=stream_network stream_vect=streams --o
 [Thread debugging using libthread_db enabled]
 Using host libthread_db library /lib/arm-linux-gnueabihf/libthread_db.so.1.
 
 Program received signal SIGILL, Illegal instruction.
 0x2c966e68 in ?? () from /usr/lib/arm-linux-gnueabihf/libcrypto.so.1.0.0
 (gdb) 
 
 
 Thanks for your help!
 
 
 Massimo.
 
 
 On Dec 6, 2013, at 4:16 PM, Rashad M mohammedrasha...@gmail.com wrote:
 
 Hello Massimo,
 
 A gdb output could be more helpful for a segfault
 
 
 On Fri, Dec 6, 2013 at 7:06 PM, epi massimodisa...@gmail.com wrote:
 Hi,
 
 i’m trying to run  “r.stream.extract” on a little linux machine, i got got 
 grass up and running on a small quad-core Arm 1gb ram
 OS :  Debian SID ArmHF.
 
 the command i’m using is :
 
 r.stream.extract elevation=elevation@PERMANENT accumulation=accum 
 threshold=40 stream_rast=stream_network stream_vect=streams --o --q
 
 location :  nc_spm_08_grass7/PERMANENT/
 
 
 i set the debug level to 5, this the segfault log :
 
 https://gist.github.com/epifanio/7829206
 
 if helpful, this is the log of make clean and make :
 
 https://gist.github.com/epifanio/7829256
 
 On other platform (same grass and r.stream.extract version it wks just fine)
 Have you any idea on what’s wrong ?
 
 Thanks,
 
 Massimo.
 
 ___
 grass-dev mailing list
 grass-dev@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-dev
 
 
 
 -- 
 Regards,
Rashad
 

___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Re: [GRASS-dev] segfault on 'r.stream.extract' - debian armh

2013-12-06 Thread epi
googling …

is it possible that in :

http://svn.osgeo.org/grass/grass/trunk/raster/r.stream.extract/

there may be some assembly code that gets executed which won't work under armhf 
 ?

Massimo.



On Dec 6, 2013, at 4:37 PM, epi massimodisa...@gmail.com wrote:

 i tried to add :
 handle SIGILL nostop
 to have more debug ..
 
 i can’t understand the log :/
 
 
 (gdb) handle SIGILL nostop
 SignalStopPrint   Pass to program Description
 SIGILLNo  Yes Yes Illegal instruction
 (gdb) 
 
 
 
 
 (gdb) exec-file r.stream.extract elevation=elevation@PERMANENT 
 accumulation=accum threshold=20 stream_rast=stream_network 
 stream_vect=streams --o
 (gdb) r
 Starting program: /home/epi/Envs/env1/grass-7.0.svn/bin/r.stream.extract 
 .stream.extract elevation=elevation50m@PERMANENT accumulation=accum 
 threshold=40 stream_rast=stream_network stream_vect=streams --o
 [Thread debugging using libthread_db enabled]
 Using host libthread_db library /lib/arm-linux-gnueabihf/libthread_db.so.1.
 
 Program received signal SIGILL, Illegal instruction.
 D2/5: G_file_name(): path = 
 /home/epi/Envs/grass7data/nc_spm_08_grass7/PERMANENT
 
 Description:
  Performs stream network extraction.
 
 Keywords:
  raster, hydrology, stream network
 
 Usage:
  r.stream.extract elevation=name [accumulation=name] [depression=name]
threshold=value [d8cut=value] [mexp=value] [stream_length=value]
[memory=value] [stream_rast=name] [stream_vect=name] [direction=name]
[--overwrite] [--verbose] [--quiet]
 
 Flags:
  --o   Allow output files to overwrite existing files
  --v   Verbose module output
  --q   Quiet module output
 
 Parameters:
   elevation   Name of input elevation raster map
accumulation   Name of input accumulation raster map
Stream extraction will use provided accumulation instead 
 of calculating it a new
  depression   Name of raster map with real depressions
Streams will not be routed out of real depressions
   threshold   Minimum flow accumulation for streams
Must be  0
   d8cut   Use SFD above this threshold
If accumulation is larger than d8cut, SFD is used instead 
 of MFD. Applies only if no accumulation map is given.
   default: infinity
mexp   Montgomery exponent for slope, disabled with 0
Montgomery: accumulation is multiplied with 
 pow(slope,mexp) and then compared with threshold.
   default: 0
   stream_length   Delete stream segments shorter than stream_length cells.
Applies only to first-order stream segments 
 (springs/stream heads).
   default: 0
  memory   Maximum memory to be used in MB
   default: 300
 stream_rast   Name for output raster map with unique stream ids
 stream_vect   Name for output vector map with unique stream ids
   direction   Name for output raster map with flow direction
 
 ERROR: Option elevation does not accept multiple answers
 [Inferior 1 (process 7565) exited with code 01]
 (gdb) 
 
 
 
 On Dec 6, 2013, at 4:32 PM, epi massimodisa...@gmail.com wrote:
 
 Hi !
 
 this the output of gdb :
 
 Starting program: /home/epi/Envs/env1/grass-7.0.svn/bin/r.stream.extract 
 .stream.extract elevation=elevation50m@PERMANENT accumulation=accum 
 threshold=40 stream_rast=stream_network stream_vect=streams --o
 [Thread debugging using libthread_db enabled]
 Using host libthread_db library /lib/arm-linux-gnueabihf/libthread_db.so.1.
 
 Program received signal SIGILL, Illegal instruction.
 0x2c966e68 in ?? () from /usr/lib/arm-linux-gnueabihf/libcrypto.so.1.0.0
 (gdb) 
 
 
 Thanks for your help!
 
 
 Massimo.
 
 
 On Dec 6, 2013, at 4:16 PM, Rashad M mohammedrasha...@gmail.com wrote:
 
 Hello Massimo,
 
 A gdb output could be more helpful for a segfault
 
 
 On Fri, Dec 6, 2013 at 7:06 PM, epi massimodisa...@gmail.com wrote:
 Hi,
 
 i’m trying to run  “r.stream.extract” on a little linux machine, i got got 
 grass up and running on a small quad-core Arm 1gb ram
 OS :  Debian SID ArmHF.
 
 the command i’m using is :
 
 r.stream.extract elevation=elevation@PERMANENT accumulation=accum 
 threshold=40 stream_rast=stream_network stream_vect=streams --o --q
 
 location :  nc_spm_08_grass7/PERMANENT/
 
 
 i set the debug level to 5, this the segfault log :
 
 https://gist.github.com/epifanio/7829206
 
 if helpful, this is the log of make clean and make :
 
 https://gist.github.com/epifanio/7829256
 
 On other platform (same grass and r.stream.extract version it wks just fine)
 Have you any idea on what’s wrong ?
 
 Thanks,
 
 Massimo.
 
 ___
 grass-dev mailing list
 grass-dev@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-dev
 
 
 
 -- 
 Regards,
Rashad
 
 

___
grass-dev mailing list
grass-dev@lists.osgeo.org

Re: [GRASS-dev] segfault on 'r.stream.extract' - debian armh

2013-12-06 Thread Maris Nartiss
Any warnings during compilation?
Try to compile with -g and -O0 to get slow but debugging-friendly version.

When program gets stopped by SIGILL in gdb, issue bt command to get backtrace.

No help form me, but more info never hurts ;)

Maris.


2013/12/6 epi massimodisa...@gmail.com:
 Hi !

 this the output of gdb :

 Starting program: /home/epi/Envs/env1/grass-7.0.svn/bin/r.stream.extract
 .stream.extract elevation=elevation50m@PERMANENT accumulation=accum
 threshold=40 stream_rast=stream_network stream_vect=streams --o
 [Thread debugging using libthread_db enabled]
 Using host libthread_db library
 /lib/arm-linux-gnueabihf/libthread_db.so.1.

 Program received signal SIGILL, Illegal instruction.
 0x2c966e68 in ?? () from /usr/lib/arm-linux-gnueabihf/libcrypto.so.1.0.0
 (gdb)


 Thanks for your help!


 Massimo.


 On Dec 6, 2013, at 4:16 PM, Rashad M mohammedrasha...@gmail.com wrote:

 Hello Massimo,

 A gdb output could be more helpful for a segfault


 On Fri, Dec 6, 2013 at 7:06 PM, epi massimodisa...@gmail.com wrote:

 Hi,

 i’m trying to run  “r.stream.extract” on a little linux machine, i got got
 grass up and running on a small quad-core Arm 1gb ram
 OS :  Debian SID ArmHF.

 the command i’m using is :

 r.stream.extract elevation=elevation@PERMANENT accumulation=accum
 threshold=40 stream_rast=stream_network stream_vect=streams --o --q

 location :  nc_spm_08_grass7/PERMANENT/


 i set the debug level to 5, this the segfault log :

 https://gist.github.com/epifanio/7829206

 if helpful, this is the log of make clean and make :

 https://gist.github.com/epifanio/7829256

 On other platform (same grass and r.stream.extract version it wks just
 fine)
 Have you any idea on what’s wrong ?

 Thanks,

 Massimo.

 ___
 grass-dev mailing list
 grass-dev@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-dev




 --
 Regards,
Rashad



 ___
 grass-dev mailing list
 grass-dev@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-dev
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Re: [GRASS-dev] r.univar, ERROR G_realloc:

2013-12-06 Thread Glynn Clements

Pietro wrote:

 I tried to use r,univar in grass7 (r58395 on linux 64bit with 24Gb of
 ram and 8Gb of swap), before to run the command I did a distclean, and
 I configured with:
 
 CFLAGS=-ggdb -Wall -Werror-implicit-function-declaration ./configure ...
 
 when I run:
 
 r.univar map=Combabula_Nearmap.red zones=seg005_64_zones
 percentile=90. output=red.csv -e --o
 
 The module termanate with:
 
 D1/1: G_find_raster(): name=MASK mapset=pietro
 Current region rows: 28545, cols: 27645
 ERROR: G_realloc: unable to allocate 68000 bytes of memory at
r.univar_main.c:327

Don't use r.univar -e, particularly for large maps. r.quantile and
r.statistics3 are far more efficient.

 Maybe is a stupid idea but since I had some problems also with
 v.build, do you think that could be possible that the problem is due
 not to the GRASS code, but to my physical memory that maybe is damaged
 in some sector?

That's unlikely. Is the OS recognising that you have 24 GiB of RAM? 
Are you allowed to use all (or most) of it for a single process (check
the output from ulimit -a)?

-- 
Glynn Clements gl...@gclements.plus.com
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev


Re: [GRASS-dev] segfault on 'r.stream.extract' - debian armh

2013-12-06 Thread Glynn Clements

epi wrote:

 googling �
 
 is it possible that in :
 
 http://svn.osgeo.org/grass/grass/trunk/raster/r.stream.extract/
 
 there may be some assembly code that gets executed which won't work under 
 armhf  ?

GRASS doesn't use assembly. And the SIGILL is reported as occurring in
libcrypto, not in the GRASS code.

The libcrypto dependency typically exists because GDAL links to libpq
(PostgreSQL client library) which uses libcrypto for certain
authentication methods.

libcrypto probably isn't even being used in this situation, so I
suspect that a bug is causing either a function pointer or a return
address to be corrupted, resulting in a jump to a random memory
location which just happens to be inside libcrypto.

As r.stream.extract is relatively new, it's possible that it hasn't
seen significant testing on platforms other than x86 and x86-64. Apart
from anything else, alignment bugs won't show up on those platforms
(x86 supports unaligned access, ARM doesn't AFAIK).

-- 
Glynn Clements gl...@gclements.plus.com
___
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev