This paper presents an adaptation of a method of "importance sampling" to simulate large ensembles of extreme heat waves (i.e., the most extreme heat waves that could be), given a fixed returned period. We illustrate how this algorithm works for European heat waves and investigate the atmospheric features of such ensembles of events. We argue that such an algorithm can be used to simulate other types of events, including cold spells or prolonged episodes of precipitation.