Rewrite Unit.__rmul__ and __rdiv__ to speed up array [*/] unit

This partially addresses #3174, by making `<array> * <unit>` faster by a factor 2.5.  See https://github.com/astropy/astropy/issues/3174#issuecomment-65345656