def METHOD_NAME(mol_h2o):
    basis = psi4.core.BasisSet.build(mol_h2o, 'orbital', 'cc-pvdz')
    mints = psi4.core.MintsHelper(basis)
    order = 6
    M = mints.ao_multipoles(order=order, origin=[0.0, 0.0, 0.0])
    with pytest.raises(Exception):
        # overlap integrals not accessible via multipole interface
        mints.ao_multipoles(order=0, origin=[0.0, 0.0, 0.0])
    with pytest.raises(Exception):
        # wrong origin specification
        mints.ao_multipoles(order=0, origin=[0.0, 0.0, 0.0, 0.0])
    Mnp = matlist_to_ndarray(M)
    # reference from l2
    dips = matlist_to_ndarray(mints.ao_dipole())
    quads = matlist_to_ndarray(mints.ao_quadrupole())
    np.testing.assert_allclose(dips, Mnp[:3], atol=1e-14)
    np.testing.assert_allclose(quads, Mnp[3:9], atol=1e-14)