help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

precision loss in besselh, bessely for n > 16


From: jessica . alfonsi
Subject: precision loss in besselh, bessely for n > 16
Date: Fri, 5 Jun 2009 11:23:05 +0200 (CEST)
User-agent: SquirrelMail/1.4.17

Hello,
I'm running some MatLab codes for Mie theory found in a book, which I
ported very easily to Octave 3.0.1 on Ubuntu Jaunty 9.04. These codes rely
on the computation of special Bessel functions, in particular I
experienced problems (NaN, Inf + Infi) with besselh and bessely due to
precision loss for order higher than 16. With lower orders the codes
produced the correct results as with MatLab.

See for example the following command in Octave:

besselh(50,[0.1:0.1:2.0]')

which gives an array with all elements
Inf + Infi

Note also that the same code on MatLab 2009a (now expired) trial run
without issues for order up to 50. It would be very important for my
research to be able to run these codes in Octave for orders of besselh as
close as those obtained with MatLab. Is there any workaround to improve
this limit ?

Thank you for your attention

Best

Jessica



reply via email to

[Prev in Thread] Current Thread [Next in Thread]