節点-要素対応表の作成

各節点がどの要素に所属しているか調べます。全数検査です。

int i,j,k,l;//index
int **nbe=malloc(sizeof(int*)*nn);//節点要素対応表(nodes by elements)
for(i=0;i<nn;i++)
 nbe[i]=malloc(sizeof(int)*8);//高々8個です
int *noean=malloc(sizeof(int)*nn);//number of elements at the node
for(i=0;i<nn;i++){
 l=0;
 for(j=0;j<noe;j++)
  for(k=0;k<3;k++)
   if(i==en[j][k]){
    ne[i][l]=j;
    l++;
   }
  noean[i]=l;
}