Left: empirical fits to g and r-band LCs used for interpolation. The LCs were fit with both a purely empirical Contardo et al. (2000) model and with smoothing splines. Both these models were modified by a Gaussian process that allowed for both random noise and smooth deviations from the fits. The resulting models were used for interpolation. Right: restframe g − r color evolution of SN 2018bcc using the interpolations from the left panel. Both methods produced similar results. The black data points represent spline interpolation of the g and r bands to each other at the epochs of photometry with a linear interpolation of the photometric error added in quadrature.