A few questions:

Is the size of the matrices fixed or variable?

What limitations are you encountering with your current approach?

Is performance an issue beyond: "It would be nice if it went quicker"?

And an unnecessary one, just because I'm intrigued - what's the problem your solving?