Should with_flux_unit also updates uncertainty unit?

[pllim]

specutils/specutils/spectra/spectrum_mixin.py

Line 88 in cadfe73

def with_flux_unit(self, unit, equivalencies=None, suppress_conversion=False):

I find that while flux unit is updated, the uncertainty stays as old unit, when I tried to use this method for spacetelescope/jdaviz#3168 . In the end, I did not use this method (also the copy is unnecessary for me too but addressing that is irrelevant here).

Should uncertainty unit also be updated in this method?