This paper focuses on the design of a dedicated P1 function space to model elliptic boundary value problem on a manifold embedded in a space of higher dimension. Using the traces of the linear P1 shape functions, it introduces an algorithm to reduce the function space into an equivalent space having the same properties than a P1 Lagrange approximation. Convergence