comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Memory Leakage......(I think)
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Memory Leakage......(I think) [message #30479 is a reply to message #30477] Wed, 01 May 2002 05:03 Go to previous messageGo to previous message
a.fielding is currently offline  a.fielding
Messages: 4
Registered: May 2002
Junior Member
Thanks for the comments regarding my problem. Someone suggested I post
the code to clarify the problem so here it is.
This is a loop within a sub-procedure that is called (256*256) times
from the main procedure. The count variable controlling the loop
varies from a minimum of 0 up to 510 depending on the call. I am
fairly certain the problem is with the arithmetic within the loop as
removing it and replacing it with a simple statement results in the
main program processing in around 2 minutes! Any hints/suggestions
gratefully received.

Andrew

; Weighted Interpolation voxel density calculation

fval=0.0
sum=0.0

for i=0L,count-1 do begin
index=ishft(X1,-14) + row[ishft(Y1,-14)] + slice[ishft(Z1,-14)]
val=ctcube[index]
if (val gt config.air) then begin
d0= X1/SCALE - fix(X1/SCALE)
d1= Y1/SCALE - fix(Y1/SCALE)
d2= Z1/SCALE - fix(Z1/SCALE)
fval=(1.0-d0)*(1.0-d1)*(1.0-d2)*float(val)
fval=temporary(fval) + d0*(1.0-d1)*(1.0-d2)*float(ctcube[index+1])
fval=temporary(fval) +
(1.0-d0)*d1*(1.0-d2)*float(ctcube[index+xdim])
fval=temporary(fval) + d0*d1*(1.0-d2)*float(ctcube[index+1+xdim])
index=temporary(index) + xdim*ydim
fval=temporary(fval) + (1.0-d0)*(1.0-d1)*d2*float(ctcube[index])
fval=temporary(fval) + d0*(1.0-d1)*d2*float(ctcube[index+1])
fval=temporary(fval) + (1.0-d0)*d1*d2*float(ctcube[index+xdim])
fval=temporary(fval) + d0*d1*d2*float(ctcube[index+1+xdim])

val= fix(fval)
sum=temporary(sum) + convert_ct2mu[val-(Hmin-1024)]
endif

X1=X1+dX
Y1=Y1+dY
Z1=Z1+dZ
endfor
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Direct->Object Graphics Help!
Next Topic: color_quan(...., Cube=6) makes white white, but ...

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sun Oct 12 04:43:28 PDT 2025

Total time taken to generate the page: 0.24232 seconds