Re: [Rd] Speeding up matrix multiplies

[James Cloos]

I just tried out the effects of using R's NA value¹ with C arithmetic
on an amd64 linux box.

I always got a NAN result for which R's R_IsNA() would return true.

At least on this platform, NAN's propagate w/o a change in their
lower 32 bits.

If NA is supposed to propagate the way NaN is spec'ed to in IEEE754,
then on some platforms it should be OK to skip the NA/NaN checks and
let the optimised routines handle all of the matrix multiplies.

A configure test could be used to determine whether the current
platform is one where NA/NaN checks are required.

-JimC

1] given:
 union rd { double d; uint64_t l; }; union rd u;
   then:
 u.d is an R NA iff both of:
 u.l & 0x7ff0 == 0x7ff0
 u.l & 0x == 1954;
     are true.
-- 
James Cloos  OpenPGP: 1024D/ED7DAEA6

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel

Re: [Rd] Development version of R: Improved nchar(), nzchar() but changed API

[James Cloos]

>>>>> "GC" == Gábor Csárdi  writes:

GC> You can get an RSS/Atom feed, however, if that's good:
GC> https://github.com/wch/r-source/commits/master.atom

That is available in gwene/gmane as:

 gwene.com.github.wch.r-source.commits.trunk

-JimC
-- 
James Cloos  OpenPGP: 0x997A9F17ED7DAEA6

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel