def METHOD_NAME(XYZ: NDArrayFloat, Ljg: NDArrayFloat) -> NDArrayFloat:
    """Error function."""
    # Error must be computed in "reference" domain and range.
    with domain_range_scale("ignore"):
        error = np.linalg.norm(XYZ_to_OSA_UCS(XYZ) - Ljg)
    return as_float(error)