Add BLAS bindings for matrix multiplication

[AtheMathmo]

trafficstars

Following requirements:

• Should be behind a feature flag (to avoid new users running into dependency hell by default).
• Hopefully will fail at compile time if bindings are not found. (Prefer this to silently continuing without failure).

It would also be good to include benchmarks for matrix multiplication across different sizes.