1 // GetDP - Copyright (C) 1997-2021 P. Dular and C. Geuzaine, University of Liege
2 //
3 // See the LICENSE.txt file for license information. Please report all
4 // issues on https://gitlab.onelab.info/getdp/getdp/issues.
5 
6 /* 1 integration point */
7 
8 static double xq1[1] = {0.};
9 static double yq1[1] = {0.};
10 static double pq1[1] = {4.};
11 
12 /* 3 integration points */
13 
14 static double xq3[3] = {0.816496580928,-0.408248290464,-0.408248290464};
15 static double yq3[3] = {0.0,0.840896415255,-0.840896415255};
16 static double pq3[3] = {1.33333333333,1.33333333333,1.33333333333};
17 
18 /* 4 integration points */
19 
20 static double xq4[4] = {0.577350269189,-0.577350269189,0.577350269189,-0.577350269189};
21 static double yq4[4] = {0.577350269189,0.577350269189,-0.577350269189,-0.577350269189};
22 static double pq4[4] = {1.,1.,1.,1.};
23 
24 /* 7 integration points */
25 
26 static double xq7[7] = {0., 0., 0.,
27                         0.7745966692414834, 0.7745966692414834,-0.7745966692414834,
28                         -0.7745966692414834};
29 static double yq7[7] = {0., 0.9660917830792959, -0.9660917830792959,
30                         0.7745966692414834, -0.7745966692414834, 0.7745966692414834,
31                         -0.7745966692414834};
32 static double pq7[7] = {1.1428571428571428, 0.31746031746031744,0.31746031746031744,
33                         0.5555555555555556, 0.5555555555555556, 0.5555555555555556,
34                         0.5555555555555556};
35 
36 /*
37    GAUSS QUADRANGLE WITH 1/R SINGULARITY OVER NODE (-1,-1,0)
38 
39    ref.: H. L. G. Pina, J. L. M. Fernandes, C. A. Brebbia,
40          Some numerical integration formulae over triangles and
41 	 squares with a 1/R singularity,
42 	 Appl. Math. Modelling, Vol 5, June 1981, pp 209--211
43 */
44 
45 /* 1 integration point */
46 
47 static double xqs1[1] = {-0.26501817};
48 static double yqs1[1] = {-0.26501817};
49 static double pqs1[1] = {3.52549435};
50 
51 /* 3 integration points */
52 
53 static double xqs3[3] = {-0.58105530,1.0,-0.21877566};
54 static double yqs3[3] = {-0.58105530,-0.21877566,1.0};
55 static double pqs3[3] = {2.37881900,0.57333767,0.57333767};
56 
57 /* 4 integration points */
58 
59 static double xqs4[4] = {-0.37512304,0.69629093,-0.92928746,-0.15602536};
60 static double yqs4[4] = {-0.92928746,-0.15602536,-0.37512304,0.69629093};
61 static double pqs4[4] = {1.02276580,0.73998134,1.02276580,0.73998134};
62