def METHOD_NAME(self):
    """
    Test :func:`colour.adaptation.zhai2018.chromatic_adaptation_Zhai2018`
    definition n-dimensional arrays support.
    """
    XYZ_b = np.array([48.900, 43.620, 6.250])
    XYZ_wb = np.array([109.850, 100, 35.585])
    XYZ_wd = np.array([95.047, 100, 108.883])
    D_b = 0.9407
    D_d = 0.9800
    XYZ_d = chromatic_adaptation_Zhai2018(XYZ_b, XYZ_wb, XYZ_wd, D_b, D_d)
    XYZ_b = np.tile(XYZ_b, (6, 1))
    XYZ_d = np.tile(XYZ_d, (6, 1))
    np.testing.assert_array_almost_equal(
        chromatic_adaptation_Zhai2018(XYZ_b, XYZ_wb, XYZ_wd, D_b, D_d),
        XYZ_d,
        decimal=7,
    )
    XYZ_wb = np.tile(XYZ_wb, (6, 1))
    XYZ_wd = np.tile(XYZ_wd, (6, 1))
    D_b = np.tile(D_b, (6, 1))
    D_d = np.tile(D_d, (6, 1))
    np.testing.assert_array_almost_equal(
        chromatic_adaptation_Zhai2018(XYZ_b, XYZ_wb, XYZ_wd, D_b, D_d),
        XYZ_d,
        decimal=7,
    )
    XYZ_b = np.reshape(XYZ_b, (2, 3, 3))
    XYZ_wb = np.reshape(XYZ_wb, (2, 3, 3))
    XYZ_wd = np.reshape(XYZ_wd, (2, 3, 3))
    D_b = np.reshape(D_b, (2, 3, 1))
    D_d = np.reshape(D_d, (2, 3, 1))
    XYZ_d = np.reshape(XYZ_d, (2, 3, 3))
    np.testing.assert_array_almost_equal(
        chromatic_adaptation_Zhai2018(XYZ_b, XYZ_wb, XYZ_wd, D_b, D_d),
        XYZ_d,
        decimal=7,
    )