On the Efficacy and High-Performance Implementation of Quaternion Matrix Multiplication