The algorithm only adapts to semiregular meshes whose valence of each vertice is 3,so the arbitrary initial control mesh needs to be semiregularized by preprocessing.