fixed handling of input vectors in factorization update routines

The input vectors were not being densified.  Function vector_process was added
to check vector length and return a densified copy.

LUSOL does not support sparse input vectors in the solve, multiply, or update
routines.
the readme was updated

Merge branch 'master' of github.com:nwh/lusol