comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » Convolution

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: convolution [message #60415 is a reply to message #26609]

Mon, 12 May 2008 19:31

David Fanning
Messages: 11724
Registered: August 2001

Senior Member

sarah writes:

> x=make_array(1024)
> sigma=15
> mu =15
> const=1/(sigma*sqrt(2*!pi))
> for i = 0,1024 do x[i]= array[0,*]
> f= const*( EXP(-1.0*(x - mu)^2/(2*sigma^2)))
>
> z = convol(array,array2,/center)
> z = z*2
> print,f
> end
>
> here is the message I get:% Out of range subscript encountered: X.
> % Execution halted at: CONV1 29
> /Users/Dave/Desktop/conv1.pro
> % $MAIN$
>
> I don't see why this doesn't work? I am very frustrated

The problem is on this line:

for I = 0,1024 do x[I]= array[0,*]

0 to 1024 is 1025 numbers. (Count them if you
don't believe me.) But X is only big enough to
hold 1024 numbers, so you are, uh, going out
of its subscript range, as the error message
suggests.

But this line of code is completely unnecessary.
Simply typing this is enough:

x = Reform(array[0,*])

Cheers,

David

--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")

Report message to a moderator

[Message index]

		Convolution By: Kay Bente on Tue, 11 September 2001 05:22
		Re: convolution By: Chris[1] on Fri, 28 March 2003 08:17
		Re: convolution By: Vince Hradil on Fri, 25 April 2008 06:29
		Re: convolution By: sarah[1] on Mon, 12 May 2008 18:56
		Re: convolution By: d.poreh on Fri, 25 April 2008 01:15
		Re: convolution By: Chris[5] on Tue, 20 May 2008 04:31
		Re: convolution By: Chris[5] on Tue, 20 May 2008 03:49
		Re: convolution By: sarah[1] on Mon, 19 May 2008 15:32
		Re: convolution By: David Fanning on Mon, 12 May 2008 19:31

Previous Topic:	How to export programmatically a volume with IVOLUME to images
Next Topic:	Re: Least square fitting

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Nov 30 04:30:24 PST 2025

Total time taken to generate the page: 0.96390 seconds