Evolution to reduce inbreeding can favor disassortative (inter-morph) over assortative (intramorph) mating in hermaphroditic sexually polymorphic plant species. Heterostyly enhances disassortative pollination through reciprocal placement of stigmas and anthers of morphs and appropriate pollinators. Stylar dimorphism in which there is not reciprocal anther placement may compromise disassortative mating, particularly when there is not intra-morph incompatibility. Variable rates of disassortative mating along with differential female fecundity or siring success among floral morphs could lead to variation in morph ratio. We investigated mating patterns, female fecundity and siring success of style-length morphs in Narcissus papyraceus, a self-incompatible but morph-compatible species with dimorphic (long- and short-styled) and monomorphic (long-styled) populations in central and north regions of its range respectively. We established experimental populations in both regions and exposed them to ambient pollinators. Using paternity analysis, we found similar siring success of morphs and high disassortative mating in most populations. Female fecundity of morphs was similar in all populations. Although these results could not completely explain the loss of dimorphism in the species’ northern range, they provided evidence for the evolutionary stability of stylar dimorphism in N. papyraceus in at least some populations. Our findings support the hypothesis that prevailing inter-morph mating is key for the maintenance of stylar dimorphism.