Feature/use matrix

Simple factorization using the internal library seems to work (the objective decreases and the relative error decreases on artificial matrices).