A student of mine has noticed some strange behavior with some code we have implemented in Rcpp Armadillo. Specifically, when we use matrix multiplication with large sparse m
