Nevertheless, the FEM is flexible. Since changing inclusion shape does not affect the boundary conditions employed but only the complexity of the mesh used, arbitrary shapes could be used to model inclusions where appropriate. However, the fitting for inclusion and interfacial permittivity must be done separately for each type of shape used. Aside from ellipsoids, we discuss analysis done with rounded prism inclusions in Sec. S9.