Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]


Groups > comp.lang.java.programmer > #16614

Re: What is the command to do a power of a value

Date 2012-07-30 07:48 -0700
From Patricia Shanahan <pats@acm.org>
Newsgroups comp.lang.java.programmer
Subject Re: What is the command to do a power of a value
References <05854398-bcfe-4de2-9eda-7b748239f025@googlegroups.com> <lfkc18hncb6tnv6q66937v5re36onqfl3q@4ax.com> <slrnk1d2km.u9l.avl@gamma.logic.tuwien.ac.at>
Message-ID <xcSdnUmqfrZJA4vNnZ2dnUVZ_uWdnZ2d@earthlink.com> (permalink)

Show all headers | View raw


On 7/30/2012 6:22 AM, Andreas Leitgeb wrote:
> On Sun, 29 Jul 2012 21:28:11 -0700 (PDT), xvictoryeohx@gmail.com wrote:
>> C=L(1+i/100)power of n
>
>   x ^ n  =  exp ( log(x) * n )    |  x = (1 + i/100)
>          = exp ( log( 1 + i/100 ) * n)
>          = exp ( log1p ( i/100 ) * n)
>
> If you're doing more calculations with same interest-rate but
> different periods, then you may want to calculate
>    double logBase = Math.log1p( i / 100 );
> once, and use that for the individual calculations:
>    C = L * Math.exp( logBase * n )
>
> The gist of this response is, that for the kind of base (1+i/100),
> you better separate the pow operation out into log and exp, and
> actually use log1p on (i/100) instead of log on (1+i/100) for
> efficiency's and precision's sake.
>
> For "Math.log1p" see:
> http://docs.oracle.com/javase/6/docs/api/java/lang/Math.html#log1p%28double%29
>

I am curious about why you expect this to be more precise than: "The
computed result must be within 1 ulp of the exact result. Results must
be semi-monotonic." (From the pow description). Or do you know of cases
where Math.pow gets an over-large rounding error?

Note that I am not disagreeing with your method for calculating a power
of a number slightly greater than 1, just questioning whether doing it
explicitly gets more precise answers than using Math.pow.

Patricia

Back to comp.lang.java.programmer | Previous | NextPrevious in thread | Next in thread | Find similar | Unroll thread


Thread

What is the command to do a power of a value xvictoryeohx@gmail.com - 2012-07-29 21:28 -0700
  Re: What is the command to do a power of a value glen herrmannsfeldt <gah@ugcs.caltech.edu> - 2012-07-30 04:31 +0000
    Re: What is the command to do a power of a value xvictoryeohx@gmail.com - 2012-07-29 21:57 -0700
  Re: What is the command to do a power of a value Roedy Green <see_website@mindprod.com.invalid> - 2012-07-30 02:22 -0700
    Re: What is the command to do a power of a value Andreas Leitgeb <avl@gamma.logic.tuwien.ac.at> - 2012-07-30 13:22 +0000
      Re: What is the command to do a power of a value Patricia Shanahan <pats@acm.org> - 2012-07-30 07:48 -0700
        Re: What is the command to do a power of a value Andreas Leitgeb <avl@gamma.logic.tuwien.ac.at> - 2012-07-30 23:06 +0000
  Re: What is the command to do a power of a value bugbear <bugbear@trim_papermule.co.uk_trim> - 2012-07-30 10:44 +0100

csiph-web