An algorithm for the perception of constitutional symmetry in molecules (graphs) is presented, which partitions not only atoms (vertices) but also all pairwise relations among skeleton atoms into equivalence classes. The method works without canonical numbering, essentially by raising the connectivity matrix of the arbitrarily numbered molecule (graph) to its second, third, etc. power and evaluating the entries in these higher order matrices.