Merge pull request #46 from paranoya/feature/tabularCEM_interpolator

Change TabularCEM interpolator
paranoya · Oct 28, 2024 · bab22e6 · bab22e6
2 parents 1e3df60 + 64aa593
commit bab22e6
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/pst/models.py b/src/pst/models.py
@@ -489,7 +489,7 @@ def stellar_mass_formed(self, times: u.Gyr):
         integral : astropy.Quantity
             The cumulative stellar mass formed at each input time.
         """
-        interpolator = interpolate.Akima1DInterpolator(
+        interpolator = interpolate.PchipInterpolator(
            self.table_t, self.table_mass)
         integral = interpolator(times) << self.table_mass.unit
         integral[times > self.table_t[-1]] = self.table_mass[-1]
@@ -515,7 +515,7 @@ def ism_metallicity(self, times: u.Gyr):
         z_t: astropy.units.Quantity
             Vector with the ISM metallicity at each input time.
         """
-        interpolator = interpolate.Akima1DInterpolator(
+        interpolator = interpolate.PchipInterpolator(
            self.table_t, self.table_metallicity)
         integral = interpolator(times)
         integral[times > self.table_t[-1]] = self.table_metallicity[-1]