sparse/dense flexible use of solve instead of inv

[TaylorOshan]

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.