NodeConnectivityx.h File Reference

header file for node connectivity computation More...

Go to the source code of this file.

Functions
void	NodeConnectivityx (int **pconnectivity, int *pwidth, int *elements, int nels, int nods)

Detailed Description

header file for node connectivity computation

Definition in file NodeConnectivityx.h.

Function Documentation

NodeConnectivityx()

void NodeConnectivityx	(	int **	pconnectivity,
		int *	pwidth,
		int *	elements,
		int	nels,
		int	nods
	)

Definition at line 18 of file NodeConnectivityx.cpp.

                                                                                          {
 
   int i,j,n;
   const int maxels=50;
   const int width=maxels+1;
 
   /*intermediary: */
   int     node;
   int     index;
   int     num_elements;
   int     already_plugged=0;
   int element;
 
   /*Allocate connectivity: */
   int* connectivity=xNewZeroInit<int>(nods*width);
 
   /*Go through all elements, and for each elements, plug into the connectivity, all the nodes. 
    * If nodes are already plugged into the connectivity, skip them.: */
   for(n=0;n<nels;n++){
 
      element=n+1; //matlab indexing
 
      for(i=0;i<3;i++){
 
         node=elements[n*3+i]; //already matlab indexed, elements comes directly from the workspace.
         index=node-1;
 
         num_elements=connectivity[width*index+maxels]; //retrieve number of elements already  plugged into the connectivity of this node.
 
         already_plugged=0;
         for(j=0;j<num_elements;j++){
            if (element==*(connectivity+width*index+j)){
               already_plugged=1;
               break;
            }
         }
         if(already_plugged)break;
 
         /*this elements is not yet plugged  into the connectivity for this node, do it, and increase counter: */
         connectivity[width*index+num_elements]=element;
         connectivity[width*index+maxels]=num_elements+1;
 
      }
   }
 
   /*Last check: is the number of elements on last column of the connectivity superior to maxels? If so, then error out and 
    * warn the user to increase the connectivity width: */
   for(i=0;i<nods;i++){
      if (*(connectivity+width*i+maxels)>maxels)
       _error_("max connectivity width reached (" << *(connectivity+width*i+maxels) << ")! increase width of connectivity table");
   }
 
   /*Assign output pointers: */
   *pconnectivity=connectivity;
   *pwidth=width;
}