If we have the estimator $(X^{T}WX)^{-1}X^{T}Wy$ where the diagonal of W contains the inverted weights, what is the fastest/most efficient way to solve it?
I know with OLS, the estimator is typically solved using QR decomposition, but does this method work just as well with WLS?