Added a matrix class and brought lapacke in as a dependency. For now, lapacke will be the linear solver until Petsc is needed.