GetFEM に新しい近似積分法を組み込む

perlスクリプトは記述ファイルから新しい立方体法を自動的に組み込みます．このような記述ファイル（拡張子は .IM ）はディレクトリ cubature にあります．たとえば， IM_TETRAHEDRON(5) の場合，次のファイルは積分法を記述しています．

NAME = IM_TETRAHEDRON(5)
N = 3
GEOTRANS = GT_PK(3,1)
NBPT = 4
0, 0.25, 0.25, 0.25, 0.008818342151675485
1, 0.31979362782962991, 0.31979362782962991, 0.31979362782962991, 0.011511367871045398
1, 0.091971078052723033, 0.091971078052723033, 0.091971078052723033, 0.01198951396316977
1, 0.056350832689629156, 0.056350832689629156, 0.44364916731037084, 0.008818342151675485
NBF = 4 IM_TRIANGLE(5)
IM_TRIANGLE(5)
IM_TRIANGLE(5)
IM_TRIANGLE(5)

ここで， NAME は， GetFEM （定数整数パラメータが許される）の手法名です．N は次元， GEOTRANS は GetFEM の有効な幾何変換を記述します．この幾何学変換は，積分法が記述されている基準要素を定義するだけです． NBPT は積分節点定義の数です．積分節点定義には積分節点の総数が NBPT より大きい対称定義が含まれています．

積分節点定義の構成は次の通りです．

• 整数：0 = 対称性なし，1 = 完全対称（3角形の場合はx6，4角形の場合はx4，4面体の場合はx24 ...），

• 積分節点の N 座標，

• 荷重．

NBF は参照要素の面の数です（ GEOTRANS に対応する必要があります）．次に，既存の積分法を各面（各面に1つずつ）に設定します．これは，境界を積分するために必要です．

ファイルフォーマットは [EncyclopCubature] からインスピレーションを受けています．