def METHOD_NAME(
    J: ArrayLike,
    in_int: bool = False,
    constants: Structure = CONSTANTS_DICOMGSDF,
) -> NDArrayFloat:
    """
    Define the *DICOM - Grayscale Standard Display Function* electro-optical
    transfer function (EOTF).
    Parameters
    ----------
    J
        Just-Noticeable Difference (JND) Index, :math:`j`.
    in_int
        Whether to treat the input value as int code value or float
        equivalent of a code value at a given bit-depth.
    constants
        *DICOM - Grayscale Standard Display Function* constants.
    Returns
    -------
    :class:`numpy.ndarray`
        Corresponding *luminance* :math:`L`.
    Notes
    -----
    +------------+-----------------------+---------------+
    | **Domain** | **Scale - Reference** | **Scale - 1** |
    +============+=======================+===============+
    | ``J``      | [0, 1]                | [0, 1]        |
    +------------+-----------------------+---------------+
    +------------+-----------------------+---------------+
    | **Range**  | **Scale - Reference** | **Scale - 1** |
    +============+=======================+===============+
    | ``L``      | [0, 1]                | [0, 1]        |
    +------------+-----------------------+---------------+
    References
    ----------
    :cite:`NationalElectricalManufacturersAssociation2004b`
    Examples
    --------
    >>> eotf_DICOMGSDF(0.500486263438448)  # doctest: +ELLIPSIS
    130.0628647...
    >>> eotf_DICOMGSDF(512, in_int=True)  # doctest: +ELLIPSIS
    130.0652840...
    """
    J = to_domain_1(J)
    if not in_int:
        J = J * 1023
    a = constants.a
    b = constants.b
    c = constants.c
    d = constants.d
    e = constants.e
    f = constants.f
    g = constants.g
    h = constants.h
    k = constants.k
    m = constants.m
    J_ln = np.log(J)
    J_ln2 = J_ln**2
    J_ln3 = J_ln**3
    J_ln4 = J_ln**4
    J_ln5 = J_ln**5
    L = (a + c * J_ln + e * J_ln2 + g * J_ln3 + m * J_ln4) / (
        1 + b * J_ln + d * J_ln2 + f * J_ln3 + h * J_ln4 + k * J_ln5
    )
    L = 10**L
    return as_float(from_range_1(L))