In terms of smooth surface generation over arbitrary topological mesh, simplicity of geometric calculation, and easiness of LOD, subdivision surface has the attractive features as a surface to be fit in reverse engineering (RE). Recently, several fitting algorithms have been proposed for RE using subdivision surfaces. However, these are not suitable for RE in the industrial design process where the physical mock-up models have very complex product geometries. To solve this problem, we propose a new method for fitting a Loop subdivision surface to a given high resolution point cloud by using a quasi-interpolation and a non-uniform subdivision algorithm. The resulting initial control mesh that was fit to about 1,000,000 points had 12,000 faces with 10^<-6>6 normalized fitting error.