/*
* Zoltan_Oct_subtree_dref(oct)
*
* Coarsen octree so that leaf octants do not have less than MINOCTREGIONS
* regions. Refinement takes precedence, so coarsening will not take place
* unless all subtrees agree on it.
*/
static int Zoltan_Oct_subtree_dref(ZZ *zz, OCT_Global_Info *OCT_info,
pOctant oct)
{
pOctant child; /* child of an octant */
int coarsen; /* flag to indicate need to coarsen */
int i; /* index counter */
int nregions; /* number of regions */
int total; /* total number of regions */
/* if terminal octant, cannot coarsen on own */
if (Zoltan_Oct_isTerminal(oct)) {
nregions=Zoltan_Oct_nRegions(oct);
if (nregions > (MAXOCTREGIONS*2) ) {
fprintf(stderr,
"OCT Zoltan_Oct_subtree_dref: warning: too many (%d) regions "
"in oct (id=%d)\n",Zoltan_Oct_nRegions(oct),Zoltan_Oct_id(oct));
return(-1);
}
else
if (nregions < MINOCTREGIONS)
return(nregions);
else
return(-1);
}
coarsen=1; /* assume to be coarsen */
total=0;
/* look at each child, see if they need to be coarsened */
for (i=0; coarsen && i<8; i++) {
child = Zoltan_Oct_child(oct,i);
/* if child is off processor cannot coarsen */
/* if (!Zoltan_Oct_local(child) ) */
if(!Zoltan_Oct_POct_local(OCT_info, oct, i))
coarsen=0;
else {
/* get the number of region of the child */
nregions=Zoltan_Oct_subtree_dref(zz, OCT_info,child);
if (nregions<0)
coarsen=0;
else
total+=nregions; /* total the region count */
}
}
/* check if octant can be coarsened */
if (coarsen && total<MAXOCTREGIONS) {
Zoltan_Oct_terminal_coarsen(zz, OCT_info,oct);
if (total < MINOCTREGIONS)
return(total);
else
return(-1);
}
return(-1);
}
/*
* int Zoltan_Oct_dfs_SetIds(pOctant octant, int number_of_previous_octants)
*
* sets the ids of all the octants so that there is a global numbering
*/
static int Zoltan_Oct_dfs_SetIds(OCT_Global_Info *OCT_info, pOctant oct,
int nprevoct) {
int id, /* the id number of an octant */
i; /* index counter */
int pid;
pOctant child; /* ith child of an octant */
if(Zoltan_Oct_isTerminal(oct)) {
id = Zoltan_Oct_id(oct);
Zoltan_Oct_setID(oct, id+nprevoct); /* now have global id's */
}
else {
for(i=0; i<8; i++) {
child = Zoltan_Oct_child(oct, i);
pid = Zoltan_Oct_Cpid(oct,i);
if ((pid == OCT_info->OCT_localpid) && child != NULL)
Zoltan_Oct_dfs_SetIds(OCT_info,child, nprevoct);
}
id = Zoltan_Oct_id(oct);
Zoltan_Oct_setID(oct, id+nprevoct); /* now have global id's */
}
return 0;
}