I learnt this from looking at code for old 3d engines, where you have to first make up your matrix functions, but then also unroll them for optimising stuff like rotations round a single axis.