Makie.jl 中混合精度浮点向量颜色循环的Bug分析

问题描述

在Makie.jl数据可视化库中，当用户同时绘制包含Float32和Float64类型数据的曲线时，出现了颜色循环不正确的问题。具体表现为：当多条曲线中同时包含不同精度的浮点数据时，图例中的颜色分配会出现异常，无法按照预期的循环顺序显示。

问题复现

该问题可以通过以下代码复现：

using CairoMakie
fig = Figure()
ax = Axis(fig[1,1])
xs = Float32.(10:100)  # x轴使用Float32数据
# 绘制三条曲线，其中两条使用Float32，一条使用Float64
lines!(xs, cumsum(randn(Float32, length(xs))); label = "1")
lines!(xs, cumsum(randn(Float32, length(xs))); label = "2")
lines!(xs, cumsum(randn(Float64, length(xs))); label = "3") 
axislegend()
fig

问题根源

经过分析，问题出在Makie.jl的Cycler类型实现上。当前版本中，颜色循环的索引计算过于严格地依赖于绘图类型的完整类型参数，包括数据的具体精度类型（Float32/Float64）。这导致当数据精度不同时，系统会将其视为完全不同的绘图类型，从而重新开始颜色循环。

技术分析

在Makie.jl的源码中，get_cycler_index!函数的实现如下：

function get_cycler_index!(c::Cycler, T::Type{Plot{Func, Args}}) where {Func, Args}
    P = Plot{Func}
    if !haskey(c.counters, P)
        return c.counters[P] = 1
    else
        return c.counters[P] += 1
    end
end

问题在于该函数原本考虑了完整的Plot类型参数（包括Func和Args），而实际上颜色循环应该只基于绘图的基本类型（如lines、scatter等），而不应考虑数据的具体精度类型。

解决方案

修复方案是修改get_cycler_index!函数，使其只关注Plot类型的第一个参数（Func），忽略数据类型的差异。修改后的实现如下：

function get_cycler_index!(c::Cycler, T::Type{Plot{Func, Args}}) where {Func, Args}
    P = Plot{Func}  # 只关注Func部分，忽略Args
    if !haskey(c.counters, P)
        return c.counters[P] = 1
    else
        return c.counters[P] += 1
    end
end

影响范围

该问题影响所有Makie后端（GLMakie、CairoMakie等），在绘制混合精度数据时会出现颜色循环异常。对于纯Float32或纯Float64数据集，则不会触发此问题。

总结

这个Bug揭示了在实现类型相关的功能时，需要仔细考虑类型系统的粒度。过度严格的类型匹配可能会导致不符合用户预期的行为。在数据可视化场景中，颜色循环这类功能通常应该基于更高层次的抽象（如绘图种类），而不是底层的数据实现细节（如数据精度）。

该修复已经合并到Makie.jl的主分支中，用户可以通过更新到最新版本来解决这个问题。这也提醒我们，在处理数值数据时，需要考虑不同精度类型之间的互操作性，以提供更一致的用户体验。