Optimize visualization of non-Euclidean geometries and domains

The following never finishes for me, the julia process just keeps spinning at 100%:

julia> import GeoArtifacts

julia> using GeoTables

julia> using Meshes
       #using Plots

julia> using CairoMakie

julia> BRA = GeoArtifacts.GADM.get("BRA", depth = 0)
               1×3 GeoTable over 1 GeometrySet
┌─────────────┬─────────────┬───────────────────────────────┐
│    GID_0    │   COUNTRY   │           geometry            │
│ Categorical │ Categorical │         MultiPolygon          │
│  [NoUnits]  │  [NoUnits]  │ 🖈 GeodeticLatLon{WGS84Latest} │
├─────────────┼─────────────┼───────────────────────────────┤
│     BRA     │   Brazil    │     Multi(1007×PolyArea)      │
└─────────────┴─────────────┴───────────────────────────────┘

julia> BRA_domain = domain(BRA)
1 GeometrySet
└─ Multi(1007×PolyArea)

julia> viz(BRA_domain)

(jl_gIjdF8) pkg> st
Status `/tmp/jl_gIjdF8/Project.toml`
  [13f3f980] CairoMakie v0.15.11
  [2a3f1dca] GeoArtifacts v1.3.4
  [e502b557] GeoTables v1.28.5
  [eacbb407] Meshes v0.57.12

julia> VERSION
v"1.11.9"