A new composite representation method is proposed for surface approximation.
The results do not rely on basis functions with the parameterization optimization.
Non-smooth features are fitted well even with smooth polynomial basis functions.
An augmented Lagrange method based numerical algorithm is presented.