For doing computation inside Coq using matrices, the most efficient way is likely as operations on Coq's primitive arrays. However, to make the connection to abstract MathComp matrices, there needs to be a refinement from them down to primitive arrays of arrays (e.g., via the refinement to lists of lists).
Suggested by a problem faced by @spitters.