Select Git revision
Gauss_Tetrahedron.h
Gauss_Tetrahedron.h 5.03 KiB
// GetDP - Copyright (C) 1997-2018 P. Dular and C. Geuzaine, University of Liege
//
// See the LICENSE.txt file for license information. Please report all
// bugs and problems to the public mailing list <getdp@onelab.info>.
/* 1 integration point */
static double xtet1[1] = {0.25};
static double ytet1[1] = {0.25};
static double ztet1[1] = {0.25};
static double ptet1[1] = {0.166666666667};
/* 4 integration points */
static double xtet4[4] = {0.138196601125,0.138196601125,0.138196601125,0.585410196625};
static double ytet4[4] = {0.138196601125,0.138196601125,0.585410196625,0.138196601125};
static double ztet4[4] = {0.138196601125,0.585410196625,0.138196601125,0.138196601125};
static double ptet4[4] = {0.0416666666667,0.0416666666667,0.0416666666667,0.0416666666667};
/* 5 integration points */
static double xtet5[5] = {0.25,0.166666666667,0.166666666667,0.166666666667,0.5};
static double ytet5[5] = {0.25,0.166666666667,0.166666666667,0.5,0.166666666667};
static double ztet5[5] = {0.25,0.166666666667,0.5,0.166666666667,0.166666666667};
static double ptet5[5] = {-0.133333333333,0.075,0.075,0.075,0.075};
/* 15 integration points */
static double xtet15[15] = {0.25,
0.0919710780526,0.0919710780526,0.0919710780526,0.724086765841,
0.319793627829,0.319793627829,0.319793627829,0.0406191165118,
0.0563508326895,0.0563508326895,0.44364916731,0.0563508326895,
0.44364916731,0.44364916731};
static double ytet15[15] = {0.25,
0.0919710780526,0.0919710780526,0.724086765841,0.0919710780526,
0.319793627829,0.319793627829,0.0406191165118,0.319793627829,
0.0563508326895,0.44364916731,0.0563508326895,0.44364916731,
0.0563508326895,0.44364916731};
static double ztet15[15] = {0.25,
0.0919710780526,0.724086765841,0.0919710780526,0.0919710780526,
0.319793627829,0.0406191165118,0.319793627829,0.319793627829,
0.44364916731,0.0563508326895,0.0563508326895,0.44364916731,
0.44364916731,0.0563508326895};
static double ptet15[15] = {0.0197530864198,
0.0119895139632,0.0119895139632,0.0119895139632,0.0119895139632,
0.011511367871,0.011511367871,0.011511367871,0.011511367871,
0.00881834215168,0.00881834215168,0.00881834215168,0.00881834215168,
0.00881834215168,0.00881834215168};
/* 16 integration points */
#define a16 0.0503737941001228 / 6.0
#define b16 0.0665420686332923 / 6.0
#define c16 0.7716429020672371
#define d16 0.0761190326442543
#define e16 0.1197005277978019
#define f16 0.0718316452676693
#define g16 0.4042339134672644
static double xtet16[16] = {c16,d16,d16,d16,e16,f16,e16,f16,
g16,g16,g16,g16,e16,f16,g16,g16};
static double ytet16[16] = {d16,c16,d16,d16,f16,e16,g16,g16,
g16,g16,e16,f16,g16,g16,e16,f16};
static double ztet16[16] = {d16,d16,c16,d16,g16,g16,g16,g16,
e16,f16,f16,e16,f16,e16,g16,g16};
static double ptet16[16] = {a16,a16,a16,a16,b16,b16,b16,b16,
b16,b16,b16,b16,b16,b16,b16,b16};
#undef a16
#undef b16
#undef c16