Astronomical surveys of celestial sources produce streams of noisy time series measuring flux versus time (light curves'). Unlike in many other physical domains, however, large (and source- specific) temporal gaps in data arise naturally due to intranight cadence choices as well as diurnal and seasonal constraints$^1-5$. With nightly observations of millions of variable stars and transients from upcoming surveys$^4,6$, efficient and accurate discovery and classification techniques on noisy, irregularly sampled data must be employed with minimal human-in-the-loop involvement. Machine learning for inference tasks on such data traditionally requires the laborious hand- coding of domain-specific numerical summaries of raw data (
features’)