R has a crossprod function.
"Given matrices x and y as arguments, return a matrix cross-product. This is formally equivalent to (but usually slightly faster than) the call t(x) %% y (crossprod) or x %% t(y) (tcrossprod)."
I suggest that a wrapper for gemm be created to do this.
