[R] Numerical Mathematics Consortium

Rau, Roland Rau at demogr.mpg.de
Thu Mar 9 13:58:26 CET 2006

Dear all,

I found a link to the Numerical Mathematics Consortium
Their rationale can be summarized (as far as I understood) by the
paragraphs which I copied from their homepage and pasted below.
Maybe this could be of interest also for the development of R?


And here the copy+paste from their homepage:

As modern numerical programming languages have replaced FORTRAN as the
language of choice, and as new interactive design and development tools
were introduced over the years, the standardized representation of
common functions began to disperse again. The problem is as simple as
looking at the Fast Fourier Transform (FFT) function:

t := 0..127
x := sin(2.p..01.t)
w :=fft(x)

t = 0:127
x = sin(2*pi*.01*t)
w = fft(x)

Size of w:
Mathcad = 65
MATLAB = 128

As in this case, even if the syntax is identical, the calculations yield
different results!

By definition, the FFT function of real data produces a result where the
second half of the data is the complex conjugate of the first half.
Mathcad uses this and only returns the first half of the data (i.e. N/2
+ 1). MATLAB returns all of the data (N). Both are valid

The Numerical Mathematics Consortium is proposing a semantic standard
for core numerical functions so as to help solve issues such as this.

This mail has been sent through the MPI for Demographic Rese...{{dropped}}

More information about the R-help mailing list