Mailing-List: contact dev-help@commons.apache.org; run by ezmlm
Precedence: bulk
Reply-To: "Commons Developers List" <dev@commons.apache.org>
Received-SPF: neutral (nike.apache.org: local policy)
Message-ID: <4C0E92ED.9080004@free.fr>
Date: Tue, 08 Jun 2010 20:58:53 +0200
From: Luc Maisonobe <Luc.Maisonobe@free.fr>
User-Agent: Mozilla/5.0 (X11; U; Linux x86_64; en-US;
 rv:1.9.1.9) Gecko/20100423 Thunderbird/3.0.4
MIME-Version: 1.0
To: Commons Developers List <dev@commons.apache.org>
Subject: Re: [math] elementary functions.
References: <alpine.DEB.1.10.1006080948260.17482@mantissa.rossi.com>
  <AANLkTimy5Yd-E-UuU5EW3X37FCv1nAH44DNxRec_MmHT@mail.gmail.com>
 <AANLkTimirB51kS8hSoU-tVzlzTqu_bTwibj36yd_3yNe@mail.gmail.com>
 <alpine.DEB.1.10.1006081111350.21986@mantissa.rossi.com>
In-Reply-To: <alpine.DEB.1.10.1006081111350.21986@mantissa.rossi.com>
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: 8bit

Le 08/06/2010 17:29, Bill Rossi a �crit :
> 
> I've looked at StrictMath, generally Math appears to delegates to
> StrictMath.  StrictMath is implemented by the C fdlibm library.
> 
> This is what I've implemented to date, the code will compile and run
> on JDK 1.4 and later.  It may work on older JDKs, but I don't have them
> available for testing.  Note that functions like expm1 are not available
> in earlier JDKs.

This is really impressive and interesting. Do you intend to maintain it
if is included in commons math ? The Apache Software Foundation is all
about growing and maintaining a community around code, it is not only
about code, so people are welcome to join.

Also as you developed it outside of the project, if you want to
contribute it you would have to provide a Software Grant for it (see
<http://www.apache.org/licenses/#grants> from the IP owner, which may be
you or your employer depending on how the code was developed in the
first time.

Thank you for this interesting proposal
Luc

> 
> public class FastMath extends java.lang.Object{
>     public FastMath();
>     public static double exp(double);
>     public static double expm1(double);
>     public static double log(double);
>     public static double log1p(double);
>     public static double log10(double);
>     public static double pow(double, double);
>     public static double sin(double);
>     public static double cos(double);
>     public static double tan(double);
>     public static double atan(double);
>     public static double atan2(double, double);
>     public static double toRadians(double);
>     public static double toDegrees(double);
>     public static double abs(double);
>     public static double ulp(double);
>     public static double floor(double);
>     public static double ceil(double);
>     public static double rint(double);
>     public static long round(double);
>     static {};
> }
> 
> Performance test gives these results:
> 
> Function       Time     Result                  Function     Time      
> Result
> ----------------------------------------------------------------------------------------------
> 
> StrictMath.log 967      1.5118099917827207E8    FastMath.log 553       
> 1.5118099917827207E8
> StrictMath.pow 3199     4.6455095486440872E16   FastMath.pow 1967      
> 4.645509548644088E16
> StrictMath.exp 1079     2.2025454782076317E10   FastMath.exp 562       
> 2.2025454782076317E10
> StrictMath.sin 1151     1839071.8010869955      FastMath.sin 766       
> 1839071.8010869955
> StrictMath.cos 1173     -544020.191353572       FastMath.cos 665       
> -544020.191353572
> StrictMath.tan 1568     -5.024600819552688E7    FastMath.tan 1081      
> -5.024600819552688E7
> StrictMath.atan 1079    1.2403715749052648E7    FastMath.atan 902      
> 1.2403715749052648E7
> StrictMath.expm1 727    -9899999.500018543      FastMath.expm1 773     
> -9899999.500018543
> 
> This table shows execution time for 10,000,000 calls in milliseconds.  
> The result printed is there to prevent the JIT from optimizing away the
> calculation entirely.
> 
> Note that some functions such as exp are nearly twice as fast.   I've
> seen it 3 times faster on different processors.   The preformance varies
> by the relative speed of calculation vs memory lookups.
> 
> The functions are implemented as tables of values in extra precision
> (approx 70 bits), and then interpolated with a minimax polynomial.
> 
> Typical test code:
> 
>       x = 0;
>       time = System.currentTimeMillis();
>       for (int i=0; i<10000000; i++)
>         x+=StrictMath.exp(i/1000000.0);
>       time = System.currentTimeMillis() - time;
>       System.out.print("StrictMath.exp "+time+"\t"+x+"\t");
> 
>       x = 0;
>       time = System.currentTimeMillis();
>       for (int i=0; i<10000000; i++)
>         x+=FastMath.exp(i/1000000.0);
>       time = System.currentTimeMillis() - time;
>       System.out.println("FastMath.exp "+time+"\t"+x);
> 
> To test accuracy, I'd compute results and compare them to an aribitrary
> precision math library.
> 
> On Tue, 8 Jun 2010, James Carman wrote:
> 
>> Have you tried looking at StrictMath?
>>
>> On Tue, Jun 8, 2010 at 10:44 AM, Ted Dunning <ted.dunning@gmail.com>
>> wrote:
>>> Bill,
>>>
>>> Which functions do you have?
>>>
>>> Anything more than the standard sin, cos, exp and log?
>>>
>>>
>>> On Tue, Jun 8, 2010 at 6:52 AM, Bill Rossi <bill@rossi.com> wrote:
>>>
>>>> I have developed over the past year a set of elementary functions
>>>> similar
>>>> to those in java.lang.Math, but with the following characteristics:
>>>>
>>>> * Higher performance.
>>>> * Better accuracy.  Results are accurate to slightly more that +/-
>>>> 0.5 ULP.
>>>> * Pure Java.  The standard Math class is impleneted via JNI, and
>>>> thus takes
>>>> a performance hit.
>>>>
>>>
>>
>> ---------------------------------------------------------------------
>> To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
>> For additional commands, e-mail: dev-help@commons.apache.org
>>
> 
> 
> 
> ---------------------------------------------------------------------
> To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
> For additional commands, e-mail: dev-help@commons.apache.org


---------------------------------------------------------------------
To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
For additional commands, e-mail: dev-help@commons.apache.org