/* mapGraphs:
*/
static void mapGraphs(graph_t * g, graph_t * cg, distfn dist)
{
node_t *n;
edge_t *e;
edge_t *ce;
node_t *t;
node_t *h;
nitem *tp;
nitem *hp;
int delta;
for (n = agfstnode(g); n; n = agnxtnode(g, n)) {
tp = (nitem *) ND_alg(n);
t = tp->cnode;
for (e = agfstout(g, n); e; e = agnxtout(g, e)) {
hp = (nitem *) ND_alg(aghead(e));
delta = dist(&tp->bb, &hp->bb);
h = hp->cnode;
#ifndef WITH_CGRAPH
ce = agedge(cg, t, h);
#else
ce = agedge(cg, t, h, NULL, 1);
agbindrec(ce, "Agedgeinfo_t", sizeof(Agedgeinfo_t), TRUE);
#endif
ED_weight(ce) = 1;
if (ED_minlen(ce) < delta) {
if (ED_minlen(ce) == 0.0) {
elist_append(ce, ND_out(t));
elist_append(ce, ND_in(h));
}
ED_minlen(ce) = delta;
}
}
}
}
/* mapGraphs:
*/
static void mapGraphs(graph_t * g, graph_t * cg, distfn dist)
{
node_t *n;
edge_t *e;
edge_t *ce;
node_t *t;
node_t *h;
nitem *tp;
nitem *hp;
int delta;
for (n = agfstnode(g); n; n = agnxtnode(g, n)) {
tp = (nitem *) ND_alg(n);
t = tp->cnode;
for (e = agfstout(g, n); e; e = agnxtout(g, e)) {
hp = (nitem *) ND_alg(e->head);
delta = dist(&tp->bb, &hp->bb);
h = hp->cnode;
ce = agedge(cg, t, h);
ED_weight(ce) = 1;
if (ED_minlen(ce) < delta) {
if (ED_minlen(ce) == 0.0) {
elist_append(ce, ND_out(t));
elist_append(ce, ND_in(h));
}
ED_minlen(ce) = delta;
}
}
}
}
static void
sameedge (same_t* same, node_t* n, edge_t *e, char *id)
/* register E in the SAME structure of N under ID. Uses static int N_SAME */
{
int i,sflag,eflag,flag;
for (i=0; i<n_same; i++)
if (streq(same[i].id,id)) {
elist_append(e,same[i].l);
goto set_arrow;
}
if (++n_same > MAXSAME) {
fprintf(stderr,"too many same{head,tail} groups for node %s\n", n->name);
return;
}
alloc_elist(1,same[i].l);
elist_fastapp(e,same[i].l);
same[i].id = id;
same[i].n_arr = 0;
same[i].arr_len = 0;
set_arrow:
arrow_flags (e, &sflag, &eflag);
if ((flag = e->head==n ? eflag : sflag))
same[i].arr_len =
/* only consider arrows if there's exactly one arrow */
(++same[i].n_arr==1) ? arrow_length(e,flag) : 0;
}
static void checkChain(graph_t * g)
{
node_t *t;
node_t *h;
edge_t *e;
t = GD_nlist(g);
for (h = ND_next(t); h; h = ND_next(h)) {
if (!agfindedge(g, t, h)) {
e = agedge(g, t, h);
ED_minlen(e) = 0;
elist_append(e, ND_out(t));
elist_append(e, ND_in(h));
}
t = h;
}
}
/* safe_list_append - append e to list L only if e not already a member */
void safe_list_append(edge_t *e, elist *L)
{
int i;
for (i = 0; i < L->size; i++) if (e == L->list[i]) return;
elist_append(e,(*L));
}
edge_t*
fast_edge(edge_t *e)
{
#ifdef DEBUG
int i;
edge_t *f;
for (i = 0; (f = ND_out(e->tail).list[i]); i++) {
if (e == f) {fprintf(stderr,"duplicate fast edge\n"); return;}
assert (e->head != f->head);
}
for (i = 0; (f = ND_in(e->head).list[i]); i++) {
if (e == f) {fprintf(stderr,"duplicate fast edge\n"); return;}
assert (e->tail != f->tail);
}
#endif
elist_append(e,ND_out(e->tail));
elist_append(e,ND_in(e->head));
return e;
}
/**
* Append `key' to the list.
*
* It is safe to call this routine whilst iterating.
*/
void
hash_list_append(hash_list_t *hl, const void *key)
{
struct hash_list_item *item;
hash_list_check(hl);
WALLOC(item);
item->key = key;
elist_append(&hl->list, item);
hash_list_insert_item(hl, item);
}
static void checkChain(graph_t * g)
{
node_t *t;
node_t *h;
edge_t *e;
t = GD_nlist(g);
for (h = ND_next(t); h; h = ND_next(h)) {
if (!agfindedge(g, t, h)) {
#ifdef WITH_CGRAPH
e = agedge(g, t, h, NULL, 1);
agbindrec(e, "Agedgeinfo_t", sizeof(Agedgeinfo_t), TRUE);
#else
e = agedge(g, t, h);
#endif
ED_minlen(e) = 0;
elist_append(e, ND_out(t));
elist_append(e, ND_in(h));
}
t = h;
}
}
/**
* Add value to the table.
*
* If it was already present, its lifetime is reset to the aging delay.
*
* The key argument is freed immediately if there is a free routine for
* keys and the key was present in the table.
*
* The previous value is freed and replaced by the new one if there is
* an insertion conflict and the key pointers are different.
*/
void
aging_insert(aging_table_t *ag, const void *key, void *value)
{
bool found;
void *ovalue;
time_t now = tm_time();
struct aging_value *aval;
aging_check(ag);
aging_synchronize(ag);
found = hikset_lookup_extended(ag->table, key, &ovalue);
if (found) {
aval = ovalue;
if (ag->kvfree != NULL) {
/*
* We discard the new and keep the old key instead.
* That way, we don't have to update the hash table.
*/
(*ag->kvfree)(deconstify_pointer(key), aval->value);
}
/*
* Value existed for this key, reset its lifetime by moving the
* entry to the tail of the list.
*/
aval->value = value;
aval->last_insert = now;
elist_moveto_tail(&ag->list, aval);
} else {
WALLOC(aval);
aval->value = value;
aval->key = deconstify_pointer(key);
aval->last_insert = now;
hikset_insert(ag->table, aval);
elist_append(&ag->list, aval);
}
aging_return_void(ag);
}
/**
* Perform the transformation (can be called from case or into commands)
* @param matchme
* @param destvalue
* @return
*/
static int revm_case_transform(revmexpr_t *matchme, char *destvalue)
{
u_int dstnbr;
char *curptr;
char *foundptr;
u_int curidx;
list_t *exprlist;
aspectype_t *type;
char namebuf[BUFSIZ];
char *rname;
revmexpr_t *candid;
list_t *looplist;
int ret;
PROFILER_IN(__FILE__, __FUNCTION__, __LINE__);
/* We matched : first find how many elements there is in the target (list) type */
exprlist = world.curjob->recur[world.curjob->curscope].rwrt.transformed;
dstnbr = 1;
for (curidx = world.curjob->iter[world.curjob->curloop].elmidx - 1, curptr = destvalue;
curptr && *curptr;
curptr = foundptr + 2, curidx++, dstnbr++)
{
foundptr = strstr(curptr, "::");
if (!foundptr && dstnbr == 1)
break;
if (foundptr)
*foundptr = 0x00;
type = revm_exprtype_get(curptr);
snprintf(namebuf, BUFSIZ, "%s-%u", world.curjob->iter[world.curjob->curloop].curkey, curidx);
rname = strdup(namebuf);
candid = revm_expr_create(type, rname, curptr);
if (!candid || !candid->annot)
PROFILER_ERR(__FILE__, __FUNCTION__, __LINE__,
"Invalid target type(s) for transformation", -1);
candid->annot->inhash = 1;
elist_append(exprlist, rname, candid);
if (!foundptr)
break;
}
/* XXX: Case where the rewritten element is not part of an iterated list -- unsupported */
looplist = (list_t *) world.curjob->iter[world.curjob->curloop].container;
if (!looplist)
PROFILER_ERR(__FILE__, __FUNCTION__, __LINE__,
"Rewriting of non-list element unimplemented", -1);
else if (looplist->type != ASPECT_TYPE_EXPR)
PROFILER_ERR(__FILE__, __FUNCTION__, __LINE__,
"Rewrite of non-expression variables unimplemented", -1);
/* Simply replace the current list element now */
/* The type of the list (list_t->type) does not change : it still is a list of revmexpr_t */
/* Just one element to swap */
else if (dstnbr == 1)
{
/* No transformation, keep the original expression */
if (!strcmp(destvalue, "."))
candid = matchme;
else
{
/* Case where we did not have a simple value here but a real expression */
rname = revm_tmpvar_create();
type = revm_exprtype_get(destvalue);
candid = revm_expr_create(type, rname, destvalue);
if (!candid || !candid->annot)
PROFILER_ERR(__FILE__, __FUNCTION__, __LINE__,
"Malformed destination type", -1);
candid->annot->inhash = 1;
/**
**
** XXX: Disabled -- do not remove -- to enable when SSA transform is ready
** if (revm_links_propagate(candid, matchme) < 0)
** PROFILER_ERR(__FILE__, __FUNCTION__, __LINE__,
** "Error while propagating dataflow links", -1);
**
** XXX: Also address the propagate for the case where dstnbr != 1
*/
ret = elist_set(looplist, world.curjob->iter[world.curjob->curloop].curkey, candid);
if (ret < 0)
PROFILER_ERR(__FILE__, __FUNCTION__, __LINE__,
"Failed to store transformed expr in input list", -1);
/*
** We reuse the same key for transform list, so that any change in transform list
** can be recorded back to the element of the original list. We also refresh the
** induction variable of current loop context.
**
** XXX: we should probably update context hashes too
*/
elist_append(exprlist, strdup(world.curjob->iter[world.curjob->curloop].curkey), candid);
char *name = world.curjob->iter[world.curjob->curloop].curind->label;
world.curjob->iter[world.curjob->curloop].curind = candid;
world.curjob->iter[world.curjob->curloop].curind->label = name;
}
}
/* Insert a list at a certain offset of the list and remove matched element */
else
{
elist_replace(looplist, world.curjob->iter[world.curjob->curloop].curkey,
elist_copy(exprlist, ELIST_DATA_NOCOPY));
world.curjob->iter[world.curjob->curloop].elmidx += exprlist->elmnbr - 1;
}
PROFILER_ROUT(__FILE__, __FUNCTION__, __LINE__, 0);
}
void flat_edge(graph_t * g, edge_t * e)
{
elist_append(e, ND_flat_out(e->tail));
elist_append(e, ND_flat_in(e->head));
GD_has_flat_edges(g->root) = GD_has_flat_edges(g) = TRUE;
}
void other_edge(edge_t * e)
{
elist_append(e, ND_other(e->tail));
}
/* makeGraphs:
* Generate dags modeling the row and column constraints.
* If the table has cc columns, we create the graph
* 0 -> 1 -> 2 -> ... -> cc
* and if a cell starts in column c with span cspan, with
* width w, we add the edge c -> c+cspan [minlen = w].
*
* We might simplify the graph by removing multiedges,
* using the max minlen, but will affect the balancing?
*/
void makeGraphs(htmltbl_t * tbl, graph_t * rowg, graph_t * colg)
{
htmlcell_t *cp;
htmlcell_t **cells;
node_t *t;
node_t *lastn;
node_t *h;
edge_t *e;
int i;
int* minc;
int* minr;
lastn = NULL;
for (i = 0; i <= tbl->cc; i++) {
t = agnode(colg, nToName(i));
alloc_elist(tbl->rc, ND_in(t));
alloc_elist(tbl->rc, ND_out(t));
if (lastn) {
ND_next(lastn) = t;
lastn = t;
} else {
lastn = GD_nlist(colg) = t;
}
}
lastn = NULL;
for (i = 0; i <= tbl->rc; i++) {
t = agnode(rowg, nToName(i));
alloc_elist(tbl->cc, ND_in(t));
alloc_elist(tbl->cc, ND_out(t));
if (lastn) {
ND_next(lastn) = t;
lastn = t;
} else {
lastn = GD_nlist(rowg) = t;
}
}
minr = N_NEW(tbl->rc, int);
minc = N_NEW(tbl->cc, int);
for (cells = tbl->u.n.cells; *cells; cells++) {
int x, y, c, r;
cp = *cells;
x = (cp->data.box.UR.x + (cp->cspan-1))/cp->cspan;
for (c = 0; c < cp->cspan; c++)
minc[cp->col + c] = MAX(minc[cp->col + c],x);
y = (cp->data.box.UR.y + (cp->rspan-1))/cp->rspan;
for (r = 0; r < cp->rspan; r++)
minr[cp->row + r] = MAX(minr[cp->row + r],y);
}
for (cells = tbl->u.n.cells; *cells; cells++) {
int x, y, c, r;
cp = *cells;
t = agfindnode(colg, nToName(cp->col));
h = agfindnode(colg, nToName(cp->col + cp->cspan));
e = agedge(colg, t, h);
x = 0;
for (c = 0; c < cp->cspan; c++)
x += minc[cp->col + c];
ED_minlen(e) = x;
/* ED_minlen(e) = cp->data.box.UR.x; */
#ifdef DEBUG
fprintf(stderr, "col edge %s -> %s %d\n", t->name, h->name,
ED_minlen(e));
#endif
elist_append(e, ND_out(t));
elist_append(e, ND_in(h));
t = agfindnode(rowg, nToName(cp->row));
h = agfindnode(rowg, nToName(cp->row + cp->rspan));
e = agedge(rowg, t, h);
y = 0;
for (r = 0; r < cp->rspan; r++)
y += minr[cp->row + r];
ED_minlen(e) = y;
/* ED_minlen(e) = cp->data.box.UR.y; */
#ifdef DEBUG
fprintf(stderr, "row edge %s -> %s %d\n", t->name, h->name,
ED_minlen(e));
#endif
elist_append(e, ND_out(t));
elist_append(e, ND_in(h));
}
/* Make sure that 0 <= 1 <= 2 ...k. This implies graph connected. */
checkChain(colg);
checkChain(rowg);
free (minc);
free (minr);
}
/* mkConstraintG:
*/
static graph_t *mkConstraintG(graph_t * g, Dt_t * list,
intersectfn intersect, distfn dist)
{
nitem *p;
nitem *nxt = NULL;
nitem *nxp;
graph_t *vg;
node_t *prev = NULL;
node_t *root = NULL;
node_t *n = NULL;
edge_t *e;
int lcnt, cnt;
int oldval = -INT_MAX;
#ifdef OLD
double root_val;
#endif
node_t *lastn = NULL;
#ifndef WITH_CGRAPH
graph_t *cg = agopen("cg", AGDIGRAPHSTRICT);
#else
graph_t *cg = agopen("cg", Agstrictdirected, NIL(Agdisc_t *));
agbindrec(cg, "Agraphinfo_t", sizeof(Agraphinfo_t), TRUE); // graph custom data
#endif
/* count distinct nodes */
cnt = 0;
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
if (oldval != p->val) {
oldval = p->val;
cnt++;
}
}
/* construct basic chain to enforce left to right order */
oldval = -INT_MAX;
lcnt = 0;
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
if (oldval != p->val) {
oldval = p->val;
/* n = newNode (cg); */
#ifndef WITH_CGRAPH
n = agnode(cg, agnameof(p->np)); /* FIX */
#else
n = agnode(cg, agnameof(p->np), 1); /* FIX */
agbindrec(n, "Agnodeinfo_t", sizeof(Agnodeinfo_t), TRUE); //node custom data
#endif
ND_alg(n) = p;
if (root) {
ND_next(lastn) = n;
lastn = n;
} else {
root = n;
#ifdef OLD
root_val = p->val;
#endif
lastn = GD_nlist(cg) = n;
}
alloc_elist(lcnt, ND_in(n));
if (prev) {
if (prev == root)
alloc_elist(2 * (cnt - 1), ND_out(prev));
else
alloc_elist(cnt - lcnt - 1, ND_out(prev));
#ifndef WITH_CGRAPH
e = agedge(cg, prev, n);
#else
e = agedge(cg, prev, n, NULL, 1);
agbindrec(e, "Agedgeinfo_t", sizeof(Agedgeinfo_t), TRUE); // edge custom data
#endif
ED_minlen(e) = SCALE;
ED_weight(e) = 1;
elist_append(e, ND_out(prev));
elist_append(e, ND_in(n));
}
lcnt++;
prev = n;
}
p->cnode = n;
}
alloc_elist(0, ND_out(prev));
/* add immediate right neighbor constraints
* Construct visibility graph, then perform transitive reduction.
* Remaining outedges are immediate right neighbors.
* FIX: Incremental algorithm to construct trans. reduction?
*/
#ifndef WITH_CGRAPH
vg = agopen("vg", AGDIGRAPHSTRICT);
#else
vg = agopen("vg", Agstrictdirected, NIL(Agdisc_t *));
#endif
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
#ifndef WITH_CGRAPH
n = agnode(vg, agnameof(p->np)); /* FIX */
#else
n = agnode(vg, agnameof(p->np), 1); /* FIX */
agbindrec(n, "Agnodeinfo_t", sizeof(Agnodeinfo_t), TRUE); //node custom data
#endif
p->vnode = n;
ND_alg(n) = p;
}
oldval = -INT_MAX;
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
if (oldval != p->val) { /* new pos: reset nxt */
oldval = p->val;
for (nxt = (nitem *) dtlink(link, (Dtlink_t *) p); nxt;
nxt = (nitem *) dtlink(list, (Dtlink_t *) nxt)) {
if (nxt->val != oldval)
break;
}
if (!nxt)
break;
}
for (nxp = nxt; nxp;
nxp = (nitem *) dtlink(list, (Dtlink_t *) nxp)) {
if (intersect(p, nxp))
#ifndef WITH_CGRAPH
agedge(vg, p->vnode, nxp->vnode);
#else
agedge(vg, p->vnode, nxp->vnode, NULL, 1);
#endif
}
}
/* Remove redundant constraints here. However, the cost of doing this
* may be a good deal more than the time saved in network simplex. Also,
* if the graph is changed, the ND_in and ND_out data has to be updated.
*/
mapGraphs(vg, cg, dist);
agclose(vg);
/* add dummy constraints for absolute values and initial positions */
#ifdef OLD
for (n = agfstnode(cg); n; n = agnxtnode(cg, n)) {
node_t *vn; /* slack node for absolute value */
node_t *an; /* node representing original position */
p = (nitem *) ND_alg(n);
if ((n == root) || (!p))
continue;
vn = newNode(cg);
ND_next(lastn) = vn;
lastn = vn;
alloc_elist(0, ND_out(vn));
alloc_elist(2, ND_in(vn));
an = newNode(cg);
ND_next(lastn) = an;
lastn = an;
alloc_elist(1, ND_in(an));
alloc_elist(1, ND_out(an));
#ifndef WITH_CGRAPH
e = agedge(cg, root, an);
#else
e = agedge(cg, root, an, 1);
#endif
ED_minlen(e) = p->val - root_val;
elist_append(e, ND_out(root));
elist_append(e, ND_in(an));
#ifndef WITH_CGRAPH
e = agedge(cg, an, vn);
#else
e = agedge(cg, an, vn, 1);
#endif
elist_append(e, ND_out(an));
elist_append(e, ND_in(vn));
#ifndef WITH_CGRAPH
e = agedge(cg, n, vn);
#else
e = agedge(cg, n, vn, 1);
#endif
elist_append(e, ND_out(n));
elist_append(e, ND_in(vn));
}
#endif /* OLD */
return cg;
}
/* mkNConstraintG:
* Similar to mkConstraintG, except it doesn't enforce orthogonal
* ordering. If there is overlap, as defined by intersect, the
* nodes will kept/pushed apart in the current order. If not, no
* constraint is enforced. If a constraint edge is added, and it
* corresponds to a real edge, we increase the weight in an attempt
* to keep the resulting shift short.
*/
static graph_t *mkNConstraintG(graph_t * g, Dt_t * list,
intersectfn intersect, distfn dist)
{
nitem *p;
nitem *nxp;
node_t *n;
edge_t *e;
node_t *lastn = NULL;
#ifndef WITH_CGRAPH
graph_t *cg = agopen("cg", AGDIGRAPHSTRICT);
#else
graph_t *cg = agopen("cg", Agstrictdirected, NIL(Agdisc_t *));
agbindrec(cg, "Agraphinfo_t", sizeof(Agraphinfo_t), TRUE); // graph custom data
#endif
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
#ifndef WITH_CGRAPH
n = agnode(cg, agnameof(p->np)); /* FIX */
#else
n = agnode(cg, agnameof(p->np), 1); /* FIX */
agbindrec(n, "Agnodeinfo_t", sizeof(Agnodeinfo_t), TRUE); //node custom data
#endif
ND_alg(n) = p;
p->cnode = n;
alloc_elist(0, ND_in(n));
alloc_elist(0, ND_out(n));
if (lastn) {
ND_next(lastn) = n;
lastn = n;
} else {
lastn = GD_nlist(cg) = n;
}
}
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
for (nxp = (nitem *) dtlink(link, (Dtlink_t *) p); nxp;
nxp = (nitem *) dtlink(list, (Dtlink_t *) nxp)) {
e = NULL;
if (intersect(p, nxp)) {
double delta = dist(&p->bb, &nxp->bb);
#ifndef WITH_CGRAPH
e = agedge(cg, p->cnode, nxp->cnode);
#else
e = agedge(cg, p->cnode, nxp->cnode, NULL, 1);
agbindrec(e, "Agedgeinfo_t", sizeof(Agedgeinfo_t), TRUE); // edge custom data
#endif
assert (delta <= 0xFFFF);
ED_minlen(e) = delta;
ED_weight(e) = 1;
}
if (e && agfindedge(g,p->np, nxp->np)) {
ED_weight(e) = 100;
}
#if 0
if (agfindedge(g,p->np, nxp->np)) {
if (e == NULL)
e = agedge(cg, p->cnode, nxp->cnode);
ED_weight(e) = 100;
/* If minlen < SCALE, the nodes can't conflict or there's
* an overlap but it will be removed in the orthogonal pass.
* So we just keep the node's basically where they are.
*/
if (SCALE > ED_minlen(e))
ED_minlen(e) = SCALE;
}
#endif
}
}
for (p = (nitem *) dtflatten(list); p;
p = (nitem *) dtlink(list, (Dtlink_t *) p)) {
n = p->cnode;
for (e = agfstout(cg,n); e; e = agnxtout(cg,e)) {
elist_append(e, ND_out(n));
elist_append(e, ND_in(aghead(e)));
}
}
/* We could remove redundant constraints here. However, the cost of doing
* this may be a good deal more than the time saved in network simplex.
* Also, if the graph is changed, the ND_in and ND_out data has to be
* updated.
*/
return cg;
}
