dcmri.flux_nscomp#

dcmri.flux_nscomp(J, T, t=None, dt=1.0)[source]#

Indicator flux out of a non-stationary compartment.

See section Non-stationary compartment for more detail.

Parameters:

  • J (array_like) – the indicator flux entering the compartment.

  • T (array_like) – array with the mean transit time as a function of time, with the same length as J. Only finite and strictly positive values are allowed.

  • t (array_like, optional) – the time points of the indicator flux J, in the same units as T. If t=None, the time points are assumed to be uniformly spaced with spacing dt. Defaults to None.

  • dt (float, optional) – spacing between time points for uniformly spaced time points, in the same units as T. This parameter is ignored if t is explicity provided. Defaults to 1.0.

Returns:

outflux as a 1D array.

Return type:

numpy.ndarray

See also

conc_nscomp

Example

>>> import dcmri as dc
>>> t = [0,5,15,30,60]
>>> J = [1,2,3,3,2]
>>> T = [1,2,3,4,5]
>>> dc.flux_nscomp(J, T, t)
array([0.        , 1.54942844, 2.65376974, 2.88280904, 2.05727851])