work in progress on matrix operations package