Here's a better way. Start with a mesh with 24 vertices which results in the cube with all vertices and all edges cut with planes. Each vertex of the original

Here's a better way. Start with a mesh with 24 vertices which results in the cube with all vertices and all edges cut with planes. Each vertex of the original cube becomes a triangle. Then iteratively split edges in half, projecting the split points onto the desired surface. The algorithm only needs to split edges near the vertices of the original cube, so the split+project step is like this: float3 midpoint = ( v1 + v2 ) * 0.5; float3 rv = (midpoint - sphereCenter).normalize(); midpoint = sphereCenter + rv * sphereRadius; Unlike the OP's code, will result in uniform triangle density of the mesh. That's how people are usually generating good quality spherical meshes: start with icosahedron, then a few iterations of splitting edges in half + reprojecting back to the sphere.