def METHOD_NAME(get_fn):
    traj = md.load(get_fn("tip3p_300K_1ATM.xtc"), top=get_fn("tip3p_300K_1ATM.pdb"))
    charges = np.tile(tip3p_charges, traj.n_residues)
    moments0 = md.geometry.dipole_moments(traj, charges)
    # Now we screw up the molecule wholeness referencing all distances relative to atom zero.
    # E.g. [[0, 0], [0, 1], [0, 2], [0, 3]]...
    atom_indices = np.array([np.zeros(traj.n_atoms), np.arange(traj.n_atoms)], dtype='int32').T
    # Define coordinates relative to atom 0, PBC corrected.
    xyz = md.compute_displacements(traj, atom_indices, periodic=True)
    traj.xyz = xyz
    moments1 = md.geometry.dipole_moments(traj, charges)
    eq(moments0, moments1, decimal=4)