Faster default matrix multiplication with benchmark to check it

I've updated default implementation and corresponding test accordingly to your feedback in PR #48. It now maintains current develop's behaviour.