init_data = precompute orbit decomposition

Orb_AutF4.jl

@@ -199,6 +199,23 @@ function λandP(m::JuMP.Model, data::ProblemData)
    return λ, recP
 end
 
+function init_data(name::String, N::Int; radius=2)
+
+   ex(fname) = isfile(joinpath(name, fname))
+
+   conditions = ex.(["delta.jld", "pm.jld", "U_pis.jld", "orbits.jld"])
+   if !all(conditions)
+      SOutFN = AutGroup(FreeGroup(N), special=true, outer=true)
+      info(logger, SOutFN)
+      S = generators(SOutFN);
+      S = [S; [inv(s) for s in S]]
+      info(logger, "Symmetric generating set of size $(length(S))")
+      info(logger, S)
+      compute_orbit_data(name, SOutFN, S, radius=radius)
+   end
+   return 0
+end
+
 function main()
    name = "SOutF4_E4"