[Rd] Extending %*%

[Uwe Ligges]

Duncan Murdoch wrote:
> I have lists of matrices stored in various ways.  I'd like to extend
> %*% to work on these.  Is this possible, or should I create my own new
> operator?
> 
> A simplified example would be as follows:
> 
> A <- list( A1, A2, A3)
> B <- list( B1, B2, B3)
> 
> where A1,...,B3 are all matrices, and I'd like A %*% B to return 
> 
> list( A1 %*% B1, A2 %*% B2, A3 %*% B3)
> 
> In the real case A and B are sometimes arrays, sometimes something
> else; lists haven't come up so far, but the principle should be the
> same.
> 
> Duncan Murdoch

The straightforward way without extending %*% (I feel that extending it 
is dangerous, in a way) is
   mapply("%*%", A, B, SIMPLIFY = FALSE)

Uwe Ligges