Geometric Inversion, 3D Cubes
geometric inversion of cubes on 3d grid

geoInv = ((With[{x662 = #.#}, If[ x662 < 0.00000001, #, #/x662] ]) &); cubes = Table[Cube[{x, y, z}, 0.8], {x, -3, 3}, {y, 0, 3}, {z, -3, 3}] /. Cube[{(0)..}, _] -> Nothing; Graphics3D[{Cyan, cubes /. x_Cube -> CanonicalizePolyhedron[x] /. Polyhedron[pts_, g_] :> Polyhedron[geoInv /@ pts, g]}, Boxed -> False, PlotRange -> 0.4]