diff --git a/src/stellar_evolution_wrapper.c b/src/stellar_evolution_wrapper.c
index 3a23c8d0fd80f6fbff27cc02c2b06e8fc267cf28..c81f81cb23b201ab64e39d6d9c103626b252137e 100644
--- a/src/stellar_evolution_wrapper.c
+++ b/src/stellar_evolution_wrapper.c
@@ -435,7 +435,7 @@ PyObject* pysupernovae_ii_yields(PyObject* self, PyObject* args) {
     const float log_m2 = log10(m2);
 
     /* Compute the yields */
-    supernovae_ii_get_yields(&fp.stellar_model.snii, log_m1, log_m2, yields);
+    supernovae_ii_get_yields_from_integral(&fp.stellar_model.snii, log_m1, log_m2, yields);
 
     for(int j = 0; j < GEAR_CHEMISTRY_ELEMENT_COUNT; j++) {
       float *y = (float*) PyArray_GETPTR2(yields_out, i, j);
@@ -506,7 +506,7 @@ PyObject* pysupernovae_ii_mass_ejected(PyObject* self, PyObject* args) {
     const float log_m2 = log10(m2);
     float *y = (float*) PyArray_GETPTR1(mass_ejected, i);
 
-    *y = supernovae_ii_get_ejected_mass_fraction_processed(&fp.stellar_model.snii, log_m1, log_m2);
+    *y = supernovae_ii_get_ejected_mass_fraction_processed_from_integral(&fp.stellar_model.snii, log_m1, log_m2);
 
   }
 
@@ -573,7 +573,7 @@ PyObject* pysupernovae_ii_mass_ejected_non_processed(PyObject* self, PyObject* a
     const float log_m2 = log10(m2);
     float *y = (float*) PyArray_GETPTR1(mass_ejected, i);
 
-    *y = supernovae_ii_get_ejected_mass_fraction(&fp.stellar_model.snii, log_m1, log_m2);
+    *y = supernovae_ii_get_ejected_mass_fraction_non_processed_from_integral(&fp.stellar_model.snii, log_m1, log_m2);
     
   }