HydPy 6.3.dev2__tar.gz → 6.4.dev0__tar.gz

{hydpy-6.3.dev2 → hydpy-6.4.dev0/HydPy.egg-info}/PKG-INFO

@@ -1,16 +1,15 @@
 Metadata-Version: 2.4
 Name: HydPy
-Version: 6.3.dev2
+Version: 6.4.dev0
 Summary: A framework for the development and application of hydrological models.
 Home-page: https://github.com/hydpy-dev/hydpy
 Author: HydPy Developers
 Author-email: c.tyralla@bjoernsen.de
-License: LGPL-3.0
+License-Expression: LGPL-3.0-only
 Keywords: hydrology modelling water balance rainfall runoff
 Classifier: Intended Audience :: Education
 Classifier: Intended Audience :: End Users/Desktop
 Classifier: Intended Audience :: Science/Research
-Classifier: License :: OSI Approved :: GNU Lesser General Public License v3 (LGPLv3)
 Classifier: Operating System :: POSIX :: Linux
 Classifier: Operating System :: Microsoft :: Windows
 Classifier: Programming Language :: Python :: 3.11
@@ -45,8 +44,8 @@ Dynamic: classifier
 Dynamic: description
 Dynamic: home-page
 Dynamic: keywords
-Dynamic: license
 Dynamic: license-file
+Dynamic: license-expression
 Dynamic: requires-dist
 Dynamic: requires-python
 Dynamic: summary

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/HydPy.egg-info/SOURCES.txt

@@ -1,20 +1,15 @@
-.gitignore
-.travis.yml
 LICENSE
 README.rst
-appveyor.yml
 call_installer.py
 execute_exe.py
 licenses_hydpy_installer.txt
 logo.ico
 make_and_install_dists.py
-make_hydpy_installer.cfgt
 noxfile.py
 prepare_build.py
 prepare_hydpy_installer.py
 pylintrc
 pyproject.toml
-requirements.txt
 setup.py
 HydPy.egg-info/PKG-INFO
 hydpy/__init__.py

{hydpy-6.3.dev2/HydPy.egg-info → hydpy-6.4.dev0}/PKG-INFO

@@ -1,16 +1,15 @@
 Metadata-Version: 2.4
 Name: HydPy
-Version: 6.3.dev2
+Version: 6.4.dev0
 Summary: A framework for the development and application of hydrological models.
 Home-page: https://github.com/hydpy-dev/hydpy
 Author: HydPy Developers
 Author-email: c.tyralla@bjoernsen.de
-License: LGPL-3.0
+License-Expression: LGPL-3.0-only
 Keywords: hydrology modelling water balance rainfall runoff
 Classifier: Intended Audience :: Education
 Classifier: Intended Audience :: End Users/Desktop
 Classifier: Intended Audience :: Science/Research
-Classifier: License :: OSI Approved :: GNU Lesser General Public License v3 (LGPLv3)
 Classifier: Operating System :: POSIX :: Linux
 Classifier: Operating System :: Microsoft :: Windows
 Classifier: Programming Language :: Python :: 3.11
@@ -45,8 +44,8 @@ Dynamic: classifier
 Dynamic: description
 Dynamic: home-page
 Dynamic: keywords
-Dynamic: license
 Dynamic: license-file
+Dynamic: license-expression
 Dynamic: requires-dist
 Dynamic: requires-python
 Dynamic: summary

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/execute_exe.py

@@ -1 +1,3 @@
+"""Entry point for the HydPy executable for Windows."""
+
 from hydpy.exe.commandtools import execute_scriptfunction

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/__init__.py

@@ -116,7 +116,7 @@ from hydpy.exe.replacetools import xml_replace
 from hydpy.exe.servertools import await_server, start_server
 from hydpy.exe.xmltools import XMLInterface, run_simulation, xml_validate
 
-__version__ = "6.3dev2"
+__version__ = "6.4dev0"
 
 pub.options = optiontools.Options()
 pub.indexer = indextools.Indexer()

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/anntools.py

@@ -1283,8 +1283,8 @@ of element `?` is not properly defined.
 
     def __eq__(self, other: object) -> bool:
         def _equal_array(
-            x: VectorFloat | MatrixInt | MatrixFloat,
-            y: VectorFloat | MatrixInt | MatrixFloat,
+            x: VectorFloat | MatrixInt | MatrixFloat | TensorFloat,
+            y: VectorFloat | MatrixInt | MatrixFloat | TensorFloat,
         ) -> bool:
             idxs = ~(numpy.isnan(x) * numpy.isnan(y))
             return bool(numpy.all(x[idxs] == y[idxs]))

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/armatools.py

@@ -12,6 +12,7 @@ from hydpy import config
 from hydpy.core import exceptiontools
 from hydpy.core import objecttools
 from hydpy.auxs import statstools
+from hydpy.auxs import iuhtools
 from hydpy.core.typingtools import *
 
 if TYPE_CHECKING:
@@ -215,13 +216,61 @@ check the calculated coefficients: 1.0.
     |MA.nmb_nodes_extra| is ignored if set to a smaller value than |MA.nmb_nodes|. 
     """
 
+    _iuh: iuhtools.IUH | None
     _coefs: VectorFloat | None = None
 
-    def __init__(self, iuh=None, coefs=None) -> None:
+    @overload
+    def __init__(self, *, iuh: iuhtools.IUH) -> None: ...
+
+    @overload
+    def __init__(self, *, coefs: VectorInputFloat) -> None: ...
+
+    def __init__(
+        self, *, iuh: iuhtools.IUH | None = None, coefs: VectorInputFloat | None = None
+    ) -> None:
         self.iuh = iuh
         if coefs is not None:
             self.coefs = coefs
 
+    @property
+    def iuh(self) -> iuhtools.IUH:
+        """The currently handled instantaneous unit hydrograph.
+
+        >>> from hydpy import TranslationDiffusionEquation
+        >>> iuh = TranslationDiffusionEquation(u=5.0, d=15.0, x=50.0)
+        >>> from hydpy.auxs.armatools import MA
+        >>> ma = MA(iuh=iuh)
+        >>> assert ma.iuh is iuh
+        >>> del ma.iuh
+        >>> ma.iuh
+        Traceback (most recent call last):
+        ...
+        hydpy.core.exceptiontools.AttributeNotReady: No instantaneous unit hydrograph \
+(`iuh`) available.
+
+        >>> ma.iuh = iuh
+        >>> assert ma.iuh is iuh
+        >>> ma.iuh = None
+        >>> ma.iuh
+        Traceback (most recent call last):
+        ...
+        hydpy.core.exceptiontools.AttributeNotReady: No instantaneous unit hydrograph \
+(`iuh`) available.
+        """
+        if (iuh := self._iuh) is not None:
+            return iuh
+        raise exceptiontools.AttributeNotReady(
+            "No instantaneous unit hydrograph (`iuh`) available."
+        )
+
+    @iuh.setter
+    def iuh(self, iuh: iuhtools.IUH | None) -> None:
+        self._iuh = iuh
+
+    @iuh.deleter
+    def iuh(self) -> None:
+        self._iuh = None
+
     @property
     def coefs(self) -> VectorFloat:
         """|numpy.ndarray| containing all MA coefficents."""
@@ -264,7 +313,7 @@ check the calculated coefficients: 1.0.
         for t in itertools.count(0.0, 1.0):
 
             responses[:n] = responses[n:]
-            responses[n:] = self.iuh(nodes + t)
+            responses[n:] = self.iuh(cast(VectorFloat, nodes + t))
 
             if (t == 0) and (moment1 < 1.0) and ((m := self.nmb_nodes_extra) > n):
                 fine_nodes = numpy.linspace(0.0 / m, 1.0 - 1.0 / m, m)
@@ -339,7 +388,7 @@ check the calculated coefficients: 1.0.
         moment2 = statstools.calc_mean_time_deviation(self.delays, self.coefs, moment1)
         return moment1, moment2
 
-    def plot(self, threshold=None, **kwargs) -> figure.Figure:
+    def plot(self, threshold: float | None = None, **kwargs: Any) -> figure.Figure:
         """Create a barplot of the MA coefficients."""
 
         try:
@@ -583,11 +632,29 @@ by the MA coefficients `1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0`.
     """Maximum deviation of the sum of all coefficents from one to be accepted by 
     method |ARMA.update_coefs|."""
 
+    _ma: MA | None
     _ma_coefs: VectorFloat | None = None
     _ar_coefs: VectorFloat | None = None
     _rel_rmse: float | None
 
-    def __init__(self, ma_model=None, ar_coefs=None, ma_coefs=None) -> None:
+    @overload
+    def __init__(self, *, ma_model: MA) -> None: ...
+
+    @overload
+    def __init__(
+        self,
+        *,
+        ar_coefs: VectorInputFloat | None = None,
+        ma_coefs: VectorInputFloat | None = None,
+    ) -> None: ...
+
+    def __init__(
+        self,
+        *,
+        ma_model: MA | None = None,
+        ar_coefs: VectorInputFloat | None = None,
+        ma_coefs: VectorInputFloat | None = None,
+    ) -> None:
         self.ma = ma_model
         if ar_coefs is not None:
             self.ar_coefs = ar_coefs
@@ -595,6 +662,44 @@ by the MA coefficients `1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0`.
             self.ma_coefs = ma_coefs
         self._rel_rmse = None
 
+    @property
+    def ma(self) -> MA:
+        """The currently handled moving average model.
+
+        >>> from hydpy.auxs.armatools import ARMA, MA
+        >>> ma = MA(coefs=[0.6, 0.4])
+        >>> arma = ARMA(ma_model=ma)
+        >>> assert arma.ma is ma
+        >>> del arma.ma
+        >>> arma.ma
+        Traceback (most recent call last):
+        ...
+        hydpy.core.exceptiontools.AttributeNotReady: No moving average model (`ma`) \
+available
+
+        >>> arma.ma = ma
+        >>> assert arma.ma is ma
+        >>> arma.ma = None
+        >>> arma.ma
+        Traceback (most recent call last):
+        ...
+        hydpy.core.exceptiontools.AttributeNotReady: No moving average model (`ma`) \
+available
+        """
+        if (ma := self._ma) is not None:
+            return ma
+        raise exceptiontools.AttributeNotReady(
+            "No moving average model (`ma`) available"
+        )
+
+    @ma.setter
+    def ma(self, ma: MA | None) -> None:
+        self._ma = ma
+
+    @ma.deleter
+    def ma(self) -> None:
+        self._ma = None
+
     @property
     def rel_rmse(self) -> float:
         """Relative root mean squared error the last time achieved by method
@@ -646,7 +751,7 @@ far.
         return ar_coefs
 
     @ar_coefs.setter
-    def ar_coefs(self, values) -> None:
+    def ar_coefs(self, values: VectorInputFloat) -> None:
         self._ar_coefs = numpy.array(values, ndmin=1, dtype=config.NP_FLOAT)
 
     @ar_coefs.deleter
@@ -797,16 +902,18 @@ far.
         """
         turning_idx, _ = ma_model.turningpoint
         values = ma_model.coefs[turning_idx:]
-        self.ar_coefs, residuals = numpy.linalg.lstsq(
+        result = numpy.linalg.lstsq(
             self.get_a(values, ar_order), self.get_b(values, ar_order), rcond=-1
-        )[:2]
+        )
+        self.ar_coefs = cast(VectorFloat, result[0])
+        residuals = result[1]
         if len(residuals) == 1:
             self._rel_rmse = numpy.sqrt(residuals[0]) / numpy.sum(values)
         else:
             self._rel_rmse = 0.0
 
     @staticmethod
-    def get_a(values, n):
+    def get_a(values: VectorInputFloat, n: int) -> MatrixFloat:
         """Extract the independent variables of the given values and return them as a
         matrix with n columns in a form suitable for the least squares approach applied
         in method |ARMA.update_ar_coefs|.
@@ -817,15 +924,15 @@ far.
             i0 = i - 1 if i > 0 else None
             i1 = i + n - 1
             a[i] = values[i1:i0:-1]
-        return numpy.array(a)
+        return a
 
     @staticmethod
-    def get_b(values, n):
+    def get_b(values: VectorInputFloat, n: int) -> VectorFloat:
         """Extract the dependent variables of the values in a vector with n entries in
         a form suitable for the least squares approach applied in method
         |ARMA.update_ar_coefs|.
         """
-        return numpy.array(values[n:])
+        return numpy.asarray(values[n:])
 
     def update_ma_coefs(self) -> None:
         """Determine the MA coefficients.
@@ -859,7 +966,7 @@ far.
                     f"(`{objecttools.repr_(numpy.min(self.response))}`)."
                 )
 
-    def calc_next_ma_coef(self, ma_order, ma_model) -> None:
+    def calc_next_ma_coef(self, ma_order: int, ma_model: MA) -> None:
         """Determine the MA coefficients of the ARMA model based on its predetermined
         AR coefficients and the MA ordinates of the given |MA| model.
 
@@ -894,7 +1001,7 @@ far.
                     value += ar_coef * values[zdx]
             values.append(value)
             sum_values += value
-        return numpy.array(values)
+        return numpy.asarray(values)
 
     @property
     def moments(self) -> tuple[float, float]:
@@ -906,7 +1013,7 @@ far.
         moment2 = statstools.calc_mean_time_deviation(timepoints, response, moment1)
         return moment1, moment2
 
-    def plot(self, threshold=None, **kwargs) -> figure.Figure:
+    def plot(self, threshold: float | None = None, **kwargs: Any) -> figure.Figure:
         """Barplot of the ARMA response."""
         try:
             # Works under matplotlib 3.

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/calibtools.py

@@ -2881,9 +2881,9 @@ parameterstep="1d"
         ]
         if self.keyword is not None:
             arguments.append(f'keyword="{self.keyword}"')
-        if not numpy.isinf(self.lower):
+        if not math.isinf(self.lower):
             arguments.append(f"lower={objecttools.repr_(self.lower)}")
-        if not numpy.isinf(self.upper):
+        if not math.isinf(self.upper):
             arguments.append(f"upper={objecttools.repr_(self.upper)}")
         if self.parameterstep is not None:
             arguments.append(f'parameterstep="{self.parameterstep}"')

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/interptools.py

@@ -13,12 +13,14 @@ the related methods in the Cython extension module |interputils|.
 from __future__ import annotations
 import abc
 import itertools
+import math
 
 import numpy
 
 import hydpy
 from hydpy import config
 from hydpy.core import exceptiontools
+from hydpy.core import modeltools
 from hydpy.core import objecttools
 from hydpy.core import parametertools
 from hydpy.core import propertytools
@@ -78,7 +80,7 @@ class InterpAlgorithm(_Labeled):
         """Return a string representation of the actual |InterpAlgorithm| object
         prefixed with the given string."""
 
-    def print_table(self, xs: VectorFloat | MatrixFloat) -> None:
+    def print_table(self, xs: VectorInputFloat | MatrixInputFloat) -> None:
         """Process the given input data and print the interpolated output values as
         well as all partial first-order derivatives.
 
@@ -156,8 +158,7 @@ dy1/dx3   dy2/dx3
         table[0, ni:nt] = yns
         table[0, nt:] = [f"d{yn}/d{xn}" for xn, yn in itertools.product(xns, yns)]
 
-        # Mypy problem? See issue https://github.com/python/mypy/issues/8586:
-        xs_: float | Iterable[float]
+        xs_: float | VectorInputFloat
         for ri, xs_ in enumerate(xs):
             ri += 1
             if isinstance(xs_, float):
@@ -220,9 +221,9 @@ class BaseInterpolator(_Labeled):
 
     name: str
     """Class name in lowercase letters."""
-    subvars: parametertools.SubParameters
+    subvars: parametertools.SubParameters[Any]
     """The |SubParameters| object containing the current |BaseInterpolator| object."""
-    subpars: parametertools.SubParameters
+    subpars: parametertools.SubParameters[Any]
     """The |SubParameters| object containing the current |BaseInterpolator| object."""
     fastaccess: parametertools.FastAccessParameter
     """The `fastaccess` object providing access to the interpolator functionalities
@@ -312,7 +313,7 @@ class SimpleInterpolator(BaseInterpolator):
 
     __simpleinterpolator: interputils.SimpleInterpolator | None
 
-    def __init__(self, subvars: parametertools.SubParameters) -> None:
+    def __init__(self, subvars: parametertools.SubParameters[modeltools.Model]) -> None:
         self.subvars = subvars
         self.subpars = subvars
         self.fastaccess = parametertools.FastAccessParameter()
@@ -333,7 +334,12 @@ class SimpleInterpolator(BaseInterpolator):
 
     @property
     def shape(self) -> ShapeHookGet:
-        """This property exists for type consistency; we might remove it later."""
+        """This property exists for type consistency; we might remove it later.
+
+        >>> from hydpy.auxs.interptools import SimpleInterpolator
+        >>> SimpleInterpolator(None).shape
+        ()
+        """
         return ()
 
     @property
@@ -704,7 +710,7 @@ error occurred: Value `1` of type `int` has been given, but an object of type \
     _do_refresh: bool
     __seasonalinterpolator: interputils.SeasonalInterpolator | None
 
-    def __init__(self, subvars: parametertools.SubParameters) -> None:
+    def __init__(self, subvars: parametertools.SubParameters[modeltools.Model]) -> None:
         self.subvars = subvars
         self.subpars = subvars
         self.fastaccess = parametertools.FastAccessParameter()
@@ -925,7 +931,7 @@ interpolation algorithm object, but for parameter `seasonalinterpolator` of elem
     def _prepare_shape(self) -> None:
         """Private on purpose."""
         nmb_weights = timetools.Period("366d") / hydpy.pub.options.simulationstep
-        nmb_weights = int(numpy.ceil(round(nmb_weights, 10)))
+        nmb_weights = int(math.ceil(round(nmb_weights, 10)))
         shape = (nmb_weights, self._seasonalinterpolator.nmb_algorithms)
         getattr(self.fastaccess, self.name).ratios = numpy.zeros(
             shape, dtype=config.NP_FLOAT

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/iuhtools.py

@@ -156,13 +156,13 @@ class IUH(metaclass=MetaIUH):
     _SECONDARY_PARAMETERS: dict[str, SecondaryParameterIUH] = {}
 
     def __init__(self, **kwargs: float) -> None:
-        self.ma = armatools.MA(self)
+        self.ma = armatools.MA(iuh=self)
         self.arma = armatools.ARMA(ma_model=self.ma)
         if kwargs:
             self.set_primary_parameters(**kwargs)
 
     @abc.abstractmethod
-    def __call__(self, t: float) -> float:
+    def __call__(self, t: VectorFloatFlex) -> VectorFloatFlex:
         """Must be implemented by the concrete |IUH| subclass."""
 
     def set_primary_parameters(self, **kwargs: float) -> None:
@@ -224,9 +224,9 @@ class IUH(metaclass=MetaIUH):
             sum_responses += self.dt_response * response
             if (sum_responses > 0.9) and (response < self.smallest_response):
                 break
-        return numpy.array(delays), numpy.array(responses)
+        return numpy.asarray(delays), numpy.asarray(responses)
 
-    def plot(self, threshold: float | None = None, **kwargs) -> None:
+    def plot(self, threshold: float | None = None, **kwargs: Any) -> None:
         """Plot the instanteneous unit hydrograph.
 
         The optional argument allows for defining a threshold of the cumulative sum of
@@ -423,22 +423,17 @@ keywords were given: d and u.
         self.a = self.x / (2.0 * self.d**0.5)
         self.b = self.u / (2.0 * self.d**0.5)
 
-    @overload
-    def __call__(self, t: float) -> float: ...
-
-    @overload
-    def __call__(self, t: VectorFloat) -> VectorFloat: ...
-
-    def __call__(self, t: float | VectorFloat) -> float | VectorFloat:
+    def __call__(self, t: VectorFloatFlex) -> VectorFloatFlex:
         if isinstance(t, float):
             if t < 1e-10:  # pylint: disable=consider-using-max-builtin
                 t = 1e-10
         else:
             t = numpy.clip(t, 1e-10, numpy.inf)
-        return (
+        return cast(
+            VectorFloatFlex,
             self._a
             / (t * numpy.sqrt(numpy.pi * t))
-            * numpy.exp(-t * numpy.square(self._a / t - self._b))
+            * numpy.exp(-t * numpy.square(self._a / t - self._b)),
         )
 
     @property
@@ -503,20 +498,14 @@ class LinearStorageCascade(IUH):
         """Determine the values of the secondary parameters |LinearStorageCascade.c|,
         |LinearStorageCascade.log_c|, and |LinearStorageCascade.log_k|."""
         self.c = 1.0 / (self.k * special.gamma(self.n))
-        self.log_c = -numpy.log(self.k) - special.gammaln(self.n)
-        self.log_k = numpy.log(self.k)
-
-    @overload
-    def __call__(self, t: float) -> float: ...
-
-    @overload
-    def __call__(self, t: VectorFloat) -> VectorFloat: ...
+        self.log_c = -math.log(self.k) - special.gammaln(self.n)
+        self.log_k = math.log(self.k)
 
-    def __call__(self, t: float | VectorFloat) -> float | VectorFloat:
+    def __call__(self, t: VectorFloatFlex) -> VectorFloatFlex:
         if isinstance(t, float):
             if t == 0.0:
                 return 0.0
-            return numpy.e ** (
+            return math.exp(
                 self._log_c
                 + (self._n - 1.0) * (math.log(t) - self._log_k)
                 - t / self._k

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/networktools.py

@@ -198,7 +198,7 @@ basin number.
         return f"RiverBasinNumber({str.__repr__(self)[1:-1]})"
 
 
-class RiverBasinNumbers(tuple):
+class RiverBasinNumbers(tuple[RiverBasinNumber]):
     """A sorted collection of |RiverBasinNumber| objects.
 
     >>> from hydpy import RiverBasinNumbers
@@ -229,7 +229,7 @@ class RiverBasinNumbers(tuple):
     False
     """
 
-    def __new__(cls, values: Iterable[int | str]):
+    def __new__(cls, values: Iterable[int | str]) -> RiverBasinNumbers:
         _values = tuple(RiverBasinNumber(value) for value in values)
         obj = tuple.__new__(RiverBasinNumbers, sorted(_values))
         vars(obj)["_tree"] = None
@@ -249,7 +249,7 @@ class RiverBasinNumbers(tuple):
                         subtree = subtree[digit]
                 subtree["0"] = {}
             vars(self)["_tree"] = tree
-        return vars(self)["_tree"]
+        return cast(_Tree, vars(self)["_tree"])
 
     def select(self, number: int | str) -> RiverBasinNumbers:
         """Select and return all river basin numbers starting with the given number.
@@ -276,7 +276,7 @@ class RiverBasinNumbers(tuple):
             except KeyError:
                 return RiverBasinNumbers(())
 
-        def _walk(number_: str, tree_: _Tree, numbers: set) -> set:
+        def _walk(number_: str, tree_: _Tree, numbers: set[str]) -> set[str]:
             if tree_:
                 numbers.remove(number_)
                 for digit_, subtree in tree_.items():
@@ -360,7 +360,7 @@ class RiverBasinNumbers(tuple):
         """
         return tuple(self._get_next_number(rbn) for rbn in self)
 
-    def __contains__(self, number) -> bool:
+    def __contains__(self, number: object) -> bool:
         if not isinstance(number, (int, str)):
             return False
         tree = self._tree

{hydpy-6.3.dev2 → hydpy-6.4.dev0}/hydpy/auxs/ppolytools.py

@@ -244,8 +244,8 @@ polynomial function by passing at leas one `Poly` object.
                 "polynomial function by passing at leas one `Poly` object."
             )
         nmb_ps = len(polynomials)
-        nmb_cs = numpy.array([len(p.cs) for p in polynomials], dtype=config.NP_INT)
-        x0s = numpy.array([p.x0 for p in polynomials], dtype=config.NP_FLOAT)
+        nmb_cs = numpy.asarray([len(p.cs) for p in polynomials], dtype=config.NP_INT)
+        x0s = numpy.asarray([p.x0 for p in polynomials], dtype=config.NP_FLOAT)
         cs = numpy.zeros((nmb_ps, max(nmb_cs)))
         for idx, (nmb, polynomial) in enumerate(zip(nmb_cs, polynomials)):
             cs[idx, :nmb] = polynomial.cs
@@ -411,16 +411,16 @@ vectors `x` (2) and `y` (3) must be identical.
             if (len(xs) == 2) or (method == "linear"):
                 nmb_ps = len(xs) - 1
                 nmb_cs = numpy.full((nmb_ps,), 2, dtype=config.NP_INT)
-                x0s = numpy.array(xs, dtype=config.NP_FLOAT)[:-1]
+                x0s = numpy.asarray(xs, dtype=config.NP_FLOAT)[:-1]
                 cs = numpy.zeros((nmb_ps, numpy.max(nmb_cs)), dtype=config.NP_FLOAT)
-                cs[:, 0] = numpy.array(ys, dtype=config.NP_FLOAT)[:-1]
+                cs[:, 0] = numpy.asarray(ys, dtype=config.NP_FLOAT)[:-1]
                 cs[:, 1] = numpy.diff(ys) / numpy.diff(xs)
             else:
                 interpolator = method(x=xs, y=ys)
-                x0s = numpy.array(xs, dtype=config.NP_FLOAT)[:-1]
+                x0s = numpy.asarray(xs, dtype=config.NP_FLOAT)[:-1]
                 cs = interpolator.c[::-1].T
                 nmb_ps = len(x0s)
-                nmb_cs = numpy.array(
+                nmb_cs = numpy.asarray(
                     [numpy.max(numpy.nonzero(cs_), initial=0) + 1 for cs_ in cs],
                     dtype=config.NP_INT,
                 )