def XYZ100_to_sRGB1_linear(XYZ100):
    XYZ100 = np.asarray(XYZ100, dtype=float)
    # this is broadcasting matrix * array-of-vectors, where the vector is the
    # last dim
    RGB_linear = np.einsum("...ij,...j->...i", XYZ100_to_sRGB1_matrix, XYZ100 / 100)
    return RGB_linear
Convert XYZ to linear sRGB, where XYZ is normalized so that reference
dataset/ETHPy150Open njsmith/colorspacious/colorspacious/basics.py/XYZ100_to_sRGB1_linear