Fast simulation of fully non-stationary wind fields using a new matrix factorization assisted interpolation method