mxlpy 0.21.0__py3-none-any.whl → 0.22.0__py3-none-any.whl

mxlpy/types.py

@@ -36,6 +36,7 @@ __all__ = [
     "McSteadyStates",
     "MockSurrogate",
     "Param",
+    "Parameter",
     "ProtocolByPars",
     "RateFn",
     "Reaction",
@@ -45,13 +46,14 @@ __all__ = [
     "RetType",
     "SteadyStates",
     "TimeCourseByPars",
+    "Variable",
     "unwrap",
     "unwrap2",
 ]
 
 type RateFn = Callable[..., float]
 type Array = NDArray[np.floating[Any]]
-type ArrayLike = NDArray[np.floating[Any]] | list[float]
+type ArrayLike = NDArray[np.floating[Any]] | pd.Index | list[float]
 
 
 Param = ParamSpec("Param")
@@ -59,6 +61,8 @@ RetType = TypeVar("RetType")
 
 
 if TYPE_CHECKING:
+    import sympy
+
     from mxlpy.model import Model
 
 
@@ -147,34 +151,53 @@ type IntegratorType = Callable[
 ]
 
 
+@dataclass
+class Parameter:
+    """Container for parameter meta information."""
+
+    value: float
+    unit: sympy.Expr | None = None
+    source: str | None = None
+
+
+@dataclass
+class Variable:
+    """Container for variable meta information."""
+
+    initial_value: float | Derived
+    unit: sympy.Expr | None = None
+    source: str | None = None
+
+
 @dataclass(kw_only=True, slots=True)
 class Derived:
     """Container for a derived value."""
 
     fn: RateFn
     args: list[str]
+    unit: sympy.Expr | None = None
 
-    def calculate(self, dependent: dict[str, Any]) -> float:
+    def calculate(self, args: dict[str, Any]) -> float:
         """Calculate the derived value.
 
         Args:
-            dependent: Dictionary of dependent variables.
+            args: Dictionary of args variables.
 
         Returns:
             The calculated derived value.
 
         """
-        return cast(float, self.fn(*(dependent[arg] for arg in self.args)))
+        return cast(float, self.fn(*(args[arg] for arg in self.args)))
 
-    def calculate_inpl(self, name: str, dependent: dict[str, Any]) -> None:
+    def calculate_inpl(self, name: str, args: dict[str, Any]) -> None:
         """Calculate the derived value in place.
 
         Args:
             name: Name of the derived variable.
-            dependent: Dictionary of dependent variables.
+            args: Dictionary of args variables.
 
         """
-        dependent[name] = cast(float, self.fn(*(dependent[arg] for arg in self.args)))
+        args[name] = cast(float, self.fn(*(args[arg] for arg in self.args)))
 
 
 @dataclass(kw_only=True, slots=True)
@@ -183,28 +206,29 @@ class Readout:
 
     fn: RateFn
     args: list[str]
+    unit: sympy.Expr | None = None
 
-    def calculate(self, dependent: dict[str, Any]) -> float:
+    def calculate(self, args: dict[str, Any]) -> float:
         """Calculate the derived value.
 
         Args:
-            dependent: Dictionary of dependent variables.
+            args: Dictionary of args variables.
 
         Returns:
             The calculated derived value.
 
         """
-        return cast(float, self.fn(*(dependent[arg] for arg in self.args)))
+        return cast(float, self.fn(*(args[arg] for arg in self.args)))
 
-    def calculate_inpl(self, name: str, dependent: dict[str, Any]) -> None:
+    def calculate_inpl(self, name: str, args: dict[str, Any]) -> None:
         """Calculate the reaction in place.
 
         Args:
             name: Name of the derived variable.
-            dependent: Dictionary of dependent variables.
+            args: Dictionary of args variables.
 
         """
-        dependent[name] = cast(float, self.fn(*(dependent[arg] for arg in self.args)))
+        args[name] = cast(float, self.fn(*(args[arg] for arg in self.args)))
 
 
 @dataclass(kw_only=True, slots=True)
@@ -214,35 +238,36 @@ class Reaction:
     fn: RateFn
     stoichiometry: Mapping[str, float | Derived]
     args: list[str]
+    unit: sympy.Expr | None = None
 
     def get_modifiers(self, model: Model) -> list[str]:
         """Get the modifiers of the reaction."""
-        include = set(model.variables)
+        include = set(model.get_variable_names())
         exclude = set(self.stoichiometry)
 
         return [k for k in self.args if k in include and k not in exclude]
 
-    def calculate(self, dependent: dict[str, Any]) -> float:
+    def calculate(self, args: dict[str, Any]) -> float:
         """Calculate the derived value.
 
         Args:
-            dependent: Dictionary of dependent variables.
+            args: Dictionary of args variables.
 
         Returns:
             The calculated derived value.
 
         """
-        return cast(float, self.fn(*(dependent[arg] for arg in self.args)))
+        return cast(float, self.fn(*(args[arg] for arg in self.args)))
 
-    def calculate_inpl(self, name: str, dependent: dict[str, Any]) -> None:
+    def calculate_inpl(self, name: str, args: dict[str, Any]) -> None:
         """Calculate the reaction in place.
 
         Args:
             name: Name of the derived variable.
-            dependent: Dictionary of dependent variables.
+            args: Dictionary of args variables.
 
         """
-        dependent[name] = cast(float, self.fn(*(dependent[arg] for arg in self.args)))
+        args[name] = cast(float, self.fn(*(args[arg] for arg in self.args)))
 
 
 @dataclass(kw_only=True)

mxlpy/units.py

@@ -0,0 +1,128 @@
+"""Unit definitions for MxlPy."""
+
+from sympy.physics.units import (
+    ampere,
+    becquerel,
+    candela,
+    coulomb,
+    farad,
+    gram,
+    gray,
+    henry,
+    hertz,
+    hour,
+    joule,
+    katal,
+    kelvin,
+    kilogram,
+    liter,
+    lux,
+    meter,
+    micro,
+    milli,
+    minute,
+    mol,
+    nano,
+    newton,
+    ohm,
+    pascal,
+    pico,
+    radian,
+    second,
+    siemens,
+    steradian,
+    tesla,
+    volt,
+    watt,
+    weber,
+)
+
+__all__ = [
+    "ampere",
+    "becquerel",
+    "coulomb",
+    "dimensionless",
+    "farad",
+    "gram",
+    "gray",
+    "henry",
+    "hertz",
+    "hour",
+    "item",
+    "joule",
+    "katal",
+    "kelvin",
+    "liter",
+    "lumen",
+    "lux",
+    "micro",
+    "milli",
+    "minute",
+    "mmol",
+    "mmol_g",
+    "mmol_s",
+    "mol",
+    "nano",
+    "newton",
+    "nmol",
+    "ohm",
+    "pascal",
+    "ppfd",
+    "radian",
+    "second",
+    "siemens",
+    "sievert",
+    "sqm",
+    "tesla",
+    "volt",
+    "watt",
+    "weber",
+]
+
+# time unit
+per_second = 1 / second  # type: ignore
+per_minute = 1 / minute  # type: ignore
+per_hour = 1 / hour  # type: ignore
+
+
+sqm = meter**2
+cbm = meter**3
+
+mol_s = mol / second  # type: ignore
+mol_m = mol / minute  # type: ignore
+mol_h = mol / hour  # type: ignore
+mol_g = mol / gram  # type: ignore
+
+mmol = mol * milli
+mmol_s = mmol / second
+mmol_m = mmol / minute
+mmol_h = mmol / hour
+mmol_g = mmol / gram
+
+mumol = mol * micro
+mumol_s = mumol / second
+mumol_m = mumol / minute
+mumol_h = mumol / hour
+mumol_g = mumol / gram
+
+nmol = mol * nano
+nmol_s = nmol / second
+nmol_m = nmol / minute
+nmol_h = nmol / hour
+nmol_g = nmol / gram
+
+pmol = mol * pico
+pmol_s = pmol / second
+pmol_m = pmol / minute
+pmol_h = pmol / hour
+pmol_g = pmol / gram
+
+ppfd = mumol / sqm / second
+
+
+# SBML units
+avogadro = 6.02214076e23
+sievert = joule / kilogram  # type: ignore
+lumen = candela * steradian  # type: ignore
+dimensionless = None
+item = 1  # pseudounit for one thing

{mxlpy-0.21.0.dist-info → mxlpy-0.22.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: mxlpy
-Version: 0.21.0
+Version: 0.22.0
 Summary: A package to build metabolic models
 Author-email: Marvin van Aalst <marvin.vanaalst@gmail.com>
 Maintainer-email: Marvin van Aalst <marvin.vanaalst@gmail.com>

{mxlpy-0.21.0.dist-info → mxlpy-0.22.0.dist-info}/RECORD

@@ -1,34 +1,36 @@
-mxlpy/__init__.py,sha256=4pbDeyLhQjfL76h2oXdodARzKkrkX5wESV7kEjwC3K8,4399
+mxlpy/__init__.py,sha256=GlBZ-WmSoMgVIeAyVqtEJiRe_jLar84bjbg_FlAr2vw,4477
 mxlpy/carousel.py,sha256=o72YKzfPCDhT5oHhow4oNvIShG-i3-Z0UMEQLt2iE5A,4699
-mxlpy/compare.py,sha256=6iIl6yKXP9guSVLgqqnaqILP_BF_oqyx7DTGbdpwAjM,7800
+mxlpy/compare.py,sha256=PJbb_R9GTGrkcEpGNpUfwZpZbUmZHNjsMQ_5qPAxVNo,7746
 mxlpy/distributions.py,sha256=ce6RTqn19YzMMec-u09fSIUA8A92M6rehCuHuXWcX7A,8734
-mxlpy/fit.py,sha256=vJGd_kl-MqyI7db96uM9STeBnvPVu1VGLCiGuH8mEKc,19009
+mxlpy/fit.py,sha256=3hGUqJ2tOOToZLMMaJw5M9b6_UlUJwT_MhUvfPmRBd8,22355
 mxlpy/fns.py,sha256=NLxYwa3ylS7SkISBjw_TgQSKEm7WnkZF9wPigX_ZCAM,13915
-mxlpy/identify.py,sha256=yU6ccd0yDJhtyo5gkemMuecZALzjR1KzT0vKPmlL4kg,2107
-mxlpy/label_map.py,sha256=Zla9tey-7_POTE57XNEuCSeTqdAbMWZdj_j_OwokngY,17823
-mxlpy/linear_label_map.py,sha256=5FyD0MMdSGsC3eKeBnpd1LBHyVBqIDWCDjgR8_q6XZo,10289
-mxlpy/mc.py,sha256=bt2DrMaovWO_LM3PfkVr0cvK6k_ZSjLRnudasgC5SRM,17132
-mxlpy/mca.py,sha256=B1bRb_EHim3uJ90KJkTZ5HXrVOBjrctcRCsQG4PXU-U,9351
-mxlpy/model.py,sha256=14gncyYft39rwoiJPb5AynL3whXnZrJY3N7SLExH0Qk,62056
+mxlpy/identify.py,sha256=G-Zyr_l-K2mDtIV1xGrQ52QJxoBYqRoNA5RW6GpbNjs,2213
+mxlpy/label_map.py,sha256=kNzqDVp5X6T4uod-y79d6cItOd7_9jmpojDA1TSPRoE,17872
+mxlpy/linear_label_map.py,sha256=6Ye6IziWGKkYD_5z3FmTVwnJ2T9SvVGN3U4CfXjXseo,10320
+mxlpy/mc.py,sha256=AvvnyNIEvaB9gJyd0RtaFs_uVhZ7Xh2FpljHH04YRkc,17217
+mxlpy/mca.py,sha256=BRNbisYijT2SUT6VdIpIh3Id3VgL3NTtycn0VARDWlE,9375
+mxlpy/model.py,sha256=GD1FGppxDJwvEx4ehcQeOx3qkqr8PxDTQ0TwoixBvKI,71617
 mxlpy/parallel.py,sha256=yLQLw5O4vnPVp_Zmtw1UhPWtB3483foimxQB-TwFKPg,5016
 mxlpy/parameterise.py,sha256=IgbvfEnkmaqVq_5GgFjHqApGUN9CJrsVD3Fr7pg9iFA,758
 mxlpy/paths.py,sha256=TK2wO4N9lG-UV1JGfeB64q48JVDbwqIUj63rl55MKuQ,1022
 mxlpy/plot.py,sha256=PA7tAmy2XXACxBLtdnfpxKUFRzi-lnCQjr7gw_nzxKU,32544
 mxlpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-mxlpy/report.py,sha256=6V_kH5usFtar2lUGLjG5k7WIJjUi1TD5qIO7V_6V3Gc,8773
-mxlpy/scan.py,sha256=_MwwTiL-IQPlUU9e2L1Wc-V95sElcaPw-DkxoXYanzc,19421
-mxlpy/simulator.py,sha256=tdxTUbsWo3SdYi7a1lsMah9D_Dz_FqnMjTfK-3oT2Gk,21091
-mxlpy/types.py,sha256=FXBkwHgQ3v_k4ER49hDqyFMIA6i1BQf8isPma97LJdg,12605
+mxlpy/report.py,sha256=v597yzKecgtoNmlNZ_nVhBCa3czNw0tksOK5mLtAQvE,10631
+mxlpy/scan.py,sha256=WARFQgFwOsO2PzC1d2rGDfgDjYV8cYwc196IB2rzrzY,19416
+mxlpy/simulator.py,sha256=so7Dax42EKppeBUTxWWKxYhdAFKS2jAgWzS04fXuFpY,27751
+mxlpy/types.py,sha256=Z697aAiJ3Dp2nvMvdm60Z8rEhQjKKgrYHa1riZpY-QE,12994
+mxlpy/units.py,sha256=WHtP9s90VR8rUpPLsQ_eQeiW3wMPEryMvLpsDfu0tEo,1946
 mxlpy/experimental/__init__.py,sha256=kZTE-92OErpHzNRqmgSQYH4CGXrogGJ5EL35XGZQ81M,206
-mxlpy/experimental/diff.py,sha256=MoM15rbMAHa7p9Zva8NxIc7K585kHJYKFaD1LnN5e10,9088
+mxlpy/experimental/diff.py,sha256=g5hKvFsEUdEk5OGQ_aQuCxLAnenD_jG4G__EcVfKsx4,9104
 mxlpy/integrators/__init__.py,sha256=OLdcNCDIOiD1Z2LO143YtD47cMadNJt0app41nLAx5o,456
 mxlpy/integrators/int_assimulo.py,sha256=8gLR1D4zJ-TnJ9DTkfkqA2uVE0H2w_npZhZ8RoWZOX8,5013
-mxlpy/integrators/int_scipy.py,sha256=MEwhTNhMMVQE2UFWxv5fifN6TKVjRsyDmyibuuNNiHI,4649
-mxlpy/meta/__init__.py,sha256=Z3LnN3a9qDAJTukHZs_nF_G6DrjKXOqBvOb47rSsAsM,314
-mxlpy/meta/codegen_latex.py,sha256=ocdn_mrjPXllYBwImOBQcFzjFR6LONnBs3fhRIA0yzs,22875
-mxlpy/meta/codegen_modebase.py,sha256=ziUuwod1F10ak7WTj5gcuVL7MLtK65kUhqKGCxgn3mY,3131
-mxlpy/meta/codegen_py.py,sha256=bpwXrGUaf8lO81VVcIh0cbtf4cd84CYDZrL3ngf1FHo,3587
-mxlpy/meta/source_tools.py,sha256=8kZD0_FnO2t8MTG9FvEFOhhU52uXKNpQJW6xDOGWGck,13540
+mxlpy/integrators/int_scipy.py,sha256=82nU6cN4PjPoTEXSK5GuFxG89EnnKBHbjG8A9YDPKnc,4749
+mxlpy/meta/__init__.py,sha256=_Bec5aPJ6YyAkxUXlsQtAy_2VzX0tPGVSj-eGACqrXc,404
+mxlpy/meta/codegen_latex.py,sha256=i4tPvk2-toAYqtf4TynuE9sfUSHUp21AMUgjgFEB0uo,23215
+mxlpy/meta/codegen_model.py,sha256=LT767mzniKgAAtMTbejGcl6YmivNa3J1wZjEsH6iX9M,5067
+mxlpy/meta/codegen_mxlpy.py,sha256=Yx6GtsHUup4bIrEyMw9FYE60Y35xAK81IC-r0EdZDEM,3576
+mxlpy/meta/source_tools.py,sha256=imdNbGfDjmEEI7gLWm3BFUfB1xwhW4Lt6mhKvE0PhLw,20735
+mxlpy/meta/sympy_tools.py,sha256=TJSCFK9yuxwNrUHbZrcayo_g1UQWmiUiw4Z5V_pUye4,2853
 mxlpy/nn/__init__.py,sha256=Qjr-ERsY2lbD75sFBOhCUwEasQDSJKcpBn_kReLZ6oA,633
 mxlpy/nn/_keras.py,sha256=-5zjHRu8OjSiZeuBSIZFyB63uBsNNH5o9y4kBcPnhx8,2263
 mxlpy/nn/_torch.py,sha256=GUJmLU282VU4O-fs3Sz90SEaAnfuXN2MPlBr_tHmvn4,5775
@@ -37,8 +39,8 @@ mxlpy/npe/_keras.py,sha256=ytvXMPK9KUCGOzTQm08_SgafiMb-MOIUdZQV7JjAO40,9721
 mxlpy/npe/_torch.py,sha256=v3joh6lFJJxvYJj--wzmKXL9UMTaIN3h6hPNq0uX9NU,11250
 mxlpy/sbml/__init__.py,sha256=Mt97CddpLi3wIrA1b_5cagLmDdNxAVF_S7QV57Pzw8s,226
 mxlpy/sbml/_data.py,sha256=yYli7ZQ1_pnH9kt5EmcuHM0moQoa43rrFVdrseXlG0o,1136
-mxlpy/sbml/_export.py,sha256=4tU3SVxfEvl0E1urZWHyphkiAeH5HeRO1cODvvrczAQ,20342
-mxlpy/sbml/_import.py,sha256=5odQBdpD93mQJp2bVIabmPo6NK60nxqrdSVB8fEsF_A,22099
+mxlpy/sbml/_export.py,sha256=DibzxWLsHQeI4rmvhpV0AfbDmud3lIcfcdm5BV51-PE,20538
+mxlpy/sbml/_import.py,sha256=aCxNKEO6qaZRc7XWFwzOLE1PV0hXFstJbD-SkLjef2k,22082
 mxlpy/sbml/_mathml.py,sha256=oaU9q5yb9jvDGxDJrqOkbOiurCB1Vv_P99oUwJ7v1VE,24437
 mxlpy/sbml/_name_conversion.py,sha256=93muW47M7qJoE227HKHmThWpPeWsXDN9eM8cRH2pqPs,1340
 mxlpy/sbml/_unit_conversion.py,sha256=dW_I6_Ou09ccwnp6LIdrPriIQnQUK5lJcjzM2Fawm6U,1927
@@ -49,8 +51,8 @@ mxlpy/surrogates/_qss.py,sha256=9w-hPPhdcwybkyaSX0sIfYfvcKu1U5j4HHj4SlgZcYQ,723
 mxlpy/surrogates/_torch.py,sha256=gU0secuRBYgewhNqZmSo6_Xf804dSzwWwIYmdKA7y60,6389
 mxlpy/symbolic/__init__.py,sha256=_vM5YM5I6OH0QDbFt9uGYKO8Z5Vly0wbGuvUScVrPRU,258
 mxlpy/symbolic/strikepy.py,sha256=tzo3uvPpXLDex09hWTuitVzuTNwbgl7jWGjD8g6a8iI,20033
-mxlpy/symbolic/symbolic_model.py,sha256=JFzcIdyfJihvKjef748DMXU6WI8nHjgjIk5BwUuB4HQ,2543
-mxlpy-0.21.0.dist-info/METADATA,sha256=ZKxvZTNcrNCaadkRmKSU_6bQlSsJdTtuz_idNVEVbVs,4601
-mxlpy-0.21.0.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-mxlpy-0.21.0.dist-info/licenses/LICENSE,sha256=lHX9Eu70g3Iv1aOxXTWNHa3vq9vaVYSPQx4jOLYmDpw,1096
-mxlpy-0.21.0.dist-info/RECORD,,
+mxlpy/symbolic/symbolic_model.py,sha256=cKfWoktvFmXjuo8egE7gXKrKBq2iBUiy_BcIKIvvz8A,3026
+mxlpy-0.22.0.dist-info/METADATA,sha256=LsEH4wwslk1cTqQzL_JvkHvPlk5zdfnNebzUzPPhHeY,4601
+mxlpy-0.22.0.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+mxlpy-0.22.0.dist-info/licenses/LICENSE,sha256=lHX9Eu70g3Iv1aOxXTWNHa3vq9vaVYSPQx4jOLYmDpw,1096
+mxlpy-0.22.0.dist-info/RECORD,,

mxlpy/meta/codegen_py.py

@@ -1,115 +0,0 @@
-"""Module to export models as code."""
-
-from __future__ import annotations
-
-import warnings
-from typing import TYPE_CHECKING
-
-import sympy
-
-from mxlpy.meta.source_tools import fn_to_sympy, sympy_to_inline
-from mxlpy.types import Derived
-
-if TYPE_CHECKING:
-    from collections.abc import Callable, Generator, Iterable, Iterator
-
-    from mxlpy.model import Model
-
-__all__ = [
-    "generate_model_code_py",
-]
-
-
-def _conditional_join[T](
-    iterable: Iterable[T],
-    question: Callable[[T], bool],
-    true_pat: str,
-    false_pat: str,
-) -> str:
-    """Join an iterable, applying a pattern to each element based on a condition."""
-
-    def inner(it: Iterator[T]) -> Generator[str, None, None]:
-        yield str(next(it))
-        while True:
-            try:
-                el = next(it)
-                if question(el):
-                    yield f"{true_pat}{el}"
-                else:
-                    yield f"{false_pat}{el}"
-            except StopIteration:
-                break
-
-    return "".join(inner(iter(iterable)))
-
-
-def _list_of_symbols(args: list[str]) -> list[sympy.Symbol | sympy.Expr]:
-    return [sympy.Symbol(arg) for arg in args]
-
-
-# FIXME: generate from SymbolicModel, should be easier?
-def generate_model_code_py(model: Model) -> str:
-    """Transform the model into a single function, inlining the function calls."""
-    source = [
-        "from collections.abc import Iterable\n",
-        "from mxlpy.types import Float\n",
-        "def model(t: Float, variables: Float) -> Iterable[Float]:",
-    ]
-
-    # Variables
-    variables = model.variables
-    if len(variables) > 0:
-        source.append("    {} = variables".format(", ".join(variables)))
-
-    # Parameters
-    parameters = model.parameters
-    if len(parameters) > 0:
-        source.append("\n".join(f"    {k} = {v}" for k, v in model.parameters.items()))
-
-    # Derived
-    for name, derived in model.derived.items():
-        expr = fn_to_sympy(derived.fn, model_args=_list_of_symbols(derived.args))
-        source.append(f"    {name} = {sympy_to_inline(expr)}")
-
-    # Reactions
-    for name, rxn in model.reactions.items():
-        expr = fn_to_sympy(rxn.fn, model_args=_list_of_symbols(rxn.args))
-        source.append(f"    {name} = {sympy_to_inline(expr)}")
-
-    # Stoichiometries; FIXME: do this with sympy instead as well?
-    stoich_srcs = {}
-    for rxn_name, rxn in model.reactions.items():
-        for i, (cpd_name, factor) in enumerate(rxn.stoichiometry.items()):
-            if isinstance(factor, Derived):
-                expr = fn_to_sympy(factor.fn, model_args=_list_of_symbols(factor.args))
-                src = f"{sympy_to_inline(expr)} * {rxn_name}"
-            elif factor == 1:
-                src = rxn_name
-            elif factor == -1:
-                src = f"-{rxn_name}" if i == 0 else f"- {rxn_name}"
-            else:
-                src = f"{factor} * {rxn_name}"
-            stoich_srcs.setdefault(cpd_name, []).append(src)
-    for variable, stoich in stoich_srcs.items():
-        source.append(
-            f"    d{variable}dt = {_conditional_join(stoich, lambda x: x.startswith('-'), ' ', ' + ')}"
-        )
-
-    # Surrogates
-    if len(model._surrogates) > 0:  # noqa: SLF001
-        warnings.warn(
-            "Generating code for Surrogates not yet supported.",
-            stacklevel=1,
-        )
-
-    # Return
-    if len(variables) > 0:
-        source.append(
-            "    return {}".format(
-                ", ".join(f"d{i}dt" for i in variables),
-            ),
-        )
-    else:
-        source.append("    return ()")
-
-    return "\n".join(source)