Implementation of catmull-clark subdivion for quads