Discussion
Rowan
GC: Mark and sweep

1. Set unvisited nodes to the set of environments in the heap.
1. Start at root node
	1. Fold thru slots
    	1. Call/3 on type
        2. Determine environment tag of result
        3. Add environment tag to unvisited node set

Types of reference:

1. array elements
2. block variables
3. block parents
4. compound functions

Schema:
```
-record(env,{parent,slots}).
-record(ctx,{ptr,env,stack}).
-record(prog,{root,heap}).
```
Top Answer
Marshall Lochbaum
Go VM:

```
package main
import "math"
import "os"

////////////////////////////////////////////////////////////////
// Types
type (
    Value     interface { call (x, w Value) Value }
    Modifier1 interface { call1(f    Value) Value }
    Modifier2 interface { call2(f, g Value) Value }

    Number float64; N=Number
    Character rune
    Array struct {
        data []Value
        shape []Value
    }

    Train  struct { f, g, h Value }
    RawFun struct { id uint8; f func(   x, w Value) Value }
    RawM1  struct { id uint8; m func(f, x, w Value) Value }
    DervR1 struct { f Value; m RawM1 }

    Block struct {
        typ uint8; imm bool; locals int
        bc []uint8
        objs []Value
        blocks []Block
    }
    Environment struct {
        vars []Value
        parent *Environment
    }
    BlockInst struct {
        id int
        typ uint8
        def *Block
        args []Value
        env Environment
    }

    Variable []Value
)

func assert(t bool) { if (!t) { panic("Assertion failed") } }

var ID int = 0
func (b Block) derive(env Environment) Value {
    if (b.typ==0 && b.imm) {
        return exec(b, Environment{make([]Value, b.locals), &env})
    } else {
        ID++; return BlockInst{ ID, b.typ, &b, nil, env }
    }
}
func (b BlockInst) eval(loc []Value) Value {
    return exec(*b.def, Environment{loc, &b.env})
}
func (e *Environment) get(n int) []Value {
    for ;n>0;n-- { e=e.parent }
    return e.vars
}

func call(f Value, x Value, w Value) Value {
    if x == nil { return x; } else { return f.call(x, w); }
}

func (v Number)    call(x, w Value) Value { return v; }
func (v Character) call(x, w Value) Value { return v; }
func (v Array)     call(x, w Value) Value { return v; }
func (t Train)     call(x, w Value) Value {
    r := call(t.h, x, w)
    var l Value; if t.f != nil { l = call(t.f, x, w); }
    return call(t.g, r, l)
}
func (f RawFun)    call(x, w Value) Value { return f.f(x,w); }
func (m RawM1)     call(x, w Value) Value { panic("Called 1-modifier as function"); }
func (f DervR1)    call(x, w Value) Value { return f.m.m(f.f, x, w) }
func (f BlockInst) call(x, w Value) Value {
    d := f.def
    assert(f.typ == 0)
    loc := make([]Value, 3 + len(f.args) + d.locals)
    loc[0]=f; loc[1]=x; loc[2]=w
    copy(loc[3:], f.args)
    return f.eval(loc)
}
func (v Variable)  call(x, w Value) Value { panic("Called symbol as function"); }

func (m RawM1) call1(f Value) Value { return DervR1{ f, m } }
func callBlockMod(m BlockInst, args []Value) Value {
    d := m.def
    if (d.imm) {
        loc := make([]Value, len(args) + d.locals)
        copy(loc, args)
        return m.eval(loc)
    } else {
        ID++; m.id=ID; m.args=args; m.typ=0; return m
    }
}
func (m BlockInst) call1(f Value) Value {
    assert(m.typ == 1)
    return callBlockMod(m, []Value{m, f})
}
func (m BlockInst) call2(f, g Value) Value {
    assert(m.typ == 2)
    return callBlockMod(m, []Value{m, f, g})
}


////////////////////////////////////////////////////////////////
// Virtual machine

func get(v Value) Value {
    switch i := v.(type) {
        case Variable: x:=i[0]; assert(x!=nil); return x
        case Array: return each(get, i)
        default: panic("System error: modification target is not a variable or array")
    }
}
func set(decl bool, id Value, vv Value) Value {
    switch i := id.(type) {
        case Variable: assert((i[0]==nil)==decl); i[0]=vv;
        case Array:
            v := vv.(Array); r := len(i.shape)
            assert(r==len(v.shape))
            for j:=0; j<r; j++ { assert(i.shape[j].(N)==v.shape[j].(N)) }
            for j,ij:=range i.data { set(decl, ij, v.data[j]) }
        default: panic("System error: assignment target is not a variable or array")
    }
    return vv
}

var stack = make([]Value, 0, 256)
func push(v Value) { stack = append(stack, v); }
func pop() Value { l:=len(stack); v:=stack[l-1]; stack=stack[:l-1]; return v; }
func exec(bl Block, env Environment) Value {
    bc:=bl.bc; objs:=bl.objs; blocks:=bl.blocks
    inst := func() uint8 { r:=bc[0]; bc=bc[1:]; return r }
    num := func() int {
        b:=128; n:=int(inst())
        if n<b { return n; }
        i:=1; t:=0
        for { t+=i*(n-b); i*=b; n=int(inst()); if n<b { break; } }
        return t+i*n
    }
    for { switch inst() {
    case  0: push(objs[num()])
    case 15: push(blocks[num()].derive(env))
    case  3,  4: l:=num(); i:=len(stack)-l; d:=make([]Value,l); copy(d,stack[i:]); stack=append(stack[:i],list(d))
    case  5, 16:          f:=pop();x:=pop();push(call(f,x,nil))
    case  6, 17: w:=pop();f:=pop();x:=pop();push(call(f,x,w))
    case  7    : f:=pop();m:=pop();         push(m.(Modifier1).call1(f  ))
    case  8    : f:=pop();m:=pop();g:=pop();push(m.(Modifier2).call2(f,g))
    case  9    :          g:=pop();h:=pop();push(Train{nil,g,h})
    case 10, 19: f:=pop();g:=pop();h:=pop();push(Train{f,  g,h})
    case 11    : i:=pop();         v:=pop();push(set(true ,i,v))
    case 12    : i:=pop();         v:=pop();push(set(false,i,v))
    case 13    : i:=pop();f:=pop();x:=pop();push(set(false,i,call(f,x,get(i))))
    case 14: pop()
    case 21: v:=env.get(num())[num()]; /*assert(v!=nil);*/ push(v)
    case 22: e:=env.get(num()); i:=num(); push(Variable(e[i:i+1]))
    case 25: return pop()
    } }
}

type BlockInd struct {
    typ uint8; imm int
    start int
    locals int
}
func run(bc []uint8, objs []Value, secs []BlockInd) Value {
    blocks := make([]Block, len(secs))
    for i,s:=range secs {
        blocks[i] = Block{ s.typ, s.imm!=0, s.locals,
                           bc[s.start:], objs, blocks }
    }
    return exec(blocks[0], Environment{make([]Value,blocks[0].locals),nil})
}
```
Answer #2
Rowan
starting to write gc using builtin erlang reference tracking
```erlang
ref(A,Fn) when is_function(Fn) -> % compound fn
    {env,[nil,X]} = fun_info(Fn,env),
    {env,[Env,_,_,_,_,_,_]} = fun_info(X,env),
    sets:del_element(Env,A);
ref(A,#m1{f=Fn} = M1) when is_record(M1,m1) -> % 1-modifier
    {env,[X]} = fun_info(Fn,env),
    {env,[Y]} = fun_info(X,env),
    {env,[Env,_,_,_,_,_,_]} = fun_info(Y,env),
    sets:del_element(Env,A);
ref(A,#v{r=R} = V) when is_record(V,v) ->
    sets:intersection(A,foldl(fun(I,X,B) -> ref(B,X) end,A,R));
ref(A,{R,I}) when is_reference(R) ->
    sets:del_element(R,A);
ref(A,Value) when Value =:= undefined; is_number(Value) ->
    A.
mark(Root,Heap,An,E,Stack) ->
    % init list of unvisited nodes
    Keys = sets:from_list(fetch_keys(Heap)),
    Unvisited = sets:del_element(E,Keys),
    Slots = fetch(Root,Heap),
    Refs = foldl(fun(_I,X,A) -> ref(A,X) end,Unvisited,Slots),
    StackRefs = case queue:is_queue(Stack) of
        true -> lists:foldl(fun(X,A) -> ref(A,X) end,Refs,queue:to_list(Stack));
        false -> ref(Refs,Stack)
    end,
    io:format("~p~n",[{mark,sets:size(Refs),sets:size(Keys),sets:size(StackRefs),sets:to_list(StackRefs)}]),
    ok.
gc() ->
    ok.
```
Answer #3
Rowan
```
-module(ebqn).

-export([test/0,fixed/1,concat/2,load_block/1]).
-import(array,[fix/1,from_list/1,resize/2,foldl/3,set/3]).
-import(queue,[cons/2,tail/1,liat/1,head/1,len/1]).
-import(dict,[store/3,fetch/2,fetch_keys/1]).

%-record(prog,{root,an,heap}).
%-record(e,{p,s}). % environment (parent, slots)

-record(vm,{f}). % vm (function) % HOF rtn'ed from load_vm
-record(v,{sh,r}). % value (shape, ravel)
-record(m1,{f}).
-record(m2,{f}).

test([B,O,S]) ->
    run(list_to_binary(B),list_to_tuple(O),list_to_tuple(lists:map(fun list_to_tuple/1,S))).
test() ->
    5 = test([[0,0,25],[5],[[0,1,0,0]]]), % 5
    3 = test([[0,0,14,0,1,25],[4,3],[[0,1,0,0]]]), % 4⋄3
    5 = test([[0,0,22,0,0,11,25],[5],[[0,1,0,1]]]), % a←5
    4 = test([[0,0,22,0,0,11,14,0,1,22,0,0,12,25],[5,4],[[0,1,0,1]]]), % a←5⋄a↩4
    2 = test([[0,0,22,0,0,11,14,0,1,22,0,1,11,14,21,0,0,25],[2,3],[[0,1,0,2]]]), % a←2⋄b←3⋄a
    1 = test([[0,0,22,0,0,11,14,0,1,21,0,0,16,25],[1,4],[[0,1,0,1]]]), % a←1⋄A 4
    2 = test([[0,0,22,0,0,11,14,0,2,21,0,0,0,1,17,25],[2,3,4],[[0,1,0,1]]]), % a←2⋄3 A 4
    6 = test([[0,0,15,1,16,25,21,0,1,25],[6],[[0,1,0,0],[0,0,6,3]]]), % {𝕩}6
    3 = test([[15,1,22,0,0,11,14,0,1,21,0,0,0,0,17,25,21,0,2,25],[3,4],[[0,1,0,1],[0,0,16,3]]]), % A←{𝕨}⋄3 A 4
    7 = test([[0,0,0,1,3,2,22,0,0,22,0,1,4,2,11,14,21,0,0,25],[7,2],[[0,1,0,2]]]), % a‿b←7‿2⋄a
    4 = test([[0,1,15,1,0,0,7,16,25,21,0,1,25],[4,6],[[0,1,0,0],[1,1,9,2]]]), % "4{𝔽}6"
    6 = test([[0,1,15,1,0,0,7,16,25,21,0,4,14,21,0,1,25],[4,6],[[0,1,0,0],[1,0,9,5]]]), % "4{𝔽⋄𝕩}6"
    1 = test([[0,1,15,1,15,2,0,0,8,16,25,21,0,1,25,21,0,2,25],[3,1],[[0,1,0,0],[0,0,11,3],[2,1,15,3]]]), % "3{𝔾}{𝕩} 1 "
    2 = test([[0,1,15,1,15,2,0,0,7,9,16,25,21,0,1,25,21,0,1,25],[2,3],[[0,1,0,0],[0,0,12,3],[1,1,16,2]]]), % "(2{𝔽}{𝕩})3 "
    3 = test([[0,1,15,1,15,2,9,0,0,17,25,21,0,2,21,0,1,3,2,25,21,0,1,22,0,3,22,0,4,4,2,11,14,21,0,3,25],[3,4],[[0,1,0,0],[0,0,11,3],[0,0,20,5]]]), % "3({a‿b←𝕩⋄a}{𝕨‿𝕩})4 "
    4 = test([[0,1,15,1,15,2,15,3,19,0,0,17,25,21,0,2,25,21,0,1,25,21,0,2,21,0,1,3,2,25],[4,5],[[0,1,0,0],[0,0,13,3],[0,0,17,3],[0,0,21,3]]]), % "4({𝕨‿𝕩}{𝕩}{𝕨})5"
    2 = test([[0,1,15,1,15,2,0,0,19,16,22,0,0,22,0,1,4,2,11,14,21,0,0,25,21,0,1,25,21,0,2,21,0,1,3,2,25],[2,5],[[0,1,0,2],[0,0,24,3],[0,0,28,3]]]), % "a‿b←(2{𝕨‿𝕩}{𝕩})5⋄a "
    2 = test([[0,2,22,0,0,11,15,1,15,2,15,3,19,16,25,0,1,22,1,0,12,14,21,0,1,25,21,0,1,14,21,1,0,25,0,0,22,1,0,12,14,21,0,1,25],[2,3,4],[[0,1,0,1],[0,0,15,3],[0,0,26,3],[0,0,34,3]]]), % "({a↩2⋄𝕩}{𝕩⋄a}{a↩3⋄𝕩})a←4 "
    8 = test([[0,0,22,0,0,11,14,0,1,15,1,22,0,0,13,14,21,0,0,25,21,0,1,25],[3,8],[[0,1,0,1],[0,0,20,3]]]), % "a←3⋄a{𝕩}↩8⋄a  "
    4 = test([[0,0,0,1,3,2,22,0,0,22,0,1,4,2,11,14,0,2,15,1,22,0,0,22,0,1,4,2,13,14,21,0,0,25,21,0,1,21,0,2,3,2,25],[2,1,4],[[0,1,0,2],[0,0,34,3]]]), % "a‿b←2‿1⋄a‿b{𝕩‿𝕨}↩4⋄a "
    1 = test([[0,0,22,0,0,11,14,15,1,14,21,0,0,25,0,1,22,0,0,11,25],[1,2],[[0,1,0,1],[0,1,14,1]]]), % "a←1⋄{a←2}⋄a"
    2 = test([[0,0,22,0,0,11,14,15,1,14,21,0,0,25,0,1,22,1,0,12,25],[1,2],[[0,1,0,1],[0,1,14,0]]]), % "a←1⋄{a↩2}⋄a"
    6 = test([[15,1,22,0,0,22,0,1,4,2,11,14,0,1,21,0,0,16,14,0,2,21,0,1,16,25,0,0,22,0,0,11,14,15,2,15,3,3,2,25,21,0,1,22,1,0,12,25,21,0,1,14,21,1,0,25],[2,6,0],[[0,1,0,2],[0,1,26,1],[0,0,40,3],[0,0,48,3]]]), % "f‿g←{a←2⋄{a↩𝕩}‿{𝕩⋄a}}⋄F 6⋄G 0"
    5 = test([[15,1,22,0,0,11,14,0,0,21,0,0,16,21,0,0,16,21,0,0,16,15,2,16,25,15,3,21,0,1,7,25,21,0,0,21,0,1,16,25,21,0,4,21,0,1,16,25],[5],[[0,1,0,1],[0,0,25,3],[0,0,32,3],[1,0,40,5]]]), % "L←{𝕩{𝕏𝕗}}⋄{𝕏𝕤}L L L 5"
    3 = test([[15,1,22,0,0,11,14,0,1,21,0,0,0,0,7,16,21,0,0,15,2,7,16,15,3,16,25,21,0,4,15,4,21,0,1,7,9,25,21,0,1,25,21,0,0,21,0,1,16,25,21,0,4,21,0,1,16,25],[3,5],[[0,1,0,1],[1,0,27,5],[0,0,38,3],[0,0,42,3],[1,0,50,5]]]), % "_l←{𝕩{𝕏𝕗} 𝔽}⋄{𝕏𝕤} {𝕩}_l 3 _l 5"
    1 = test([[0,1,15,1,15,2,15,3,8,0,0,17,25,21,0,1,25,21,0,1,21,0,2,15,4,21,0,1,7,19,25,21,0,2,25,21,0,2,21,0,4,21,0,1,17,25],[1,0],[[0,1,0,0],[0,0,13,3],[2,1,17,3],[0,0,31,3],[1,0,35,5]]]), % "1{𝕨}{𝔽{𝕩𝔽𝕨}𝔾𝔽}{𝕩}0 "
    2 = test([[0,0,0,1,0,2,0,3,0,4,15,1,3,2,3,2,3,2,3,2,3,2,15,2,0,0,0,0,15,3,3,2,3,2,7,16,25,21,0,1,25,21,0,1,15,4,16,25,21,0,1,25,21,0,1,22,0,3,22,0,4,4,2,11,14,21,0,0,22,0,5,11,14,0,0,15,5,15,6,7,16,14,15,7,21,0,4,21,0,5,16,7,25,21,0,1,21,1,4,16,25,21,0,0,14,15,8,22,1,5,12,25,21,0,1,22,0,5,22,0,6,4,2,11,14,21,0,6,21,0,4,16,25,21,0,0,14,15,9,25,21,0,1,22,0,3,22,0,4,4,2,11,14,21,0,3,25],[0,1,2,3,4],[[0,1,0,0],[0,0,37,3],[1,1,41,2],[0,0,48,3],[0,0,52,6],[1,1,93,2],[0,0,101,3],[1,0,112,7],[0,0,133,3],[0,0,140,5]]]), % "0‿(0‿{𝕩}){{a‿b←𝕩⋄t←𝕤⋄{𝕤⋄T↩{𝕤⋄{a‿b←𝕩⋄a}}}{B𝕗}0⋄(T b){a‿b←𝕩⋄𝔽b}}𝕗} 0‿(1‿(2‿(3‿(4‿{𝕩}))))"
    true.

arr(R,Sh) ->
    #v{r=R,sh=Sh}.
list(A) ->
    arr(A,[array:size(A)]).
call(_F,undefined,_W) ->
    undefined;
call(F,_X,_W) when not is_function(F) ->
    F;
call(F,X,W) ->
    true = (not is_record(F,m1) and not is_record(F,m2)),
    F(X,W).
tr3o(H,G,undefined) ->
    fun(X,W) ->
        call(G,call(H,X,W),undefined)
    end;
tr3o(H,G,F) ->
    fun(X,W) ->
        call(G,call(H,X,W),call(F,X,W))
    end.
ge(I,E,An) when I =:= 0 ->
    E;
ge(I,E,An) when I =/= 0 ->
    #{E := Parent} = An,
    ge(I-1,Parent,An).
hset(Heap,D,#v{r=Id},#v{r=V}) ->
    foldl(fun(J,N,A) -> hset(A,D,N,array:get(J,V)) end,Heap,Id);
hset(Heap,D,{E,I},V) ->
    A = fetch(E,Heap),
    D = (array:get(I,A) =:= undefined),
    store(E,set(I,V,A),Heap).
hget1(I) ->
    bad.
hget(Heap,{T,I}) when is_reference(T) ->
    Slots = fetch(T,Heap),
    Z = array:get(I,Slots),
    true = (null =/= Z),
    Z;
hget(Heap,#v{sh=S,r=R} = I) when is_record(I,v) ->
    arr(array:map(fun(_J,E) -> hget(Heap,E) end,R),S).
num(Binary,Ptr) ->
    {Size,Bitstring} = num(Binary,Ptr,0,<<>>),
    <<Value:Size/unsigned-integer>> = Bitstring,
    {Value,Ptr+trunc(Size/7)}.
num(Binary,Ptr,Size,Acc) ->
    <<H:1,Chunk:7/bitstring>> = binary_part(Binary,Ptr,1),
    num(Binary,Ptr,Size,Chunk,Acc,H).
num(_Binary,_Ptr,Size,Chunk,Acc,0) ->
    {Size+7,<<Chunk/bitstring,Acc/bitstring>>};
num(Binary,Ptr,Size,Chunk,Acc,1) ->
    num(Binary,Ptr+1,Size+7,<<Chunk/bitstring,Acc/bitstring>>).
fixed(X) when X =:= nil;X =:= [nil] ->
    nil;
fixed(X) when is_list(X) ->
    fix(resize(length(X),from_list(X))).
concat(nil,nil) ->
    nil;
concat(X,nil) when X =/= nil ->
    X;
concat(nil,W) when W =/= nil ->
    W;
concat(X,W) ->
    Xs = array:size(X),
    Z = array:resize(Xs+array:size(W),X),
    array:foldl(fun(I,V,A) -> array:set(Xs+I,V,A) end,Z,W).
tail(L,A,S) when L =:= -1 ->
    {A,S};
tail(L,A,S) ->
    tail(L-1,set(L,head(S),A),tail(S)).

args(B,P,Op) when Op =:= 7; Op =:= 8; Op =:= 9; Op =:= 11; Op =:= 12; Op =:= 13; Op =:= 14; Op =:= 16; Op =:= 17; Op =:= 19; Op =:= 25 ->
    {undefined,P};
args(B,P,Op) when Op =:= 0; Op =:= 3; Op =:= 4; Op =:= 15 ->
    num(B,P);
args(B,P,Op) when Op =:= 21; Op =:= 22 ->
    {X,Xp} = num(B,P),
    {Y,Yp} = num(B,Xp),
    {{X,Y},Yp}.

stack(B,O,S,Root,Heap,An,E,Stack,X,0) ->
    cons(element(1+X,O),Stack);
stack(B,O,S,Root,Heap,An,E,Stack,X,Op) when Op =:= 3; Op =:= 4 ->
    {T,Si} = case X of
        0 -> {list(fixed([])),Stack};
        _ -> tail(X-1,array:new(X),Stack)
    end,
    cons(list(T),Si);
stack(B,O,S,Root,Heap,An,E,Stack,undefined,7) ->
    F = head(Stack),
    #m1{f=M} = head(tail(Stack)),
    cons(M(F),tail(tail(Stack)));
stack(B,O,S,Root,Heap,An,E,Stack,undefined,8) ->
    F = head(Stack),
    #m2{f=M} = head(tail(Stack)),
    G = head(tail(tail(Stack))),
    cons(M(F,G),tail(tail(tail(Stack))));
stack(B,O,S,Root,Heap,An,E,Stack,undefined,9) ->
    G = head(Stack),
    J = head(tail(Stack)),
    cons(fun(X,W) -> call(G,call(J,X,W),undefined) end,tail(tail(Stack)));
stack(B,O,S,Root,Heap,An,E,Stack,X,Op) when Op =:= 11; Op =:= 12 ->
    tail(Stack);
stack(B,O,S,Root,Heap,An,E,Stack,X,13) ->
    tail(tail(Stack));
stack(B,O,S,Root,Heap,An,E,Stack,X,14) ->
    liat(Stack);
stack(B,O,S,Root,Heap,An,E,Stack,X,15) ->
    #vm{f=Block} = load_block(element(1+X,S)),
    Fn = Block(B,O,S,make_ref(),E),
    cons(Fn,Stack);
stack(B,O,S,Root,Heap,An,E,Stack,undefined,16) ->
    F = head(Stack),
    X = head(tail(Stack)),
    cons(call(F,X,undefined),tail(tail(Stack)));
stack(B,O,S,Root,Heap,An,E,Stack,undefined,17) ->
    W = head(Stack),
    F = head(tail(Stack)),
    X = head(tail(tail(Stack))),
    cons(call(F,X,W),tail(tail(tail(Stack))));
stack(B,O,S,Root,Heap,An,E,Stack,undefined,19) ->
    F = head(Stack),
    G = head(tail(Stack)),
    H = head(tail(tail(Stack))),
    cons(tr3o(H,G,F),tail(tail(tail(Stack))));
stack(B,O,S,Root,Heap,An,E,Stack,{X,Y},21) ->
    T = ge(X,E,An),
    Slots = fetch(T,Heap),
    Z = array:get(Y,Slots),
    true = (null =/= Z),
    cons(Z,Stack);
stack(B,O,S,Root,Heap,An,E,Stack,{X,Y},22) ->
    T = ge(X,E,An),
    cons({T,Y},Stack);
stack(B,O,S,Root,Heap,An,E,Stack,X,25) ->
    1 = len(Stack),
    head(Stack).

heap(Root,Heap,Stack,Op) when Op =:= 0; Op =:= 3; Op =:= 4; Op =:= 7; Op =:= 8; Op =:= 9; Op =:= 14; Op =:= 15; Op =:= 16; Op =:= 17; Op =:= 19; Op =:= 21; Op =:= 22; Op =:= 25 ->
    Heap;
heap(Root,Heap,Stack,Op) when Op =:= 11 ->
    I = head(Stack),
    V = head(tail(Stack)),
    hset(Heap,true,I,V);
heap(Root,Heap,Stack,Op) when Op =:= 12 ->
    I = head(Stack),
    V = head(tail(Stack)),
    hset(Heap,false,I,V);
heap(Root,Heap,Stack,Op) when Op =:= 13 ->
    I = head(Stack),
    F = head(tail(Stack)),
    X = head(tail(tail(Stack))),
    hset(Heap,false,I,call(F,X,hget(Heap,I))).

ctrl(Op) when Op =:= 0; Op =:= 3; Op =:= 4; Op =:= 7; Op =:= 8; Op =:= 9; Op =:= 11; Op =:= 12; Op =:= 13; Op =:= 14; Op =:= 15; Op =:= 16; Op =:= 17; Op =:= 19; Op =:= 21; Op =:= 22 ->
    cont;
ctrl(Op) when Op =:= 25 ->
    rtn.

ref(A,#v{r=R} = Value) when is_record(Value,v) ->
    foldl(fun(I,X,B) -> ref(B,X) end,A,R);
ref(A,Value) when Value =:= undefined; is_number(Value); is_function(Value) ->
    A;
ref(A,{R,I}) when is_reference(R) ->
    io:format("~p~n",[{ref,R,I}]),
    exit(self(),exit),A.
mark(Root,Heap,An,E,Stack) ->
    % init list of unvisited nodes
    Unvisited = fetch_keys(Heap),
    Slots = fetch(Root,Heap),
    Refs = foldl(fun(I,X,A) -> ref(A,X) end,sets:new(),Slots),
    % start at root
    % go thru stack, visit nodes (delay doing by only running when stack is empty)
    % go thru heap, starting at root node, fold over values, collect references to child environments
    % visit child environments, repeat fold over values...
    io:format("~p~n",[{mark,Unvisited,Slots,An,E,Stack}]),
    Refs.

vm(B,O,S,E,P,Stack,rtn) ->
    Stack;
vm(B,O,S,E,P,Stack,cont) ->
    Pi = P+1,
    {Op,Pi} = num(B,P),
    {Arg,Pn} = args(B,Pi,Op), % advances the ptr and reads the args
    Sn = stack(B,O,S,get(root),get(heap),get(an),E,Stack,Arg,Op), % mutates the stack
    put(heap,heap(get(root),get(heap),Stack,Op)), % mutates the heap
    Ctrl = ctrl(Op), % set ctrl atom
    mark(get(root),get(heap),get(an),E,Sn),
    vm(B,O,S,E,Pn,Sn,Ctrl). % call itself with new state

load_vm(B,O,S,E,Parent,V,ST) ->
    put(heap,store(E,V,get(heap))), % alloc
    An = get(an),
    put(an,An#{E => Parent}), % alloc
    vm(B,O,S,E,ST,queue:new(),cont). % run vm

load_block({T,I,ST,L}) -> % lexically scoped block
    Program = fun (B,O,S,E,Parent) ->
        C = fun(SV) -> load_vm(B,O,S,E,Parent,concat(SV,array:new(L)),ST) end,
        F =
            case T of
                0 -> fun(N) -> N(nil) end;
                1 -> fun(N) -> R = fun R(F  ) -> N(fixed([R,F  ])) end,#m1{f=R} end;
                2 -> fun(N) -> R = fun R(F,G) -> N(fixed([R,F,G])) end,#m2{f=R} end
            end,
        G =
            case I of
                0 -> fun(V) -> fun R(X,W) -> C(concat(fixed([R,X,W]),V)) end end;
                1 -> C
            end,
        F(G)
    end,
    #vm{f=Program}.

run(B,O,S) ->
    Root = make_ref(),
    Heap = dict:new(),
    An = #{}, % ancestors
    put(heap,Heap), % init the proc_dict
    put(root,Root),
    put(an,An),
    #vm{f=Block} = load_block(element(1,S)),
    Block(B,O,S,Root,Root). % set the root environment, and root as its own parent.
```

Enter question or answer id or url (and optionally further answer ids/urls from the same question) from

Separate each id/url with a space. No need to list your own answers; they will be imported automatically.