As we discussed on the call, this would enable making any spherical potential a triaxial model (in the potential). API to be worked out, but the idea would be to not have to implement flattened and triaxial versions of all potential models. Instead, we want something like (names of arguments and wrappers TBD):

pot = Triaxial(HernquistPotential(...), axis_ratio1=..., axis_ratio2=...)
pot2 = Flattened(HernquistPotential(...), axis_ratio=...)