Improve package efficiency using more efficient matrix products

Check if crossprod or equivalent tricks can not be used to improve computation speed.