The topology of the elements of a metamaterial can be engineered from its desired electromagnetic constitutive tensor using an inverse group theory method. Therefore, given a desired electromagnetic response and a generic metamaterial elemental design, group theory is applied to predict the various ways that the element can be arranged in three dimensions to produce the desired functionality. An optimizer can then be applied to an electromagnetic modeling tool to fine tune the values of the electromagnetic properties of the resulting metamaterial topology.