def _phonons_allclose(self, ph, fc333): ph333 = Phonopy(ph.unitcell, supercell_matrix=[3, 3, 3], primitive_matrix=ph.primitive_matrix) ph333.force_constants = fc333 ph333.nac_params = ph.nac_params comm_points = self._get_comm_points(ph) ph.run_qpoints(comm_points) ph333.run_qpoints(comm_points) np.testing.assert_allclose(ph.get_qpoints_dict()['frequencies'], ph333.get_qpoints_dict()['frequencies'], atol=1e-5)