[GRASS-user] r.terraflow - WORKS

[Francesco Mirabella]

Hi Markus,
thanks for you reply, it seems that the patch works.... below I attach 
the command result
Now I'll try to play with the maps .....
many thanks
Francesco

============================================
GRASS 6.4.0RC5 (WGS84_UTM33):~ > r.terraflow elevation=dem_tagliato 
filled=flood direction=flow swatershed=sink accumulation=accumulation 
tci=tci d8cut=infinity memory=1600 STREAM_DIR=/tmp stats=stats.out
STREAM temporary files in /tmp  (THESE INTERMEDIATE STREAMS WILL NOT BE 
DELETED IN CASE OF ABNORMAL TERMINATION OF THE PROGRAM. TO SAVE SPACE 
PLEASE DELETE THESE FILES MANUALLY!)
MFD flow direction
D8CUT=999999986991104.000000
Memory size: 1.56G (1677721600) bytes
Memory manager registering memory in MM_IGNORE_MEMORY_EXCEEDED mode.
total elements=67071935, nodata elements=61864898
largest temporary files:
FILL: 3.50G (3756028360) [67071935 elements, 56B each]
FLOW: 397.27M (416562960) [5207037 elements, 80B each]
Will need at least 7.00G (7512056720) space available in /tmp
------------------------------
COMPUTING FLOW DIRECTIONS
classifying nodata (inner & boundary)
EMPQUEUEADAPTIVE: starting in-memory pqueue
EMPQUEUEADAPTIVE: available memory: 1597.93MB
EMPQUEUEADAPTIVE: desired memory: 1597.93MB
sz_stream: 270388 buf_arity: 200 mm_overhead: 8665728 mm_avail: 1675549602.
EMPQUEUEADAPTIVE: memory overhead set to 8.26428MB
EMPQUEUEADAPTIVE: pqsize set to 208360484
assigning preliminary directions
finding flat areas (plateaus and depressions)
EMPQUEUEADAPTIVE: starting in-memory pqueue
EMPQUEUEADAPTIVE: available memory: 1597.41MB
EMPQUEUEADAPTIVE: desired memory: 1597.41MB
sz_stream: 270388 buf_arity: 200 mm_overhead: 8665728 mm_avail: 1675008754.
EMPQUEUEADAPTIVE: memory overhead set to 8.26428MB
EMPQUEUEADAPTIVE: pqsize set to 208292878
assigning directions on plateaus
generating watersheds and watershed graph
EMPQUEUEADAPTIVE: starting in-memory pqueue
EMPQUEUEADAPTIVE: available memory: 1598.96MB
EMPQUEUEADAPTIVE: desired memory: 1598.96MB
sz_stream: 270396 buf_arity: 200 mm_overhead: 8678784 mm_avail: 1676631202.
EMPQUEUEADAPTIVE: memory overhead set to 8.27673MB
EMPQUEUEADAPTIVE: pqsize set to 104247026
flooding depressions
available memory: 1599MB (1676901678B)
UnionFind::makeSet: reallocate double 2000
warning: watershed 1 (R=1) not done
------------------------------
REASSIGNING DIRECTIONS
finding flat areas (plateaus and depressions)
EMPQUEUEADAPTIVE: starting in-memory pqueue
EMPQUEUEADAPTIVE: available memory: 1597.15MB
EMPQUEUEADAPTIVE: desired memory: 1597.15MB
sz_stream: 270388 buf_arity: 200 mm_overhead: 8665728 mm_avail: 1674738378.
EMPQUEUEADAPTIVE: memory overhead set to 8.26428MB
EMPQUEUEADAPTIVE: pqsize set to 208259081
EMPQUEUEADAPTIVE: starting in-memory pqueue
EMPQUEUEADAPTIVE: available memory: 1596.38MB
EMPQUEUEADAPTIVE: desired memory: 1596.38MB
sz_stream: 270388 buf_arity: 200 mm_overhead: 8665728 mm_avail: 1673927126.
EMPQUEUEADAPTIVE: memory overhead set to 8.26428MB
EMPQUEUEADAPTIVE: pqsize set to 208157674
assigning directions on plateaus
creating flowStream: [AMI_STREAM /tmp/flowStream 0]
compute flow directions done.
writing stream /tmp/STREAM_fE0Ows to raster map flow:  100%
writing stream /tmp/STREAM_G82bHw to raster map flood:  100%
writing stream /tmp/STREAM_tSF3ML to raster map sink:  100%
------------------------------
COMPUTING FLOW ACCUMULATION
creating sweep stream from fill output stream
sorting sweep stream
sweeping: EMPQUEUEADAPTIVE: starting in-memory pqueue
EMPQUEUEADAPTIVE: available memory: 1599.73MB
EMPQUEUEADAPTIVE: desired memory: 1599.73MB
sz_stream: 270396 buf_arity: 200 mm_overhead: 8678784 mm_avail: 1677442317.
EMPQUEUEADAPTIVE: memory overhead set to 8.27673MB
EMPQUEUEADAPTIVE: pqsize set to 104297720
  100%
sorting sweep output stream
writing stream /tmp/STREAM_a2UrtD to raster maps accumulation, tci:  100%
r.terraflow complete.
GRASS 6.4.0RC5 (WGS84_UTM33):~ >
============================================




Markus Neteler wrote:
> Hi Francesco,
> 
> On Wed, Nov 11, 2009 at 12:53 PM, Francesco Mirabella <mirabell at unipg.it> wrote:
>> Hi Markus,
>> I am new to code patching, yes there is space in the /tmp directory
>> (about 14Gb) :
> ..
> 
>> I have now re-configured with both 64bit and large file support.
> 
> well, did you are recompile?
> 
>> I get the same error however I need e clarification:
>>
>> There exist two ami_sort_impl.h files, one in
>> /../../../grass640_rc5/dist.i686-pc-linux-gnu/include/grass/iostream/
> 
> ... that's not relevant.
> 
>> and the other in
>> /../../../grass640_rc5/include/iostream/
> 
> ... that's relevant.
> 
>> I have tried to apply the patch to the first, then to second and to both.
>> However, I get very similar errors to the original error when using the
>> patch for the first and for the second and for both.
> 
> Unless you don't recompile it won't work.
> The "first" instance above is a copy created during compilation
> of the "second" instance which is the real code.
> 
> 1. Patch
>  grass640_rc5/include/iostream/
> 
> 2. configure
> 
> 3. compile
> 
> 4. install or run directly from GRASS from
>    grass640_rc5/bin.i686-pc-linux-gnu/grass64
> 
> All clear like mud? :)
> 
> Markus


-- 
**********************************************
Francesco Mirabella,
Geologia Strutturale e Geofisica
Universita' di Perugia,
Dipartimento di Scienze della Terra,
Piazza Universita' 1, 06100 Perugia (Italy)
tel: ++39.(0)75.584.7948
fax: ++39.(0)75.585.2603
skype: francesco.mirabella
web: http://www.unipg.it/~mirabell/
**********************************************