I think that’s fixed in the meantime:
using GR ρ = LinRange(0., 7, 200) θ = LinRange(0., 2π, 360) funp(ρ,θ) = sin(2ρ) * cos(θ) polarheatmap(ρ,θ,funp.(ρ,θ')')