def METHOD_NAME(xp, scp, lmb, data):
    result = scp.stats.boxcox_llf(lmb, data)
    if data.ndim == 1:
        if data.dtype.kind == 'c':
            assert result.dtype == xp.complex128
        else:
            assert result.dtype == xp.float64
    else:
        if data.dtype.kind in 'cf':
            assert result.dtype == data.dtype
        elif lmb == 0:
            for dtype1 in [xp.float16, xp.float32, xp.float64]:
                if xp.can_cast(data.dtype, dtype1):
                    break
            assert result.dtype == dtype1
        else:
            assert result.dtype == xp.float64
    return result, _dtype(data.dtype, xp)