Comment by fn-mote

Isn’t the linear algebra conventional wisdom that matrix ops are ALWAYS the bottleneck?

I’m sure this is true in scientific computing.

In R maybe a bunch of resampling would be expected to dominate?