6 static void D_Q_P7_2D_Funct(
double xi,
double eta,
double *values)
12 double t13 = -1.0+3.0*t1;
14 double t17 = -1.0+3.0*t5;
17 double t30 = xi*(5.0*t1-3.0);
18 double t36 = 5.0*t5-3.0;
21 double t50 = 3.0+35.0*t24-30.0*t1;
23 double t59 = 3.0+35.0*t39-30.0*t5;
26 double t75 = xi*(63.0*t24-70.0*t1+15.0);
27 double t87 = 63.0*t39-70.0*t5+15.0;
33 values[3] = -1.0/2.0+3.0/2.0*t1;
35 values[5] = -1.0/2.0+3.0/2.0*t5;
36 values[6] = 5.0/2.0*t8-3.0/2.0*xi;
38 values[8] = xi*t17/2.0;
39 values[9] = 5.0/2.0*t20-3.0/2.0*eta;
40 values[10] = 3.0/8.0+35.0/8.0*t24-15.0/4.0*t1;
41 values[11] = t30*eta/2.0;
42 values[12] = t13*t17/4.0;
43 values[13] = t4*t36/2.0;
44 values[14] = 3.0/8.0+35.0/8.0*t39-15.0/4.0*t5;
45 values[15] = 63.0/8.0*t43-35.0/4.0*t8+15.0/8.0*xi;
47 values[17] = t30*t17/4.0;
48 values[18] = t14*t36/4.0;
49 values[19] = xi*t59/8.0;
50 values[20] = 63.0/8.0*t62-35.0/4.0*t20+15.0/8.0*eta;
51 values[21] = -5.0/16.0+231.0/16.0*t67-315.0/16.0*t24+105.0/16.0*t1;
52 values[22] = t75*eta/8.0;
53 values[23] = t50*t17/16.0;
54 values[24] = t30*eta*t36/4.0;
55 values[25] = t13*t59/16.0;
56 values[26] = t4*t87/8.0;
57 values[27] = -5.0/16.0+231.0/16.0*t90-315.0/16.0*t39+105.0/16.0*t5;
58 values[28] = 429.0/16.0*t24*t8-693.0/16.0*t43+315.0/16.0*t8-35.0/16.0*xi;
59 values[29] = (-5.0+231.0*t67-315.0*t24+105.0*t1)*eta/16.0;
60 values[30] = t75*t17/16.0;
61 values[31] = t51*t36/16.0;
62 values[32] = t30*t59/16.0;
63 values[33] = t14*t87/16.0;
64 values[34] = xi*(-5.0+231.0*t90-315.0*t39+105.0*t5)/16.0;
65 values[35] = 429.0/16.0*t39*t20-693.0/16.0*t62+315.0/16.0*t20-35.0/16.0*eta;
69 static void D_Q_P7_2D_DeriveXi(
double xi,
double eta,
double *values)
76 double t19 = -1.0+3.0*t7;
77 double t23 = 5.0*t7-3.0;
80 double t32 = xi*(7.0*t2-3.0);
84 double t72 = 63.0*t42-70.0*t7+15.0;
92 values[6] = 15.0/2.0*t2-3.0/2.0;
94 values[8] = -1.0/2.0+3.0/2.0*t7;
96 values[10] = 35.0/2.0*t10-15.0/2.0*xi;
97 values[11] = 15.0/2.0*t14-3.0/2.0*eta;
98 values[12] = 3.0/2.0*xi*t19;
101 values[15] = 315.0/8.0*t26-105.0/4.0*t2+15.0/8.0;
102 values[16] = 5.0/2.0*t32*eta;
103 values[17] = -15.0/4.0*t2+45.0/4.0*t36+3.0/4.0-9.0/4.0*t7;
104 values[18] = 3.0/2.0*t5*t23;
105 values[19] = 3.0/8.0+35.0/8.0*t42-15.0/4.0*t7;
107 values[21] = 693.0/8.0*t26*xi-315.0/4.0*t10+105.0/8.0*xi;
108 values[22] = 315.0/8.0*eta*t26-105.0/4.0*t14+15.0/8.0*eta;
109 values[23] = 5.0/4.0*t32*t19;
110 values[24] = 75.0/4.0*t2*t58-45.0/4.0*t14-15.0/4.0*t58+9.0/4.0*eta;
111 values[25] = 3.0/8.0*xi*(3.0+35.0*t42-30.0*t7);
112 values[26] = eta*t72/8.0;
114 values[28] = 3003.0/16.0*t26*t2-3465.0/16.0*t26+945.0/16.0*t2-35.0/16.0;
115 values[29] = 21.0/8.0*xi*(33.0*t26-30.0*t2+5.0)*eta;
116 values[30] = -315.0/16.0*t26+945.0/16.0*t26*t7+105.0/8.0*t2-315.0/8.0*t36-15.0/16.0+45.0/16.0*t7;
117 values[31] = 5.0/4.0*t32*t24;
118 values[32] = 45.0/16.0*t2+525.0/16.0*t2*t42-225.0/8.0*t36-9.0/16.0-105.0/16.0*t42+45.0/8.0*t7;
119 values[33] = 3.0/8.0*t5*t72;
120 values[34] = -5.0/16.0+231.0/16.0*t42*t7-315.0/16.0*t42+105.0/16.0*t7;
125 static void D_Q_P7_2D_DeriveEta(
double xi,
double eta,
double *values)
130 double t12 = xi*(5.0*t2-3.0);
131 double t16 = (-1.0+3.0*t2)*eta;
136 double t38 = 7.0*t7-3.0;
138 double t48 = xi*(63.0*t26-70.0*t2+15.0);
148 values[7] = -1.0/2.0+3.0/2.0*t2;
150 values[9] = 15.0/2.0*t7-3.0/2.0;
152 values[11] = t12/2.0;
153 values[12] = 3.0/2.0*t16;
154 values[13] = 15.0/2.0*t18-3.0/2.0*xi;
155 values[14] = 35.0/2.0*t22-15.0/2.0*eta;
157 values[16] = 3.0/8.0+35.0/8.0*t26-15.0/4.0*t2;
158 values[17] = 3.0/2.0*t12*eta;
159 values[18] = -15.0/4.0*t7+3.0/4.0+45.0/4.0*t33-9.0/4.0*t2;
160 values[19] = 5.0/2.0*t5*t38;
161 values[20] = 315.0/8.0*t41-105.0/4.0*t7+15.0/8.0;
163 values[22] = t48/8.0;
164 values[23] = 3.0/8.0*(3.0+35.0*t26-30.0*t2)*eta;
165 values[24] = 75.0/4.0*t55*t7-15.0/4.0*t55-45.0/4.0*t18+9.0/4.0*xi;
166 values[25] = 5.0/4.0*t16*t38;
167 values[26] = 315.0/8.0*xi*t41-105.0/4.0*t18+15.0/8.0*xi;
168 values[27] = 693.0/8.0*t41*eta-315.0/4.0*t22+105.0/8.0*eta;
170 values[29] = -5.0/16.0+231.0/16.0*t26*t2-315.0/16.0*t26+105.0/16.0*t2;
171 values[30] = 3.0/8.0*t48*eta;
172 values[31] = 45.0/16.0*t7-9.0/16.0+525.0/16.0*t26*t7-105.0/16.0*t26-225.0/8.0*t33+45.0/8.0*t2;
173 values[32] = 5.0/4.0*t12*eta*t38;
174 values[33] = -315.0/16.0*t41+105.0/8.0*t7-15.0/16.0+945.0/16.0*t2*t41-315.0/8.0*t33+45.0/16.0*t2;
175 values[34] = 21.0/8.0*t5*(33.0*t41-30.0*t7+5.0);
176 values[35] = 3003.0/16.0*t41*t7-3465.0/16.0*t41+945.0/16.0*t7-35.0/16.0;
180 static void D_Q_P7_2D_DeriveXiXi(
double xi,
double eta,
190 double t38 = 45.0/4.0*t8;
204 values[10] = 105.0/2.0*t3-15.0/2.0;
205 values[11] = 15.0*t6;
206 values[12] = -3.0/2.0+9.0/2.0*t8;
209 values[15] = 315.0/2.0*t11-105.0/2.0*xi;
210 values[16] = 105.0/2.0*t15-15.0/2.0*eta;
211 values[17] = -15.0/2.0*xi+45.0/2.0*t20;
212 values[18] = 3.0/2.0*eta*(5.0*t8-3.0);
215 values[21] = 3465.0/8.0*t27-945.0/4.0*t3+105.0/8.0;
216 values[22] = 315.0/2.0*eta*t11-105.0/2.0*t6;
217 values[23] = -105.0/4.0*t3+315.0/4.0*t3*t8+15.0/4.0-t38;
218 values[24] = 75.0/2.0*xi*t40-45.0/2.0*t6;
219 values[25] = 9.0/8.0+105.0/8.0*t45-t38;
222 values[28] = 9009.0/8.0*t27*xi-3465.0/4.0*t11+945.0/8.0*xi;
223 values[29] = 3465.0/8.0*eta*t27-945.0/4.0*t15+105.0/8.0*eta;
224 values[30] = -315.0/4.0*t11+945.0/4.0*t11*t8+105.0/4.0*xi-315.0/4.0*t20;
225 values[31] = 525.0/4.0*t3*t40-315.0/4.0*t15-75.0/4.0*t40+45.0/4.0*eta;
226 values[32] = 45.0/8.0*xi+525.0/8.0*xi*t45-225.0/4.0*t20;
227 values[33] = 3.0/8.0*eta*(63.0*t45-70.0*t8+15.0);
233 static void D_Q_P7_2D_DeriveXiEta(
double xi,
double eta,
239 double t13 = xi*(7.0*t3-3.0);
258 values[11] = 15.0/2.0*t3-3.0/2.0;
260 values[13] = 15.0/2.0*t8-3.0/2.0;
263 values[16] = 5.0/2.0*t13;
264 values[17] = 45.0/2.0*t15-9.0/2.0*eta;
265 values[18] = 45.0/2.0*t19-9.0/2.0*xi;
266 values[19] = 35.0/2.0*t23-15.0/2.0*eta;
269 values[22] = 315.0/8.0*t27-105.0/4.0*t3+15.0/8.0;
270 values[23] = 15.0/2.0*t13*eta;
271 values[24] = 225.0/4.0*t3*t8-45.0/4.0*t3-45.0/4.0*t8+9.0/4.0;
272 values[25] = 15.0/2.0*t6*(7.0*t8-3.0);
273 values[26] = 315.0/8.0*t42-105.0/4.0*t8+15.0/8.0;
276 values[29] = 21.0/8.0*xi*(33.0*t27-30.0*t3+5.0);
277 values[30] = 945.0/8.0*eta*t27-315.0/4.0*t15+45.0/8.0*eta;
278 values[31] = 525.0/4.0*t56*t8-105.0/4.0*t56-225.0/4.0*t19+45.0/4.0*xi;
279 values[32] = 525.0/4.0*t3*t23-225.0/4.0*t15-105.0/4.0*t23+45.0/4.0*eta;
280 values[33] = 945.0/8.0*xi*t42-315.0/4.0*t19+45.0/8.0*xi;
281 values[34] = 693.0/8.0*t42*eta-315.0/4.0*t23+105.0/8.0*eta;
286 static void D_Q_P7_2D_DeriveEtaEta(
double xi,
double eta,
296 double t29 = 45.0/4.0*t3;
309 values[9] = 15.0*eta;
312 values[12] = 9.0/2.0*t3-3.0/2.0;
313 values[13] = 15.0*t6;
314 values[14] = 105.0/2.0*t8-15.0/2.0;
317 values[17] = 3.0/2.0*xi*(5.0*t3-3.0);
318 values[18] = -15.0/2.0*eta+45.0/2.0*t16;
319 values[19] = 105.0/2.0*t19-15.0/2.0*xi;
320 values[20] = 315.0/2.0*t23-105.0/2.0*eta;
323 values[23] = 9.0/8.0+105.0/8.0*t27-t29;
324 values[24] = 75.0/2.0*eta*t31-45.0/2.0*t6;
325 values[25] = -105.0/4.0*t8+15.0/4.0+315.0/4.0*t3*t8-t29;
326 values[26] = 315.0/2.0*xi*t23-105.0/2.0*t6;
327 values[27] = 3465.0/8.0*t44-945.0/4.0*t8+105.0/8.0;
330 values[30] = 3.0/8.0*xi*(63.0*t27-70.0*t3+15.0);
331 values[31] = 45.0/8.0*eta+525.0/8.0*eta*t27-225.0/4.0*t16;
332 values[32] = 525.0/4.0*t31*t8-75.0/4.0*t31-315.0/4.0*t19+45.0/4.0*xi;
333 values[33] = -315.0/4.0*t23+105.0/4.0*eta+945.0/4.0*t3*t23-315.0/4.0*t16;
334 values[34] = 3465.0/8.0*xi*t44-945.0/4.0*t19+105.0/8.0*xi;
335 values[35] = 9009.0/8.0*t44*eta-3465.0/4.0*t23+945.0/8.0*eta;
341 (36, BF_D_Q_P7_2D, BFUnitSquare,
342 D_Q_P7_2D_Funct, D_Q_P7_2D_DeriveXi,
343 D_Q_P7_2D_DeriveEta, D_Q_P7_2D_DeriveXiXi,
344 D_Q_P7_2D_DeriveXiEta, D_Q_P7_2D_DeriveEtaEta, 7, 7,
Definition: BaseFunct2D.h:27