sparse/dense flexible use of solve instead of inv

Could be useful for ``compute_betas`` in IWLS to find a sparse-dennse flexible way of using the more robust ``scipy.linalg.solve`` method rather than directly computing the inverse of X'X, which would then match how we are computing the GWR estimators in ``compute_beta_gwr``.