I have an expression that contains both matrices and scalars:
H0 and H[k] are (non-commuting) matrices, and all the us are scalars of the underlying field (complex).
How do I inform Mathematica of this fact?
In particular, I want to compute expressions that contain commuations of various of these A matrices. When I use . or ** as the operator for matrix multiplications, Mathematica is unable to simplify the results (because, naturally, it assumes that u do not commute either).
I have come across the NCalgebra package, though I am unsure how I would apply it to my problem, as I have an symbolic number n+1 of members in the algebra.

.only in the cases you want it to—you could also maybe make it look the way you'd want via e.g. the notation package. (You could also maybe unprotectDotand add some definitions to it that perform the simplifications you're looking for...) $\endgroup$