Smoothing via Weaverbird’s Catmull-Clark Subdivision, 3 subdivisions per face (L=3), naked edges smooth (S=1).

Relaxation via Zombie solver, threshold 0.00001, tolerance 0.0001
edge length factor 0.7, strength 5
Anchor points: naked boundary edges (783 points)