def METHOD_NAME(self):
    """grouped pauli test"""
    pauli = 2 * (I ^ I) + (X ^ I) + 3 * (X ^ Y)
    grouped_pauli = PauliSumOp(pauli.primitive, grouping_type="TPB")
    converter = PauliBasisChange()
    cob = converter.convert(grouped_pauli)
    np.testing.assert_array_almost_equal(pauli.to_matrix(), cob.to_matrix())
    origin_x = reduce(np.logical_or, pauli.primitive.paulis.x)
    origin_z = reduce(np.logical_or, pauli.primitive.paulis.z)
    origin_pauli = Pauli((origin_z, origin_x))
    inst, dest = converter.get_cob_circuit(origin_pauli)
    self.assertEqual(str(dest), "ZZ")
    expected_inst = np.array(
        [
            [0.5, -0.5j, 0.5, -0.5j],
            [0.5, 0.5j, 0.5, 0.5j],
            [0.5, -0.5j, -0.5, 0.5j],
            [0.5, 0.5j, -0.5, -0.5j],
        ]
    )
    np.testing.assert_array_almost_equal(inst.to_matrix(), expected_inst)