def METHOD_NAME(L):
    """
    Define *Recommendation ITU-R BT.1361* extended color gamut system
    opto-electronic transfer function (OETF).
    Parameters
    ----------
    L
        Scene *Luminance* :math:`L`.
    Returns
    -------
    :class:`numpy.ndarray`
        Corresponding non-linear primary signal :math:`E'`.
    Notes
    -----
    +------------+-----------------------+---------------+
    | **Domain** | **Scale - Reference** | **Scale - 1** |
    +============+=======================+===============+
    | ``L``      | [0, 1]                | [0, 1]        |
    +------------+-----------------------+---------------+
    +------------+-----------------------+-------------------+
    | **Range**  | **Scale - Reference** | **Scale - 1**     |
    +============+=======================+===================+
    | ``E_p'``   | [0, 1]                | [0, 1]            |
    +------------+-----------------------+-------------------+
    References
    ----------
    :cite:`InternationalTelecommunicationUnion1998`
    Examples
    --------
    >>> oetf_BT1361(0.18)  # doctest: +ELLIPSIS
    0.4090077288641...
    >>> oetf_BT1361(-0.25)  # doctest: +ELLIPSIS
    -0.25
    >>> oetf_BT1361(1.33)  # doctest: +ELLIPSIS
    1.1504846663972...
    """
    L = to_domain_1(L)
    with domain_range_scale("ignore"):
        E_p = np.where(
            L >= 0,
            oetf_BT709(L),
            np.where(
                -0.0045 >= L,
                # L in [-0.25, -0.0045] range
                -(1.099 * spow(-4 * L, 0.45) - 0.099) / 4,
                # L in [-0.0045, 0] range
                4.500 * L,
            ),
        )
    return as_float(from_range_1(E_p))