dot-make_matrix function

Make the underlying matrix that will be modified in place by the solvers