Constructing of 16x16 Involutory Binary Matrices of Branch Number 7
I guess you can construct those matrices simply with M1h and M2h with a double and add algorithm. 

Clearly, $M_{4h} = M_{2h} \cdot M_{2h}$ and $M_{Fh} = ((((M_{2h} \oplus M_{1h}) \cdot M_{2h}) \oplus M_{1h}) \cdot M_{2h}) \oplus M_{1h} $. From this you can deduce how to construct all the remaining matrices.