algebraicjulia / combinatorialchains.jl Goto Github PK

J₀ = @acset IsingModel begin
    V1 = 9
    L1 = 12
    src1 = [1,2,5,5,5,5,7,8,3,6,1,4]
    tgt1 = [2,3,2,4,6,8,8,9,6,9,4,7]
  end

  m = homomorphisms(codom(rule(4)[1]), J₀, monic=true)[1]
  J₁ = rewrite_match(rule(4)..., m)
  to_graphviz(J₁)




  function rewrite_ising(j::IsingCats.AbstractIsingModel, T, maxtries=100)
    # choose a random rule
    for i in 1:maxtries
      l,r = rule(rand(0:4))
    if rand(Bool)
        r,l = l,r
    end
  
  #      l,r = rule(4)

      qₗ = homomorphism_query(codom(l))
      matches = query(j, qₗ)
      @show length(matches)

      αs = map(ρ -> make_homomorphism(matches[ρ], codom(l), j),
         1:length(matches))


      # quick hack for "monic on V1"
      αs = filter(αs) do α
        length(unique(collect(components(α).V1))) == length(collect(components(α).V1))
      end

      @show length(αs)
      for α in αs
        println(components(α).V1)
      end

      αs = filter(α->valid_dpo(l, α), αs)
      @show length(αs)

      if length(αs) > 0

        for i in 1:maxtries
          # compute table of matches
          # pick random match
          @show ρ = rand(1:length(αs))
          α = αs[ρ]
          @show components(α)
          @show valid_dpo(l, α)
          @show accept_rewrite(Span(l, r), T)

          if valid_dpo(l, α) && accept_rewrite(Span(l, r), T)
            α₂ = rewrite_match(l, r, α)
            return α₂
          end
        end
      end
  end 
     error("Could not find a valid match in $maxtries attempts")
end

 K = rewrite_ising(J₀, 2, 80)
 to_graphviz(K)