using DrWatson
@quickactivate
using Revise
using CairoMakie
includet("../src/main.jl")
include("../src/io.jl")

┌ Warning: Using an un-quoted Symbol on the LHS is deprecated. Write :tmax = ... instead.
└ @ DataFramesMeta /Users/zijin/.julia/packages/DataFramesMeta/Ga6VV/src/parsing.jl:387

split_results (generic function with 1 method)

For the same parameter, the left-hand transition matrix \(M_l\) has a 180-degree rotation about the right-hand matrix \(M_r\) : \(M_l[i,j] = M_r[n-i,n-j]\).

Test the different algorithms for the problem

dir = "test_alg"
df = get_result(; dir);

┌ Info: Scanning folder /Users/zijin/projects/ion_scattering_by_SWD/data/test_alg for result files.
└ @ DrWatson /Users/zijin/.julia/packages/DrWatson/8XYbJ/src/result_collection.jl:117
┌ Info: Added 54 entries.
└ @ DrWatson /Users/zijin/.julia/packages/DrWatson/8XYbJ/src/result_collection.jl:193

for subdf in groupby(df, :alg)
    pa_pair_hist(subset_leave(subdf)) |> display
    pa_pair_hist(subset_trap(subdf)) |> display
end

Test the different algorithms for the problem

pa_pair_plot(ldf, rdf) |> display
pa_pair_plot(ldf, rdf; func) |> display
pa_pair_plot(ldf, rdf; func=func_t) |> display

BoundsError: BoundsError: attempt to access 0-element Vector{Any} at index [1]
BoundsError: attempt to access 0-element Vector{Any} at index [1]



Stacktrace:

 [1] throw_boundserror(A::Vector{Any}, I::Tuple{Int64})

   @ Base ./essentials.jl:14

 [2] getindex(A::Vector{Any}, i::Int64)

   @ Base ./essentials.jl:915

 [3] pa_pair_plot(ldf::DataFrame, rdf::DataFrame; s::Symbol, func::Function, figure::@NamedTuple{size::Tuple{Int64, Int64}})

   @ Main ~/projects/ion_scattering_by_SWD/src/plot.jl:42

 [4] top-level scope

   @ ~/projects/ion_scattering_by_SWD/notebooks/jl_notebook_cell_df34fa98e69747e1a8f8a730347b8e2f_W4sZmlsZQ==.jl:3

pa_pair_plot(ldf, rdf; s = :β, func)

fig = Figure(; size=(1200, 1000,))
gridsl = pa_diff_plot!(fig[1,1], func(ldf))
gridsr = pa_diff_plot!(fig[3,1], func(rdf))
sign_label(fig)
legend!(fig[:,end+1], gridsl)
fig

tdf = @subset(func(df), :v .== 8, :θ .== 45);
l = data(tdf) * mapping(col=sign_map,row=β_map)
w1_map_plot(l; color=Δα)

# Create a layer with faceting for α, β, and v
begin
    using GLMakie
    GLMakie.activate!()
    fg = pa_pair_plot(layer; figure = figure)
end