More options
Change derived data type to arrays, add matmul_2by2 subroutine
Refactorise all matrix multiplications
Change all tests to new definitions of matrices
Add test matmul_2by2