In a category $C$ with biproducts, morphisms between finite biproducts are naturally encoded in terms of arrays of morphisms between the direct summands of the objects. The natural operations on morphisms (addition, composition) correspond to the usual matrix calculus operations on these arrays.
For the special case that $C =$ Vect this reproduces the standard matrix calculus of linear algebra.
Let $f : X \to Y$ be a morphism in a category with biproducts where the objects $X$ and $Y$ are given as direct sums
Since a biproduct is both a product as well as a coproduct, the morphism $f$ is fixed by all its compositions $f^i_j$ with the product projections $\pi^i : Y \to Y_i$ and the coproduct injections $\iota_j : X_j \to X$:
In matrix calculus one therefore writes
With this notation one has the following rules for computation:
matrix addition
matrix multiplication
where in each case the sum of morphisms is taken using the canonical enrichment of $C$ in abelian monoids (as described at biproduct).
As can be seen in the above formulas, particularly for matrix multiplication, this is a context is which the Einstein summation convention can be used, with a distinction drawn between upper and lower indices. Then repeated indices (in formulas with general applicability) will always appear once upper and once lower, summed over. However, this convention can apply only to the morphisms, not to the objects.
If the category $C$ is in addition a dagger category with an obvious compatibility condition between the dagger operation $(-)^\dagger : C \to C$ and the biproduct structure, then the usual rules of computation for matrices over complex numbers have analogs in $C$.
conjugation
Here the distinction between upper and lower indices cannot be maintained, although it is still true that repeated indices will be summed in formulas with general applicability.
Discussion in the generality of monoidal category theory is in
Formalization in terms of dependent linear type theory is in