ParMooN
 All Classes Functions Variables Friends Pages
NF_N_Q_Q5_2D.h
1 
2 static double NF_N_Q_Q5_2D_Xi[60] =
3  {
4  -0.932469514203152027812301554495,
5  -0.661209386466264513661399595021,
6  -0.238619186083196908630501721681,
7  0.238619186083196908630501721681,
8  0.661209386466264513661399595021,
9  0.932469514203152027812301554495,
10  1, 1, 1, 1, 1, 1,
11  0.932469514203152027812301554495,
12  0.661209386466264513661399595021,
13  0.238619186083196908630501721681,
14  -0.238619186083196908630501721681,
15  -0.661209386466264513661399595021,
16  -0.932469514203152027812301554495,
17  -1, -1, -1, -1, -1, -1,
18 
19  -0.932469514203152027812301554495,
20  -0.661209386466264513661399595021,
21  -0.238619186083196908630501721681,
22  0.238619186083196908630501721681,
23  0.661209386466264513661399595021,
24  0.932469514203152027812301554495,
25  -0.932469514203152027812301554495,
26  -0.661209386466264513661399595021,
27  -0.238619186083196908630501721681,
28  0.238619186083196908630501721681,
29  0.661209386466264513661399595021,
30  0.932469514203152027812301554495,
31  -0.932469514203152027812301554495,
32  -0.661209386466264513661399595021,
33  -0.238619186083196908630501721681,
34  0.238619186083196908630501721681,
35  0.661209386466264513661399595021,
36  0.932469514203152027812301554495,
37  -0.932469514203152027812301554495,
38  -0.661209386466264513661399595021,
39  -0.238619186083196908630501721681,
40  0.238619186083196908630501721681,
41  0.661209386466264513661399595021,
42  0.932469514203152027812301554495,
43  -0.932469514203152027812301554495,
44  -0.661209386466264513661399595021,
45  -0.238619186083196908630501721681,
46  0.238619186083196908630501721681,
47  0.661209386466264513661399595021,
48  0.932469514203152027812301554495,
49  -0.932469514203152027812301554495,
50  -0.661209386466264513661399595021,
51  -0.238619186083196908630501721681,
52  0.238619186083196908630501721681,
53  0.661209386466264513661399595021,
54  0.932469514203152027812301554495
55  };
56 
57 static double NF_N_Q_Q5_2D_Eta[60] =
58  {
59  -1, -1, -1, -1, -1, -1,
60  -0.932469514203152027812301554495,
61  -0.661209386466264513661399595021,
62  -0.238619186083196908630501721681,
63  0.238619186083196908630501721681,
64  0.661209386466264513661399595021,
65  0.932469514203152027812301554495,
66  1, 1, 1, 1, 1, 1,
67  0.932469514203152027812301554495,
68  0.661209386466264513661399595021,
69  0.238619186083196908630501721681,
70  -0.238619186083196908630501721681,
71  -0.661209386466264513661399595021,
72  -0.932469514203152027812301554495,
73 
74  -0.932469514203152027812301554495,
75  -0.932469514203152027812301554495,
76  -0.932469514203152027812301554495,
77  -0.932469514203152027812301554495,
78  -0.932469514203152027812301554495,
79  -0.932469514203152027812301554495,
80  -0.661209386466264513661399595021,
81  -0.661209386466264513661399595021,
82  -0.661209386466264513661399595021,
83  -0.661209386466264513661399595021,
84  -0.661209386466264513661399595021,
85  -0.661209386466264513661399595021,
86  -0.238619186083196908630501721681,
87  -0.238619186083196908630501721681,
88  -0.238619186083196908630501721681,
89  -0.238619186083196908630501721681,
90  -0.238619186083196908630501721681,
91  -0.238619186083196908630501721681,
92  0.238619186083196908630501721681,
93  0.238619186083196908630501721681,
94  0.238619186083196908630501721681,
95  0.238619186083196908630501721681,
96  0.238619186083196908630501721681,
97  0.238619186083196908630501721681,
98  0.661209386466264513661399595021,
99  0.661209386466264513661399595021,
100  0.661209386466264513661399595021,
101  0.661209386466264513661399595021,
102  0.661209386466264513661399595021,
103  0.661209386466264513661399595021,
104  0.932469514203152027812301554495,
105  0.932469514203152027812301554495,
106  0.932469514203152027812301554495,
107  0.932469514203152027812301554495,
108  0.932469514203152027812301554495,
109  0.932469514203152027812301554495
110  };
111 
112 static double NF_N_Q_Q5_2D_T[6] =
113  {
114  -0.932469514203152027812301554495,
115  -0.661209386466264513661399595021,
116  -0.238619186083196908630501721681,
117  0.238619186083196908630501721681,
118  0.661209386466264513661399595021,
119  0.932469514203152027812301554495
120  };
121 
122 
123 static double NF_N_Q_Q5_2D_EdgeWeight0[6] = {
124  0.0856622461895851725201480710875,
125  0.180380786524069303784916756918,
126  0.233956967286345523694935171995,
127  0.233956967286345523694935171995,
128  0.180380786524069303784916756918,
129  0.0856622461895851725201480710875 };
130 
131 static double NF_N_Q_Q5_2D_EdgeWeight1[6] = {
132  -0.239632299269859890105518370648,
133  -0.357808407563646294838734602300,
134  -0.167479863337082682628861535781,
135  0.167479863337082682628861535781,
136  0.357808407563646294838734602300,
137  0.239632299269859890105518370648 };
138 
139 static double NF_N_Q_Q5_2D_EdgeWeight2[6] = {
140  0.344468918744913553951585996483,
141  0.140513727783900954708003276278,
142  -0.484982646528814508659589272762,
143  -0.484982646528814508659589272762,
144  0.140513727783900954708003276278,
145  0.344468918744913553951585996483 };
146 
147 static double NF_N_Q_Q5_2D_EdgeWeight3[6] = {
148  -0.376721097993911994360121892502,
149  0.339803199487927794882556709270,
150  0.530551948742438674219510675340,
151  -0.530551948742438674219510675340,
152  -0.339803199487927794882556709270,
153  0.376721097993911994360121892502 };
154 
155 static double NF_N_Q_Q5_2D_EdgeWeight4[6] = {
156  0.32534907297640428362836948428805172,
157  -0.69522592887681074559282921930495380,
158  0.36987685590040646196445973503819405,
159  0.36987685590040646196445973503819405,
160  -0.69522592887681074559282921930495380,
161  0.32534907297640428362836948428805172 };
162 
163 static double NF_N_Q_Q5_2D_CellWeight0[36] = {
164  .00733802042224509940868926568334, .0154518233430958321543579495904,
165  .0200412793294516546970268202701, .0200412793294516546970268202701,
166  .0154518233430958321543579495904, .00733802042224509940868926568334,
167  .0154518233430958321543579495904, .0325372281470418621450407542845,
168  .0422013417718969575930596745842, .0422013417718969575930596745842,
169  .0325372281470418621450407542845, .0154518233430958321543579495904,
170  .0200412793294516546970268202701, .0422013417718969575930596745842,
171  .0547358625418241495573810911433, .0547358625418241495573810911433,
172  .0422013417718969575930596745842, .0200412793294516546970268202701,
173  .0200412793294516546970268202701, .0422013417718969575930596745842,
174  .0547358625418241495573810911433, .0547358625418241495573810911433,
175  .0422013417718969575930596745842, .0200412793294516546970268202701,
176  .0154518233430958321543579495904, .0325372281470418621450407542845,
177  .0422013417718969575930596745842, .0422013417718969575930596745842,
178  .0325372281470418621450407542845, .0154518233430958321543579495904,
179  .00733802042224509940868926568334, .0154518233430958321543579495904,
180  .0200412793294516546970268202701, .0200412793294516546970268202701,
181  .0154518233430958321543579495904, .00733802042224509940868926568334 };
182 
183 static double NF_N_Q_Q5_2D_CellWeight1[36] = {
184  -.0205274410150310890902280205411, -.0306506718974204982425775505983,
185  -.0143467012849792564661051058036, .0143467012849792564661051058036,
186  .0306506718974204982425775505983, .0205274410150310890902280205411,
187  -.0432250626188684853298413153462, -.0645417619812552667292685475461,
188  -.0302101494756866126166405168285, .0302101494756866126166405168285,
189  .0645417619812552667292685475461, .0432250626188684853298413153462,
190  -.0560636460010303706326898494374, -.0837117699031473824475212030064,
191  -.0391830809078754722316851452608, .0391830809078754722316851452608,
192  .0837117699031473824475212030064, .0560636460010303706326898494374,
193  -.0560636460010303706326898494374, -.0837117699031473824475212030064,
194  -.0391830809078754722316851452608, .0391830809078754722316851452608,
195  .0837117699031473824475212030064, .0560636460010303706326898494374,
196  -.0432250626188684853298413153462, -.0645417619812552667292685475461,
197  -.0302101494756866126166405168285, .0302101494756866126166405168285,
198  .0645417619812552667292685475461, .0432250626188684853298413153462,
199  -.0205274410150310890902280205411, -.0306506718974204982425775505983,
200  -.0143467012849792564661051058036, .0143467012849792564661051058036,
201  .0306506718974204982425775505983, .0205274410150310890902280205411 };
202 
203 static double NF_N_Q_Q5_2D_CellWeight2[36] = {
204  -.0205274410150310890902280205411, -.0432250626188684853298413153462,
205  -.0560636460010303706326898494374, -.0560636460010303706326898494374,
206  -.0432250626188684853298413153462, -.0205274410150310890902280205411,
207  -.0306506718974204982425775505983, -.0645417619812552667292685475461,
208  -.0837117699031473824475212030064, -.0837117699031473824475212030064,
209  -.0645417619812552667292685475461, -.0306506718974204982425775505983,
210  -.0143467012849792564661051058036, -.0302101494756866126166405168285,
211  -.0391830809078754722316851452608, -.0391830809078754722316851452608,
212  -.0302101494756866126166405168285, -.0143467012849792564661051058036,
213  .0143467012849792564661051058036, .0302101494756866126166405168285,
214  .0391830809078754722316851452608, .0391830809078754722316851452608,
215  .0302101494756866126166405168285, .0143467012849792564661051058036,
216  .0306506718974204982425775505983, .0645417619812552667292685475461,
217  .0837117699031473824475212030064, .0837117699031473824475212030064,
218  .0645417619812552667292685475461, .0306506718974204982425775505983,
219  .0205274410150310890902280205411, .0432250626188684853298413153462,
220  .0560636460010303706326898494374, .0560636460010303706326898494374,
221  .0432250626188684853298413153462, .0205274410150310890902280205411 };
222 
223 static double NF_N_Q_Q5_2D_CellWeight3[36] = {
224  .0295079813221869954954183356070, .0120367215424408777451918206188,
225  -.0415447028646278732406101562252, -.0415447028646278732406101562252,
226  .0120367215424408777451918206188, .0295079813221869954954183356070,
227  .0621355744963032267860303955680, .0253459767350890238483123189352,
228  -.0874815512313922506343427145015, -.0874815512313922506343427145015,
229  .0253459767350890238483123189352, .0621355744963032267860303955680,
230  .0805909035539665546943442670685, .0328741656144205757604974985854,
231  -.113465069168387130454841765652, -.113465069168387130454841765652,
232  .0328741656144205757604974985854, .0805909035539665546943442670685,
233  .0805909035539665546943442670685, .0328741656144205757604974985854,
234  -.113465069168387130454841765652, -.113465069168387130454841765652,
235  .0328741656144205757604974985854, .0805909035539665546943442670685,
236  .0621355744963032267860303955680, .0253459767350890238483123189352,
237  -.0874815512313922506343427145015, -.0874815512313922506343427145015,
238  .0253459767350890238483123189352, .0621355744963032267860303955680,
239  .0295079813221869954954183356070, .0120367215424408777451918206188,
240  -.0415447028646278732406101562252, -.0415447028646278732406101562252,
241  .0120367215424408777451918206188, .0295079813221869954954183356070 };
242 
243 static double NF_N_Q_Q5_2D_CellWeight4[36] = {
244  .0574236388533596928205697997784, .0857424514025636879989939242949,
245  .0401335847328670327064738215239, -.0401335847328670327064738215239,
246  -.0857424514025636879989939242949, -.0574236388533596928205697997784,
247  .0857424514025636879989939242952, .128026856523232415053097045576,
248  .0599257031996186631287638017479, -.0599257031996186631287638017479,
249  -.128026856523232415053097045576, -.0857424514025636879989939242952,
250  .0401335847328670327064738215239, .0599257031996186631287638017478,
251  .0280495046234078921263331546470, -.0280495046234078921263331546470,
252  -.0599257031996186631287638017478, -.0401335847328670327064738215239,
253  -.0401335847328670327064738215239, -.0599257031996186631287638017478,
254  -.0280495046234078921263331546470, .0280495046234078921263331546470,
255  .0599257031996186631287638017478, .0401335847328670327064738215239,
256  -.0857424514025636879989939242952, -.128026856523232415053097045576,
257  -.0599257031996186631287638017479, .0599257031996186631287638017479,
258  .128026856523232415053097045576, .0857424514025636879989939242952,
259  -.0574236388533596928205697997784, -.0857424514025636879989939242949,
260  -.0401335847328670327064738215239, .0401335847328670327064738215239,
261  .0857424514025636879989939242949, .0574236388533596928205697997784 };
262 
263 static double NF_N_Q_Q5_2D_CellWeight5[36] = {
264  .0295079813221869954954183356070, .0621355744963032267860303955680,
265  .0805909035539665546943442670685, .0805909035539665546943442670685,
266  .0621355744963032267860303955680, .0295079813221869954954183356070,
267  .0120367215424408777451918206188, .0253459767350890238483123189352,
268  .0328741656144205757604974985854, .0328741656144205757604974985854,
269  .0253459767350890238483123189352, .0120367215424408777451918206188,
270  -.0415447028646278732406101562252, -.0874815512313922506343427145015,
271  -.113465069168387130454841765652, -.113465069168387130454841765652,
272  -.0874815512313922506343427145015, -.0415447028646278732406101562252,
273  -.0415447028646278732406101562252, -.0874815512313922506343427145015,
274  -.113465069168387130454841765652, -.113465069168387130454841765652,
275  -.0874815512313922506343427145015, -.0415447028646278732406101562252,
276  .0120367215424408777451918206188, .0253459767350890238483123189352,
277  .0328741656144205757604974985854, .0328741656144205757604974985854,
278  .0253459767350890238483123189352, .0120367215424408777451918206188,
279  .0295079813221869954954183356070, .0621355744963032267860303955680,
280  .0805909035539665546943442670685, .0805909035539665546943442670685,
281  .0621355744963032267860303955680, .0295079813221869954954183356070 };
282 
283 static double NF_N_Q_Q5_2D_CellWeight6[36] = {
284  -.0322707754411653301183827870092, .0291083053305435929932841956179,
285  .0454482716495389550841415739731, -.0454482716495389550841415739731,
286  -.0291083053305435929932841956179, .0322707754411653301183827870092,
287  -.0679532479563528323042652136542, .0612939683870276393316878586185,
288  .0957013778060387900116939496570, -.0957013778060387900116939496570,
289  -.0612939683870276393316878586185, .0679532479563528323042652136542,
290  -.0881365255994378347574129455911, .0794993260263926651163063003988,
291  .124126324915641592013919814046, -.124126324915641592013919814046,
292  -.0794993260263926651163063003988, .0881365255994378347574129455911,
293  -.0881365255994378347574129455911, .0794993260263926651163063003988,
294  .124126324915641592013919814046, -.124126324915641592013919814046,
295  -.0794993260263926651163063003988, .0881365255994378347574129455911,
296  -.0679532479563528323042652136542, .0612939683870276393316878586185,
297  .0957013778060387900116939496570, -.0957013778060387900116939496570,
298  -.0612939683870276393316878586185, .0679532479563528323042652136542,
299  -.0322707754411653301183827870092, .0291083053305435929932841956179,
300  .0454482716495389550841415739731, -.0454482716495389550841415739731,
301  -.0291083053305435929932841956179, .0322707754411653301183827870092 };
302 
303 static double NF_N_Q_Q5_2D_CellWeight7[36] = {
304  -.0825458790258461740469481958082, -.0336716276678353801030539338124,
305  .116217506693681554150002129619, .116217506693681554150002129619,
306  -.0336716276678353801030539338124, -.0825458790258461740469481958082,
307  -.123253875271288587829801198130, -.0502769931791892828892241340110,
308  .173530868450477870719025332138, .173530868450477870719025332138,
309  -.0502769931791892828892241340110, -.123253875271288587829801198130,
310  -.0576916074352707611754535192208, -.0235332199262317698076445226032,
311  .0812248273615025309830980418228, .0812248273615025309830980418228,
312  -.0235332199262317698076445226032, -.0576916074352707611754535192208,
313  .0576916074352707611754535192208, .0235332199262317698076445226032,
314  -.0812248273615025309830980418228, -.0812248273615025309830980418228,
315  .0235332199262317698076445226032, .0576916074352707611754535192208,
316  .123253875271288587829801198130, .0502769931791892828892241340110,
317  -.173530868450477870719025332138, -.173530868450477870719025332138,
318  .0502769931791892828892241340110, .123253875271288587829801198130,
319  .0825458790258461740469481958082, .0336716276678353801030539338124,
320  -.116217506693681554150002129619, -.116217506693681554150002129619,
321  .0336716276678353801030539338124, .0825458790258461740469481958082 };
322 
323 static double NF_N_Q_Q5_2D_CellWeight8[36] = {
324  -.0825458790258461740469481958082, -.123253875271288587829801198130,
325  -.0576916074352707611754535192208, .0576916074352707611754535192208,
326  .123253875271288587829801198130, .0825458790258461740469481958082,
327  -.0336716276678353801030539338122, -.0502769931791892828892241340111,
328  -.0235332199262317698076445226033, .0235332199262317698076445226033,
329  .0502769931791892828892241340111, .0336716276678353801030539338122,
330  .116217506693681554150002129619, .173530868450477870719025332138,
331  .0812248273615025309830980418227, -.0812248273615025309830980418227,
332  -.173530868450477870719025332138, -.116217506693681554150002129619,
333  .116217506693681554150002129619, .173530868450477870719025332138,
334  .0812248273615025309830980418227, -.0812248273615025309830980418227,
335  -.173530868450477870719025332138, -.116217506693681554150002129619,
336  -.0336716276678353801030539338122, -.0502769931791892828892241340111,
337  -.0235332199262317698076445226033, .0235332199262317698076445226033,
338  .0502769931791892828892241340111, .0336716276678353801030539338122,
339  -.0825458790258461740469481958082, -.123253875271288587829801198130,
340  -.0576916074352707611754535192208, .0576916074352707611754535192208,
341  .123253875271288587829801198130, .0825458790258461740469481958082 };
342 
343 static double NF_N_Q_Q5_2D_CellWeight9[36] = {
344  -.0322707754411653301183827870092, -.0679532479563528323042652136542,
345  -.0881365255994378347574129455911, -.0881365255994378347574129455911,
346  -.0679532479563528323042652136542, -.0322707754411653301183827870092,
347  .0291083053305435929932841956179, .0612939683870276393316878586185,
348  .0794993260263926651163063003988, .0794993260263926651163063003988,
349  .0612939683870276393316878586185, .0291083053305435929932841956179,
350  .0454482716495389550841415739731, .0957013778060387900116939496570,
351  .124126324915641592013919814046, .124126324915641592013919814046,
352  .0957013778060387900116939496570, .0454482716495389550841415739731,
353  -.0454482716495389550841415739731, -.0957013778060387900116939496570,
354  -.124126324915641592013919814046, -.124126324915641592013919814046,
355  -.0957013778060387900116939496570, -.0454482716495389550841415739731,
356  -.0291083053305435929932841956179, -.0612939683870276393316878586185,
357  -.0794993260263926651163063003988, -.0794993260263926651163063003988,
358  -.0612939683870276393316878586185, -.0291083053305435929932841956179,
359  .0322707754411653301183827870092, .0679532479563528323042652136542,
360  .0881365255994378347574129455911, .0881365255994378347574129455911,
361  .0679532479563528323042652136542, .0322707754411653301183827870092 };
362 
363 
364 
365 void NF_N_Q_Q5_2D_EvalAll(TCollection *Coll, TBaseCell *Cell, double *PointValues,
366  double *Functionals)
367 {
368  int OwnNum, NeighNum;
369  TBaseCell *neigh;
370 
371  Functionals[0] =( NF_N_Q_Q5_2D_EdgeWeight0[0]*PointValues[0]
372  +NF_N_Q_Q5_2D_EdgeWeight0[1]*PointValues[1]
373  +NF_N_Q_Q5_2D_EdgeWeight0[2]*PointValues[2]
374  +NF_N_Q_Q5_2D_EdgeWeight0[3]*PointValues[3]
375  +NF_N_Q_Q5_2D_EdgeWeight0[4]*PointValues[4]
376  +NF_N_Q_Q5_2D_EdgeWeight0[5]*PointValues[5]);
377  Functionals[1] =( NF_N_Q_Q5_2D_EdgeWeight0[0]*PointValues[6]
378  +NF_N_Q_Q5_2D_EdgeWeight0[1]*PointValues[7]
379  +NF_N_Q_Q5_2D_EdgeWeight0[2]*PointValues[8]
380  +NF_N_Q_Q5_2D_EdgeWeight0[3]*PointValues[9]
381  +NF_N_Q_Q5_2D_EdgeWeight0[4]*PointValues[10]
382  +NF_N_Q_Q5_2D_EdgeWeight0[5]*PointValues[11]);
383  Functionals[2] =( NF_N_Q_Q5_2D_EdgeWeight0[0]*PointValues[12]
384  +NF_N_Q_Q5_2D_EdgeWeight0[1]*PointValues[13]
385  +NF_N_Q_Q5_2D_EdgeWeight0[2]*PointValues[14]
386  +NF_N_Q_Q5_2D_EdgeWeight0[3]*PointValues[15]
387  +NF_N_Q_Q5_2D_EdgeWeight0[4]*PointValues[16]
388  +NF_N_Q_Q5_2D_EdgeWeight0[5]*PointValues[17]);
389  Functionals[3] =( NF_N_Q_Q5_2D_EdgeWeight0[0]*PointValues[18]
390  +NF_N_Q_Q5_2D_EdgeWeight0[1]*PointValues[19]
391  +NF_N_Q_Q5_2D_EdgeWeight0[2]*PointValues[20]
392  +NF_N_Q_Q5_2D_EdgeWeight0[3]*PointValues[21]
393  +NF_N_Q_Q5_2D_EdgeWeight0[4]*PointValues[22]
394  +NF_N_Q_Q5_2D_EdgeWeight0[5]*PointValues[23]);
395 
396  Functionals[4] =( NF_N_Q_Q5_2D_EdgeWeight1[0]*PointValues[0]
397  +NF_N_Q_Q5_2D_EdgeWeight1[1]*PointValues[1]
398  +NF_N_Q_Q5_2D_EdgeWeight1[2]*PointValues[2]
399  +NF_N_Q_Q5_2D_EdgeWeight1[3]*PointValues[3]
400  +NF_N_Q_Q5_2D_EdgeWeight1[4]*PointValues[4]
401  +NF_N_Q_Q5_2D_EdgeWeight1[5]*PointValues[5]);
402  Functionals[5] =( NF_N_Q_Q5_2D_EdgeWeight1[0]*PointValues[6]
403  +NF_N_Q_Q5_2D_EdgeWeight1[1]*PointValues[7]
404  +NF_N_Q_Q5_2D_EdgeWeight1[2]*PointValues[8]
405  +NF_N_Q_Q5_2D_EdgeWeight1[3]*PointValues[9]
406  +NF_N_Q_Q5_2D_EdgeWeight1[4]*PointValues[10]
407  +NF_N_Q_Q5_2D_EdgeWeight1[5]*PointValues[11]);
408  Functionals[6] =( NF_N_Q_Q5_2D_EdgeWeight1[0]*PointValues[12]
409  +NF_N_Q_Q5_2D_EdgeWeight1[1]*PointValues[13]
410  +NF_N_Q_Q5_2D_EdgeWeight1[2]*PointValues[14]
411  +NF_N_Q_Q5_2D_EdgeWeight1[3]*PointValues[15]
412  +NF_N_Q_Q5_2D_EdgeWeight1[4]*PointValues[16]
413  +NF_N_Q_Q5_2D_EdgeWeight1[5]*PointValues[17]);
414  Functionals[7] =( NF_N_Q_Q5_2D_EdgeWeight1[0]*PointValues[18]
415  +NF_N_Q_Q5_2D_EdgeWeight1[1]*PointValues[19]
416  +NF_N_Q_Q5_2D_EdgeWeight1[2]*PointValues[20]
417  +NF_N_Q_Q5_2D_EdgeWeight1[3]*PointValues[21]
418  +NF_N_Q_Q5_2D_EdgeWeight1[4]*PointValues[22]
419  +NF_N_Q_Q5_2D_EdgeWeight1[5]*PointValues[23]);
420 
421  Functionals[8] =( NF_N_Q_Q5_2D_EdgeWeight2[0]*PointValues[0]
422  +NF_N_Q_Q5_2D_EdgeWeight2[1]*PointValues[1]
423  +NF_N_Q_Q5_2D_EdgeWeight2[2]*PointValues[2]
424  +NF_N_Q_Q5_2D_EdgeWeight2[3]*PointValues[3]
425  +NF_N_Q_Q5_2D_EdgeWeight2[4]*PointValues[4]
426  +NF_N_Q_Q5_2D_EdgeWeight2[5]*PointValues[5]);
427  Functionals[9] =( NF_N_Q_Q5_2D_EdgeWeight2[0]*PointValues[6]
428  +NF_N_Q_Q5_2D_EdgeWeight2[1]*PointValues[7]
429  +NF_N_Q_Q5_2D_EdgeWeight2[2]*PointValues[8]
430  +NF_N_Q_Q5_2D_EdgeWeight2[3]*PointValues[9]
431  +NF_N_Q_Q5_2D_EdgeWeight2[4]*PointValues[10]
432  +NF_N_Q_Q5_2D_EdgeWeight2[5]*PointValues[11]);
433  Functionals[10]=( NF_N_Q_Q5_2D_EdgeWeight2[0]*PointValues[12]
434  +NF_N_Q_Q5_2D_EdgeWeight2[1]*PointValues[13]
435  +NF_N_Q_Q5_2D_EdgeWeight2[2]*PointValues[14]
436  +NF_N_Q_Q5_2D_EdgeWeight2[3]*PointValues[15]
437  +NF_N_Q_Q5_2D_EdgeWeight2[4]*PointValues[16]
438  +NF_N_Q_Q5_2D_EdgeWeight2[5]*PointValues[17]);
439  Functionals[11]=( NF_N_Q_Q5_2D_EdgeWeight2[0]*PointValues[18]
440  +NF_N_Q_Q5_2D_EdgeWeight2[1]*PointValues[19]
441  +NF_N_Q_Q5_2D_EdgeWeight2[2]*PointValues[20]
442  +NF_N_Q_Q5_2D_EdgeWeight2[3]*PointValues[21]
443  +NF_N_Q_Q5_2D_EdgeWeight2[4]*PointValues[22]
444  +NF_N_Q_Q5_2D_EdgeWeight2[5]*PointValues[23]);
445 
446  Functionals[12]=( NF_N_Q_Q5_2D_EdgeWeight3[0]*PointValues[0]
447  +NF_N_Q_Q5_2D_EdgeWeight3[1]*PointValues[1]
448  +NF_N_Q_Q5_2D_EdgeWeight3[2]*PointValues[2]
449  +NF_N_Q_Q5_2D_EdgeWeight3[3]*PointValues[3]
450  +NF_N_Q_Q5_2D_EdgeWeight3[4]*PointValues[4]
451  +NF_N_Q_Q5_2D_EdgeWeight3[5]*PointValues[5]);
452  Functionals[13]=( NF_N_Q_Q5_2D_EdgeWeight3[0]*PointValues[6]
453  +NF_N_Q_Q5_2D_EdgeWeight3[1]*PointValues[7]
454  +NF_N_Q_Q5_2D_EdgeWeight3[2]*PointValues[8]
455  +NF_N_Q_Q5_2D_EdgeWeight3[3]*PointValues[9]
456  +NF_N_Q_Q5_2D_EdgeWeight3[4]*PointValues[10]
457  +NF_N_Q_Q5_2D_EdgeWeight3[5]*PointValues[11]);
458  Functionals[14]=( NF_N_Q_Q5_2D_EdgeWeight3[0]*PointValues[12]
459  +NF_N_Q_Q5_2D_EdgeWeight3[1]*PointValues[13]
460  +NF_N_Q_Q5_2D_EdgeWeight3[2]*PointValues[14]
461  +NF_N_Q_Q5_2D_EdgeWeight3[3]*PointValues[15]
462  +NF_N_Q_Q5_2D_EdgeWeight3[4]*PointValues[16]
463  +NF_N_Q_Q5_2D_EdgeWeight3[5]*PointValues[17]);
464  Functionals[15]=( NF_N_Q_Q5_2D_EdgeWeight3[0]*PointValues[18]
465  +NF_N_Q_Q5_2D_EdgeWeight3[1]*PointValues[19]
466  +NF_N_Q_Q5_2D_EdgeWeight3[2]*PointValues[20]
467  +NF_N_Q_Q5_2D_EdgeWeight3[3]*PointValues[21]
468  +NF_N_Q_Q5_2D_EdgeWeight3[4]*PointValues[22]
469  +NF_N_Q_Q5_2D_EdgeWeight3[5]*PointValues[23]);
470 
471  Functionals[16]=( NF_N_Q_Q5_2D_EdgeWeight4[0]*PointValues[0]
472  +NF_N_Q_Q5_2D_EdgeWeight4[1]*PointValues[1]
473  +NF_N_Q_Q5_2D_EdgeWeight4[2]*PointValues[2]
474  +NF_N_Q_Q5_2D_EdgeWeight4[3]*PointValues[3]
475  +NF_N_Q_Q5_2D_EdgeWeight4[4]*PointValues[4]
476  +NF_N_Q_Q5_2D_EdgeWeight4[5]*PointValues[5]);
477  Functionals[17]=( NF_N_Q_Q5_2D_EdgeWeight4[0]*PointValues[6]
478  +NF_N_Q_Q5_2D_EdgeWeight4[1]*PointValues[7]
479  +NF_N_Q_Q5_2D_EdgeWeight4[2]*PointValues[8]
480  +NF_N_Q_Q5_2D_EdgeWeight4[3]*PointValues[9]
481  +NF_N_Q_Q5_2D_EdgeWeight4[4]*PointValues[10]
482  +NF_N_Q_Q5_2D_EdgeWeight4[5]*PointValues[11]);
483  Functionals[18]=( NF_N_Q_Q5_2D_EdgeWeight4[0]*PointValues[12]
484  +NF_N_Q_Q5_2D_EdgeWeight4[1]*PointValues[13]
485  +NF_N_Q_Q5_2D_EdgeWeight4[2]*PointValues[14]
486  +NF_N_Q_Q5_2D_EdgeWeight4[3]*PointValues[15]
487  +NF_N_Q_Q5_2D_EdgeWeight4[4]*PointValues[16]
488  +NF_N_Q_Q5_2D_EdgeWeight4[5]*PointValues[17]);
489  Functionals[19]=( NF_N_Q_Q5_2D_EdgeWeight4[0]*PointValues[18]
490  +NF_N_Q_Q5_2D_EdgeWeight4[1]*PointValues[19]
491  +NF_N_Q_Q5_2D_EdgeWeight4[2]*PointValues[20]
492  +NF_N_Q_Q5_2D_EdgeWeight4[3]*PointValues[21]
493  +NF_N_Q_Q5_2D_EdgeWeight4[4]*PointValues[22]
494  +NF_N_Q_Q5_2D_EdgeWeight4[5]*PointValues[23]);
495 
496  Functionals[20] =( NF_N_Q_Q5_2D_CellWeight0[ 0]*PointValues[24]
497  +NF_N_Q_Q5_2D_CellWeight0[ 1]*PointValues[25]
498  +NF_N_Q_Q5_2D_CellWeight0[ 2]*PointValues[26]
499  +NF_N_Q_Q5_2D_CellWeight0[ 3]*PointValues[27]
500  +NF_N_Q_Q5_2D_CellWeight0[ 4]*PointValues[28]
501  +NF_N_Q_Q5_2D_CellWeight0[ 5]*PointValues[29]
502  +NF_N_Q_Q5_2D_CellWeight0[ 6]*PointValues[30]
503  +NF_N_Q_Q5_2D_CellWeight0[ 7]*PointValues[31]
504  +NF_N_Q_Q5_2D_CellWeight0[ 8]*PointValues[32]
505  +NF_N_Q_Q5_2D_CellWeight0[ 9]*PointValues[33]
506  +NF_N_Q_Q5_2D_CellWeight0[10]*PointValues[34]
507  +NF_N_Q_Q5_2D_CellWeight0[11]*PointValues[35]
508  +NF_N_Q_Q5_2D_CellWeight0[12]*PointValues[36]
509  +NF_N_Q_Q5_2D_CellWeight0[13]*PointValues[37]
510  +NF_N_Q_Q5_2D_CellWeight0[14]*PointValues[38]
511  +NF_N_Q_Q5_2D_CellWeight0[15]*PointValues[39]
512  +NF_N_Q_Q5_2D_CellWeight0[16]*PointValues[40]
513  +NF_N_Q_Q5_2D_CellWeight0[17]*PointValues[41]
514  +NF_N_Q_Q5_2D_CellWeight0[18]*PointValues[42]
515  +NF_N_Q_Q5_2D_CellWeight0[19]*PointValues[43]
516  +NF_N_Q_Q5_2D_CellWeight0[20]*PointValues[44]
517  +NF_N_Q_Q5_2D_CellWeight0[21]*PointValues[45]
518  +NF_N_Q_Q5_2D_CellWeight0[22]*PointValues[46]
519  +NF_N_Q_Q5_2D_CellWeight0[23]*PointValues[47]
520  +NF_N_Q_Q5_2D_CellWeight0[24]*PointValues[48]
521  +NF_N_Q_Q5_2D_CellWeight0[25]*PointValues[49]
522  +NF_N_Q_Q5_2D_CellWeight0[26]*PointValues[50]
523  +NF_N_Q_Q5_2D_CellWeight0[27]*PointValues[51]
524  +NF_N_Q_Q5_2D_CellWeight0[28]*PointValues[52]
525  +NF_N_Q_Q5_2D_CellWeight0[29]*PointValues[53]
526  +NF_N_Q_Q5_2D_CellWeight0[30]*PointValues[54]
527  +NF_N_Q_Q5_2D_CellWeight0[31]*PointValues[55]
528  +NF_N_Q_Q5_2D_CellWeight0[32]*PointValues[56]
529  +NF_N_Q_Q5_2D_CellWeight0[33]*PointValues[57]
530  +NF_N_Q_Q5_2D_CellWeight0[34]*PointValues[58]
531  +NF_N_Q_Q5_2D_CellWeight0[35]*PointValues[59] );
532  Functionals[21]= ( NF_N_Q_Q5_2D_CellWeight1[ 0]*PointValues[24]
533  +NF_N_Q_Q5_2D_CellWeight1[ 1]*PointValues[25]
534  +NF_N_Q_Q5_2D_CellWeight1[ 2]*PointValues[26]
535  +NF_N_Q_Q5_2D_CellWeight1[ 3]*PointValues[27]
536  +NF_N_Q_Q5_2D_CellWeight1[ 4]*PointValues[28]
537  +NF_N_Q_Q5_2D_CellWeight1[ 5]*PointValues[29]
538  +NF_N_Q_Q5_2D_CellWeight1[ 6]*PointValues[30]
539  +NF_N_Q_Q5_2D_CellWeight1[ 7]*PointValues[31]
540  +NF_N_Q_Q5_2D_CellWeight1[ 8]*PointValues[32]
541  +NF_N_Q_Q5_2D_CellWeight1[ 9]*PointValues[33]
542  +NF_N_Q_Q5_2D_CellWeight1[10]*PointValues[34]
543  +NF_N_Q_Q5_2D_CellWeight1[11]*PointValues[35]
544  +NF_N_Q_Q5_2D_CellWeight1[12]*PointValues[36]
545  +NF_N_Q_Q5_2D_CellWeight1[13]*PointValues[37]
546  +NF_N_Q_Q5_2D_CellWeight1[14]*PointValues[38]
547  +NF_N_Q_Q5_2D_CellWeight1[15]*PointValues[39]
548  +NF_N_Q_Q5_2D_CellWeight1[16]*PointValues[40]
549  +NF_N_Q_Q5_2D_CellWeight1[17]*PointValues[41]
550  +NF_N_Q_Q5_2D_CellWeight1[18]*PointValues[42]
551  +NF_N_Q_Q5_2D_CellWeight1[19]*PointValues[43]
552  +NF_N_Q_Q5_2D_CellWeight1[20]*PointValues[44]
553  +NF_N_Q_Q5_2D_CellWeight1[21]*PointValues[45]
554  +NF_N_Q_Q5_2D_CellWeight1[22]*PointValues[46]
555  +NF_N_Q_Q5_2D_CellWeight1[23]*PointValues[47]
556  +NF_N_Q_Q5_2D_CellWeight1[24]*PointValues[48]
557  +NF_N_Q_Q5_2D_CellWeight1[25]*PointValues[49]
558  +NF_N_Q_Q5_2D_CellWeight1[26]*PointValues[50]
559  +NF_N_Q_Q5_2D_CellWeight1[27]*PointValues[51]
560  +NF_N_Q_Q5_2D_CellWeight1[28]*PointValues[52]
561  +NF_N_Q_Q5_2D_CellWeight1[29]*PointValues[53]
562  +NF_N_Q_Q5_2D_CellWeight1[30]*PointValues[54]
563  +NF_N_Q_Q5_2D_CellWeight1[31]*PointValues[55]
564  +NF_N_Q_Q5_2D_CellWeight1[32]*PointValues[56]
565  +NF_N_Q_Q5_2D_CellWeight1[33]*PointValues[57]
566  +NF_N_Q_Q5_2D_CellWeight1[34]*PointValues[58]
567  +NF_N_Q_Q5_2D_CellWeight1[35]*PointValues[59] );
568  Functionals[22] =( NF_N_Q_Q5_2D_CellWeight2[ 0]*PointValues[24]
569  +NF_N_Q_Q5_2D_CellWeight2[ 1]*PointValues[25]
570  +NF_N_Q_Q5_2D_CellWeight2[ 2]*PointValues[26]
571  +NF_N_Q_Q5_2D_CellWeight2[ 3]*PointValues[27]
572  +NF_N_Q_Q5_2D_CellWeight2[ 4]*PointValues[28]
573  +NF_N_Q_Q5_2D_CellWeight2[ 5]*PointValues[29]
574  +NF_N_Q_Q5_2D_CellWeight2[ 6]*PointValues[30]
575  +NF_N_Q_Q5_2D_CellWeight2[ 7]*PointValues[31]
576  +NF_N_Q_Q5_2D_CellWeight2[ 8]*PointValues[32]
577  +NF_N_Q_Q5_2D_CellWeight2[ 9]*PointValues[33]
578  +NF_N_Q_Q5_2D_CellWeight2[10]*PointValues[34]
579  +NF_N_Q_Q5_2D_CellWeight2[11]*PointValues[35]
580  +NF_N_Q_Q5_2D_CellWeight2[12]*PointValues[36]
581  +NF_N_Q_Q5_2D_CellWeight2[13]*PointValues[37]
582  +NF_N_Q_Q5_2D_CellWeight2[14]*PointValues[38]
583  +NF_N_Q_Q5_2D_CellWeight2[15]*PointValues[39]
584  +NF_N_Q_Q5_2D_CellWeight2[16]*PointValues[40]
585  +NF_N_Q_Q5_2D_CellWeight2[17]*PointValues[41]
586  +NF_N_Q_Q5_2D_CellWeight2[18]*PointValues[42]
587  +NF_N_Q_Q5_2D_CellWeight2[19]*PointValues[43]
588  +NF_N_Q_Q5_2D_CellWeight2[20]*PointValues[44]
589  +NF_N_Q_Q5_2D_CellWeight2[21]*PointValues[45]
590  +NF_N_Q_Q5_2D_CellWeight2[22]*PointValues[46]
591  +NF_N_Q_Q5_2D_CellWeight2[23]*PointValues[47]
592  +NF_N_Q_Q5_2D_CellWeight2[24]*PointValues[48]
593  +NF_N_Q_Q5_2D_CellWeight2[25]*PointValues[49]
594  +NF_N_Q_Q5_2D_CellWeight2[26]*PointValues[50]
595  +NF_N_Q_Q5_2D_CellWeight2[27]*PointValues[51]
596  +NF_N_Q_Q5_2D_CellWeight2[28]*PointValues[52]
597  +NF_N_Q_Q5_2D_CellWeight2[29]*PointValues[53]
598  +NF_N_Q_Q5_2D_CellWeight2[30]*PointValues[54]
599  +NF_N_Q_Q5_2D_CellWeight2[31]*PointValues[55]
600  +NF_N_Q_Q5_2D_CellWeight2[32]*PointValues[56]
601  +NF_N_Q_Q5_2D_CellWeight2[33]*PointValues[57]
602  +NF_N_Q_Q5_2D_CellWeight2[34]*PointValues[58]
603  +NF_N_Q_Q5_2D_CellWeight2[35]*PointValues[59] );
604  Functionals[23] =( NF_N_Q_Q5_2D_CellWeight3[ 0]*PointValues[24]
605  +NF_N_Q_Q5_2D_CellWeight3[ 1]*PointValues[25]
606  +NF_N_Q_Q5_2D_CellWeight3[ 2]*PointValues[26]
607  +NF_N_Q_Q5_2D_CellWeight3[ 3]*PointValues[27]
608  +NF_N_Q_Q5_2D_CellWeight3[ 4]*PointValues[28]
609  +NF_N_Q_Q5_2D_CellWeight3[ 5]*PointValues[29]
610  +NF_N_Q_Q5_2D_CellWeight3[ 6]*PointValues[30]
611  +NF_N_Q_Q5_2D_CellWeight3[ 7]*PointValues[31]
612  +NF_N_Q_Q5_2D_CellWeight3[ 8]*PointValues[32]
613  +NF_N_Q_Q5_2D_CellWeight3[ 9]*PointValues[33]
614  +NF_N_Q_Q5_2D_CellWeight3[10]*PointValues[34]
615  +NF_N_Q_Q5_2D_CellWeight3[11]*PointValues[35]
616  +NF_N_Q_Q5_2D_CellWeight3[12]*PointValues[36]
617  +NF_N_Q_Q5_2D_CellWeight3[13]*PointValues[37]
618  +NF_N_Q_Q5_2D_CellWeight3[14]*PointValues[38]
619  +NF_N_Q_Q5_2D_CellWeight3[15]*PointValues[39]
620  +NF_N_Q_Q5_2D_CellWeight3[16]*PointValues[40]
621  +NF_N_Q_Q5_2D_CellWeight3[17]*PointValues[41]
622  +NF_N_Q_Q5_2D_CellWeight3[18]*PointValues[42]
623  +NF_N_Q_Q5_2D_CellWeight3[19]*PointValues[43]
624  +NF_N_Q_Q5_2D_CellWeight3[20]*PointValues[44]
625  +NF_N_Q_Q5_2D_CellWeight3[21]*PointValues[45]
626  +NF_N_Q_Q5_2D_CellWeight3[22]*PointValues[46]
627  +NF_N_Q_Q5_2D_CellWeight3[23]*PointValues[47]
628  +NF_N_Q_Q5_2D_CellWeight3[24]*PointValues[48]
629  +NF_N_Q_Q5_2D_CellWeight3[25]*PointValues[49]
630  +NF_N_Q_Q5_2D_CellWeight3[26]*PointValues[50]
631  +NF_N_Q_Q5_2D_CellWeight3[27]*PointValues[51]
632  +NF_N_Q_Q5_2D_CellWeight3[28]*PointValues[52]
633  +NF_N_Q_Q5_2D_CellWeight3[29]*PointValues[53]
634  +NF_N_Q_Q5_2D_CellWeight3[30]*PointValues[54]
635  +NF_N_Q_Q5_2D_CellWeight3[31]*PointValues[55]
636  +NF_N_Q_Q5_2D_CellWeight3[32]*PointValues[56]
637  +NF_N_Q_Q5_2D_CellWeight3[33]*PointValues[57]
638  +NF_N_Q_Q5_2D_CellWeight3[34]*PointValues[58]
639  +NF_N_Q_Q5_2D_CellWeight3[35]*PointValues[59] );
640  Functionals[24]= ( NF_N_Q_Q5_2D_CellWeight4[ 0]*PointValues[24]
641  +NF_N_Q_Q5_2D_CellWeight4[ 1]*PointValues[25]
642  +NF_N_Q_Q5_2D_CellWeight4[ 2]*PointValues[26]
643  +NF_N_Q_Q5_2D_CellWeight4[ 3]*PointValues[27]
644  +NF_N_Q_Q5_2D_CellWeight4[ 4]*PointValues[28]
645  +NF_N_Q_Q5_2D_CellWeight4[ 5]*PointValues[29]
646  +NF_N_Q_Q5_2D_CellWeight4[ 6]*PointValues[30]
647  +NF_N_Q_Q5_2D_CellWeight4[ 7]*PointValues[31]
648  +NF_N_Q_Q5_2D_CellWeight4[ 8]*PointValues[32]
649  +NF_N_Q_Q5_2D_CellWeight4[ 9]*PointValues[33]
650  +NF_N_Q_Q5_2D_CellWeight4[10]*PointValues[34]
651  +NF_N_Q_Q5_2D_CellWeight4[11]*PointValues[35]
652  +NF_N_Q_Q5_2D_CellWeight4[12]*PointValues[36]
653  +NF_N_Q_Q5_2D_CellWeight4[13]*PointValues[37]
654  +NF_N_Q_Q5_2D_CellWeight4[14]*PointValues[38]
655  +NF_N_Q_Q5_2D_CellWeight4[15]*PointValues[39]
656  +NF_N_Q_Q5_2D_CellWeight4[16]*PointValues[40]
657  +NF_N_Q_Q5_2D_CellWeight4[17]*PointValues[41]
658  +NF_N_Q_Q5_2D_CellWeight4[18]*PointValues[42]
659  +NF_N_Q_Q5_2D_CellWeight4[19]*PointValues[43]
660  +NF_N_Q_Q5_2D_CellWeight4[20]*PointValues[44]
661  +NF_N_Q_Q5_2D_CellWeight4[21]*PointValues[45]
662  +NF_N_Q_Q5_2D_CellWeight4[22]*PointValues[46]
663  +NF_N_Q_Q5_2D_CellWeight4[23]*PointValues[47]
664  +NF_N_Q_Q5_2D_CellWeight4[24]*PointValues[48]
665  +NF_N_Q_Q5_2D_CellWeight4[25]*PointValues[49]
666  +NF_N_Q_Q5_2D_CellWeight4[26]*PointValues[50]
667  +NF_N_Q_Q5_2D_CellWeight4[27]*PointValues[51]
668  +NF_N_Q_Q5_2D_CellWeight4[28]*PointValues[52]
669  +NF_N_Q_Q5_2D_CellWeight4[29]*PointValues[53]
670  +NF_N_Q_Q5_2D_CellWeight4[30]*PointValues[54]
671  +NF_N_Q_Q5_2D_CellWeight4[31]*PointValues[55]
672  +NF_N_Q_Q5_2D_CellWeight4[32]*PointValues[56]
673  +NF_N_Q_Q5_2D_CellWeight4[33]*PointValues[57]
674  +NF_N_Q_Q5_2D_CellWeight4[34]*PointValues[58]
675  +NF_N_Q_Q5_2D_CellWeight4[35]*PointValues[59] );
676  Functionals[25] =( NF_N_Q_Q5_2D_CellWeight5[ 0]*PointValues[24]
677  +NF_N_Q_Q5_2D_CellWeight5[ 1]*PointValues[25]
678  +NF_N_Q_Q5_2D_CellWeight5[ 2]*PointValues[26]
679  +NF_N_Q_Q5_2D_CellWeight5[ 3]*PointValues[27]
680  +NF_N_Q_Q5_2D_CellWeight5[ 4]*PointValues[28]
681  +NF_N_Q_Q5_2D_CellWeight5[ 5]*PointValues[29]
682  +NF_N_Q_Q5_2D_CellWeight5[ 6]*PointValues[30]
683  +NF_N_Q_Q5_2D_CellWeight5[ 7]*PointValues[31]
684  +NF_N_Q_Q5_2D_CellWeight5[ 8]*PointValues[32]
685  +NF_N_Q_Q5_2D_CellWeight5[ 9]*PointValues[33]
686  +NF_N_Q_Q5_2D_CellWeight5[10]*PointValues[34]
687  +NF_N_Q_Q5_2D_CellWeight5[11]*PointValues[35]
688  +NF_N_Q_Q5_2D_CellWeight5[12]*PointValues[36]
689  +NF_N_Q_Q5_2D_CellWeight5[13]*PointValues[37]
690  +NF_N_Q_Q5_2D_CellWeight5[14]*PointValues[38]
691  +NF_N_Q_Q5_2D_CellWeight5[15]*PointValues[39]
692  +NF_N_Q_Q5_2D_CellWeight5[16]*PointValues[40]
693  +NF_N_Q_Q5_2D_CellWeight5[17]*PointValues[41]
694  +NF_N_Q_Q5_2D_CellWeight5[18]*PointValues[42]
695  +NF_N_Q_Q5_2D_CellWeight5[19]*PointValues[43]
696  +NF_N_Q_Q5_2D_CellWeight5[20]*PointValues[44]
697  +NF_N_Q_Q5_2D_CellWeight5[21]*PointValues[45]
698  +NF_N_Q_Q5_2D_CellWeight5[22]*PointValues[46]
699  +NF_N_Q_Q5_2D_CellWeight5[23]*PointValues[47]
700  +NF_N_Q_Q5_2D_CellWeight5[24]*PointValues[48]
701  +NF_N_Q_Q5_2D_CellWeight5[25]*PointValues[49]
702  +NF_N_Q_Q5_2D_CellWeight5[26]*PointValues[50]
703  +NF_N_Q_Q5_2D_CellWeight5[27]*PointValues[51]
704  +NF_N_Q_Q5_2D_CellWeight5[28]*PointValues[52]
705  +NF_N_Q_Q5_2D_CellWeight5[29]*PointValues[53]
706  +NF_N_Q_Q5_2D_CellWeight5[30]*PointValues[54]
707  +NF_N_Q_Q5_2D_CellWeight5[31]*PointValues[55]
708  +NF_N_Q_Q5_2D_CellWeight5[32]*PointValues[56]
709  +NF_N_Q_Q5_2D_CellWeight5[33]*PointValues[57]
710  +NF_N_Q_Q5_2D_CellWeight5[34]*PointValues[58]
711  +NF_N_Q_Q5_2D_CellWeight5[35]*PointValues[59] );
712  Functionals[26]= ( NF_N_Q_Q5_2D_CellWeight6[ 0]*PointValues[24]
713  +NF_N_Q_Q5_2D_CellWeight6[ 1]*PointValues[25]
714  +NF_N_Q_Q5_2D_CellWeight6[ 2]*PointValues[26]
715  +NF_N_Q_Q5_2D_CellWeight6[ 3]*PointValues[27]
716  +NF_N_Q_Q5_2D_CellWeight6[ 4]*PointValues[28]
717  +NF_N_Q_Q5_2D_CellWeight6[ 5]*PointValues[29]
718  +NF_N_Q_Q5_2D_CellWeight6[ 6]*PointValues[30]
719  +NF_N_Q_Q5_2D_CellWeight6[ 7]*PointValues[31]
720  +NF_N_Q_Q5_2D_CellWeight6[ 8]*PointValues[32]
721  +NF_N_Q_Q5_2D_CellWeight6[ 9]*PointValues[33]
722  +NF_N_Q_Q5_2D_CellWeight6[10]*PointValues[34]
723  +NF_N_Q_Q5_2D_CellWeight6[11]*PointValues[35]
724  +NF_N_Q_Q5_2D_CellWeight6[12]*PointValues[36]
725  +NF_N_Q_Q5_2D_CellWeight6[13]*PointValues[37]
726  +NF_N_Q_Q5_2D_CellWeight6[14]*PointValues[38]
727  +NF_N_Q_Q5_2D_CellWeight6[15]*PointValues[39]
728  +NF_N_Q_Q5_2D_CellWeight6[16]*PointValues[40]
729  +NF_N_Q_Q5_2D_CellWeight6[17]*PointValues[41]
730  +NF_N_Q_Q5_2D_CellWeight6[18]*PointValues[42]
731  +NF_N_Q_Q5_2D_CellWeight6[19]*PointValues[43]
732  +NF_N_Q_Q5_2D_CellWeight6[20]*PointValues[44]
733  +NF_N_Q_Q5_2D_CellWeight6[21]*PointValues[45]
734  +NF_N_Q_Q5_2D_CellWeight6[22]*PointValues[46]
735  +NF_N_Q_Q5_2D_CellWeight6[23]*PointValues[47]
736  +NF_N_Q_Q5_2D_CellWeight6[24]*PointValues[48]
737  +NF_N_Q_Q5_2D_CellWeight6[25]*PointValues[49]
738  +NF_N_Q_Q5_2D_CellWeight6[26]*PointValues[50]
739  +NF_N_Q_Q5_2D_CellWeight6[27]*PointValues[51]
740  +NF_N_Q_Q5_2D_CellWeight6[28]*PointValues[52]
741  +NF_N_Q_Q5_2D_CellWeight6[29]*PointValues[53]
742  +NF_N_Q_Q5_2D_CellWeight6[30]*PointValues[54]
743  +NF_N_Q_Q5_2D_CellWeight6[31]*PointValues[55]
744  +NF_N_Q_Q5_2D_CellWeight6[32]*PointValues[56]
745  +NF_N_Q_Q5_2D_CellWeight6[33]*PointValues[57]
746  +NF_N_Q_Q5_2D_CellWeight6[34]*PointValues[58]
747  +NF_N_Q_Q5_2D_CellWeight6[35]*PointValues[59] );
748  Functionals[27] =( NF_N_Q_Q5_2D_CellWeight7[ 0]*PointValues[24]
749  +NF_N_Q_Q5_2D_CellWeight7[ 1]*PointValues[25]
750  +NF_N_Q_Q5_2D_CellWeight7[ 2]*PointValues[26]
751  +NF_N_Q_Q5_2D_CellWeight7[ 3]*PointValues[27]
752  +NF_N_Q_Q5_2D_CellWeight7[ 4]*PointValues[28]
753  +NF_N_Q_Q5_2D_CellWeight7[ 5]*PointValues[29]
754  +NF_N_Q_Q5_2D_CellWeight7[ 6]*PointValues[30]
755  +NF_N_Q_Q5_2D_CellWeight7[ 7]*PointValues[31]
756  +NF_N_Q_Q5_2D_CellWeight7[ 8]*PointValues[32]
757  +NF_N_Q_Q5_2D_CellWeight7[ 9]*PointValues[33]
758  +NF_N_Q_Q5_2D_CellWeight7[10]*PointValues[34]
759  +NF_N_Q_Q5_2D_CellWeight7[11]*PointValues[35]
760  +NF_N_Q_Q5_2D_CellWeight7[12]*PointValues[36]
761  +NF_N_Q_Q5_2D_CellWeight7[13]*PointValues[37]
762  +NF_N_Q_Q5_2D_CellWeight7[14]*PointValues[38]
763  +NF_N_Q_Q5_2D_CellWeight7[15]*PointValues[39]
764  +NF_N_Q_Q5_2D_CellWeight7[16]*PointValues[40]
765  +NF_N_Q_Q5_2D_CellWeight7[17]*PointValues[41]
766  +NF_N_Q_Q5_2D_CellWeight7[18]*PointValues[42]
767  +NF_N_Q_Q5_2D_CellWeight7[19]*PointValues[43]
768  +NF_N_Q_Q5_2D_CellWeight7[20]*PointValues[44]
769  +NF_N_Q_Q5_2D_CellWeight7[21]*PointValues[45]
770  +NF_N_Q_Q5_2D_CellWeight7[22]*PointValues[46]
771  +NF_N_Q_Q5_2D_CellWeight7[23]*PointValues[47]
772  +NF_N_Q_Q5_2D_CellWeight7[24]*PointValues[48]
773  +NF_N_Q_Q5_2D_CellWeight7[25]*PointValues[49]
774  +NF_N_Q_Q5_2D_CellWeight7[26]*PointValues[50]
775  +NF_N_Q_Q5_2D_CellWeight7[27]*PointValues[51]
776  +NF_N_Q_Q5_2D_CellWeight7[28]*PointValues[52]
777  +NF_N_Q_Q5_2D_CellWeight7[29]*PointValues[53]
778  +NF_N_Q_Q5_2D_CellWeight7[30]*PointValues[54]
779  +NF_N_Q_Q5_2D_CellWeight7[31]*PointValues[55]
780  +NF_N_Q_Q5_2D_CellWeight7[32]*PointValues[56]
781  +NF_N_Q_Q5_2D_CellWeight7[33]*PointValues[57]
782  +NF_N_Q_Q5_2D_CellWeight7[34]*PointValues[58]
783  +NF_N_Q_Q5_2D_CellWeight7[35]*PointValues[59] );
784  Functionals[28]= ( NF_N_Q_Q5_2D_CellWeight8[ 0]*PointValues[24]
785  +NF_N_Q_Q5_2D_CellWeight8[ 1]*PointValues[25]
786  +NF_N_Q_Q5_2D_CellWeight8[ 2]*PointValues[26]
787  +NF_N_Q_Q5_2D_CellWeight8[ 3]*PointValues[27]
788  +NF_N_Q_Q5_2D_CellWeight8[ 4]*PointValues[28]
789  +NF_N_Q_Q5_2D_CellWeight8[ 5]*PointValues[29]
790  +NF_N_Q_Q5_2D_CellWeight8[ 6]*PointValues[30]
791  +NF_N_Q_Q5_2D_CellWeight8[ 7]*PointValues[31]
792  +NF_N_Q_Q5_2D_CellWeight8[ 8]*PointValues[32]
793  +NF_N_Q_Q5_2D_CellWeight8[ 9]*PointValues[33]
794  +NF_N_Q_Q5_2D_CellWeight8[10]*PointValues[34]
795  +NF_N_Q_Q5_2D_CellWeight8[11]*PointValues[35]
796  +NF_N_Q_Q5_2D_CellWeight8[12]*PointValues[36]
797  +NF_N_Q_Q5_2D_CellWeight8[13]*PointValues[37]
798  +NF_N_Q_Q5_2D_CellWeight8[14]*PointValues[38]
799  +NF_N_Q_Q5_2D_CellWeight8[15]*PointValues[39]
800  +NF_N_Q_Q5_2D_CellWeight8[16]*PointValues[40]
801  +NF_N_Q_Q5_2D_CellWeight8[17]*PointValues[41]
802  +NF_N_Q_Q5_2D_CellWeight8[18]*PointValues[42]
803  +NF_N_Q_Q5_2D_CellWeight8[19]*PointValues[43]
804  +NF_N_Q_Q5_2D_CellWeight8[20]*PointValues[44]
805  +NF_N_Q_Q5_2D_CellWeight8[21]*PointValues[45]
806  +NF_N_Q_Q5_2D_CellWeight8[22]*PointValues[46]
807  +NF_N_Q_Q5_2D_CellWeight8[23]*PointValues[47]
808  +NF_N_Q_Q5_2D_CellWeight8[24]*PointValues[48]
809  +NF_N_Q_Q5_2D_CellWeight8[25]*PointValues[49]
810  +NF_N_Q_Q5_2D_CellWeight8[26]*PointValues[50]
811  +NF_N_Q_Q5_2D_CellWeight8[27]*PointValues[51]
812  +NF_N_Q_Q5_2D_CellWeight8[28]*PointValues[52]
813  +NF_N_Q_Q5_2D_CellWeight8[29]*PointValues[53]
814  +NF_N_Q_Q5_2D_CellWeight8[30]*PointValues[54]
815  +NF_N_Q_Q5_2D_CellWeight8[31]*PointValues[55]
816  +NF_N_Q_Q5_2D_CellWeight8[32]*PointValues[56]
817  +NF_N_Q_Q5_2D_CellWeight8[33]*PointValues[57]
818  +NF_N_Q_Q5_2D_CellWeight8[34]*PointValues[58]
819  +NF_N_Q_Q5_2D_CellWeight8[35]*PointValues[59] );
820  Functionals[29] =( NF_N_Q_Q5_2D_CellWeight9[ 0]*PointValues[24]
821  +NF_N_Q_Q5_2D_CellWeight9[ 1]*PointValues[25]
822  +NF_N_Q_Q5_2D_CellWeight9[ 2]*PointValues[26]
823  +NF_N_Q_Q5_2D_CellWeight9[ 3]*PointValues[27]
824  +NF_N_Q_Q5_2D_CellWeight9[ 4]*PointValues[28]
825  +NF_N_Q_Q5_2D_CellWeight9[ 5]*PointValues[29]
826  +NF_N_Q_Q5_2D_CellWeight9[ 6]*PointValues[30]
827  +NF_N_Q_Q5_2D_CellWeight9[ 7]*PointValues[31]
828  +NF_N_Q_Q5_2D_CellWeight9[ 8]*PointValues[32]
829  +NF_N_Q_Q5_2D_CellWeight9[ 9]*PointValues[33]
830  +NF_N_Q_Q5_2D_CellWeight9[10]*PointValues[34]
831  +NF_N_Q_Q5_2D_CellWeight9[11]*PointValues[35]
832  +NF_N_Q_Q5_2D_CellWeight9[12]*PointValues[36]
833  +NF_N_Q_Q5_2D_CellWeight9[13]*PointValues[37]
834  +NF_N_Q_Q5_2D_CellWeight9[14]*PointValues[38]
835  +NF_N_Q_Q5_2D_CellWeight9[15]*PointValues[39]
836  +NF_N_Q_Q5_2D_CellWeight9[16]*PointValues[40]
837  +NF_N_Q_Q5_2D_CellWeight9[17]*PointValues[41]
838  +NF_N_Q_Q5_2D_CellWeight9[18]*PointValues[42]
839  +NF_N_Q_Q5_2D_CellWeight9[19]*PointValues[43]
840  +NF_N_Q_Q5_2D_CellWeight9[20]*PointValues[44]
841  +NF_N_Q_Q5_2D_CellWeight9[21]*PointValues[45]
842  +NF_N_Q_Q5_2D_CellWeight9[22]*PointValues[46]
843  +NF_N_Q_Q5_2D_CellWeight9[23]*PointValues[47]
844  +NF_N_Q_Q5_2D_CellWeight9[24]*PointValues[48]
845  +NF_N_Q_Q5_2D_CellWeight9[25]*PointValues[49]
846  +NF_N_Q_Q5_2D_CellWeight9[26]*PointValues[50]
847  +NF_N_Q_Q5_2D_CellWeight9[27]*PointValues[51]
848  +NF_N_Q_Q5_2D_CellWeight9[28]*PointValues[52]
849  +NF_N_Q_Q5_2D_CellWeight9[29]*PointValues[53]
850  +NF_N_Q_Q5_2D_CellWeight9[30]*PointValues[54]
851  +NF_N_Q_Q5_2D_CellWeight9[31]*PointValues[55]
852  +NF_N_Q_Q5_2D_CellWeight9[32]*PointValues[56]
853  +NF_N_Q_Q5_2D_CellWeight9[33]*PointValues[57]
854  +NF_N_Q_Q5_2D_CellWeight9[34]*PointValues[58]
855  +NF_N_Q_Q5_2D_CellWeight9[35]*PointValues[59] );
856  /*
857  if(Cell)
858  {
859  if(Cell->GetVertex(0) > Cell->GetVertex(1))
860  {
861  Functionals[4] = -Functionals[4];
862  Functionals[12] = -Functionals[12];
863  }
864  if(Cell->GetVertex(1) > Cell->GetVertex(2))
865  {
866  Functionals[5] = -Functionals[5];
867  Functionals[13] = -Functionals[13];
868  }
869  if(Cell->GetVertex(2) > Cell->GetVertex(3))
870  {
871  Functionals[6] = -Functionals[6];
872  Functionals[14] = -Functionals[14];
873  }
874  if(Cell->GetVertex(3) > Cell->GetVertex(0))
875  {
876  Functionals[7] = -Functionals[7];
877  Functionals[15] = -Functionals[15];
878  }
879  }
880  */
881 
882  if(Cell)
883  {
884  OwnNum = Coll->GetIndex(Cell);
885 
886  neigh = Cell->GetJoint(0)->GetNeighbour(Cell);
887  if(neigh)
888  {
889  NeighNum = Coll->GetIndex(neigh);
890  if(NeighNum < OwnNum)
891  {
892  Functionals[ 4] = -Functionals[ 4];
893  Functionals[12] = -Functionals[12];
894  }
895  } // endif neigh
896 
897  neigh = Cell->GetJoint(1)->GetNeighbour(Cell);
898  if(neigh)
899  {
900  NeighNum = Coll->GetIndex(neigh);
901  if(NeighNum < OwnNum)
902  {
903  Functionals[ 5] = -Functionals[ 5];
904  Functionals[13] = -Functionals[13];
905  }
906  } // endif neigh
907 
908  neigh = Cell->GetJoint(2)->GetNeighbour(Cell);
909  if(neigh)
910  {
911  NeighNum = Coll->GetIndex(neigh);
912  if(NeighNum < OwnNum)
913  {
914  Functionals[ 6] = -Functionals[ 6];
915  Functionals[14] = -Functionals[14];
916  }
917  } // endif neigh
918 
919  neigh = Cell->GetJoint(3)->GetNeighbour(Cell);
920  if(neigh)
921  {
922  NeighNum = Coll->GetIndex(neigh);
923  if(NeighNum < OwnNum)
924  {
925  Functionals[ 7] = -Functionals[ 7];
926  Functionals[15] = -Functionals[15];
927  }
928  } // endif neigh
929  } // endif Cell
930 }
931 
932 void NF_N_Q_Q5_2D_EvalEdge(TCollection *Coll, TBaseCell *Cell, int Joint,
933  double *PointValues, double *Functionals)
934 {
935  int OwnNum, NeighNum;
936  TBaseCell *neigh;
937 
938  Functionals[0] =( NF_N_Q_Q5_2D_EdgeWeight0[0]*PointValues[0]
939  +NF_N_Q_Q5_2D_EdgeWeight0[1]*PointValues[1]
940  +NF_N_Q_Q5_2D_EdgeWeight0[2]*PointValues[2]
941  +NF_N_Q_Q5_2D_EdgeWeight0[3]*PointValues[3]
942  +NF_N_Q_Q5_2D_EdgeWeight0[4]*PointValues[4]
943  +NF_N_Q_Q5_2D_EdgeWeight0[5]*PointValues[5]);
944  Functionals[1] =( NF_N_Q_Q5_2D_EdgeWeight1[0]*PointValues[0]
945  +NF_N_Q_Q5_2D_EdgeWeight1[1]*PointValues[1]
946  +NF_N_Q_Q5_2D_EdgeWeight1[2]*PointValues[2]
947  +NF_N_Q_Q5_2D_EdgeWeight1[3]*PointValues[3]
948  +NF_N_Q_Q5_2D_EdgeWeight1[4]*PointValues[4]
949  +NF_N_Q_Q5_2D_EdgeWeight1[5]*PointValues[5]);
950  Functionals[2] =( NF_N_Q_Q5_2D_EdgeWeight2[0]*PointValues[0]
951  +NF_N_Q_Q5_2D_EdgeWeight2[1]*PointValues[1]
952  +NF_N_Q_Q5_2D_EdgeWeight2[2]*PointValues[2]
953  +NF_N_Q_Q5_2D_EdgeWeight2[3]*PointValues[3]
954  +NF_N_Q_Q5_2D_EdgeWeight2[4]*PointValues[4]
955  +NF_N_Q_Q5_2D_EdgeWeight2[5]*PointValues[5]);
956  Functionals[3] =( NF_N_Q_Q5_2D_EdgeWeight3[0]*PointValues[0]
957  +NF_N_Q_Q5_2D_EdgeWeight3[1]*PointValues[1]
958  +NF_N_Q_Q5_2D_EdgeWeight3[2]*PointValues[2]
959  +NF_N_Q_Q5_2D_EdgeWeight3[3]*PointValues[3]
960  +NF_N_Q_Q5_2D_EdgeWeight3[4]*PointValues[4]
961  +NF_N_Q_Q5_2D_EdgeWeight3[5]*PointValues[5]);
962  Functionals[4] =( NF_N_Q_Q5_2D_EdgeWeight4[0]*PointValues[0]
963  +NF_N_Q_Q5_2D_EdgeWeight4[1]*PointValues[1]
964  +NF_N_Q_Q5_2D_EdgeWeight4[2]*PointValues[2]
965  +NF_N_Q_Q5_2D_EdgeWeight4[3]*PointValues[3]
966  +NF_N_Q_Q5_2D_EdgeWeight4[4]*PointValues[4]
967  +NF_N_Q_Q5_2D_EdgeWeight4[5]*PointValues[5]);
968 
969  if(Joint != -1)
970  {
971  // if(Cell->GetVertex(Joint) > Cell->GetVertex((Joint+1)%4))
972  // {
973  // Functionals[1] = -Functionals[1];
974  // Functionals[3] = -Functionals[3];
975  // }
976  neigh = Cell->GetJoint(Joint)->GetNeighbour(Cell);
977  if(neigh)
978  {
979  OwnNum = Coll->GetIndex(Cell);
980  NeighNum = Coll->GetIndex(neigh);
981  if(NeighNum < OwnNum)
982  {
983  Functionals[1] = -Functionals[1];
984  Functionals[3] = -Functionals[3];
985  }
986  } // endif neigh
987  // */
988  }
989 }
990 
991 /*
992  TNodalFunctional2D(NodalFunctional2D id,
993  int n_allfunctionals, int n_edgefunctionals,
994  int n_pointsall, int n_pointsedge,
995  double *xi, double *eta, double *t,
996  DoubleFunctVect *evalall,
997  DoubleFunctVect *evaledge);
998 */
999 
1000 TNodalFunctional2D *NF_N_Q_Q5_2D_Obj = new TNodalFunctional2D
1001  (NF_N_Q_Q5_2D, 30, 5, 60, 6, NF_N_Q_Q5_2D_Xi, NF_N_Q_Q5_2D_Eta,
1002  NF_N_Q_Q5_2D_T, NF_N_Q_Q5_2D_EvalAll, NF_N_Q_Q5_2D_EvalEdge);
store cells in an array, used by cell iterators
Definition: Collection.h:18
Definition: NodalFunctional2D.h:20
int GetIndex(TBaseCell *cell)
return Index of cell in Cells-array
Definition: Collection.C:115
TJoint * GetJoint(int J_i)
return the pointer to face with number i
Definition: BaseCell.h:175
TBaseCell * GetNeighbour(TBaseCell *Me) const
Definition: Joint.C:57
information for finite element data structure
Definition: BaseCell.h:25