Here are some schemes:

- octree (current thrust)
- tetrahedron
- Sliver
- Easy Sliver (2nd best just now)
- Balls
- Solve Poisson’s equation
- FFT

Spherical Harmonics may facilitate some of these schemes.

As we move mass on the surface to make the potential at the surface more uniform, the vector angular momentum will change.
We compute the new covariance tensor and modify the angular velocity to preserve the old vector angular momentum.
This preserves the invariant that the spin axis remains the principle axis of the tensor inertia; the axis of rotation is aligned with the angular momentum.

I can use Jacoby’s routine to compute an eigenvector for the inertia tensor of the blob, and spin about that.
The covariance tensor has the same eigenvectors as the inertia tensor.

notes on that code, Algebra with computer help, Initial configuration

earlier notes

Solution of new spin axis: Compute biggest eigenvector after moving mass.