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

Home » Public Forums » archive » fractional powers of negative numbers?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
fractional powers of negative numbers? [message #8837] Thu, 24 April 1997 00:00
Med Bennett is currently offline  Med Bennett
Messages: 109
Registered: April 1997
Senior Member
Can someone explain this to me? Is it possible to compute the fractional
power of a negative number? IDL does it on single negative values but
not on arrays containing negative values.

IDL> print,findgen(11)-5
-5.00000 -4.00000 -3.00000 -2.00000
-1.00000 0.000000 1.00000 2.00000
3.00000 4.00000 5.00000
IDL> print,(findgen(11)-5)^1.51
-1.#IND0 -1.#IND0 -1.#IND0 -1.#IND0
-1.#IND0 0.000000 1.00000 2.84810
5.25355 8.11168 11.3617
% Program caused arithmetic error: Floating illegal operand
IDL> print,-2.^1.51
-2.84810
IDL> print,-5.^1.51
-11.3617

Med Bennett
bennettm@boulder.pti-enviro.com
Re: fractional powers of negative numbers? [message #8838 is a reply to message #8837] Thu, 24 April 1997 00:00 Go to previous message
Ewan A. Macpherson is currently offline  Ewan A. Macpherson
Messages: 14
Registered: December 1996
Junior Member
Med Bennett wrote:

> Can someone explain this to me? Is it possible to compute the fractional
> power of a negative number? IDL does it on single negative values but
> not on arrays containing negative values.

> IDL> print,-2.^1.51
> -2.84810
> IDL> print,-5.^1.51
> -11.3617

These may not give an error, but they're wrong. The results are going to
be complex. IDL is giving you -(2^1.51), not (-2)^1.51, which gives me
an error in PV-WAVE. When you do it on the array it is not treating the
negative sign as a negation operator, but as part of the value, hence
the error. You'll need to do something like this:

print, exp(1.51 * alog(complex(-2,0)))

which gives: ( 0.0894611, -2.84669)

--
Ewan Macpherson <macpherson@earvax.waisman.wisc.edu>
Hearing Development Research Lab, Waisman Center
University of Wisconsin - Madison
http://www.waisman.wisc.edu/~macpherson/
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: IDL Homepage
Next Topic: C structure definition and IDL

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

Current Time: Fri Oct 10 23:11:33 PDT 2025

Total time taken to generate the page: 0.48234 seconds