def METHOD_NAME(variant_scalar_spectral, np_rng):
    """d65 emitters should evaluate to the product of D65 and sRGB spectra,
    with intensity factored out when evaluating the sRGB model."""
    import numpy as np
    wavelengths = np.linspace(300, 800, mi.MI_WAVELENGTH_SAMPLES)
    wavelengths += (10 * np_rng.uniform(size=wavelengths.shape)).astype(int)
    d65 = mi.load_dict({ "type" : "d65" })
    ps = mi.PositionSample3f()
    d65_eval = d65.eval(mi.SurfaceInteraction3f(ps, wavelengths))
    for color in [[1, 1, 1], [0.1, 0.2, 0.3], [34, 0.1, 62],
                  [0.001, 0.02, 11.4]]:
        intensity  = (np.max(color) * 2.0) or 1.0
        normalized = np.array(color) / intensity
        d65 = mi.load_dict({
            "type" : "d65",
            "color" : mi.Color3f(color)
        })
        srgb = mi.load_dict({
            "type" : "srgb",
            "color" : mi.Color3f(normalized),
        })
        assert dr.allclose(d65.eval(mi.SurfaceInteraction3f(ps, wavelengths)),
                           d65_eval * intensity * srgb.eval(mi.SurfaceInteraction3f(ps, wavelengths)), atol=1e-5)