要素-節点対応表の作成

先ほどのメッシュ図を見ていると次の関係が見えてきます。

実装は以下の通りです。

int i,j;//index
int n=3;//分割数
int noe=2*n*n;//要素数(number of elements)
int **ebn=malloc(sizeof(int*)*noe);//要素節点対応表(elements by nodes)
for(i=0;i<noe;i++)
 ebn[i]=malloc(sizeof(int)*3);//ひとつの要素に節点は3つ
int tri[2][3]={{0,n+1,n+2},{0,n+2,1}};
for(i=0;i<noe;i++)
 for(j=0;j<3;j++)
  ebn[i][j]=i/2+i/2/n+tri[i%2][j];//i/2+i/2/n:小正方形の左下の節点番号、tri[i%2][j]:左下以外への番号増分