6 static void D_Q_P6_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;
20 double t50 = 3.0+35.0*t24-30.0*t1;
21 double t59 = 3.0+35.0*t39-30.0*t5;
26 values[3] = -1.0/2.0+3.0/2.0*t1;
28 values[5] = -1.0/2.0+3.0/2.0*t5;
29 values[6] = 5.0/2.0*t8-3.0/2.0*xi;
31 values[8] = xi*t17/2.0;
32 values[9] = 5.0/2.0*t20-3.0/2.0*eta;
33 values[10] = 3.0/8.0+35.0/8.0*t24-15.0/4.0*t1;
34 values[11] = t30*eta/2.0;
35 values[12] = t13*t17/4.0;
36 values[13] = t4*t36/2.0;
37 values[14] = 3.0/8.0+35.0/8.0*t39-15.0/4.0*t5;
38 values[15] = 63.0/8.0*t24*xi-35.0/4.0*t8+15.0/8.0*xi;
39 values[16] = t50*eta/8.0;
40 values[17] = t30*t17/4.0;
41 values[18] = t14*t36/4.0;
42 values[19] = xi*t59/8.0;
43 values[20] = 63.0/8.0*t39*eta-35.0/4.0*t20+15.0/8.0*eta;
44 values[21] = -5.0/16.0+231.0/16.0*t24*t1-315.0/16.0*t24+105.0/16.0*t1;
45 values[22] = xi*(63.0*t24-70.0*t1+15.0)*eta/8.0;
46 values[23] = t50*t17/16.0;
47 values[24] = t30*eta*t36/4.0;
48 values[25] = t13*t59/16.0;
49 values[26] = t4*(63.0*t39-70.0*t5+15.0)/8.0;
50 values[27] = -5.0/16.0+231.0/16.0*t39*t5-315.0/16.0*t39+105.0/16.0*t5;
54 static void D_Q_P6_2D_DeriveXi(
double xi,
double eta,
double *values)
61 double t19 = -1.0+3.0*t7;
62 double t23 = 5.0*t7-3.0;
64 double t32 = xi*(7.0*t2-3.0);
74 values[6] = 15.0/2.0*t2-3.0/2.0;
76 values[8] = -1.0/2.0+3.0/2.0*t7;
78 values[10] = 35.0/2.0*t10-15.0/2.0*xi;
79 values[11] = 15.0/2.0*t14-3.0/2.0*eta;
80 values[12] = 3.0/2.0*xi*t19;
81 values[13] = eta*t23/2.0;
83 values[15] = 315.0/8.0*t26-105.0/4.0*t2+15.0/8.0;
84 values[16] = 5.0/2.0*t32*eta;
85 values[17] = -15.0/4.0*t2+45.0/4.0*t2*t7+3.0/4.0-9.0/4.0*t7;
86 values[18] = 3.0/2.0*t5*t23;
87 values[19] = 3.0/8.0+35.0/8.0*t42-15.0/4.0*t7;
89 values[21] = 693.0/8.0*t26*xi-315.0/4.0*t10+105.0/8.0*xi;
90 values[22] = 315.0/8.0*eta*t26-105.0/4.0*t14+15.0/8.0*eta;
91 values[23] = 5.0/4.0*t32*t19;
92 values[24] = 75.0/4.0*t2*t58-45.0/4.0*t14-15.0/4.0*t58+9.0/4.0*eta;
93 values[25] = 3.0/8.0*xi*(3.0+35.0*t42-30.0*t7);
94 values[26] = eta*(63.0*t42-70.0*t7+15.0)/8.0;
99 static void D_Q_P6_2D_DeriveEta(
double xi,
double eta,
double *values)
104 double t12 = xi*(5.0*t2-3.0);
105 double t16 = (-1.0+3.0*t2)*eta;
109 double t38 = 7.0*t7-3.0;
120 values[7] = -1.0/2.0+3.0/2.0*t2;
122 values[9] = 15.0/2.0*t7-3.0/2.0;
124 values[11] = t12/2.0;
125 values[12] = 3.0/2.0*t16;
126 values[13] = 15.0/2.0*t18-3.0/2.0*xi;
127 values[14] = 35.0/2.0*t22-15.0/2.0*eta;
129 values[16] = 3.0/8.0+35.0/8.0*t26-15.0/4.0*t2;
130 values[17] = 3.0/2.0*t12*eta;
131 values[18] = -15.0/4.0*t7+3.0/4.0+45.0/4.0*t2*t7-9.0/4.0*t2;
132 values[19] = 5.0/2.0*t5*t38;
133 values[20] = 315.0/8.0*t41-105.0/4.0*t7+15.0/8.0;
135 values[22] = xi*(63.0*t26-70.0*t2+15.0)/8.0;
136 values[23] = 3.0/8.0*(3.0+35.0*t26-30.0*t2)*eta;
137 values[24] = 75.0/4.0*t55*t7-15.0/4.0*t55-45.0/4.0*t18+9.0/4.0*xi;
138 values[25] = 5.0/4.0*t16*t38;
139 values[26] = 315.0/8.0*xi*t41-105.0/4.0*t18+15.0/8.0*xi;
140 values[27] = 693.0/8.0*t41*eta-315.0/4.0*t22+105.0/8.0*eta;
144 static void D_Q_P6_2D_DeriveXiXi(
double xi,
double eta,
152 double t38 = 45.0/4.0*t8;
165 values[10] = 105.0/2.0*t3-15.0/2.0;
166 values[11] = 15.0*t6;
167 values[12] = -3.0/2.0+9.0/2.0*t8;
170 values[15] = 315.0/2.0*t11-105.0/2.0*xi;
171 values[16] = 105.0/2.0*t3*eta-15.0/2.0*eta;
172 values[17] = -15.0/2.0*xi+45.0/2.0*xi*t8;
173 values[18] = 3.0/2.0*eta*(5.0*t8-3.0);
176 values[21] = 3465.0/8.0*t27-945.0/4.0*t3+105.0/8.0;
177 values[22] = 315.0/2.0*eta*t11-105.0/2.0*t6;
178 values[23] = -105.0/4.0*t3+315.0/4.0*t3*t8+15.0/4.0-t38;
179 values[24] = 75.0/2.0*xi*t8*eta-45.0/2.0*t6;
180 values[25] = 9.0/8.0+105.0/8.0*t45-t38;
186 static void D_Q_P6_2D_DeriveXiEta(
double xi,
double eta,
192 double t13 = xi*(7.0*t3-3.0);
207 values[11] = 15.0/2.0*t3-3.0/2.0;
209 values[13] = -3.0/2.0+15.0/2.0*t8;
212 values[16] = 5.0/2.0*t13;
213 values[17] = 45.0/2.0*t3*eta-9.0/2.0*eta;
214 values[18] = 45.0/2.0*xi*t8-9.0/2.0*xi;
215 values[19] = 35.0/2.0*t8*eta-15.0/2.0*eta;
218 values[22] = 315.0/8.0*t27-105.0/4.0*t3+15.0/8.0;
219 values[23] = 15.0/2.0*t13*eta;
220 values[24] = 225.0/4.0*t3*t8-45.0/4.0*t3-45.0/4.0*t8+9.0/4.0;
221 values[25] = 15.0/2.0*t6*(7.0*t8-3.0);
222 values[26] = 315.0/8.0*t42-105.0/4.0*t8+15.0/8.0;
227 static void D_Q_P6_2D_DeriveEtaEta(
double xi,
double eta,
235 double t29 = 45.0/4.0*t3;
247 values[9] = 15.0*eta;
250 values[12] = 9.0/2.0*t3-3.0/2.0;
251 values[13] = 15.0*t6;
252 values[14] = 105.0/2.0*t8-15.0/2.0;
255 values[17] = 3.0/2.0*xi*(5.0*t3-3.0);
256 values[18] = -15.0/2.0*eta+45.0/2.0*t3*eta;
257 values[19] = 105.0/2.0*xi*t8-15.0/2.0*xi;
258 values[20] = 315.0/2.0*t23-105.0/2.0*eta;
261 values[23] = 9.0/8.0+105.0/8.0*t27-t29;
262 values[24] = 75.0/2.0*eta*t3*xi-45.0/2.0*t6;
263 values[25] = -105.0/4.0*t8+15.0/4.0+315.0/4.0*t3*t8-t29;
264 values[26] = 315.0/2.0*xi*t23-105.0/2.0*t6;
265 values[27] = 3465.0/8.0*t44-945.0/4.0*t8+105.0/8.0;
271 (28, BF_D_Q_P6_2D, BFUnitSquare,
272 D_Q_P6_2D_Funct, D_Q_P6_2D_DeriveXi,
273 D_Q_P6_2D_DeriveEta, D_Q_P6_2D_DeriveXiXi,
274 D_Q_P6_2D_DeriveXiEta, D_Q_P6_2D_DeriveEtaEta, 6, 6,
Definition: BaseFunct2D.h:27