We consider the solution of elliptic problems in general regions by embedding and least squares approximation of overdetermined collocated tensor product of basis functions. The resulting least squares problem will generally be ill-conditioned or even singular, and thus, regularization...