def METHOD_NAME():
    """Test dev_head_t computation for Artemis123."""
    # test a random selected point
    raw = read_raw_artemis123(short_hpi_1kz_fname, preload=True, add_head_trans=False)
    meg_picks = pick_types(raw.info, meg=True, eeg=False)
    # checked against matlab reader.
    assert_allclose(raw[meg_picks[12]][0][0][123], 1.08239606023e-11)
    dev_head_t_1 = np.array(
        [
            [9.713e-01, 2.340e-01, -4.164e-02, 1.302e-04],
            [-2.371e-01, 9.664e-01, -9.890e-02, 1.977e-03],
            [1.710e-02, 1.059e-01, 9.942e-01, -8.159e-03],
            [0.0, 0.0, 0.0, 1.0],
        ]
    )
    dev_head_t_2 = np.array(
        [
            [9.890e-01, 1.475e-01, -8.090e-03, 4.997e-04],
            [-1.476e-01, 9.846e-01, -9.389e-02, 1.962e-03],
            [-5.888e-03, 9.406e-02, 9.955e-01, -1.610e-02],
            [0.0, 0.0, 0.0, 1.0],
        ]
    )
    expected_dev_hpi_rr = np.array(
        [
            [-0.01579644, 0.06527367, 0.00152648],
            [0.06666813, 0.0148956, 0.00545488],
            [-0.06699212, -0.01732376, 0.0112027],
        ]
    )
    # test with head loc no digitization
    raw = read_raw_artemis123(short_HPI_dip_fname, add_head_trans=True)
    _assert_trans(raw.info["dev_head_t"]["trans"], dev_head_t_1)
    assert_equal(raw.info["sfreq"], 5000.0)
    # test with head loc and digitization
    with pytest.warns(RuntimeWarning, match="Large difference"):
        raw = read_raw_artemis123(
            short_HPI_dip_fname, add_head_trans=True, pos_fname=dig_fname
        )
    _assert_trans(raw.info["dev_head_t"]["trans"], dev_head_t_1)
    # test cHPI localization..
    dev_hpi_rr = np.array(
        [p["r"] for p in raw.info["dig"] if p["coord_frame"] == FIFF.FIFFV_COORD_DEVICE]
    )
    # points should be within 0.1 mm (1e-4m) and within 1%
    assert_allclose(dev_hpi_rr, expected_dev_hpi_rr, atol=1e-4, rtol=0.01)
    # test 1kz hpi head loc (different freq)
    raw = read_raw_artemis123(short_hpi_1kz_fname, add_head_trans=True)
    _assert_trans(raw.info["dev_head_t"]["trans"], dev_head_t_2)
    assert_equal(raw.info["sfreq"], 1000.0)