brainstate 0.1.7__py2.py3-none-any.whl → 0.1.9__py2.py3-none-any.whl

brainstate/{functional → nn}/_activations.py

@@ -1,813 +1,808 @@
-# Copyright 2024 BDP Ecosystem Limited. All Rights Reserved.
-#
-# Licensed under the Apache License, Version 2.0 (the "License");
-# you may not use this file except in compliance with the License.
-# You may obtain a copy of the License at
-#
-#     http://www.apache.org/licenses/LICENSE-2.0
-#
-# Unless required by applicable law or agreed to in writing, software
-# distributed under the License is distributed on an "AS IS" BASIS,
-# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-# See the License for the specific language governing permissions and
-# limitations under the License.
-# ==============================================================================
-
-
-"""
-Shared neural network activations and other functions.
-"""
-
-from typing import Any, Union, Sequence
-
-import brainunit as u
-import jax
-from jax.scipy.special import logsumexp
-
-from brainstate import random
-from brainstate.typing import ArrayLike
-
-__all__ = [
-    "tanh",
-    "relu",
-    "squareplus",
-    "softplus",
-    "soft_sign",
-    "sigmoid",
-    "silu",
-    "swish",
-    "log_sigmoid",
-    "elu",
-    "leaky_relu",
-    "hard_tanh",
-    "celu",
-    "selu",
-    "gelu",
-    "glu",
-    "logsumexp",
-    "log_softmax",
-    "softmax",
-    "standardize",
-    "one_hot",
-    "relu6",
-    "hard_sigmoid",
-    "hard_silu",
-    "hard_swish",
-    'hard_shrink',
-    'rrelu',
-    'mish',
-    'soft_shrink',
-    'prelu',
-    'tanh_shrink',
-    'softmin',
-    'sparse_plus',
-    'sparse_sigmoid',
-]
-
-
-def tanh(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Hyperbolic tangent activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{tanh}(x) = \frac{e^x - e^{-x}}{e^x + e^{-x}}
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-    """
-    return u.math.tanh(x)
-
-
-def softmin(x, axis=-1):
-    r"""
-    Applies the Softmin function to an n-dimensional input Tensor
-    rescaling them so that the elements of the n-dimensional output Tensor
-    lie in the range `[0, 1]` and sum to 1.
-
-    Softmin is defined as:
-
-    .. math::
-        \text{Softmin}(x_{i}) = \frac{\exp(-x_i)}{\sum_j \exp(-x_j)}
-
-    Shape:
-        - Input: :math:`(*)` where `*` means, any number of additional
-          dimensions
-        - Output: :math:`(*)`, same shape as the input
-
-    Args:
-        axis (int): A dimension along which Softmin will be computed (so every slice
-            along dim will sum to 1).
-    """
-    unnormalized = u.math.exp(-x)
-    return unnormalized / unnormalized.sum(axis, keepdims=True)
-
-
-def tanh_shrink(x):
-    r"""
-    Applies the element-wise function:
-
-    .. math::
-        \text{Tanhshrink}(x) = x - \tanh(x)
-    """
-    return x - u.math.tanh(x)
-
-
-def prelu(x, a=0.25):
-    r"""
-    Applies the element-wise function:
-
-    .. math::
-        \text{PReLU}(x) = \max(0,x) + a * \min(0,x)
-
-    or
-
-    .. math::
-        \text{PReLU}(x) =
-        \begin{cases}
-        x, & \text{ if } x \geq 0 \\
-        ax, & \text{ otherwise }
-        \end{cases}
-
-    Here :math:`a` is a learnable parameter. When called without arguments, `nn.PReLU()` uses a single
-    parameter :math:`a` across all input channels. If called with `nn.PReLU(nChannels)`,
-    a separate :math:`a` is used for each input channel.
-    """
-    return u.math.where(x >= 0., x, a * x)
-
-
-def soft_shrink(x, lambd=0.5):
-    r"""
-    Applies the soft shrinkage function elementwise:
-
-    .. math::
-        \text{SoftShrinkage}(x) =
-        \begin{cases}
-        x - \lambda, & \text{ if } x > \lambda \\
-        x + \lambda, & \text{ if } x < -\lambda \\
-        0, & \text{ otherwise }
-        \end{cases}
-
-    Args:
-        lambd: the :math:`\lambda` (must be no less than zero) value for the Softshrink formulation. Default: 0.5
-
-    Shape:
-        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
-        - Output: :math:`(*)`, same shape as the input.
-    """
-    return u.math.where(x > lambd,
-                        x - lambd,
-                        u.math.where(x < -lambd,
-                                     x + lambd,
-                                     u.Quantity(0., unit=u.get_unit(lambd))))
-
-
-def mish(x):
-    r"""Applies the Mish function, element-wise.
-
-    Mish: A Self Regularized Non-Monotonic Neural Activation Function.
-
-    .. math::
-        \text{Mish}(x) = x * \text{Tanh}(\text{Softplus}(x))
-
-    .. note::
-        See `Mish: A Self Regularized Non-Monotonic Neural Activation Function <https://arxiv.org/abs/1908.08681>`_
-
-    Shape:
-        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
-        - Output: :math:`(*)`, same shape as the input.
-    """
-    return x * u.math.tanh(softplus(x))
-
-
-def rrelu(x, lower=0.125, upper=0.3333333333333333):
-    r"""Applies the randomized leaky rectified liner unit function, element-wise,
-    as described in the paper:
-
-    `Empirical Evaluation of Rectified Activations in Convolutional Network`_.
-
-    The function is defined as:
-
-    .. math::
-        \text{RReLU}(x) =
-        \begin{cases}
-            x & \text{if } x \geq 0 \\
-            ax & \text{ otherwise }
-        \end{cases}
-
-    where :math:`a` is randomly sampled from uniform distribution
-    :math:`\mathcal{U}(\text{lower}, \text{upper})`.
-
-     See: https://arxiv.org/pdf/1505.00853.pdf
-
-    Args:
-        lower: lower bound of the uniform distribution. Default: :math:`\frac{1}{8}`
-        upper: upper bound of the uniform distribution. Default: :math:`\frac{1}{3}`
-
-    Shape:
-        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
-        - Output: :math:`(*)`, same shape as the input.
-
-    .. _`Empirical Evaluation of Rectified Activations in Convolutional Network`:
-        https://arxiv.org/abs/1505.00853
-    """
-    a = random.uniform(lower, upper, size=u.math.shape(x), dtype=x.dtype)
-    return u.math.where(u.get_mantissa(x) >= 0., x, a * x)
-
-
-def hard_shrink(x, lambd=0.5):
-    r"""Applies the Hard Shrinkage (Hardshrink) function element-wise.
-
-    Hardshrink is defined as:
-
-    .. math::
-        \text{HardShrink}(x) =
-        \begin{cases}
-        x, & \text{ if } x > \lambda \\
-        x, & \text{ if } x < -\lambda \\
-        0, & \text{ otherwise }
-        \end{cases}
-
-    Args:
-        lambd: the :math:`\lambda` value for the Hardshrink formulation. Default: 0.5
-
-    Shape:
-        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
-        - Output: :math:`(*)`, same shape as the input.
-
-    """
-    return u.math.where(x > lambd,
-                        x,
-                        u.math.where(x < -lambd,
-                                     x,
-                                     u.Quantity(0., unit=u.get_unit(x))))
-
-
-def _keep_unit(fun, x, **kwargs):
-    unit = u.get_unit(x)
-    x = fun(u.get_mantissa(x), **kwargs)
-    return x if unit.is_unitless else u.Quantity(x, unit=unit)
-
-
-def relu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Rectified linear unit activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{relu}(x) = \max(x, 0)
-
-    except under differentiation, we take:
-
-    .. math::
-      \nabla \mathrm{relu}(0) = 0
-
-    For more information see
-    `Numerical influence of ReLU’(0) on backpropagation
-    <https://openreview.net/forum?id=urrcVI-_jRm>`_.
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    Example:
-      >>> jax.nn.relu(jax.numpy.array([-2., -1., -0.5, 0, 0.5, 1., 2.]))
-      Array([0. , 0. , 0. , 0. , 0.5, 1. , 2. ], dtype=float32)
-
-    See also:
-      :func:`relu6`
-
-    """
-    return _keep_unit(jax.nn.relu, x)
-
-
-def squareplus(x: ArrayLike, b: ArrayLike = 4) -> Union[jax.Array, u.Quantity]:
-    r"""Squareplus activation function.
-
-    Computes the element-wise function
-
-    .. math::
-      \mathrm{squareplus}(x) = \frac{x + \sqrt{x^2 + b}}{2}
-
-    as described in https://arxiv.org/abs/2112.11687.
-
-    Args:
-      x : input array
-      b : smoothness parameter
-    """
-    return _keep_unit(jax.nn.squareplus, x, b=b)
-
-
-def softplus(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Softplus activation function.
-
-    Computes the element-wise function
-
-    .. math::
-      \mathrm{softplus}(x) = \log(1 + e^x)
-
-    Args:
-      x : input array
-    """
-    return _keep_unit(jax.nn.softplus, x)
-
-
-def soft_sign(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Soft-sign activation function.
-
-    Computes the element-wise function
-
-    .. math::
-      \mathrm{soft\_sign}(x) = \frac{x}{|x| + 1}
-
-    Args:
-      x : input array
-    """
-    return _keep_unit(jax.nn.soft_sign, x)
-
-
-def sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Sigmoid activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{sigmoid}(x) = \frac{1}{1 + e^{-x}}
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`log_sigmoid`
-
-    """
-    return _keep_unit(jax.nn.sigmoid, x)
-
-
-def silu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""SiLU (a.k.a. swish) activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{silu}(x) = x \cdot \mathrm{sigmoid}(x) = \frac{x}{1 + e^{-x}}
-
-    :func:`swish` and :func:`silu` are both aliases for the same function.
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`sigmoid`
-    """
-    return _keep_unit(jax.nn.silu, x)
-
-
-swish = silu
-
-
-def log_sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Log-sigmoid activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{log\_sigmoid}(x) = \log(\mathrm{sigmoid}(x)) = -\log(1 + e^{-x})
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`sigmoid`
-    """
-    return _keep_unit(jax.nn.log_sigmoid, x)
-
-
-def elu(x: ArrayLike, alpha: ArrayLike = 1.0) -> Union[jax.Array, u.Quantity]:
-    r"""Exponential linear unit activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{elu}(x) = \begin{cases}
-        x, & x > 0\\
-        \alpha \left(\exp(x) - 1\right), & x \le 0
-      \end{cases}
-
-    Args:
-      x : input array
-      alpha : scalar or array of alpha values (default: 1.0)
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`selu`
-    """
-    return _keep_unit(jax.nn.elu, x)
-
-
-def leaky_relu(x: ArrayLike, negative_slope: ArrayLike = 1e-2) -> Union[jax.Array, u.Quantity]:
-    r"""Leaky rectified linear unit activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{leaky\_relu}(x) = \begin{cases}
-        x, & x \ge 0\\
-        \alpha x, & x < 0
-      \end{cases}
-
-    where :math:`\alpha` = :code:`negative_slope`.
-
-    Args:
-      x : input array
-      negative_slope : array or scalar specifying the negative slope (default: 0.01)
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`relu`
-    """
-    return _keep_unit(jax.nn.leaky_relu, x, negative_slope=negative_slope)
-
-
-def _hard_tanh(x, min_val=- 1.0, max_val=1.0):
-    return jax.numpy.where(x > max_val, max_val, jax.numpy.where(x < min_val, min_val, x))
-
-
-def hard_tanh(
-    x: ArrayLike,
-    min_val: float = - 1.0,
-    max_val: float = 1.0
-) -> Union[jax.Array, u.Quantity]:
-    r"""Hard :math:`\mathrm{tanh}` activation function.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{hard\_tanh}(x) = \begin{cases}
-        -1, & x < -1\\
-        x, & -1 \le x \le 1\\
-        1, & 1 < x
-      \end{cases}
-
-    Args:
-      x : input array
-      min_val: float. minimum value of the linear region range. Default: -1
-      max_val: float. maximum value of the linear region range. Default: 1
-
-    Returns:
-      An array.
-    """
-    return _keep_unit(_hard_tanh, x, min_val=min_val, max_val=max_val)
-
-
-def celu(x: ArrayLike, alpha: ArrayLike = 1.0) -> Union[jax.Array, u.Quantity]:
-    r"""Continuously-differentiable exponential linear unit activation.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{celu}(x) = \begin{cases}
-        x, & x > 0\\
-        \alpha \left(\exp(\frac{x}{\alpha}) - 1\right), & x \le 0
-      \end{cases}
-
-    For more information, see
-    `Continuously Differentiable Exponential Linear Units
-    <https://arxiv.org/pdf/1704.07483.pdf>`_.
-
-    Args:
-      x : input array
-      alpha : array or scalar (default: 1.0)
-
-    Returns:
-      An array.
-    """
-    return _keep_unit(jax.nn.celu, x, alpha=alpha)
-
-
-def selu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Scaled exponential linear unit activation.
-
-    Computes the element-wise function:
-
-    .. math::
-      \mathrm{selu}(x) = \lambda \begin{cases}
-        x, & x > 0\\
-        \alpha e^x - \alpha, & x \le 0
-      \end{cases}
-
-    where :math:`\lambda = 1.0507009873554804934193349852946` and
-    :math:`\alpha = 1.6732632423543772848170429916717`.
-
-    For more information, see
-    `Self-Normalizing Neural Networks
-    <https://papers.nips.cc/paper/6698-self-normalizing-neural-networks.pdf>`_.
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`elu`
-    """
-    return _keep_unit(jax.nn.selu, x)
-
-
-def gelu(x: ArrayLike, approximate: bool = True) -> Union[jax.Array, u.Quantity]:
-    r"""Gaussian error linear unit activation function.
-
-    If ``approximate=False``, computes the element-wise function:
-
-    .. math::
-      \mathrm{gelu}(x) = \frac{x}{2} \left(1 + \mathrm{erf} \left(
-        \frac{x}{\sqrt{2}} \right) \right)
-
-    If ``approximate=True``, uses the approximate formulation of GELU:
-
-    .. math::
-      \mathrm{gelu}(x) = \frac{x}{2} \left(1 + \mathrm{tanh} \left(
-        \sqrt{\frac{2}{\pi}} \left(x + 0.044715 x^3 \right) \right) \right)
-
-    For more information, see `Gaussian Error Linear Units (GELUs)
-    <https://arxiv.org/abs/1606.08415>`_, section 2.
-
-    Args:
-      x : input array
-      approximate: whether to use the approximate or exact formulation.
-    """
-    return _keep_unit(jax.nn.gelu, x, approximate=approximate)
-
-
-def glu(x: ArrayLike, axis: int = -1) -> Union[jax.Array, u.Quantity]:
-    r"""Gated linear unit activation function.
-
-    Computes the function:
-
-    .. math::
-      \mathrm{glu}(x) =  x\left[\ldots, 0:\frac{n}{2}, \ldots\right] \cdot
-        \mathrm{sigmoid} \left( x\left[\ldots, \frac{n}{2}:n, \ldots\right]
-          \right)
-
-    where the array is split into two along ``axis``. The size of the ``axis``
-    dimension must be divisible by two.
-
-    Args:
-      x : input array
-      axis: the axis along which the split should be computed (default: -1)
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`sigmoid`
-    """
-    return _keep_unit(jax.nn.glu, x, axis=axis)
-
-
-def log_softmax(x: ArrayLike,
-                axis: int | tuple[int, ...] | None = -1,
-                where: ArrayLike | None = None) -> Union[jax.Array, u.Quantity]:
-    r"""Log-Softmax function.
-
-    Computes the logarithm of the :code:`softmax` function, which rescales
-    elements to the range :math:`[-\infty, 0)`.
-
-    .. math ::
-      \mathrm{log\_softmax}(x)_i = \log \left( \frac{\exp(x_i)}{\sum_j \exp(x_j)}
-      \right)
-
-    Args:
-      x : input array
-      axis: the axis or axes along which the :code:`log_softmax` should be
-        computed. Either an integer or a tuple of integers.
-      where: Elements to include in the :code:`log_softmax`.
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`softmax`
-    """
-    return _keep_unit(jax.nn.log_softmax, x, axis=axis, where=where)
-
-
-def softmax(x: ArrayLike,
-            axis: int | tuple[int, ...] | None = -1,
-            where: ArrayLike | None = None) -> Union[jax.Array, u.Quantity]:
-    r"""Softmax function.
-
-    Computes the function which rescales elements to the range :math:`[0, 1]`
-    such that the elements along :code:`axis` sum to :math:`1`.
-
-    .. math ::
-      \mathrm{softmax}(x) = \frac{\exp(x_i)}{\sum_j \exp(x_j)}
-
-    Args:
-      x : input array
-      axis: the axis or axes along which the softmax should be computed. The
-        softmax output summed across these dimensions should sum to :math:`1`.
-        Either an integer or a tuple of integers.
-      where: Elements to include in the :code:`softmax`.
-      initial: The minimum value used to shift the input array. Must be present
-        when :code:`where` is not None.
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`log_softmax`
-    """
-    return _keep_unit(jax.nn.softmax, x, axis=axis, where=where)
-
-
-def standardize(x: ArrayLike,
-                axis: int | tuple[int, ...] | None = -1,
-                variance: ArrayLike | None = None,
-                epsilon: ArrayLike = 1e-5,
-                where: ArrayLike | None = None) -> Union[jax.Array, u.Quantity]:
-    r"""Normalizes an array by subtracting ``mean`` and dividing by :math:`\sqrt{\mathrm{variance}}`."""
-    return _keep_unit(jax.nn.standardize, x, axis=axis, where=where, variance=variance, epsilon=epsilon)
-
-
-def one_hot(x: Any,
-            num_classes: int, *,
-            dtype: Any = jax.numpy.float_,
-            axis: Union[int, Sequence[int]] = -1) -> Union[jax.Array, u.Quantity]:
-    """One-hot encodes the given indices.
-
-    Each index in the input ``x`` is encoded as a vector of zeros of length
-    ``num_classes`` with the element at ``index`` set to one::
-
-      >>> one_hot(jnp.array([0, 1, 2]), 3)
-      Array([[1., 0., 0.],
-             [0., 1., 0.],
-             [0., 0., 1.]], dtype=float32)
-
-    Indices outside the range [0, num_classes) will be encoded as zeros::
-
-      >>> one_hot(jnp.array([-1, 3]), 3)
-      Array([[0., 0., 0.],
-             [0., 0., 0.]], dtype=float32)
-
-    Args:
-      x: A tensor of indices.
-      num_classes: Number of classes in the one-hot dimension.
-      dtype: optional, a float dtype for the returned values (default :obj:`jnp.float_`).
-      axis: the axis or axes along which the function should be
-        computed.
-    """
-    return _keep_unit(jax.nn.one_hot, x, axis=axis, num_classes=num_classes, dtype=dtype)
-
-
-def relu6(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Rectified Linear Unit 6 activation function.
-
-    Computes the element-wise function
-
-    .. math::
-      \mathrm{relu6}(x) = \min(\max(x, 0), 6)
-
-    except under differentiation, we take:
-
-    .. math::
-      \nabla \mathrm{relu}(0) = 0
-
-    and
-
-    .. math::
-      \nabla \mathrm{relu}(6) = 0
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`relu`
-    """
-    return _keep_unit(jax.nn.relu6, x)
-
-
-def hard_sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Hard Sigmoid activation function.
-
-    Computes the element-wise function
-
-    .. math::
-      \mathrm{hard\_sigmoid}(x) = \frac{\mathrm{relu6}(x + 3)}{6}
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`relu6`
-    """
-    return _keep_unit(jax.nn.hard_sigmoid, x)
-
-
-def hard_silu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Hard SiLU (swish) activation function
-
-    Computes the element-wise function
-
-    .. math::
-      \mathrm{hard\_silu}(x) = x \cdot \mathrm{hard\_sigmoid}(x)
-
-    Both :func:`hard_silu` and :func:`hard_swish` are aliases for the same
-    function.
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`hard_sigmoid`
-    """
-    return _keep_unit(jax.nn.hard_silu, x)
-
-    return jax.nn.hard_silu(x)
-
-
-hard_swish = hard_silu
-
-
-def sparse_plus(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Sparse plus function.
-
-    Computes the function:
-
-    .. math::
-
-      \mathrm{sparse\_plus}(x) = \begin{cases}
-        0, & x \leq -1\\
-        \frac{1}{4}(x+1)^2, & -1 < x < 1 \\
-        x, & 1 \leq x
-      \end{cases}
-
-    This is the twin function of the softplus activation ensuring a zero output
-    for inputs less than -1 and a linear output for inputs greater than 1,
-    while remaining smooth, convex, monotonic by an adequate definition between
-    -1 and 1.
-
-    Args:
-      x: input (float)
-    """
-    return _keep_unit(jax.nn.sparse_plus, x)
-
-
-def sparse_sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
-    r"""Sparse sigmoid activation function.
-
-    Computes the function:
-
-    .. math::
-
-      \mathrm{sparse\_sigmoid}(x) = \begin{cases}
-        0, & x \leq -1\\
-        \frac{1}{2}(x+1), & -1 < x < 1 \\
-        1, & 1 \leq x
-      \end{cases}
-
-    This is the twin function of the ``sigmoid`` activation ensuring a zero output
-    for inputs less than -1, a 1 output for inputs greater than 1, and a linear
-    output for inputs between -1 and 1. It is the derivative of ``sparse_plus``.
-
-    For more information, see `Learning with Fenchel-Young Losses (section 6.2)
-    <https://arxiv.org/abs/1901.02324>`_.
-
-    Args:
-      x : input array
-
-    Returns:
-      An array.
-
-    See also:
-      :func:`sigmoid`
-    """
-    return _keep_unit(jax.nn.sparse_sigmoid, x)
+# Copyright 2024 BDP Ecosystem Limited. All Rights Reserved.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#     http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+# ==============================================================================
+
+
+"""
+Shared neural network activations and other functions.
+"""
+
+from typing import Any, Union, Sequence
+
+import brainunit as u
+import jax
+from jax.scipy.special import logsumexp
+
+from brainstate import random
+from brainstate.typing import ArrayLike
+
+__all__ = [
+    "tanh",
+    "relu",
+    "squareplus",
+    "softplus",
+    "soft_sign",
+    "sigmoid",
+    "silu",
+    "swish",
+    "log_sigmoid",
+    "elu",
+    "leaky_relu",
+    "hard_tanh",
+    "celu",
+    "selu",
+    "gelu",
+    "glu",
+    "logsumexp",
+    "log_softmax",
+    "softmax",
+    "standardize",
+    "one_hot",
+    "relu6",
+    "hard_sigmoid",
+    "hard_silu",
+    "hard_swish",
+    'hard_shrink',
+    'rrelu',
+    'mish',
+    'soft_shrink',
+    'prelu',
+    'tanh_shrink',
+    'softmin',
+    'sparse_plus',
+    'sparse_sigmoid',
+]
+
+
+def tanh(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Hyperbolic tangent activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{tanh}(x) = \frac{e^x - e^{-x}}{e^x + e^{-x}}
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+    """
+    return u.math.tanh(x)
+
+
+def softmin(x, axis=-1):
+    r"""
+    Applies the Softmin function to an n-dimensional input Tensor
+    rescaling them so that the elements of the n-dimensional output Tensor
+    lie in the range `[0, 1]` and sum to 1.
+
+    Softmin is defined as:
+
+    .. math::
+        \text{Softmin}(x_{i}) = \frac{\exp(-x_i)}{\sum_j \exp(-x_j)}
+
+    Shape:
+        - Input: :math:`(*)` where `*` means, any number of additional
+          dimensions
+        - Output: :math:`(*)`, same shape as the input
+
+    Args:
+        axis (int): A dimension along which Softmin will be computed (so every slice
+            along dim will sum to 1).
+    """
+    unnormalized = u.math.exp(-x)
+    return unnormalized / unnormalized.sum(axis, keepdims=True)
+
+
+def tanh_shrink(x):
+    r"""
+    Applies the element-wise function:
+
+    .. math::
+        \text{Tanhshrink}(x) = x - \tanh(x)
+    """
+    return x - u.math.tanh(x)
+
+
+def prelu(x, a=0.25):
+    r"""
+    Applies the element-wise function:
+
+    .. math::
+        \text{PReLU}(x) = \max(0,x) + a * \min(0,x)
+
+    or
+
+    .. math::
+        \text{PReLU}(x) =
+        \begin{cases}
+        x, & \text{ if } x \geq 0 \\
+        ax, & \text{ otherwise }
+        \end{cases}
+
+    Here :math:`a` is a learnable parameter. When called without arguments, `nn.PReLU()` uses a single
+    parameter :math:`a` across all input channels. If called with `nn.PReLU(nChannels)`,
+    a separate :math:`a` is used for each input channel.
+    """
+    return u.math.where(x >= 0., x, a * x)
+
+
+def soft_shrink(x, lambd=0.5):
+    r"""
+    Applies the soft shrinkage function elementwise:
+
+    .. math::
+        \text{SoftShrinkage}(x) =
+        \begin{cases}
+        x - \lambda, & \text{ if } x > \lambda \\
+        x + \lambda, & \text{ if } x < -\lambda \\
+        0, & \text{ otherwise }
+        \end{cases}
+
+    Args:
+        lambd: the :math:`\lambda` (must be no less than zero) value for the Softshrink formulation. Default: 0.5
+
+    Shape:
+        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
+        - Output: :math:`(*)`, same shape as the input.
+    """
+    return u.math.where(x > lambd,
+                        x - lambd,
+                        u.math.where(x < -lambd,
+                                     x + lambd,
+                                     u.Quantity(0., unit=u.get_unit(lambd))))
+
+
+def mish(x):
+    r"""Applies the Mish function, element-wise.
+
+    Mish: A Self Regularized Non-Monotonic Neural Activation Function.
+
+    .. math::
+        \text{Mish}(x) = x * \text{Tanh}(\text{Softplus}(x))
+
+    .. note::
+        See `Mish: A Self Regularized Non-Monotonic Neural Activation Function <https://arxiv.org/abs/1908.08681>`_
+
+    Shape:
+        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
+        - Output: :math:`(*)`, same shape as the input.
+    """
+    return x * u.math.tanh(softplus(x))
+
+
+def rrelu(x, lower=0.125, upper=0.3333333333333333):
+    r"""Applies the randomized leaky rectified liner unit function, element-wise,
+    as described in the paper:
+
+    `Empirical Evaluation of Rectified Activations in Convolutional Network`_.
+
+    The function is defined as:
+
+    .. math::
+        \text{RReLU}(x) =
+        \begin{cases}
+            x & \text{if } x \geq 0 \\
+            ax & \text{ otherwise }
+        \end{cases}
+
+    where :math:`a` is randomly sampled from uniform distribution
+    :math:`\mathcal{U}(\text{lower}, \text{upper})`.
+
+     See: https://arxiv.org/pdf/1505.00853.pdf
+
+    Args:
+        lower: lower bound of the uniform distribution. Default: :math:`\frac{1}{8}`
+        upper: upper bound of the uniform distribution. Default: :math:`\frac{1}{3}`
+
+    Shape:
+        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
+        - Output: :math:`(*)`, same shape as the input.
+
+    .. _`Empirical Evaluation of Rectified Activations in Convolutional Network`:
+        https://arxiv.org/abs/1505.00853
+    """
+    a = random.uniform(lower, upper, size=u.math.shape(x), dtype=x.dtype)
+    return u.math.where(u.get_mantissa(x) >= 0., x, a * x)
+
+
+def hard_shrink(x, lambd=0.5):
+    r"""Applies the Hard Shrinkage (Hardshrink) function element-wise.
+
+    Hardshrink is defined as:
+
+    .. math::
+        \text{HardShrink}(x) =
+        \begin{cases}
+        x, & \text{ if } x > \lambda \\
+        x, & \text{ if } x < -\lambda \\
+        0, & \text{ otherwise }
+        \end{cases}
+
+    Args:
+        lambd: the :math:`\lambda` value for the Hardshrink formulation. Default: 0.5
+
+    Shape:
+        - Input: :math:`(*)`, where :math:`*` means any number of dimensions.
+        - Output: :math:`(*)`, same shape as the input.
+
+    """
+    return u.math.where(x > lambd,
+                        x,
+                        u.math.where(x < -lambd,
+                                     x,
+                                     u.Quantity(0., unit=u.get_unit(x))))
+
+
+def relu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Rectified linear unit activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{relu}(x) = \max(x, 0)
+
+    except under differentiation, we take:
+
+    .. math::
+      \nabla \mathrm{relu}(0) = 0
+
+    For more information see
+    `Numerical influence of ReLU’(0) on backpropagation
+    <https://openreview.net/forum?id=urrcVI-_jRm>`_.
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    Example:
+      >>> jax.nn.relu(jax.numpy.array([-2., -1., -0.5, 0, 0.5, 1., 2.]))
+      Array([0. , 0. , 0. , 0. , 0.5, 1. , 2. ], dtype=float32)
+
+    See also:
+      :func:`relu6`
+
+    """
+    return u.math.relu(x)
+
+
+def squareplus(x: ArrayLike, b: ArrayLike = 4) -> Union[jax.Array, u.Quantity]:
+    r"""Squareplus activation function.
+
+    Computes the element-wise function
+
+    .. math::
+      \mathrm{squareplus}(x) = \frac{x + \sqrt{x^2 + b}}{2}
+
+    as described in https://arxiv.org/abs/2112.11687.
+
+    Args:
+      x : input array
+      b : smoothness parameter
+    """
+    return u.math.squareplus(x, b=b)
+
+
+def softplus(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Softplus activation function.
+
+    Computes the element-wise function
+
+    .. math::
+      \mathrm{softplus}(x) = \log(1 + e^x)
+
+    Args:
+      x : input array
+    """
+    return u.math.softplus(x)
+
+
+def soft_sign(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Soft-sign activation function.
+
+    Computes the element-wise function
+
+    .. math::
+      \mathrm{soft\_sign}(x) = \frac{x}{|x| + 1}
+
+    Args:
+      x : input array
+    """
+    return u.math.soft_sign(x)
+
+
+def sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Sigmoid activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{sigmoid}(x) = \frac{1}{1 + e^{-x}}
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`log_sigmoid`
+
+    """
+    return u.math.sigmoid(x)
+
+
+def silu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""SiLU (a.k.a. swish) activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{silu}(x) = x \cdot \mathrm{sigmoid}(x) = \frac{x}{1 + e^{-x}}
+
+    :func:`swish` and :func:`silu` are both aliases for the same function.
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`sigmoid`
+    """
+    return u.math.silu(x)
+
+
+swish = silu
+
+
+def log_sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Log-sigmoid activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{log\_sigmoid}(x) = \log(\mathrm{sigmoid}(x)) = -\log(1 + e^{-x})
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`sigmoid`
+    """
+    return u.math.log_sigmoid(x)
+
+
+def elu(x: ArrayLike, alpha: ArrayLike = 1.0) -> Union[jax.Array, u.Quantity]:
+    r"""Exponential linear unit activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{elu}(x) = \begin{cases}
+        x, & x > 0\\
+        \alpha \left(\exp(x) - 1\right), & x \le 0
+      \end{cases}
+
+    Args:
+      x : input array
+      alpha : scalar or array of alpha values (default: 1.0)
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`selu`
+    """
+    return u.math.elu(x, alpha=alpha)
+
+
+def leaky_relu(x: ArrayLike, negative_slope: ArrayLike = 1e-2) -> Union[jax.Array, u.Quantity]:
+    r"""Leaky rectified linear unit activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{leaky\_relu}(x) = \begin{cases}
+        x, & x \ge 0\\
+        \alpha x, & x < 0
+      \end{cases}
+
+    where :math:`\alpha` = :code:`negative_slope`.
+
+    Args:
+      x : input array
+      negative_slope : array or scalar specifying the negative slope (default: 0.01)
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`relu`
+    """
+    return u.math.leaky_relu(x, negative_slope=negative_slope)
+
+
+def _hard_tanh(x, min_val=- 1.0, max_val=1.0):
+    return jax.numpy.where(x > max_val, max_val, jax.numpy.where(x < min_val, min_val, x))
+
+
+def hard_tanh(
+    x: ArrayLike,
+    min_val: float = - 1.0,
+    max_val: float = 1.0
+) -> Union[jax.Array, u.Quantity]:
+    r"""Hard :math:`\mathrm{tanh}` activation function.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{hard\_tanh}(x) = \begin{cases}
+        -1, & x < -1\\
+        x, & -1 \le x \le 1\\
+        1, & 1 < x
+      \end{cases}
+
+    Args:
+      x : input array
+      min_val: float. minimum value of the linear region range. Default: -1
+      max_val: float. maximum value of the linear region range. Default: 1
+
+    Returns:
+      An array.
+    """
+    x = u.Quantity(x)
+    min_val = u.Quantity(min_val).to(x.unit).mantissa
+    max_val = u.Quantity(max_val).to(x.unit).mantissa
+    return u.maybe_decimal(_hard_tanh(x.mantissa, min_val=min_val, max_val=max_val) * x.unit)
+
+
+def celu(x: ArrayLike, alpha: ArrayLike = 1.0) -> Union[jax.Array, u.Quantity]:
+    r"""Continuously-differentiable exponential linear unit activation.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{celu}(x) = \begin{cases}
+        x, & x > 0\\
+        \alpha \left(\exp(\frac{x}{\alpha}) - 1\right), & x \le 0
+      \end{cases}
+
+    For more information, see
+    `Continuously Differentiable Exponential Linear Units
+    <https://arxiv.org/pdf/1704.07483.pdf>`_.
+
+    Args:
+      x : input array
+      alpha : array or scalar (default: 1.0)
+
+    Returns:
+      An array.
+    """
+    return u.math.celu(x, alpha=alpha)
+
+
+def selu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Scaled exponential linear unit activation.
+
+    Computes the element-wise function:
+
+    .. math::
+      \mathrm{selu}(x) = \lambda \begin{cases}
+        x, & x > 0\\
+        \alpha e^x - \alpha, & x \le 0
+      \end{cases}
+
+    where :math:`\lambda = 1.0507009873554804934193349852946` and
+    :math:`\alpha = 1.6732632423543772848170429916717`.
+
+    For more information, see
+    `Self-Normalizing Neural Networks
+    <https://papers.nips.cc/paper/6698-self-normalizing-neural-networks.pdf>`_.
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`elu`
+    """
+    return u.math.selu(x)
+
+
+def gelu(x: ArrayLike, approximate: bool = True) -> Union[jax.Array, u.Quantity]:
+    r"""Gaussian error linear unit activation function.
+
+    If ``approximate=False``, computes the element-wise function:
+
+    .. math::
+      \mathrm{gelu}(x) = \frac{x}{2} \left(1 + \mathrm{erf} \left(
+        \frac{x}{\sqrt{2}} \right) \right)
+
+    If ``approximate=True``, uses the approximate formulation of GELU:
+
+    .. math::
+      \mathrm{gelu}(x) = \frac{x}{2} \left(1 + \mathrm{tanh} \left(
+        \sqrt{\frac{2}{\pi}} \left(x + 0.044715 x^3 \right) \right) \right)
+
+    For more information, see `Gaussian Error Linear Units (GELUs)
+    <https://arxiv.org/abs/1606.08415>`_, section 2.
+
+    Args:
+      x : input array
+      approximate: whether to use the approximate or exact formulation.
+    """
+    return u.math.gelu(x, approximate=approximate)
+
+
+def glu(x: ArrayLike, axis: int = -1) -> Union[jax.Array, u.Quantity]:
+    r"""Gated linear unit activation function.
+
+    Computes the function:
+
+    .. math::
+      \mathrm{glu}(x) =  x\left[\ldots, 0:\frac{n}{2}, \ldots\right] \cdot
+        \mathrm{sigmoid} \left( x\left[\ldots, \frac{n}{2}:n, \ldots\right]
+          \right)
+
+    where the array is split into two along ``axis``. The size of the ``axis``
+    dimension must be divisible by two.
+
+    Args:
+      x : input array
+      axis: the axis along which the split should be computed (default: -1)
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`sigmoid`
+    """
+    return u.math.glu(x, axis=axis)
+
+
+def log_softmax(x: ArrayLike,
+                axis: int | tuple[int, ...] | None = -1,
+                where: ArrayLike | None = None) -> Union[jax.Array, u.Quantity]:
+    r"""Log-Softmax function.
+
+    Computes the logarithm of the :code:`softmax` function, which rescales
+    elements to the range :math:`[-\infty, 0)`.
+
+    .. math ::
+      \mathrm{log\_softmax}(x)_i = \log \left( \frac{\exp(x_i)}{\sum_j \exp(x_j)}
+      \right)
+
+    Args:
+      x : input array
+      axis: the axis or axes along which the :code:`log_softmax` should be
+        computed. Either an integer or a tuple of integers.
+      where: Elements to include in the :code:`log_softmax`.
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`softmax`
+    """
+    return jax.nn.log_softmax(x, axis=axis, where=where)
+
+
+def softmax(x: ArrayLike,
+            axis: int | tuple[int, ...] | None = -1,
+            where: ArrayLike | None = None) -> Union[jax.Array, u.Quantity]:
+    r"""Softmax function.
+
+    Computes the function which rescales elements to the range :math:`[0, 1]`
+    such that the elements along :code:`axis` sum to :math:`1`.
+
+    .. math ::
+      \mathrm{softmax}(x) = \frac{\exp(x_i)}{\sum_j \exp(x_j)}
+
+    Args:
+      x : input array
+      axis: the axis or axes along which the softmax should be computed. The
+        softmax output summed across these dimensions should sum to :math:`1`.
+        Either an integer or a tuple of integers.
+      where: Elements to include in the :code:`softmax`.
+      initial: The minimum value used to shift the input array. Must be present
+        when :code:`where` is not None.
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`log_softmax`
+    """
+    return jax.nn.softmax(x, axis=axis, where=where)
+
+
+def standardize(x: ArrayLike,
+                axis: int | tuple[int, ...] | None = -1,
+                variance: ArrayLike | None = None,
+                epsilon: ArrayLike = 1e-5,
+                where: ArrayLike | None = None) -> Union[jax.Array, u.Quantity]:
+    r"""Normalizes an array by subtracting ``mean`` and dividing by :math:`\sqrt{\mathrm{variance}}`."""
+    return jax.nn.standardize(x, axis=axis, where=where, variance=variance, epsilon=epsilon)
+
+
+def one_hot(x: Any,
+            num_classes: int, *,
+            dtype: Any = jax.numpy.float_,
+            axis: Union[int, Sequence[int]] = -1) -> Union[jax.Array, u.Quantity]:
+    """One-hot encodes the given indices.
+
+    Each index in the input ``x`` is encoded as a vector of zeros of length
+    ``num_classes`` with the element at ``index`` set to one::
+
+      >>> one_hot(jnp.array([0, 1, 2]), 3)
+      Array([[1., 0., 0.],
+             [0., 1., 0.],
+             [0., 0., 1.]], dtype=float32)
+
+    Indices outside the range [0, num_classes) will be encoded as zeros::
+
+      >>> one_hot(jnp.array([-1, 3]), 3)
+      Array([[0., 0., 0.],
+             [0., 0., 0.]], dtype=float32)
+
+    Args:
+      x: A tensor of indices.
+      num_classes: Number of classes in the one-hot dimension.
+      dtype: optional, a float dtype for the returned values (default :obj:`jnp.float_`).
+      axis: the axis or axes along which the function should be
+        computed.
+    """
+    return jax.nn.one_hot(x, axis=axis, num_classes=num_classes, dtype=dtype)
+
+
+def relu6(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Rectified Linear Unit 6 activation function.
+
+    Computes the element-wise function
+
+    .. math::
+      \mathrm{relu6}(x) = \min(\max(x, 0), 6)
+
+    except under differentiation, we take:
+
+    .. math::
+      \nabla \mathrm{relu}(0) = 0
+
+    and
+
+    .. math::
+      \nabla \mathrm{relu}(6) = 0
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`relu`
+    """
+    return u.math.relu6(x)
+
+
+def hard_sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Hard Sigmoid activation function.
+
+    Computes the element-wise function
+
+    .. math::
+      \mathrm{hard\_sigmoid}(x) = \frac{\mathrm{relu6}(x + 3)}{6}
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`relu6`
+    """
+    return u.math.hard_sigmoid(x)
+
+
+def hard_silu(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Hard SiLU (swish) activation function
+
+    Computes the element-wise function
+
+    .. math::
+      \mathrm{hard\_silu}(x) = x \cdot \mathrm{hard\_sigmoid}(x)
+
+    Both :func:`hard_silu` and :func:`hard_swish` are aliases for the same
+    function.
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`hard_sigmoid`
+    """
+    return u.math.hard_silu(x)
+
+
+hard_swish = hard_silu
+
+
+def sparse_plus(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Sparse plus function.
+
+    Computes the function:
+
+    .. math::
+
+      \mathrm{sparse\_plus}(x) = \begin{cases}
+        0, & x \leq -1\\
+        \frac{1}{4}(x+1)^2, & -1 < x < 1 \\
+        x, & 1 \leq x
+      \end{cases}
+
+    This is the twin function of the softplus activation ensuring a zero output
+    for inputs less than -1 and a linear output for inputs greater than 1,
+    while remaining smooth, convex, monotonic by an adequate definition between
+    -1 and 1.
+
+    Args:
+      x: input (float)
+    """
+    return u.math.sparse_plus(x)
+
+
+def sparse_sigmoid(x: ArrayLike) -> Union[jax.Array, u.Quantity]:
+    r"""Sparse sigmoid activation function.
+
+    Computes the function:
+
+    .. math::
+
+      \mathrm{sparse\_sigmoid}(x) = \begin{cases}
+        0, & x \leq -1\\
+        \frac{1}{2}(x+1), & -1 < x < 1 \\
+        1, & 1 \leq x
+      \end{cases}
+
+    This is the twin function of the ``sigmoid`` activation ensuring a zero output
+    for inputs less than -1, a 1 output for inputs greater than 1, and a linear
+    output for inputs between -1 and 1. It is the derivative of ``sparse_plus``.
+
+    For more information, see `Learning with Fenchel-Young Losses (section 6.2)
+    <https://arxiv.org/abs/1901.02324>`_.
+
+    Args:
+      x : input array
+
+    Returns:
+      An array.
+
+    See also:
+      :func:`sigmoid`
+    """
+    return u.math.sparse_sigmoid(x)