One of the attractive benefits of slotless machines is low losses at high speeds, which could be emphasized by a careful stator core loss assessment, potentially available already at the pre-design stage. Unfortunately, mainstream iron loss estimation methods are typically implemented in the finite element analysis (FEA) environment with a constant-coefficients dummy model, leading to weak extrapolations with huge errors. In this paper, an analytical method for iron loss prediction in the stator core of slotless PM machines is derived. It is based on the extension of the 2-D field solution over the entire machine geometry. Then, the analytical solution is combined with variable- or constant-coefficient loss models (i.e., VARCO or CCM), which can be efficiently computed by vectorized post-processing. VARCO loss models are shown to be preferred at a general level.Moreover, the paper proposes a lookup-table-based (LUT) solution as an alternative approach. The main contribution lies in the numerical link between the analytical field solution and the iron loss estimate, with the aid of a code implementation of the proposed methodology. First, the models are compared against a sufficiently dense dataset available from laminations manufacturer for validation purposes. Then, all the methods are compared for the slotless machine case. Finally, the models are applied to a real case study and validated experimentally.