We develop a new paradigm for thin-shell finite-element analysis based on the use of subdivision surfaces for: i) describing the geometry of the shell in its undeformed configuration, and ii) generating smooth interpolated displacement fields possessing bounded energy within the strict framework of the Kirchhoff-Love theory of thin shells. The particular subdivision strategy adopted here is Loop's scheme, with extensions such as required to account for creases and displacement boundary conditions. The displacement fields obtained by subdivision are H 2 and, consequently, have a finite Kirchhoff-Love energy. The resulting finite elements contain three nodes and element integrals are computed by a onepoint quadrature. The displacement field of the shell is interpolated from nodal displacements only. In particular, no nodal rotations are used in the interpolation. The interpolation scheme induced by subdivision is nonlocal, i. e., the displacement field over one element depend on the nodal displacements of the element nodes and all nodes of immediately neighboring elements. However, the use of subdivision surfaces ensures that all the local displacement fields thus constructed combine conformingly to define one single limit surface. Numerical tests, including the Belytschko et al. [7] obstacle course of benchmark problems, demonstrate the high accuracy and optimal convergence of the method.