ParMooN
 All Classes Functions Variables Friends Pages
NF_D_H_Q1_3D.h
1 /*
2  TNodalFunctional3D(NodalFunctional3D id,
3  int n_allfunctionals, int *n_facefunctionals,
4  int n_pointsall, int *n_pointsface,
5  double *xi, double *eta, double *zeta,
6  double **xiarray, double **etaarray,
7  double **zetaarray,
8  DoubleFunctVect *evalall,
9  DoubleFunctVect *evalface);
10 */
11 
12 /* for all functionals */
13 #include <stdlib.h>
14 
15 
16 static double NF_D_H_Q1_3D_Xi[] = {
17  -0.7745966692414833770358530, 0,
18  0.7745966692414833770358530,
19  -0.7745966692414833770358530, 0,
20  0.7745966692414833770358530,
21  -0.7745966692414833770358530, 0,
22  0.7745966692414833770358530,
23  -0.7745966692414833770358530, 0,
24  0.7745966692414833770358530,
25  -0.7745966692414833770358530, 0,
26  0.7745966692414833770358530,
27  -0.7745966692414833770358530, 0,
28  0.7745966692414833770358530,
29  -0.7745966692414833770358530, 0,
30  0.7745966692414833770358530,
31  -0.7745966692414833770358530, 0,
32  0.7745966692414833770358530,
33  -0.7745966692414833770358530, 0,
34  0.7745966692414833770358530 } ;
35 static double NF_D_H_Q1_3D_Eta[] = {
36  -0.7745966692414833770358530, -0.7745966692414833770358530,
37  -0.7745966692414833770358530, 0, 0, 0,
38  0.7745966692414833770358530,
39  0.7745966692414833770358530, 0.7745966692414833770358530,
40  -0.7745966692414833770358530, -0.7745966692414833770358530,
41  -0.7745966692414833770358530, 0, 0, 0,
42  0.7745966692414833770358530,
43  0.7745966692414833770358530, 0.7745966692414833770358530,
44  -0.7745966692414833770358530, -0.7745966692414833770358530,
45  -0.7745966692414833770358530, 0, 0, 0,
46  0.7745966692414833770358530,
47  0.7745966692414833770358530, 0.7745966692414833770358530 };
48 static double NF_D_H_Q1_3D_Zeta[] = {
49  -0.7745966692414833770358530, -0.7745966692414833770358530,
50  -0.7745966692414833770358530, -0.7745966692414833770358530,
51  -0.7745966692414833770358530, -0.7745966692414833770358530,
52  -0.7745966692414833770358530, -0.7745966692414833770358530,
53  -0.7745966692414833770358530, 0, 0, 0, 0, 0, 0, 0, 0, 0,
54  0.7745966692414833770358530, 0.7745966692414833770358530,
55  0.7745966692414833770358530, 0.7745966692414833770358530,
56  0.7745966692414833770358530, 0.7745966692414833770358530,
57  0.7745966692414833770358530, 0.7745966692414833770358530,
58  0.7745966692414833770358530 };
59 
60 /* face 0 0 */
61 static double *NF_D_H_Q1_3D_F0_Xi = NULL;
62 static double *NF_D_H_Q1_3D_F0_Eta = NULL;
63 static double *NF_D_H_Q1_3D_F0_Zeta = NULL;
64 
65 /* face 1 1 */
66 static double *NF_D_H_Q1_3D_F1_Xi = NULL;
67 static double *NF_D_H_Q1_3D_F1_Eta = NULL;
68 static double *NF_D_H_Q1_3D_F1_Zeta = NULL;
69 
70 /* face 2 2 */
71 static double *NF_D_H_Q1_3D_F2_Xi = NULL;
72 static double *NF_D_H_Q1_3D_F2_Eta = NULL;
73 static double *NF_D_H_Q1_3D_F2_Zeta = NULL;
74 
75 /* face 3 3 */
76 static double *NF_D_H_Q1_3D_F3_Xi = NULL;
77 static double *NF_D_H_Q1_3D_F3_Eta = NULL;
78 static double *NF_D_H_Q1_3D_F3_Zeta = NULL;
79 
80 /* face 4 4 */
81 static double *NF_D_H_Q1_3D_F4_Xi = NULL;
82 static double *NF_D_H_Q1_3D_F4_Eta = NULL;
83 static double *NF_D_H_Q1_3D_F4_Zeta = NULL;
84 
85 /* face 5 5 */
86 static double *NF_D_H_Q1_3D_F5_Xi = NULL;
87 static double *NF_D_H_Q1_3D_F5_Eta = NULL;
88 static double *NF_D_H_Q1_3D_F5_Zeta = NULL;
89 
90 static double *NF_D_H_Q1_3D_XiArray[6] = {
91  NF_D_H_Q1_3D_F0_Xi,
92  NF_D_H_Q1_3D_F1_Xi,
93  NF_D_H_Q1_3D_F2_Xi,
94  NF_D_H_Q1_3D_F3_Xi,
95  NF_D_H_Q1_3D_F4_Xi,
96  NF_D_H_Q1_3D_F5_Xi };
97 
98 static double *NF_D_H_Q1_3D_EtaArray[6] = {
99  NF_D_H_Q1_3D_F0_Eta,
100  NF_D_H_Q1_3D_F1_Eta,
101  NF_D_H_Q1_3D_F2_Eta,
102  NF_D_H_Q1_3D_F3_Eta,
103  NF_D_H_Q1_3D_F4_Eta,
104  NF_D_H_Q1_3D_F5_Eta };
105 
106 static double *NF_D_H_Q1_3D_ZetaArray[6] = {
107  NF_D_H_Q1_3D_F0_Zeta,
108  NF_D_H_Q1_3D_F1_Zeta,
109  NF_D_H_Q1_3D_F2_Zeta,
110  NF_D_H_Q1_3D_F3_Zeta,
111  NF_D_H_Q1_3D_F4_Zeta,
112  NF_D_H_Q1_3D_F5_Zeta };
113 
114 static double *NF_D_H_Q1_3D_T = NULL;
115 static double *NF_D_H_Q1_3D_S = NULL;
116 
117 void NF_D_H_Q1_3D_EvalAll(TCollection *Coll, TBaseCell *Cell,
118  double *PointValues, double *Functionals)
119 {
120  Functionals[0] =
121  ( +PointValues[0]*NF_D_H_P2_3D_Weight[0]
122  +PointValues[1]*NF_D_H_P2_3D_Weight[1]
123  +PointValues[2]*NF_D_H_P2_3D_Weight[2]
124  +PointValues[3]*NF_D_H_P2_3D_Weight[3]
125  +PointValues[4]*NF_D_H_P2_3D_Weight[4]
126  +PointValues[5]*NF_D_H_P2_3D_Weight[5]
127  +PointValues[6]*NF_D_H_P2_3D_Weight[6]
128  +PointValues[7]*NF_D_H_P2_3D_Weight[7]
129  +PointValues[8]*NF_D_H_P2_3D_Weight[8]
130  +PointValues[9]*NF_D_H_P2_3D_Weight[9]
131  +PointValues[10]*NF_D_H_P2_3D_Weight[10]
132  +PointValues[11]*NF_D_H_P2_3D_Weight[11]
133  +PointValues[12]*NF_D_H_P2_3D_Weight[12]
134  +PointValues[13]*NF_D_H_P2_3D_Weight[13]
135  +PointValues[14]*NF_D_H_P2_3D_Weight[14]
136  +PointValues[15]*NF_D_H_P2_3D_Weight[15]
137  +PointValues[16]*NF_D_H_P2_3D_Weight[16]
138  +PointValues[17]*NF_D_H_P2_3D_Weight[17]
139  +PointValues[18]*NF_D_H_P2_3D_Weight[18]
140  +PointValues[19]*NF_D_H_P2_3D_Weight[19]
141  +PointValues[20]*NF_D_H_P2_3D_Weight[20]
142  +PointValues[21]*NF_D_H_P2_3D_Weight[21]
143  +PointValues[22]*NF_D_H_P2_3D_Weight[22]
144  +PointValues[23]*NF_D_H_P2_3D_Weight[23]
145  +PointValues[24]*NF_D_H_P2_3D_Weight[24]
146  +PointValues[25]*NF_D_H_P2_3D_Weight[25]
147  +PointValues[26]*NF_D_H_P2_3D_Weight[26] ) * 0.125;
148  Functionals[1] =
149  ( +PointValues[0]*NF_D_H_Q1_3D_Xi[0]*NF_D_H_P2_3D_Weight[0]
150  +PointValues[1]*NF_D_H_Q1_3D_Xi[1]*NF_D_H_P2_3D_Weight[1]
151  +PointValues[2]*NF_D_H_Q1_3D_Xi[2]*NF_D_H_P2_3D_Weight[2]
152  +PointValues[3]*NF_D_H_Q1_3D_Xi[3]*NF_D_H_P2_3D_Weight[3]
153  +PointValues[4]*NF_D_H_Q1_3D_Xi[4]*NF_D_H_P2_3D_Weight[4]
154  +PointValues[5]*NF_D_H_Q1_3D_Xi[5]*NF_D_H_P2_3D_Weight[5]
155  +PointValues[6]*NF_D_H_Q1_3D_Xi[6]*NF_D_H_P2_3D_Weight[6]
156  +PointValues[7]*NF_D_H_Q1_3D_Xi[7]*NF_D_H_P2_3D_Weight[7]
157  +PointValues[8]*NF_D_H_Q1_3D_Xi[8]*NF_D_H_P2_3D_Weight[8]
158  +PointValues[9]*NF_D_H_Q1_3D_Xi[9]*NF_D_H_P2_3D_Weight[9]
159  +PointValues[10]*NF_D_H_Q1_3D_Xi[10]*NF_D_H_P2_3D_Weight[10]
160  +PointValues[11]*NF_D_H_Q1_3D_Xi[11]*NF_D_H_P2_3D_Weight[11]
161  +PointValues[12]*NF_D_H_Q1_3D_Xi[12]*NF_D_H_P2_3D_Weight[12]
162  +PointValues[13]*NF_D_H_Q1_3D_Xi[13]*NF_D_H_P2_3D_Weight[13]
163  +PointValues[14]*NF_D_H_Q1_3D_Xi[14]*NF_D_H_P2_3D_Weight[14]
164  +PointValues[15]*NF_D_H_Q1_3D_Xi[15]*NF_D_H_P2_3D_Weight[15]
165  +PointValues[16]*NF_D_H_Q1_3D_Xi[16]*NF_D_H_P2_3D_Weight[16]
166  +PointValues[17]*NF_D_H_Q1_3D_Xi[17]*NF_D_H_P2_3D_Weight[17]
167  +PointValues[18]*NF_D_H_Q1_3D_Xi[18]*NF_D_H_P2_3D_Weight[18]
168  +PointValues[19]*NF_D_H_Q1_3D_Xi[19]*NF_D_H_P2_3D_Weight[19]
169  +PointValues[20]*NF_D_H_Q1_3D_Xi[20]*NF_D_H_P2_3D_Weight[20]
170  +PointValues[21]*NF_D_H_Q1_3D_Xi[21]*NF_D_H_P2_3D_Weight[21]
171  +PointValues[22]*NF_D_H_Q1_3D_Xi[22]*NF_D_H_P2_3D_Weight[22]
172  +PointValues[23]*NF_D_H_Q1_3D_Xi[23]*NF_D_H_P2_3D_Weight[23]
173  +PointValues[24]*NF_D_H_Q1_3D_Xi[24]*NF_D_H_P2_3D_Weight[24]
174  +PointValues[25]*NF_D_H_Q1_3D_Xi[25]*NF_D_H_P2_3D_Weight[25]
175  +PointValues[26]*NF_D_H_Q1_3D_Xi[26]*NF_D_H_P2_3D_Weight[26] ) * 0.375;
176  Functionals[2] =
177  ( +PointValues[0] *NF_D_H_Q1_3D_Eta[0] *NF_D_H_P2_3D_Weight[0]
178  +PointValues[1] *NF_D_H_Q1_3D_Eta[1] *NF_D_H_P2_3D_Weight[1]
179  +PointValues[2] *NF_D_H_Q1_3D_Eta[2] *NF_D_H_P2_3D_Weight[2]
180  +PointValues[3] *NF_D_H_Q1_3D_Eta[3] *NF_D_H_P2_3D_Weight[3]
181  +PointValues[4] *NF_D_H_Q1_3D_Eta[4] *NF_D_H_P2_3D_Weight[4]
182  +PointValues[5] *NF_D_H_Q1_3D_Eta[5] *NF_D_H_P2_3D_Weight[5]
183  +PointValues[6] *NF_D_H_Q1_3D_Eta[6] *NF_D_H_P2_3D_Weight[6]
184  +PointValues[7] *NF_D_H_Q1_3D_Eta[7] *NF_D_H_P2_3D_Weight[7]
185  +PointValues[8] *NF_D_H_Q1_3D_Eta[8] *NF_D_H_P2_3D_Weight[8]
186  +PointValues[9] *NF_D_H_Q1_3D_Eta[9] *NF_D_H_P2_3D_Weight[9]
187  +PointValues[10]*NF_D_H_Q1_3D_Eta[10]*NF_D_H_P2_3D_Weight[10]
188  +PointValues[11]*NF_D_H_Q1_3D_Eta[11]*NF_D_H_P2_3D_Weight[11]
189  +PointValues[12]*NF_D_H_Q1_3D_Eta[12]*NF_D_H_P2_3D_Weight[12]
190  +PointValues[13]*NF_D_H_Q1_3D_Eta[13]*NF_D_H_P2_3D_Weight[13]
191  +PointValues[14]*NF_D_H_Q1_3D_Eta[14]*NF_D_H_P2_3D_Weight[14]
192  +PointValues[15]*NF_D_H_Q1_3D_Eta[15]*NF_D_H_P2_3D_Weight[15]
193  +PointValues[16]*NF_D_H_Q1_3D_Eta[16]*NF_D_H_P2_3D_Weight[16]
194  +PointValues[17]*NF_D_H_Q1_3D_Eta[17]*NF_D_H_P2_3D_Weight[17]
195  +PointValues[18]*NF_D_H_Q1_3D_Eta[18]*NF_D_H_P2_3D_Weight[18]
196  +PointValues[19]*NF_D_H_Q1_3D_Eta[19]*NF_D_H_P2_3D_Weight[19]
197  +PointValues[20]*NF_D_H_Q1_3D_Eta[20]*NF_D_H_P2_3D_Weight[20]
198  +PointValues[21]*NF_D_H_Q1_3D_Eta[21]*NF_D_H_P2_3D_Weight[21]
199  +PointValues[22]*NF_D_H_Q1_3D_Eta[22]*NF_D_H_P2_3D_Weight[22]
200  +PointValues[23]*NF_D_H_Q1_3D_Eta[23]*NF_D_H_P2_3D_Weight[23]
201  +PointValues[24]*NF_D_H_Q1_3D_Eta[24]*NF_D_H_P2_3D_Weight[24]
202  +PointValues[25]*NF_D_H_Q1_3D_Eta[25]*NF_D_H_P2_3D_Weight[25]
203  +PointValues[26]*NF_D_H_Q1_3D_Eta[26]*NF_D_H_P2_3D_Weight[26] ) * 0.375;
204  Functionals[3] =
205  ( +PointValues[0]*NF_D_H_Q1_3D_Zeta[0]*NF_D_H_P2_3D_Weight[0]
206  +PointValues[1]*NF_D_H_Q1_3D_Zeta[1]*NF_D_H_P2_3D_Weight[1]
207  +PointValues[2]*NF_D_H_Q1_3D_Zeta[2]*NF_D_H_P2_3D_Weight[2]
208  +PointValues[3]*NF_D_H_Q1_3D_Zeta[3]*NF_D_H_P2_3D_Weight[3]
209  +PointValues[4]*NF_D_H_Q1_3D_Zeta[4]*NF_D_H_P2_3D_Weight[4]
210  +PointValues[5]*NF_D_H_Q1_3D_Zeta[5]*NF_D_H_P2_3D_Weight[5]
211  +PointValues[6]*NF_D_H_Q1_3D_Zeta[6]*NF_D_H_P2_3D_Weight[6]
212  +PointValues[7]*NF_D_H_Q1_3D_Zeta[7]*NF_D_H_P2_3D_Weight[7]
213  +PointValues[8]*NF_D_H_Q1_3D_Zeta[8]*NF_D_H_P2_3D_Weight[8]
214  +PointValues[9]*NF_D_H_Q1_3D_Zeta[9]*NF_D_H_P2_3D_Weight[9]
215  +PointValues[10]*NF_D_H_Q1_3D_Zeta[10]*NF_D_H_P2_3D_Weight[10]
216  +PointValues[11]*NF_D_H_Q1_3D_Zeta[11]*NF_D_H_P2_3D_Weight[11]
217  +PointValues[12]*NF_D_H_Q1_3D_Zeta[12]*NF_D_H_P2_3D_Weight[12]
218  +PointValues[13]*NF_D_H_Q1_3D_Zeta[13]*NF_D_H_P2_3D_Weight[13]
219  +PointValues[14]*NF_D_H_Q1_3D_Zeta[14]*NF_D_H_P2_3D_Weight[14]
220  +PointValues[15]*NF_D_H_Q1_3D_Zeta[15]*NF_D_H_P2_3D_Weight[15]
221  +PointValues[16]*NF_D_H_Q1_3D_Zeta[16]*NF_D_H_P2_3D_Weight[16]
222  +PointValues[17]*NF_D_H_Q1_3D_Zeta[17]*NF_D_H_P2_3D_Weight[17]
223  +PointValues[18]*NF_D_H_Q1_3D_Zeta[18]*NF_D_H_P2_3D_Weight[18]
224  +PointValues[19]*NF_D_H_Q1_3D_Zeta[19]*NF_D_H_P2_3D_Weight[19]
225  +PointValues[20]*NF_D_H_Q1_3D_Zeta[20]*NF_D_H_P2_3D_Weight[20]
226  +PointValues[21]*NF_D_H_Q1_3D_Zeta[21]*NF_D_H_P2_3D_Weight[21]
227  +PointValues[22]*NF_D_H_Q1_3D_Zeta[22]*NF_D_H_P2_3D_Weight[22]
228  +PointValues[23]*NF_D_H_Q1_3D_Zeta[23]*NF_D_H_P2_3D_Weight[23]
229  +PointValues[24]*NF_D_H_Q1_3D_Zeta[24]*NF_D_H_P2_3D_Weight[24]
230  +PointValues[25]*NF_D_H_Q1_3D_Zeta[25]*NF_D_H_P2_3D_Weight[25]
231  +PointValues[26]*NF_D_H_Q1_3D_Zeta[26]*NF_D_H_P2_3D_Weight[26] ) * 0.375;
232  Functionals[4] =
233  ( +PointValues[0] *NF_D_H_Q1_3D_Xi[0] *NF_D_H_Q1_3D_Eta[0] *NF_D_H_P2_3D_Weight[0]
234  +PointValues[1] *NF_D_H_Q1_3D_Xi[1] *NF_D_H_Q1_3D_Eta[1] *NF_D_H_P2_3D_Weight[1]
235  +PointValues[2] *NF_D_H_Q1_3D_Xi[2] *NF_D_H_Q1_3D_Eta[2] *NF_D_H_P2_3D_Weight[2]
236  +PointValues[3] *NF_D_H_Q1_3D_Xi[3] *NF_D_H_Q1_3D_Eta[3] *NF_D_H_P2_3D_Weight[3]
237  +PointValues[4] *NF_D_H_Q1_3D_Xi[4] *NF_D_H_Q1_3D_Eta[4] *NF_D_H_P2_3D_Weight[4]
238  +PointValues[5] *NF_D_H_Q1_3D_Xi[5] *NF_D_H_Q1_3D_Eta[5] *NF_D_H_P2_3D_Weight[5]
239  +PointValues[6] *NF_D_H_Q1_3D_Xi[6] *NF_D_H_Q1_3D_Eta[6] *NF_D_H_P2_3D_Weight[6]
240  +PointValues[7] *NF_D_H_Q1_3D_Xi[7] *NF_D_H_Q1_3D_Eta[7] *NF_D_H_P2_3D_Weight[7]
241  +PointValues[8] *NF_D_H_Q1_3D_Xi[8] *NF_D_H_Q1_3D_Eta[8] *NF_D_H_P2_3D_Weight[8]
242  +PointValues[9] *NF_D_H_Q1_3D_Xi[9] *NF_D_H_Q1_3D_Eta[9] *NF_D_H_P2_3D_Weight[9]
243  +PointValues[10]*NF_D_H_Q1_3D_Xi[10]*NF_D_H_Q1_3D_Eta[10]*NF_D_H_P2_3D_Weight[10]
244  +PointValues[11]*NF_D_H_Q1_3D_Xi[11]*NF_D_H_Q1_3D_Eta[11]*NF_D_H_P2_3D_Weight[11]
245  +PointValues[12]*NF_D_H_Q1_3D_Xi[12]*NF_D_H_Q1_3D_Eta[12]*NF_D_H_P2_3D_Weight[12]
246  +PointValues[13]*NF_D_H_Q1_3D_Xi[13]*NF_D_H_Q1_3D_Eta[13]*NF_D_H_P2_3D_Weight[13]
247  +PointValues[14]*NF_D_H_Q1_3D_Xi[14]*NF_D_H_Q1_3D_Eta[14]*NF_D_H_P2_3D_Weight[14]
248  +PointValues[15]*NF_D_H_Q1_3D_Xi[15]*NF_D_H_Q1_3D_Eta[15]*NF_D_H_P2_3D_Weight[15]
249  +PointValues[16]*NF_D_H_Q1_3D_Xi[16]*NF_D_H_Q1_3D_Eta[16]*NF_D_H_P2_3D_Weight[16]
250  +PointValues[17]*NF_D_H_Q1_3D_Xi[17]*NF_D_H_Q1_3D_Eta[17]*NF_D_H_P2_3D_Weight[17]
251  +PointValues[18]*NF_D_H_Q1_3D_Xi[18]*NF_D_H_Q1_3D_Eta[18]*NF_D_H_P2_3D_Weight[18]
252  +PointValues[19]*NF_D_H_Q1_3D_Xi[19]*NF_D_H_Q1_3D_Eta[19]*NF_D_H_P2_3D_Weight[19]
253  +PointValues[20]*NF_D_H_Q1_3D_Xi[20]*NF_D_H_Q1_3D_Eta[20]*NF_D_H_P2_3D_Weight[20]
254  +PointValues[21]*NF_D_H_Q1_3D_Xi[21]*NF_D_H_Q1_3D_Eta[21]*NF_D_H_P2_3D_Weight[21]
255  +PointValues[22]*NF_D_H_Q1_3D_Xi[22]*NF_D_H_Q1_3D_Eta[22]*NF_D_H_P2_3D_Weight[22]
256  +PointValues[23]*NF_D_H_Q1_3D_Xi[23]*NF_D_H_Q1_3D_Eta[23]*NF_D_H_P2_3D_Weight[23]
257  +PointValues[24]*NF_D_H_Q1_3D_Xi[24]*NF_D_H_Q1_3D_Eta[24]*NF_D_H_P2_3D_Weight[24]
258  +PointValues[25]*NF_D_H_Q1_3D_Xi[25]*NF_D_H_Q1_3D_Eta[25]*NF_D_H_P2_3D_Weight[25]
259  +PointValues[26]*NF_D_H_Q1_3D_Xi[26]*NF_D_H_Q1_3D_Eta[26]*NF_D_H_P2_3D_Weight[26] ) * 1.125;
260  Functionals[5] =
261  ( +PointValues[0] *NF_D_H_Q1_3D_Xi[0] *NF_D_H_Q1_3D_Zeta[0] *NF_D_H_P2_3D_Weight[0]
262  +PointValues[1] *NF_D_H_Q1_3D_Xi[1] *NF_D_H_Q1_3D_Zeta[1] *NF_D_H_P2_3D_Weight[1]
263  +PointValues[2] *NF_D_H_Q1_3D_Xi[2] *NF_D_H_Q1_3D_Zeta[2] *NF_D_H_P2_3D_Weight[2]
264  +PointValues[3] *NF_D_H_Q1_3D_Xi[3] *NF_D_H_Q1_3D_Zeta[3] *NF_D_H_P2_3D_Weight[3]
265  +PointValues[4] *NF_D_H_Q1_3D_Xi[4] *NF_D_H_Q1_3D_Zeta[4] *NF_D_H_P2_3D_Weight[4]
266  +PointValues[5] *NF_D_H_Q1_3D_Xi[5] *NF_D_H_Q1_3D_Zeta[5] *NF_D_H_P2_3D_Weight[5]
267  +PointValues[6] *NF_D_H_Q1_3D_Xi[6] *NF_D_H_Q1_3D_Zeta[6] *NF_D_H_P2_3D_Weight[6]
268  +PointValues[7] *NF_D_H_Q1_3D_Xi[7] *NF_D_H_Q1_3D_Zeta[7] *NF_D_H_P2_3D_Weight[7]
269  +PointValues[8] *NF_D_H_Q1_3D_Xi[8] *NF_D_H_Q1_3D_Zeta[8] *NF_D_H_P2_3D_Weight[8]
270  +PointValues[9] *NF_D_H_Q1_3D_Xi[9] *NF_D_H_Q1_3D_Zeta[9] *NF_D_H_P2_3D_Weight[9]
271  +PointValues[10]*NF_D_H_Q1_3D_Xi[10]*NF_D_H_Q1_3D_Zeta[10]*NF_D_H_P2_3D_Weight[10]
272  +PointValues[11]*NF_D_H_Q1_3D_Xi[11]*NF_D_H_Q1_3D_Zeta[11]*NF_D_H_P2_3D_Weight[11]
273  +PointValues[12]*NF_D_H_Q1_3D_Xi[12]*NF_D_H_Q1_3D_Zeta[12]*NF_D_H_P2_3D_Weight[12]
274  +PointValues[13]*NF_D_H_Q1_3D_Xi[13]*NF_D_H_Q1_3D_Zeta[13]*NF_D_H_P2_3D_Weight[13]
275  +PointValues[14]*NF_D_H_Q1_3D_Xi[14]*NF_D_H_Q1_3D_Zeta[14]*NF_D_H_P2_3D_Weight[14]
276  +PointValues[15]*NF_D_H_Q1_3D_Xi[15]*NF_D_H_Q1_3D_Zeta[15]*NF_D_H_P2_3D_Weight[15]
277  +PointValues[16]*NF_D_H_Q1_3D_Xi[16]*NF_D_H_Q1_3D_Zeta[16]*NF_D_H_P2_3D_Weight[16]
278  +PointValues[17]*NF_D_H_Q1_3D_Xi[17]*NF_D_H_Q1_3D_Zeta[17]*NF_D_H_P2_3D_Weight[17]
279  +PointValues[18]*NF_D_H_Q1_3D_Xi[18]*NF_D_H_Q1_3D_Zeta[18]*NF_D_H_P2_3D_Weight[18]
280  +PointValues[19]*NF_D_H_Q1_3D_Xi[19]*NF_D_H_Q1_3D_Zeta[19]*NF_D_H_P2_3D_Weight[19]
281  +PointValues[20]*NF_D_H_Q1_3D_Xi[20]*NF_D_H_Q1_3D_Zeta[20]*NF_D_H_P2_3D_Weight[20]
282  +PointValues[21]*NF_D_H_Q1_3D_Xi[21]*NF_D_H_Q1_3D_Zeta[21]*NF_D_H_P2_3D_Weight[21]
283  +PointValues[22]*NF_D_H_Q1_3D_Xi[22]*NF_D_H_Q1_3D_Zeta[22]*NF_D_H_P2_3D_Weight[22]
284  +PointValues[23]*NF_D_H_Q1_3D_Xi[23]*NF_D_H_Q1_3D_Zeta[23]*NF_D_H_P2_3D_Weight[23]
285  +PointValues[24]*NF_D_H_Q1_3D_Xi[24]*NF_D_H_Q1_3D_Zeta[24]*NF_D_H_P2_3D_Weight[24]
286  +PointValues[25]*NF_D_H_Q1_3D_Xi[25]*NF_D_H_Q1_3D_Zeta[25]*NF_D_H_P2_3D_Weight[25]
287  +PointValues[26]*NF_D_H_Q1_3D_Xi[26]*NF_D_H_Q1_3D_Zeta[26]*NF_D_H_P2_3D_Weight[26] ) * 1.125;
288  Functionals[6] =
289  ( +PointValues[0] *NF_D_H_Q1_3D_Eta[0] *NF_D_H_Q1_3D_Zeta[0] *NF_D_H_P2_3D_Weight[0]
290  +PointValues[1] *NF_D_H_Q1_3D_Eta[1] *NF_D_H_Q1_3D_Zeta[1] *NF_D_H_P2_3D_Weight[1]
291  +PointValues[2] *NF_D_H_Q1_3D_Eta[2] *NF_D_H_Q1_3D_Zeta[2] *NF_D_H_P2_3D_Weight[2]
292  +PointValues[3] *NF_D_H_Q1_3D_Eta[3] *NF_D_H_Q1_3D_Zeta[3] *NF_D_H_P2_3D_Weight[3]
293  +PointValues[4] *NF_D_H_Q1_3D_Eta[4] *NF_D_H_Q1_3D_Zeta[4] *NF_D_H_P2_3D_Weight[4]
294  +PointValues[5] *NF_D_H_Q1_3D_Eta[5] *NF_D_H_Q1_3D_Zeta[5] *NF_D_H_P2_3D_Weight[5]
295  +PointValues[6] *NF_D_H_Q1_3D_Eta[6] *NF_D_H_Q1_3D_Zeta[6] *NF_D_H_P2_3D_Weight[6]
296  +PointValues[7] *NF_D_H_Q1_3D_Eta[7] *NF_D_H_Q1_3D_Zeta[7] *NF_D_H_P2_3D_Weight[7]
297  +PointValues[8] *NF_D_H_Q1_3D_Eta[8] *NF_D_H_Q1_3D_Zeta[8] *NF_D_H_P2_3D_Weight[8]
298  +PointValues[9] *NF_D_H_Q1_3D_Eta[9] *NF_D_H_Q1_3D_Zeta[9] *NF_D_H_P2_3D_Weight[9]
299  +PointValues[10]*NF_D_H_Q1_3D_Eta[10]*NF_D_H_Q1_3D_Zeta[10]*NF_D_H_P2_3D_Weight[10]
300  +PointValues[11]*NF_D_H_Q1_3D_Eta[11]*NF_D_H_Q1_3D_Zeta[11]*NF_D_H_P2_3D_Weight[11]
301  +PointValues[12]*NF_D_H_Q1_3D_Eta[12]*NF_D_H_Q1_3D_Zeta[12]*NF_D_H_P2_3D_Weight[12]
302  +PointValues[13]*NF_D_H_Q1_3D_Eta[13]*NF_D_H_Q1_3D_Zeta[13]*NF_D_H_P2_3D_Weight[13]
303  +PointValues[14]*NF_D_H_Q1_3D_Eta[14]*NF_D_H_Q1_3D_Zeta[14]*NF_D_H_P2_3D_Weight[14]
304  +PointValues[15]*NF_D_H_Q1_3D_Eta[15]*NF_D_H_Q1_3D_Zeta[15]*NF_D_H_P2_3D_Weight[15]
305  +PointValues[16]*NF_D_H_Q1_3D_Eta[16]*NF_D_H_Q1_3D_Zeta[16]*NF_D_H_P2_3D_Weight[16]
306  +PointValues[17]*NF_D_H_Q1_3D_Eta[17]*NF_D_H_Q1_3D_Zeta[17]*NF_D_H_P2_3D_Weight[17]
307  +PointValues[18]*NF_D_H_Q1_3D_Eta[18]*NF_D_H_Q1_3D_Zeta[18]*NF_D_H_P2_3D_Weight[18]
308  +PointValues[19]*NF_D_H_Q1_3D_Eta[19]*NF_D_H_Q1_3D_Zeta[19]*NF_D_H_P2_3D_Weight[19]
309  +PointValues[20]*NF_D_H_Q1_3D_Eta[20]*NF_D_H_Q1_3D_Zeta[20]*NF_D_H_P2_3D_Weight[20]
310  +PointValues[21]*NF_D_H_Q1_3D_Eta[21]*NF_D_H_Q1_3D_Zeta[21]*NF_D_H_P2_3D_Weight[21]
311  +PointValues[22]*NF_D_H_Q1_3D_Eta[22]*NF_D_H_Q1_3D_Zeta[22]*NF_D_H_P2_3D_Weight[22]
312  +PointValues[23]*NF_D_H_Q1_3D_Eta[23]*NF_D_H_Q1_3D_Zeta[23]*NF_D_H_P2_3D_Weight[23]
313  +PointValues[24]*NF_D_H_Q1_3D_Eta[24]*NF_D_H_Q1_3D_Zeta[24]*NF_D_H_P2_3D_Weight[24]
314  +PointValues[25]*NF_D_H_Q1_3D_Eta[25]*NF_D_H_Q1_3D_Zeta[25]*NF_D_H_P2_3D_Weight[25]
315  +PointValues[26]*NF_D_H_Q1_3D_Eta[26]*NF_D_H_Q1_3D_Zeta[26]*NF_D_H_P2_3D_Weight[26] ) * 1.125;
316  Functionals[7] =
317  ( +PointValues[0] *NF_D_H_Q1_3D_Xi[0] *NF_D_H_Q1_3D_Eta[0] *NF_D_H_Q1_3D_Zeta[0] *NF_D_H_P2_3D_Weight[0]
318  +PointValues[1] *NF_D_H_Q1_3D_Xi[1] *NF_D_H_Q1_3D_Eta[1] *NF_D_H_Q1_3D_Zeta[1] *NF_D_H_P2_3D_Weight[1]
319  +PointValues[2] *NF_D_H_Q1_3D_Xi[2] *NF_D_H_Q1_3D_Eta[2] *NF_D_H_Q1_3D_Zeta[2] *NF_D_H_P2_3D_Weight[2]
320  +PointValues[3] *NF_D_H_Q1_3D_Xi[3] *NF_D_H_Q1_3D_Eta[3] *NF_D_H_Q1_3D_Zeta[3] *NF_D_H_P2_3D_Weight[3]
321  +PointValues[4] *NF_D_H_Q1_3D_Xi[4] *NF_D_H_Q1_3D_Eta[4] *NF_D_H_Q1_3D_Zeta[4] *NF_D_H_P2_3D_Weight[4]
322  +PointValues[5] *NF_D_H_Q1_3D_Xi[5] *NF_D_H_Q1_3D_Eta[5] *NF_D_H_Q1_3D_Zeta[5] *NF_D_H_P2_3D_Weight[5]
323  +PointValues[6] *NF_D_H_Q1_3D_Xi[6] *NF_D_H_Q1_3D_Eta[6] *NF_D_H_Q1_3D_Zeta[6] *NF_D_H_P2_3D_Weight[6]
324  +PointValues[7] *NF_D_H_Q1_3D_Xi[7] *NF_D_H_Q1_3D_Eta[7] *NF_D_H_Q1_3D_Zeta[7] *NF_D_H_P2_3D_Weight[7]
325  +PointValues[8] *NF_D_H_Q1_3D_Xi[8] *NF_D_H_Q1_3D_Eta[8] *NF_D_H_Q1_3D_Zeta[8] *NF_D_H_P2_3D_Weight[8]
326  +PointValues[9] *NF_D_H_Q1_3D_Xi[9] *NF_D_H_Q1_3D_Eta[9] *NF_D_H_Q1_3D_Zeta[9] *NF_D_H_P2_3D_Weight[9]
327  +PointValues[10]*NF_D_H_Q1_3D_Xi[10]*NF_D_H_Q1_3D_Eta[10]*NF_D_H_Q1_3D_Zeta[10]*NF_D_H_P2_3D_Weight[10]
328  +PointValues[11]*NF_D_H_Q1_3D_Xi[11]*NF_D_H_Q1_3D_Eta[11]*NF_D_H_Q1_3D_Zeta[11]*NF_D_H_P2_3D_Weight[11]
329  +PointValues[12]*NF_D_H_Q1_3D_Xi[12]*NF_D_H_Q1_3D_Eta[12]*NF_D_H_Q1_3D_Zeta[12]*NF_D_H_P2_3D_Weight[12]
330  +PointValues[13]*NF_D_H_Q1_3D_Xi[13]*NF_D_H_Q1_3D_Eta[13]*NF_D_H_Q1_3D_Zeta[13]*NF_D_H_P2_3D_Weight[13]
331  +PointValues[14]*NF_D_H_Q1_3D_Xi[14]*NF_D_H_Q1_3D_Eta[14]*NF_D_H_Q1_3D_Zeta[14]*NF_D_H_P2_3D_Weight[14]
332  +PointValues[15]*NF_D_H_Q1_3D_Xi[15]*NF_D_H_Q1_3D_Eta[15]*NF_D_H_Q1_3D_Zeta[15]*NF_D_H_P2_3D_Weight[15]
333  +PointValues[16]*NF_D_H_Q1_3D_Xi[16]*NF_D_H_Q1_3D_Eta[16]*NF_D_H_Q1_3D_Zeta[16]*NF_D_H_P2_3D_Weight[16]
334  +PointValues[17]*NF_D_H_Q1_3D_Xi[17]*NF_D_H_Q1_3D_Eta[17]*NF_D_H_Q1_3D_Zeta[17]*NF_D_H_P2_3D_Weight[17]
335  +PointValues[18]*NF_D_H_Q1_3D_Xi[18]*NF_D_H_Q1_3D_Eta[18]*NF_D_H_Q1_3D_Zeta[18]*NF_D_H_P2_3D_Weight[18]
336  +PointValues[19]*NF_D_H_Q1_3D_Xi[19]*NF_D_H_Q1_3D_Eta[19]*NF_D_H_Q1_3D_Zeta[19]*NF_D_H_P2_3D_Weight[19]
337  +PointValues[20]*NF_D_H_Q1_3D_Xi[20]*NF_D_H_Q1_3D_Eta[20]*NF_D_H_Q1_3D_Zeta[20]*NF_D_H_P2_3D_Weight[20]
338  +PointValues[21]*NF_D_H_Q1_3D_Xi[21]*NF_D_H_Q1_3D_Eta[21]*NF_D_H_Q1_3D_Zeta[21]*NF_D_H_P2_3D_Weight[21]
339  +PointValues[22]*NF_D_H_Q1_3D_Xi[22]*NF_D_H_Q1_3D_Eta[22]*NF_D_H_Q1_3D_Zeta[22]*NF_D_H_P2_3D_Weight[22]
340  +PointValues[23]*NF_D_H_Q1_3D_Xi[23]*NF_D_H_Q1_3D_Eta[23]*NF_D_H_Q1_3D_Zeta[23]*NF_D_H_P2_3D_Weight[23]
341  +PointValues[24]*NF_D_H_Q1_3D_Xi[24]*NF_D_H_Q1_3D_Eta[24]*NF_D_H_Q1_3D_Zeta[24]*NF_D_H_P2_3D_Weight[24]
342  +PointValues[25]*NF_D_H_Q1_3D_Xi[25]*NF_D_H_Q1_3D_Eta[25]*NF_D_H_Q1_3D_Zeta[25]*NF_D_H_P2_3D_Weight[25]
343  +PointValues[26]*NF_D_H_Q1_3D_Xi[26]*NF_D_H_Q1_3D_Eta[26]*NF_D_H_Q1_3D_Zeta[26]*NF_D_H_P2_3D_Weight[26] ) * 3.375;
344 }
345 
346 void NF_D_H_Q1_3D_EvalFace(TCollection *Coll, TBaseCell *Cell, int Joint,
347  double *PointValues, double *Functionals)
348 {
349  ErrMsg("EvalFace not implemented for discontinuous Q1-elements on hexahedra");
350  exit(0);
351 }
352 
353 static int NF_D_H_Q1_3D_N_AllFunctionals = 8;
354 static int NF_D_H_Q1_3D_N_PointsAll = 27;
355 static int NF_D_H_Q1_3D_N_FaceFunctionals[] = { 0, 0, 0, 0, 0, 0 };
356 static int NF_D_H_Q1_3D_N_PointsFace[] = { 0, 0, 0, 0, 0, 0 };
357 
358 TNodalFunctional3D *NF_D_H_Q1_3D_Obj = new TNodalFunctional3D
359  (NF_D_H_Q1_3D, NF_D_H_Q1_3D_N_AllFunctionals,
360  NF_D_H_Q1_3D_N_FaceFunctionals, NF_D_H_Q1_3D_N_PointsAll,
361  NF_D_H_Q1_3D_N_PointsFace,
362  NF_D_H_Q1_3D_Xi, NF_D_H_Q1_3D_Eta, NF_D_H_Q1_3D_Zeta,
363  NF_D_H_Q1_3D_XiArray, NF_D_H_Q1_3D_EtaArray,
364  NF_D_H_Q1_3D_ZetaArray,
365  NF_D_H_Q1_3D_T, NF_D_H_Q1_3D_S,
366  NF_D_H_Q1_3D_EvalAll, NF_D_H_Q1_3D_EvalFace);
TNodalFunctional3D
Definition: NodalFunctional3D.h:21
TCollection
store cells in an array, used by cell iterators
Definition: Collection.h:18
TBaseCell
information for finite element data structure
Definition: BaseCell.h:25
Generated on Tue May 10 2016 09:08:00 for ParMooN by   doxygen 1.8.5