def mesh_surface_area(verts, faces):
 actual_verts = verts[faces]
 a = (actual_verts[:, 0, :] - actual_verts[:, 1, :])
 b = (actual_verts[:, 0, :] - actual_verts[:, 2, :])
 del actual_verts
 return (((np.cross(a, b) ** 2).sum(axis=1) ** 0.5).sum() / 2.0)

def mesh_surface_area ( verts , faces ) : actual_verts = verts [ faces ] a = ( actual_verts [ : , 0 , : ] - actual_verts [ : , 1 , : ] ) b = ( actual_verts [ : , 0 , : ] - actual_verts [ : , 2 , : ] ) del actual_verts return ( ( ( np . cross ( a , b ) ** 2 ) . sum ( axis = 1 ) ** 0.5 ) . sum ( ) / 2.0 )
compute surface area .