An algorithmic GB decomposition looks quite similar, using a finite number of 𝑑GB ! -at approximately 6,000 per atom for this demonstration-bounded by Gs intersecting at and along the nodes and edges of a triangulated sphere mesh. Gaussian quadrature tetrahedral numerical integration is performed for each GB by decomposing it, first into convex polyhedra defined by the GB intersections with lower and upper 𝜌 isosurfaces, and again into tetrahedra.