Implement SED models and GalaxySED model generator

src/pst/SSP.py

@@ -230,8 +230,12 @@ def interpolate_sed(self, new_wl, verbose=True, log=False, **interp_kwargs):
         ----------
         - new_wl: bin centers of the new interpolated points.
         """
-        if not isinstance(new_wl, units.Quantity):
-            new_wl = new_wl << self.wavelength.unit
+        new_wl = check_unit(new_wl, self.wavelength.unit)
+
+        if new_wl.size == self.wavelength.size and (
+            self.wavelength.value == new_wl.value).all():
+            print("[SSP] SSP already sampled in target wavelength grid")
+            return
 
         if verbose:
             print('[SSP] Interpolating SSP SEDs')

src/pst/__init__.py

@@ -1,6 +1,8 @@
 from . import SSP
 from . import observables
-from . import models
+from . import cem
+#from . import models
 from ._version import get_versions
+
 __version__ = get_versions()["version"]