Gauss_Quadrangle.h (getdp-3.4.0-source.tgz) | : | Gauss_Quadrangle.h (getdp-3.5.0-source.tgz) | ||
---|---|---|---|---|
// GetDP - Copyright (C) 1997-2021 P. Dular and C. Geuzaine, University of Liege | // GetDP - Copyright (C) 1997-2022 P. Dular and C. Geuzaine, University of Liege | |||
// | // | |||
// See the LICENSE.txt file for license information. Please report all | // See the LICENSE.txt file for license information. Please report all | |||
// issues on https://gitlab.onelab.info/getdp/getdp/issues. | // issues on https://gitlab.onelab.info/getdp/getdp/issues. | |||
/* 1 integration point */ | /* 1 integration point */ | |||
static double xq1[1] = {0.}; | static double xq1[1] = {0.}; | |||
static double yq1[1] = {0.}; | static double yq1[1] = {0.}; | |||
static double pq1[1] = {4.}; | static double pq1[1] = {4.}; | |||
/* 3 integration points */ | /* 3 integration points */ | |||
static double xq3[3] = {0.816496580928,-0.408248290464,-0.408248290464}; | static double xq3[3] = {0.816496580928, -0.408248290464, -0.408248290464}; | |||
static double yq3[3] = {0.0,0.840896415255,-0.840896415255}; | static double yq3[3] = {0.0, 0.840896415255, -0.840896415255}; | |||
static double pq3[3] = {1.33333333333,1.33333333333,1.33333333333}; | static double pq3[3] = {1.33333333333, 1.33333333333, 1.33333333333}; | |||
/* 4 integration points */ | /* 4 integration points */ | |||
static double xq4[4] = {0.577350269189,-0.577350269189,0.577350269189,-0.5773502 | static double xq4[4] = {0.577350269189, -0.577350269189, 0.577350269189, | |||
69189}; | -0.577350269189}; | |||
static double yq4[4] = {0.577350269189,0.577350269189,-0.577350269189,-0.5773502 | static double yq4[4] = {0.577350269189, 0.577350269189, -0.577350269189, | |||
69189}; | -0.577350269189}; | |||
static double pq4[4] = {1.,1.,1.,1.}; | static double pq4[4] = {1., 1., 1., 1.}; | |||
/* 7 integration points */ | /* 7 integration points */ | |||
static double xq7[7] = {0., 0., 0., | static double xq7[7] = {0., | |||
0.7745966692414834, 0.7745966692414834,-0.77459666924148 | 0., | |||
34, | 0., | |||
0.7745966692414834, | ||||
0.7745966692414834, | ||||
-0.7745966692414834, | ||||
-0.7745966692414834}; | -0.7745966692414834}; | |||
static double yq7[7] = {0., 0.9660917830792959, -0.9660917830792959, | static double yq7[7] = {0., | |||
0.7745966692414834, -0.7745966692414834, 0.7745966692414 | 0.9660917830792959, | |||
834, | -0.9660917830792959, | |||
0.7745966692414834, | ||||
-0.7745966692414834, | ||||
0.7745966692414834, | ||||
-0.7745966692414834}; | -0.7745966692414834}; | |||
static double pq7[7] = {1.1428571428571428, 0.31746031746031744,0.31746031746031 | static double pq7[7] = {1.1428571428571428, 0.31746031746031744, | |||
744, | 0.31746031746031744, 0.5555555555555556, | |||
0.5555555555555556, 0.5555555555555556, 0.55555555555555 | 0.5555555555555556, 0.5555555555555556, | |||
56, | ||||
0.5555555555555556}; | 0.5555555555555556}; | |||
/* | /* | |||
GAUSS QUADRANGLE WITH 1/R SINGULARITY OVER NODE (-1,-1,0) | GAUSS QUADRANGLE WITH 1/R SINGULARITY OVER NODE (-1,-1,0) | |||
ref.: H. L. G. Pina, J. L. M. Fernandes, C. A. Brebbia, | ref.: H. L. G. Pina, J. L. M. Fernandes, C. A. Brebbia, | |||
Some numerical integration formulae over triangles and | Some numerical integration formulae over triangles and | |||
squares with a 1/R singularity, | squares with a 1/R singularity, | |||
Appl. Math. Modelling, Vol 5, June 1981, pp 209--211 | Appl. Math. Modelling, Vol 5, June 1981, pp 209--211 | |||
*/ | */ | |||
/* 1 integration point */ | /* 1 integration point */ | |||
static double xqs1[1] = {-0.26501817}; | static double xqs1[1] = {-0.26501817}; | |||
static double yqs1[1] = {-0.26501817}; | static double yqs1[1] = {-0.26501817}; | |||
static double pqs1[1] = {3.52549435}; | static double pqs1[1] = {3.52549435}; | |||
/* 3 integration points */ | /* 3 integration points */ | |||
static double xqs3[3] = {-0.58105530,1.0,-0.21877566}; | static double xqs3[3] = {-0.58105530, 1.0, -0.21877566}; | |||
static double yqs3[3] = {-0.58105530,-0.21877566,1.0}; | static double yqs3[3] = {-0.58105530, -0.21877566, 1.0}; | |||
static double pqs3[3] = {2.37881900,0.57333767,0.57333767}; | static double pqs3[3] = {2.37881900, 0.57333767, 0.57333767}; | |||
/* 4 integration points */ | /* 4 integration points */ | |||
static double xqs4[4] = {-0.37512304,0.69629093,-0.92928746,-0.15602536}; | static double xqs4[4] = {-0.37512304, 0.69629093, -0.92928746, -0.15602536}; | |||
static double yqs4[4] = {-0.92928746,-0.15602536,-0.37512304,0.69629093}; | static double yqs4[4] = {-0.92928746, -0.15602536, -0.37512304, 0.69629093}; | |||
static double pqs4[4] = {1.02276580,0.73998134,1.02276580,0.73998134}; | static double pqs4[4] = {1.02276580, 0.73998134, 1.02276580, 0.73998134}; | |||
End of changes. 9 change blocks. | ||||
24 lines changed or deleted | 29 lines changed or added |