File: Gauss_Triangle.h

package info (click to toggle)
getdp 2.9.2+dfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 6,384 kB
  • ctags: 8,206
  • sloc: cpp: 55,135; fortran: 13,955; yacc: 8,493; lex: 746; sh: 56; ansic: 34; awk: 33; makefile: 24
file content (123 lines) | stat: -rw-r--r-- 6,304 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
// GetDP - Copyright (C) 1997-2016 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 xt1[1] = {0.333333333333333};
static double yt1[1] = {0.333333333333333};
static double pt1[1] = {0.5};

/* 3 integration points */

static double xt3[3] = {0.16666666666666,0.66666666666666,0.16666666666666};
static double yt3[3] = {0.16666666666666,0.16666666666666,0.66666666666666};
static double pt3[3] = {0.16666666666666,0.16666666666666,0.16666666666666};

/* 4 integration points */

static double xt4[4] = {0.333333333333333,0.6,0.2,0.2};
static double yt4[4] = {0.3333333333333333,0.2,0.6,0.2};
static double pt4[4] = {-0.28125,.260416666666,.260416666666,.260416666666};

/* 6 integration points */

static double xt6[6] = {0.816847572980459,0.091576213509771,0.091576213509771,
			0.108103018168070,0.445948490915965,0.445948490915965};
static double yt6[6] = {0.091576213509771,0.816847572980459,0.091576213509771,
			0.445948490915965,0.108103018168070,0.445948490915965};
static double pt6[6] = {0.054975871827661,0.054975871827661,0.054975871827661,
			0.111690794839,0.111690794839,0.111690794839};

/* 7 integration points */

static double xt7[7] = {0.333333333333333,0.797426985353087,0.101286507323456,
                        0.101286507323456,0.470142064105115,0.059715871789770,
                        0.470142064105115};
static double yt7[7] = {0.333333333333333,0.101286507323456,0.797426985353087,
                        0.101286507323456,0.059715871789770,0.470142064105115,
                        0.470142064105115};
static double pt7[7] = {0.112500000000000,0.062969590272414,0.062969590272414,
                        0.062969590272414,0.066197076394253,0.066197076394253,
                        0.066197076394253};
/* 12 integration points */

static double xt12[12] = {0.873821971016996,0.063089014491502,0.063089014491502,
                          0.501426509658179,0.249286745170910,0.249286745170910,
                          0.636502499121399,0.310352451033785,0.636502499121399,
                          0.310352451033785,0.053145049844816,0.053145049844816};
static double yt12[12] = {0.063089014491502,0.873821971016996,0.063089014491502,
                          0.249286745170910,0.501426509658179,0.249286745170910,
                          0.310352451033785,0.636502499121399,0.053145049844816,
                          0.053145049844816,0.310352451033785,0.636502499121399};
static double pt12[12] = {0.025422453185104,0.025422453185104,0.025422453185104,
                          0.058393137863189,0.058393137863189,0.058393137863189,
                          0.041425537809187,0.041425537809187,0.041425537809187,
                          0.041425537809187,0.041425537809187,0.041425537809187};

/* 13 integration points */

static double xt13[13] = {0.333333333333333,0.479308067841920,0.260345966079040,
                          0.260345966079040,0.869739794195568,0.065130102902216,
                          0.065130102902216,0.048690315425316,0.312865496004874,
                          0.638444188569810,0.048690315425316,0.312865496004874,
                          0.638444188569810};
static double yt13[13] = {0.333333333333333,0.260345966079040,0.479308067841920,
                          0.260345966079040,0.065130102902216,0.869739794195568,
                          0.065130102902216,0.312865496004874,0.048690315425316,
                          0.048690315425316,0.638444188569810,0.638444188569810,
                          0.312865496004874};
static double pt13[13] = {-0.074785022233841,0.087807628716604,0.087807628716604,
                          0.087807628716604,0.026673617804419,0.026673617804419,
                          0.026673617804419,0.038556880445128,0.038556880445128,
                          0.038556880445128,0.038556880445128,0.038556880445128,
                          0.038556880445128};

/* 16 integration points */

static double xt16[16] = {0.333333333333333,0.081414823414554,0.459292588292723,
                          0.459292588292723,0.658861384496480,0.170569307751760,
                          0.170569307751760,0.898905543365938,0.050547228317031,
                          0.050547228317031,0.008394777409958,0.728492392955404,
                          0.263112829634638,0.008394777409958,0.263112829634638,
                          0.728492392955404};
static double yt16[16] = {0.333333333333333,0.459292588292723,0.081414823414554,
                          0.459292588292723,0.170569307751760,0.658861384496480,
                          0.170569307751760,0.050547228317031,0.898905543365938,
                          0.050547228317031,0.728492392955404,0.008394777409958,
                          0.008394777409958,0.263112829634638,0.728492392955404,
                          0.263112829634638};
static double pt16[16] = {0.072157803838894,0.047545817133643,0.047545817133643,
                          0.047545817133643,0.051608685267359,0.051608685267359,
                          0.051608685267359,0.016229248811599,0.016229248811599,
                          0.016229248811599,0.013615157087217,0.013615157087217,
                          0.013615157087217,0.013615157087217,0.013615157087217,
                          0.013615157087217};

/*
   GAUSS TRIANGLE WITH 1/R SINGULARITY OVER NODE (0,0,0)

   ref.: H. L. G. Pina, J. L. M. Fernandes, C. A. Brebbia,
         Some numerical integration formulae over triangles and
	 squares with a 1/R singularity,
	 Appl. Math. Modelling, Vol 5, June 1981, pp 209--211
*/

/* 1 integration point */

static double xts1[1] = {0.25};
static double yts1[1] = {0.25};
static double pts1[1] = {1.24645048};

/* 3 integration points */

static double xts3[3] = {0.16666667,0.81742619,0.18257381};
static double yts3[3] = {0.16666667,0.18257381,0.81742619};
static double pts3[3] = {0.93483790,0.15580629,0.15580629};

/* 4 integration points */

static double xts4[4] = {0.16385495,0.61114353,0.04756957,0.17753138};
static double yts4[4] = {0.04756957,0.17753138,0.16385495,0.61114353};
static double pts4[4] = {0.31161231,0.31161293,0.31161231,0.31161293};