Better initialization for most factorizations

[ChrisRackauckas]

LU factorizations use lu_instance, but the other factorizations do not init a cache. This means that blank A, such as a just initialized Jacobian, can call the factorization which costs an extra factorization and can be singular due to having blank values.

The fix might be https://github.com/SciML/LinearSolve.jl/issues/61