rename matrix_mul to vector_mul
add matrix multiplication