/* This software was developed by Bruce Hendrickson and Robert Leland   *
 * at Sandia National Laboratories under US Department of Energy        *
 * contract DE-AC04-76DP00789 and is copyrighted by Sandia Corporation. */

#include	<stdio.h>
#include	"structs.h"
#include	"defs.h"
#include	"refine_map.h"


void      compute_cube_vdata(vdata, comm_graph, vtx, mask, vtx2node)
struct refine_vdata *vdata;	/* preference data for a vertex */
struct vtx_data **comm_graph;	/* communication graph data structure */
int       vtx;			/* current vertex */
int       mask;			/* bit set in current hypercube dimension */
short    *vtx2node;		/* maps graph vtxs to mesh nodes */
{
    float     same;		/* my preference to stay where I am */
    float     change;		/* my preference to change this bit */
    float     ewgt;		/* weight of an edge */
    int       neighbor;		/* neighboring vtx in comm_graph */
    int       my_side;		/* which side of hypercube I'm on */
    int       neighb_side;	/* which side of hypercube neighbor's on */
    int       j;		/* loop counter */

    my_side = (vtx2node[vtx] & mask);

    change = 0;
    same = 0;
    for (j = 1; j < comm_graph[vtx]->nedges; j++) {
	neighbor = comm_graph[vtx]->edges[j];
	ewgt = comm_graph[vtx]->ewgts[j];

	neighb_side = (vtx2node[neighbor] & mask);

	if (neighb_side != my_side) {
	    change += ewgt;
	}
	else {
	    same += ewgt;
	}
    }
    vdata->same = same;
    vdata->above = change;
}


syntax highlighted by Code2HTML, v. 0.9.1