You can have sharp edges and vertices with subdiv surfaces.

You can apply the algorithm on a subset of vertices. Exclude all the sharp edges. An extra falloff algorithm can make it even better