[Rd] Thought on crossprod

Thomas Lumley tlumley@u.washington.edu
Fri, 15 Mar 2002 08:42:58 -0800 (PST)

On Fri, 15 Mar 2002, Martin Maechler wrote:

> >>>>> "JonR" == Jonathan Rougier <J.C.Rougier@durham.ac.uk> writes:
>     JonR> does not recognise the case where y is the same as x as special.  I'd
>     JonR> like to extend crossprod to something like
>     >> crossprod
>     JonR> function(x, y = NULL)
>     JonR> if (is.null(y))
>     JonR>    .Internal(crossprod.symm(x))
>     JonR> else
>     JonR>    .Internal(crossprod(x, y))
>     JonR> Before I start mesing up my R-1.4.1 installation, does anyone have any
>     JonR> thoughts on this?
> Just for clarification:
> The main point of your suggestion is to pass only half as much
> data down to C, right?

If this can be done with the BLAS functions then halving the necessary
memory bandwidth will also help -- that's one of the limiting factors in
ATLAS, for example.


r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch