static int calc_depth
( mclx* m_transient
)
{ mclx* m_inverse = mclxTranspose(m_transient)
;
dim c, depth = 0
;
if(0)puts("")
;
for (c=0; c<N_COLS(m_inverse); c++)
{ dim this_depth = 0
;
if (!m_inverse->cols[c].n_ivps) /* no incoming nodes */
{ mclv* next = mclxGetVector(m_transient, m_inverse->cols[c].vid, RETURN_ON_FAIL, NULL)
;
if (!next)
continue
;
mclgUnionvInitList(m_transient, next)
;
do
{ mclv* next2 = mclgUnionv(m_transient, next, NULL, SCRATCH_UPDATE, NULL)
;
if (0 && next->ivps)
fprintf(stdout, "chain %d ->\n", (int) m_inverse->cols[c].vid)
, mclvaDump(next, stdout, -1, " ", 0)
;
if (this_depth) /* otherwise starting vector in matrix */
mclvFree(&next)
;
next = next2
;
this_depth++
;
}
while (next->n_ivps)
;
mclvFree(&next) /* did loop at least once, so not the starting vector */
;
mclgUnionvReset(m_transient)
;
}
if (this_depth > depth)
depth = this_depth
;
}
mclxFree(&m_inverse)
;
return depth
;
}
mclMatrix* mclInterpret
( mclMatrix* dag
)
{ mclv* v_attr = mclvCopy(NULL, dag->dom_cols)
; mclx* m_attr = NULL, *m_cls = NULL, *m_clst = NULL
; dim d
; mclvMakeCharacteristic(v_attr)
; for (d=0;d<N_COLS(dag);d++)
{ mclv* col = dag->cols+d
; if (mclvGetIvp(col, col->vid, NULL)) /* deemed attractor */
mclvInsertIdx(v_attr, col->vid, 2.0)
; }
mclvSelectGqBar(v_attr, 1.5)
; m_attr = mclxSub(dag, v_attr, v_attr)
; mclxAddTranspose(m_attr, 1.0)
; m_cls = clmUGraphComponents(m_attr, NULL) /* attractor systems as clusters */
; mclvCopy(m_cls->dom_rows, dag->dom_cols) /* add all nodes to this cluster matrix */
; m_clst = mclxTranspose(m_cls) /* nodes(columns) with zero neighbours need to be classified */
; mclgUnionvReset(dag) /* make mx->dom-rows characteristic */
; mclxFree(&m_cls)
; for (d=0;d<N_COLS(dag);d++)
{ mclv* closure, *clsids
; if (mclvGetIvp(v_attr, dag->cols[d].vid, NULL))
continue /* attractor already classified */
; closure = get_closure(dag, dag->cols+d) /* take all [neighbours of [neighbours of [..]]] */
; clsids = mclgUnionv(m_clst, closure, NULL, SCRATCH_READY, NULL)
; mclvAdd(m_clst->cols+d, clsids, m_clst->cols+d)
; mclvFree(&clsids)
; mclvFree(&closure)
; }
m_cls = mclxTranspose(m_clst)
; mclxFree(&m_attr)
; mclxFree(&m_clst)
; mclvFree(&v_attr)
; return m_cls
; }
