comp.lang.idl-pvwave archive: archive » Re: smooth function and rounding error

Home » Public Forums » archive » Re: smooth function and rounding error

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

Re: smooth function and rounding error [message #69200 is a reply to message #69199]

Fri, 18 December 2009 03:19

simona bellavista
Messages: 56
Registered: December 2009

Member

thank you for your reply

> y will be of the same type as A, so int, float etc...

yes, anyway I am using double

>> y[i]=total(A[i-width/2:i+width/2-1])/width
>
> I believe you want to divide by width+1, because you are considering

ok, I try also this, but according to the previous link it should be
just width, I am trying to follow closely the prescription they give
for smooth function. can you see any discrepancy between their recipe
and my implementation?

> A[i] in your total. A;so, be sure that width is a float!

actually it was a long(it is an even long),because I thought idl would
cast it automatically to double in division.

now I am casting it explicitly:

y[i]=total(A[i-width/2:i+width/2-1])/double(width)

or

y[i]=total(A[i-width/2:i+width/2-1])/double(width+1)

in both cases I cannot see any difference, relative errors always are
~10^-1

any suggestions/hints/comments?

I forgot to say I am on linux and using idl 7.0.3

Report message to a moderator

[Message index]

		Re: smooth function and rounding error By: jeanh on Fri, 18 December 2009 09:18
		Re: smooth function and rounding error By: simona bellavista on Fri, 18 December 2009 07:02
		Re: smooth function and rounding error By: lbusett@yahoo.it on Fri, 18 December 2009 05:05
		Re: smooth function and rounding error By: simona bellavista on Fri, 18 December 2009 03:19
		Re: smooth function and rounding error By: jeanh on Thu, 17 December 2009 16:13

Previous Topic:	array index summations
Next Topic:	ROI points don't seem to get stored properly

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Nov 29 09:25:55 PST 2025

Total time taken to generate the page: 0.96433 seconds