Numpy implementation of the polyhedral gravity algorithm from Werner - Scheeres (1996) [1]. Given a constant-density polyhedron, the Werner - Scheeres formulas allow for the efficient and accurate calculation of the grativational attraction and potential at a point which is exterior to the polyhedron.
[1] R. A. Werner, D. J. Scheeres, "Exterior Gravitation of a Polyhedron Derived and Compared with Harmonic and Mascon Gravitation Representations of Asteroid 4769 Castalia" (1996), Celestial Mechanics and Dynamical Astronomy 65(3):313-344 DOI: 10.1007/BF00053511