ParMooN
 All Classes Functions Variables Friends Pages
BF_C_Q_Q9_2D.h
1 
2 // ***********************************************************************
3 // Q9 element, conforming, 2D
4 // ***********************************************************************
5 
6 // base function values
7 static void C_Q_Q9_2D_Funct(double xi, double eta, double *values)
8 {
9 
10  double xi0= -0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi+0.216243896484375E1*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi-0.627374267578125*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi+0.5067836216517857E-1*xi*xi*xi-0.5067836216517857E-1*xi*xi-0.5340576171875E-3*xi+0.5340576171875E-3;
11  double xi1= 0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.1459646301269531E2*xi*xi*xi*xi*xi*xi*xi*xi-0.2687602713448661E2*xi*xi*xi*xi*xi*xi*xi+0.2090357666015625E2*xi*xi*xi*xi*xi*xi+0.8849981689453125E1*xi*xi*xi*xi*xi-0.6883319091796875E1*xi*xi*xi*xi-0.7487810407366071*xi*xi*xi+0.58238525390625*xi*xi+0.7945469447544643E-2*xi-0.61798095703125E-2;
12  double xi2= -0.7506752406529018E2*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.4170418003627232E2*xi*xi*xi*xi*xi*xi*xi*xi+0.129746337890625E3*xi*xi*xi*xi*xi*xi*xi-0.72081298828125E2*xi*xi*xi*xi*xi*xi-0.603881103515625E2*xi*xi*xi*xi*xi+0.335489501953125E2*xi*xi*xi*xi+0.5771589006696429E1*xi*xi*xi-0.3206438337053571E1*xi*xi-0.6229248046875E-1*xi+0.3460693359375E-1;
13  double xi3= 0.1751575561523438E3*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.5838585205078125E2*xi*xi*xi*xi*xi*xi*xi*xi-0.337340478515625E3*xi*xi*xi*xi*xi*xi*xi+0.112446826171875E3*xi*xi*xi*xi*xi*xi+0.1968620361328125E3*xi*xi*xi*xi*xi-0.656206787109375E2*xi*xi*xi*xi-0.35082861328125E2*xi*xi*xi+0.11694287109375E2*xi*xi+0.40374755859375*xi-0.13458251953125;
14  double xi4= -0.2627363342285156E3*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2919292602539063E2*xi*xi*xi*xi*xi*xi*xi*xi+0.5319599853515625E3*xi*xi*xi*xi*xi*xi*xi-0.591066650390625E2*xi*xi*xi*xi*xi*xi-0.3449490600585938E3*xi*xi*xi*xi*xi+0.3832767333984375E2*xi*xi*xi*xi+0.811760009765625E2*xi*xi*xi-0.90195556640625E1*xi*xi-0.5450592041015625E1*xi+0.605621337890625;
15  double xi5= 0.2627363342285156E3*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2919292602539063E2*xi*xi*xi*xi*xi*xi*xi*xi-0.5319599853515625E3*xi*xi*xi*xi*xi*xi*xi-0.591066650390625E2*xi*xi*xi*xi*xi*xi+0.3449490600585938E3*xi*xi*xi*xi*xi+0.3832767333984375E2*xi*xi*xi*xi-0.811760009765625E2*xi*xi*xi-0.90195556640625E1*xi*xi+0.5450592041015625E1*xi+0.605621337890625;
16  double xi6= -0.1751575561523438E3*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.5838585205078125E2*xi*xi*xi*xi*xi*xi*xi*xi+0.337340478515625E3*xi*xi*xi*xi*xi*xi*xi+0.112446826171875E3*xi*xi*xi*xi*xi*xi-0.1968620361328125E3*xi*xi*xi*xi*xi-0.656206787109375E2*xi*xi*xi*xi+0.35082861328125E2*xi*xi*xi+0.11694287109375E2*xi*xi-0.40374755859375*xi-0.13458251953125;
17  double xi7= 0.7506752406529018E2*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.4170418003627232E2*xi*xi*xi*xi*xi*xi*xi*xi-0.129746337890625E3*xi*xi*xi*xi*xi*xi*xi-0.72081298828125E2*xi*xi*xi*xi*xi*xi+0.603881103515625E2*xi*xi*xi*xi*xi+0.335489501953125E2*xi*xi*xi*xi-0.5771589006696429E1*xi*xi*xi-0.3206438337053571E1*xi*xi+0.6229248046875E-1*xi+0.3460693359375E-1;
18  double xi8= -0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.1459646301269531E2*xi*xi*xi*xi*xi*xi*xi*xi+0.2687602713448661E2*xi*xi*xi*xi*xi*xi*xi+0.2090357666015625E2*xi*xi*xi*xi*xi*xi-0.8849981689453125E1*xi*xi*xi*xi*xi-0.6883319091796875E1*xi*xi*xi*xi+0.7487810407366071*xi*xi*xi+0.58238525390625*xi*xi-0.7945469447544643E-2*xi-0.61798095703125E-2;
19  double xi9= 0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi-0.5067836216517857E-1*xi*xi*xi-0.5067836216517857E-1*xi*xi+0.5340576171875E-3*xi+0.5340576171875E-3;
20 
21  double eta0= -0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta+0.216243896484375E1*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta-0.627374267578125*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta+0.5067836216517857E-1*eta*eta*eta-0.5067836216517857E-1*eta*eta-0.5340576171875E-3*eta+0.5340576171875E-3;
22  double eta1= 0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.1459646301269531E2*eta*eta*eta*eta*eta*eta*eta*eta-0.2687602713448661E2*eta*eta*eta*eta*eta*eta*eta+0.2090357666015625E2*eta*eta*eta*eta*eta*eta+0.8849981689453125E1*eta*eta*eta*eta*eta-0.6883319091796875E1*eta*eta*eta*eta-0.7487810407366071*eta*eta*eta+0.58238525390625*eta*eta+0.7945469447544643E-2*eta-0.61798095703125E-2;
23  double eta2= -0.7506752406529018E2*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.4170418003627232E2*eta*eta*eta*eta*eta*eta*eta*eta+0.129746337890625E3*eta*eta*eta*eta*eta*eta*eta-0.72081298828125E2*eta*eta*eta*eta*eta*eta-0.603881103515625E2*eta*eta*eta*eta*eta+0.335489501953125E2*eta*eta*eta*eta+0.5771589006696429E1*eta*eta*eta-0.3206438337053571E1*eta*eta-0.6229248046875E-1*eta+0.3460693359375E-1;
24  double eta3= 0.1751575561523438E3*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.5838585205078125E2*eta*eta*eta*eta*eta*eta*eta*eta-0.337340478515625E3*eta*eta*eta*eta*eta*eta*eta+0.112446826171875E3*eta*eta*eta*eta*eta*eta+0.1968620361328125E3*eta*eta*eta*eta*eta-0.656206787109375E2*eta*eta*eta*eta-0.35082861328125E2*eta*eta*eta+0.11694287109375E2*eta*eta+0.40374755859375*eta-0.13458251953125;
25  double eta4= -0.2627363342285156E3*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2919292602539063E2*eta*eta*eta*eta*eta*eta*eta*eta+0.5319599853515625E3*eta*eta*eta*eta*eta*eta*eta-0.591066650390625E2*eta*eta*eta*eta*eta*eta-0.3449490600585938E3*eta*eta*eta*eta*eta+0.3832767333984375E2*eta*eta*eta*eta+0.811760009765625E2*eta*eta*eta-0.90195556640625E1*eta*eta-0.5450592041015625E1*eta+0.605621337890625;
26  double eta5= 0.2627363342285156E3*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2919292602539063E2*eta*eta*eta*eta*eta*eta*eta*eta-0.5319599853515625E3*eta*eta*eta*eta*eta*eta*eta-0.591066650390625E2*eta*eta*eta*eta*eta*eta+0.3449490600585938E3*eta*eta*eta*eta*eta+0.3832767333984375E2*eta*eta*eta*eta-0.811760009765625E2*eta*eta*eta-0.90195556640625E1*eta*eta+0.5450592041015625E1*eta+0.605621337890625;
27  double eta6= -0.1751575561523438E3*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.5838585205078125E2*eta*eta*eta*eta*eta*eta*eta*eta+0.337340478515625E3*eta*eta*eta*eta*eta*eta*eta+0.112446826171875E3*eta*eta*eta*eta*eta*eta-0.1968620361328125E3*eta*eta*eta*eta*eta-0.656206787109375E2*eta*eta*eta*eta+0.35082861328125E2*eta*eta*eta+0.11694287109375E2*eta*eta-0.40374755859375*eta-0.13458251953125;
28  double eta7= 0.7506752406529018E2*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.4170418003627232E2*eta*eta*eta*eta*eta*eta*eta*eta-0.129746337890625E3*eta*eta*eta*eta*eta*eta*eta-0.72081298828125E2*eta*eta*eta*eta*eta*eta+0.603881103515625E2*eta*eta*eta*eta*eta+0.335489501953125E2*eta*eta*eta*eta-0.5771589006696429E1*eta*eta*eta-0.3206438337053571E1*eta*eta+0.6229248046875E-1*eta+0.3460693359375E-1;
29  double eta8= -0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.1459646301269531E2*eta*eta*eta*eta*eta*eta*eta*eta+0.2687602713448661E2*eta*eta*eta*eta*eta*eta*eta+0.2090357666015625E2*eta*eta*eta*eta*eta*eta-0.8849981689453125E1*eta*eta*eta*eta*eta-0.6883319091796875E1*eta*eta*eta*eta+0.7487810407366071*eta*eta*eta+0.58238525390625*eta*eta-0.7945469447544643E-2*eta-0.61798095703125E-2;
30  double eta9= 0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta-0.5067836216517857E-1*eta*eta*eta-0.5067836216517857E-1*eta*eta+0.5340576171875E-3*eta+0.5340576171875E-3;
31 
32 
33  values[0] = xi0*eta0;
34  values[1] = xi1*eta0;
35  values[2] = xi2*eta0;
36  values[3] = xi3*eta0;
37  values[4] = xi4*eta0;
38  values[5] = xi5*eta0;
39  values[6] = xi6*eta0;
40  values[7] = xi7*eta0;
41  values[8] = xi8*eta0;
42  values[9] = xi9*eta0;
43  values[10] = xi0*eta1;
44  values[11] = xi1*eta1;
45  values[12] = xi2*eta1;
46  values[13] = xi3*eta1;
47  values[14] = xi4*eta1;
48  values[15] = xi5*eta1;
49  values[16] = xi6*eta1;
50  values[17] = xi7*eta1;
51  values[18] = xi8*eta1;
52  values[19] = xi9*eta1;
53  values[20] = xi0*eta2;
54  values[21] = xi1*eta2;
55  values[22] = xi2*eta2;
56  values[23] = xi3*eta2;
57  values[24] = xi4*eta2;
58  values[25] = xi5*eta2;
59  values[26] = xi6*eta2;
60  values[27] = xi7*eta2;
61  values[28] = xi8*eta2;
62  values[29] = xi9*eta2;
63  values[30] = xi0*eta3;
64  values[31] = xi1*eta3;
65  values[32] = xi2*eta3;
66  values[33] = xi3*eta3;
67  values[34] = xi4*eta3;
68  values[35] = xi5*eta3;
69  values[36] = xi6*eta3;
70  values[37] = xi7*eta3;
71  values[38] = xi8*eta3;
72  values[39] = xi9*eta3;
73  values[40] = xi0*eta4;
74  values[41] = xi1*eta4;
75  values[42] = xi2*eta4;
76  values[43] = xi3*eta4;
77  values[44] = xi4*eta4;
78  values[45] = xi5*eta4;
79  values[46] = xi6*eta4;
80  values[47] = xi7*eta4;
81  values[48] = xi8*eta4;
82  values[49] = xi9*eta4;
83  values[50] = xi0*eta5;
84  values[51] = xi1*eta5;
85  values[52] = xi2*eta5;
86  values[53] = xi3*eta5;
87  values[54] = xi4*eta5;
88  values[55] = xi5*eta5;
89  values[56] = xi6*eta5;
90  values[57] = xi7*eta5;
91  values[58] = xi8*eta5;
92  values[59] = xi9*eta5;
93  values[60] = xi0*eta6;
94  values[61] = xi1*eta6;
95  values[62] = xi2*eta6;
96  values[63] = xi3*eta6;
97  values[64] = xi4*eta6;
98  values[65] = xi5*eta6;
99  values[66] = xi6*eta6;
100  values[67] = xi7*eta6;
101  values[68] = xi8*eta6;
102  values[69] = xi9*eta6;
103  values[70] = xi0*eta7;
104  values[71] = xi1*eta7;
105  values[72] = xi2*eta7;
106  values[73] = xi3*eta7;
107  values[74] = xi4*eta7;
108  values[75] = xi5*eta7;
109  values[76] = xi6*eta7;
110  values[77] = xi7*eta7;
111  values[78] = xi8*eta7;
112  values[79] = xi9*eta7;
113  values[80] = xi0*eta8;
114  values[81] = xi1*eta8;
115  values[82] = xi2*eta8;
116  values[83] = xi3*eta8;
117  values[84] = xi4*eta8;
118  values[85] = xi5*eta8;
119  values[86] = xi6*eta8;
120  values[87] = xi7*eta8;
121  values[88] = xi8*eta8;
122  values[89] = xi9*eta8;
123  values[90] = xi0*eta9;
124  values[91] = xi1*eta9;
125  values[92] = xi2*eta9;
126  values[93] = xi3*eta9;
127  values[94] = xi4*eta9;
128  values[95] = xi5*eta9;
129  values[96] = xi6*eta9;
130  values[97] = xi7*eta9;
131  values[98] = xi8*eta9;
132  values[99] = xi9*eta9;
133 }
134 
135 
136 // values of the derivatives in xi direction
137 static void C_Q_Q9_2D_DeriveXi(double xi, double eta, double *values)
138 {
139 
140  double xi0= -0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi+0.1668167201450893E2*xi*xi*xi*xi*xi*xi*xi+0.1513707275390625E2*xi*xi*xi*xi*xi*xi-0.129746337890625E2*xi*xi*xi*xi*xi-0.3136871337890625E1*xi*xi*xi*xi+0.25094970703125E1*xi*xi*xi+0.1520350864955357*xi*xi-0.1013567243303571*xi-0.5340576171875E-3;
141  double xi1= 0.1689019291469029E3*xi*xi*xi*xi*xi*xi*xi*xi-0.1167717041015625E3*xi*xi*xi*xi*xi*xi*xi-0.1881321899414063E3*xi*xi*xi*xi*xi*xi+0.1254214599609375E3*xi*xi*xi*xi*xi+0.4424990844726563E2*xi*xi*xi*xi-0.275332763671875E2*xi*xi*xi-0.2246343122209821E1*xi*xi+0.11647705078125E1*xi+0.7945469447544643E-2;
142  double xi2= -0.6756077165876116E3*xi*xi*xi*xi*xi*xi*xi*xi+0.3336334402901786E3*xi*xi*xi*xi*xi*xi*xi+0.908224365234375E3*xi*xi*xi*xi*xi*xi-0.43248779296875E3*xi*xi*xi*xi*xi-0.3019405517578125E3*xi*xi*xi*xi+0.13419580078125E3*xi*xi*xi+0.1731476702008929E2*xi*xi-0.6412876674107143E1*xi-0.6229248046875E-1;
143  double xi3= 0.1576418005371094E4*xi*xi*xi*xi*xi*xi*xi*xi-0.46708681640625E3*xi*xi*xi*xi*xi*xi*xi-0.2361383349609375E4*xi*xi*xi*xi*xi*xi+0.67468095703125E3*xi*xi*xi*xi*xi+0.9843101806640625E3*xi*xi*xi*xi-0.26248271484375E3*xi*xi*xi-0.105248583984375E3*xi*xi+0.2338857421875E2*xi+0.40374755859375;
144  double xi4= -0.2364627008056641E4*xi*xi*xi*xi*xi*xi*xi*xi+0.233543408203125E3*xi*xi*xi*xi*xi*xi*xi+0.3723719897460938E4*xi*xi*xi*xi*xi*xi-0.354639990234375E3*xi*xi*xi*xi*xi-0.1724745300292969E4*xi*xi*xi*xi+0.153310693359375E3*xi*xi*xi+0.2435280029296875E3*xi*xi-0.18039111328125E2*xi-0.5450592041015625E1;
145  double xi5= 0.2364627008056641E4*xi*xi*xi*xi*xi*xi*xi*xi+0.233543408203125E3*xi*xi*xi*xi*xi*xi*xi-0.3723719897460938E4*xi*xi*xi*xi*xi*xi-0.354639990234375E3*xi*xi*xi*xi*xi+0.1724745300292969E4*xi*xi*xi*xi+0.153310693359375E3*xi*xi*xi-0.2435280029296875E3*xi*xi-0.18039111328125E2*xi+0.5450592041015625E1;
146  double xi6= -0.1576418005371094E4*xi*xi*xi*xi*xi*xi*xi*xi-0.46708681640625E3*xi*xi*xi*xi*xi*xi*xi+0.2361383349609375E4*xi*xi*xi*xi*xi*xi+0.67468095703125E3*xi*xi*xi*xi*xi-0.9843101806640625E3*xi*xi*xi*xi-0.26248271484375E3*xi*xi*xi+0.105248583984375E3*xi*xi+0.2338857421875E2*xi-0.40374755859375;
147  double xi7= 0.6756077165876116E3*xi*xi*xi*xi*xi*xi*xi*xi+0.3336334402901786E3*xi*xi*xi*xi*xi*xi*xi-0.908224365234375E3*xi*xi*xi*xi*xi*xi-0.43248779296875E3*xi*xi*xi*xi*xi+0.3019405517578125E3*xi*xi*xi*xi+0.13419580078125E3*xi*xi*xi-0.1731476702008929E2*xi*xi-0.6412876674107143E1*xi+0.6229248046875E-1;
148  double xi8= -0.1689019291469029E3*xi*xi*xi*xi*xi*xi*xi*xi-0.1167717041015625E3*xi*xi*xi*xi*xi*xi*xi+0.1881321899414063E3*xi*xi*xi*xi*xi*xi+0.1254214599609375E3*xi*xi*xi*xi*xi-0.4424990844726563E2*xi*xi*xi*xi-0.275332763671875E2*xi*xi*xi+0.2246343122209821E1*xi*xi+0.11647705078125E1*xi-0.7945469447544643E-2;
149  double xi9= 0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi+0.1668167201450893E2*xi*xi*xi*xi*xi*xi*xi-0.1513707275390625E2*xi*xi*xi*xi*xi*xi-0.129746337890625E2*xi*xi*xi*xi*xi+0.3136871337890625E1*xi*xi*xi*xi+0.25094970703125E1*xi*xi*xi-0.1520350864955357*xi*xi-0.1013567243303571*xi+0.5340576171875E-3;
150 
151  double eta0= -0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta+0.216243896484375E1*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta-0.627374267578125*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta+0.5067836216517857E-1*eta*eta*eta-0.5067836216517857E-1*eta*eta-0.5340576171875E-3*eta+0.5340576171875E-3;
152  double eta1= 0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.1459646301269531E2*eta*eta*eta*eta*eta*eta*eta*eta-0.2687602713448661E2*eta*eta*eta*eta*eta*eta*eta+0.2090357666015625E2*eta*eta*eta*eta*eta*eta+0.8849981689453125E1*eta*eta*eta*eta*eta-0.6883319091796875E1*eta*eta*eta*eta-0.7487810407366071*eta*eta*eta+0.58238525390625*eta*eta+0.7945469447544643E-2*eta-0.61798095703125E-2;
153  double eta2= -0.7506752406529018E2*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.4170418003627232E2*eta*eta*eta*eta*eta*eta*eta*eta+0.129746337890625E3*eta*eta*eta*eta*eta*eta*eta-0.72081298828125E2*eta*eta*eta*eta*eta*eta-0.603881103515625E2*eta*eta*eta*eta*eta+0.335489501953125E2*eta*eta*eta*eta+0.5771589006696429E1*eta*eta*eta-0.3206438337053571E1*eta*eta-0.6229248046875E-1*eta+0.3460693359375E-1;
154  double eta3= 0.1751575561523438E3*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.5838585205078125E2*eta*eta*eta*eta*eta*eta*eta*eta-0.337340478515625E3*eta*eta*eta*eta*eta*eta*eta+0.112446826171875E3*eta*eta*eta*eta*eta*eta+0.1968620361328125E3*eta*eta*eta*eta*eta-0.656206787109375E2*eta*eta*eta*eta-0.35082861328125E2*eta*eta*eta+0.11694287109375E2*eta*eta+0.40374755859375*eta-0.13458251953125;
155  double eta4= -0.2627363342285156E3*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2919292602539063E2*eta*eta*eta*eta*eta*eta*eta*eta+0.5319599853515625E3*eta*eta*eta*eta*eta*eta*eta-0.591066650390625E2*eta*eta*eta*eta*eta*eta-0.3449490600585938E3*eta*eta*eta*eta*eta+0.3832767333984375E2*eta*eta*eta*eta+0.811760009765625E2*eta*eta*eta-0.90195556640625E1*eta*eta-0.5450592041015625E1*eta+0.605621337890625;
156  double eta5= 0.2627363342285156E3*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2919292602539063E2*eta*eta*eta*eta*eta*eta*eta*eta-0.5319599853515625E3*eta*eta*eta*eta*eta*eta*eta-0.591066650390625E2*eta*eta*eta*eta*eta*eta+0.3449490600585938E3*eta*eta*eta*eta*eta+0.3832767333984375E2*eta*eta*eta*eta-0.811760009765625E2*eta*eta*eta-0.90195556640625E1*eta*eta+0.5450592041015625E1*eta+0.605621337890625;
157  double eta6= -0.1751575561523438E3*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.5838585205078125E2*eta*eta*eta*eta*eta*eta*eta*eta+0.337340478515625E3*eta*eta*eta*eta*eta*eta*eta+0.112446826171875E3*eta*eta*eta*eta*eta*eta-0.1968620361328125E3*eta*eta*eta*eta*eta-0.656206787109375E2*eta*eta*eta*eta+0.35082861328125E2*eta*eta*eta+0.11694287109375E2*eta*eta-0.40374755859375*eta-0.13458251953125;
158  double eta7= 0.7506752406529018E2*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.4170418003627232E2*eta*eta*eta*eta*eta*eta*eta*eta-0.129746337890625E3*eta*eta*eta*eta*eta*eta*eta-0.72081298828125E2*eta*eta*eta*eta*eta*eta+0.603881103515625E2*eta*eta*eta*eta*eta+0.335489501953125E2*eta*eta*eta*eta-0.5771589006696429E1*eta*eta*eta-0.3206438337053571E1*eta*eta+0.6229248046875E-1*eta+0.3460693359375E-1;
159  double eta8= -0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.1459646301269531E2*eta*eta*eta*eta*eta*eta*eta*eta+0.2687602713448661E2*eta*eta*eta*eta*eta*eta*eta+0.2090357666015625E2*eta*eta*eta*eta*eta*eta-0.8849981689453125E1*eta*eta*eta*eta*eta-0.6883319091796875E1*eta*eta*eta*eta+0.7487810407366071*eta*eta*eta+0.58238525390625*eta*eta-0.7945469447544643E-2*eta-0.61798095703125E-2;
160  double eta9= 0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta-0.5067836216517857E-1*eta*eta*eta-0.5067836216517857E-1*eta*eta+0.5340576171875E-3*eta+0.5340576171875E-3;
161 
162 
163  values[0] = xi0*eta0;
164  values[1] = xi1*eta0;
165  values[2] = xi2*eta0;
166  values[3] = xi3*eta0;
167  values[4] = xi4*eta0;
168  values[5] = xi5*eta0;
169  values[6] = xi6*eta0;
170  values[7] = xi7*eta0;
171  values[8] = xi8*eta0;
172  values[9] = xi9*eta0;
173  values[10] = xi0*eta1;
174  values[11] = xi1*eta1;
175  values[12] = xi2*eta1;
176  values[13] = xi3*eta1;
177  values[14] = xi4*eta1;
178  values[15] = xi5*eta1;
179  values[16] = xi6*eta1;
180  values[17] = xi7*eta1;
181  values[18] = xi8*eta1;
182  values[19] = xi9*eta1;
183  values[20] = xi0*eta2;
184  values[21] = xi1*eta2;
185  values[22] = xi2*eta2;
186  values[23] = xi3*eta2;
187  values[24] = xi4*eta2;
188  values[25] = xi5*eta2;
189  values[26] = xi6*eta2;
190  values[27] = xi7*eta2;
191  values[28] = xi8*eta2;
192  values[29] = xi9*eta2;
193  values[30] = xi0*eta3;
194  values[31] = xi1*eta3;
195  values[32] = xi2*eta3;
196  values[33] = xi3*eta3;
197  values[34] = xi4*eta3;
198  values[35] = xi5*eta3;
199  values[36] = xi6*eta3;
200  values[37] = xi7*eta3;
201  values[38] = xi8*eta3;
202  values[39] = xi9*eta3;
203  values[40] = xi0*eta4;
204  values[41] = xi1*eta4;
205  values[42] = xi2*eta4;
206  values[43] = xi3*eta4;
207  values[44] = xi4*eta4;
208  values[45] = xi5*eta4;
209  values[46] = xi6*eta4;
210  values[47] = xi7*eta4;
211  values[48] = xi8*eta4;
212  values[49] = xi9*eta4;
213  values[50] = xi0*eta5;
214  values[51] = xi1*eta5;
215  values[52] = xi2*eta5;
216  values[53] = xi3*eta5;
217  values[54] = xi4*eta5;
218  values[55] = xi5*eta5;
219  values[56] = xi6*eta5;
220  values[57] = xi7*eta5;
221  values[58] = xi8*eta5;
222  values[59] = xi9*eta5;
223  values[60] = xi0*eta6;
224  values[61] = xi1*eta6;
225  values[62] = xi2*eta6;
226  values[63] = xi3*eta6;
227  values[64] = xi4*eta6;
228  values[65] = xi5*eta6;
229  values[66] = xi6*eta6;
230  values[67] = xi7*eta6;
231  values[68] = xi8*eta6;
232  values[69] = xi9*eta6;
233  values[70] = xi0*eta7;
234  values[71] = xi1*eta7;
235  values[72] = xi2*eta7;
236  values[73] = xi3*eta7;
237  values[74] = xi4*eta7;
238  values[75] = xi5*eta7;
239  values[76] = xi6*eta7;
240  values[77] = xi7*eta7;
241  values[78] = xi8*eta7;
242  values[79] = xi9*eta7;
243  values[80] = xi0*eta8;
244  values[81] = xi1*eta8;
245  values[82] = xi2*eta8;
246  values[83] = xi3*eta8;
247  values[84] = xi4*eta8;
248  values[85] = xi5*eta8;
249  values[86] = xi6*eta8;
250  values[87] = xi7*eta8;
251  values[88] = xi8*eta8;
252  values[89] = xi9*eta8;
253  values[90] = xi0*eta9;
254  values[91] = xi1*eta9;
255  values[92] = xi2*eta9;
256  values[93] = xi3*eta9;
257  values[94] = xi4*eta9;
258  values[95] = xi5*eta9;
259  values[96] = xi6*eta9;
260  values[97] = xi7*eta9;
261  values[98] = xi8*eta9;
262  values[99] = xi9*eta9;
263 }
264 
265 // values of the derivatives in eta direction
266 static void C_Q_Q9_2D_DeriveEta(double xi, double eta, double *values)
267 {
268  double xi0= -0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi+0.216243896484375E1*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi-0.627374267578125*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi+0.5067836216517857E-1*xi*xi*xi-0.5067836216517857E-1*xi*xi-0.5340576171875E-3*xi+0.5340576171875E-3;
269  double xi1= 0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.1459646301269531E2*xi*xi*xi*xi*xi*xi*xi*xi-0.2687602713448661E2*xi*xi*xi*xi*xi*xi*xi+0.2090357666015625E2*xi*xi*xi*xi*xi*xi+0.8849981689453125E1*xi*xi*xi*xi*xi-0.6883319091796875E1*xi*xi*xi*xi-0.7487810407366071*xi*xi*xi+0.58238525390625*xi*xi+0.7945469447544643E-2*xi-0.61798095703125E-2;
270  double xi2= -0.7506752406529018E2*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.4170418003627232E2*xi*xi*xi*xi*xi*xi*xi*xi+0.129746337890625E3*xi*xi*xi*xi*xi*xi*xi-0.72081298828125E2*xi*xi*xi*xi*xi*xi-0.603881103515625E2*xi*xi*xi*xi*xi+0.335489501953125E2*xi*xi*xi*xi+0.5771589006696429E1*xi*xi*xi-0.3206438337053571E1*xi*xi-0.6229248046875E-1*xi+0.3460693359375E-1;
271  double xi3= 0.1751575561523438E3*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.5838585205078125E2*xi*xi*xi*xi*xi*xi*xi*xi-0.337340478515625E3*xi*xi*xi*xi*xi*xi*xi+0.112446826171875E3*xi*xi*xi*xi*xi*xi+0.1968620361328125E3*xi*xi*xi*xi*xi-0.656206787109375E2*xi*xi*xi*xi-0.35082861328125E2*xi*xi*xi+0.11694287109375E2*xi*xi+0.40374755859375*xi-0.13458251953125;
272  double xi4= -0.2627363342285156E3*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2919292602539063E2*xi*xi*xi*xi*xi*xi*xi*xi+0.5319599853515625E3*xi*xi*xi*xi*xi*xi*xi-0.591066650390625E2*xi*xi*xi*xi*xi*xi-0.3449490600585938E3*xi*xi*xi*xi*xi+0.3832767333984375E2*xi*xi*xi*xi+0.811760009765625E2*xi*xi*xi-0.90195556640625E1*xi*xi-0.5450592041015625E1*xi+0.605621337890625;
273  double xi5= 0.2627363342285156E3*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2919292602539063E2*xi*xi*xi*xi*xi*xi*xi*xi-0.5319599853515625E3*xi*xi*xi*xi*xi*xi*xi-0.591066650390625E2*xi*xi*xi*xi*xi*xi+0.3449490600585938E3*xi*xi*xi*xi*xi+0.3832767333984375E2*xi*xi*xi*xi-0.811760009765625E2*xi*xi*xi-0.90195556640625E1*xi*xi+0.5450592041015625E1*xi+0.605621337890625;
274  double xi6= -0.1751575561523438E3*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.5838585205078125E2*xi*xi*xi*xi*xi*xi*xi*xi+0.337340478515625E3*xi*xi*xi*xi*xi*xi*xi+0.112446826171875E3*xi*xi*xi*xi*xi*xi-0.1968620361328125E3*xi*xi*xi*xi*xi-0.656206787109375E2*xi*xi*xi*xi+0.35082861328125E2*xi*xi*xi+0.11694287109375E2*xi*xi-0.40374755859375*xi-0.13458251953125;
275  double xi7= 0.7506752406529018E2*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.4170418003627232E2*xi*xi*xi*xi*xi*xi*xi*xi-0.129746337890625E3*xi*xi*xi*xi*xi*xi*xi-0.72081298828125E2*xi*xi*xi*xi*xi*xi+0.603881103515625E2*xi*xi*xi*xi*xi+0.335489501953125E2*xi*xi*xi*xi-0.5771589006696429E1*xi*xi*xi-0.3206438337053571E1*xi*xi+0.6229248046875E-1*xi+0.3460693359375E-1;
276  double xi8= -0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.1459646301269531E2*xi*xi*xi*xi*xi*xi*xi*xi+0.2687602713448661E2*xi*xi*xi*xi*xi*xi*xi+0.2090357666015625E2*xi*xi*xi*xi*xi*xi-0.8849981689453125E1*xi*xi*xi*xi*xi-0.6883319091796875E1*xi*xi*xi*xi+0.7487810407366071*xi*xi*xi+0.58238525390625*xi*xi-0.7945469447544643E-2*xi-0.61798095703125E-2;
277  double xi9= 0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi-0.5067836216517857E-1*xi*xi*xi-0.5067836216517857E-1*xi*xi+0.5340576171875E-3*xi+0.5340576171875E-3;
278 
279  double eta0= -0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta+0.1668167201450893E2*eta*eta*eta*eta*eta*eta*eta+0.1513707275390625E2*eta*eta*eta*eta*eta*eta-0.129746337890625E2*eta*eta*eta*eta*eta-0.3136871337890625E1*eta*eta*eta*eta+0.25094970703125E1*eta*eta*eta+0.1520350864955357*eta*eta-0.1013567243303571*eta-0.5340576171875E-3;
280  double eta1= 0.1689019291469029E3*eta*eta*eta*eta*eta*eta*eta*eta-0.1167717041015625E3*eta*eta*eta*eta*eta*eta*eta-0.1881321899414063E3*eta*eta*eta*eta*eta*eta+0.1254214599609375E3*eta*eta*eta*eta*eta+0.4424990844726563E2*eta*eta*eta*eta-0.275332763671875E2*eta*eta*eta-0.2246343122209821E1*eta*eta+0.11647705078125E1*eta+0.7945469447544643E-2;
281  double eta2= -0.6756077165876116E3*eta*eta*eta*eta*eta*eta*eta*eta+0.3336334402901786E3*eta*eta*eta*eta*eta*eta*eta+0.908224365234375E3*eta*eta*eta*eta*eta*eta-0.43248779296875E3*eta*eta*eta*eta*eta-0.3019405517578125E3*eta*eta*eta*eta+0.13419580078125E3*eta*eta*eta+0.1731476702008929E2*eta*eta-0.6412876674107143E1*eta-0.6229248046875E-1;
282  double eta3= 0.1576418005371094E4*eta*eta*eta*eta*eta*eta*eta*eta-0.46708681640625E3*eta*eta*eta*eta*eta*eta*eta-0.2361383349609375E4*eta*eta*eta*eta*eta*eta+0.67468095703125E3*eta*eta*eta*eta*eta+0.9843101806640625E3*eta*eta*eta*eta-0.26248271484375E3*eta*eta*eta-0.105248583984375E3*eta*eta+0.2338857421875E2*eta+0.40374755859375;
283  double eta4= -0.2364627008056641E4*eta*eta*eta*eta*eta*eta*eta*eta+0.233543408203125E3*eta*eta*eta*eta*eta*eta*eta+0.3723719897460938E4*eta*eta*eta*eta*eta*eta-0.354639990234375E3*eta*eta*eta*eta*eta-0.1724745300292969E4*eta*eta*eta*eta+0.153310693359375E3*eta*eta*eta+0.2435280029296875E3*eta*eta-0.18039111328125E2*eta-0.5450592041015625E1;
284  double eta5= 0.2364627008056641E4*eta*eta*eta*eta*eta*eta*eta*eta+0.233543408203125E3*eta*eta*eta*eta*eta*eta*eta-0.3723719897460938E4*eta*eta*eta*eta*eta*eta-0.354639990234375E3*eta*eta*eta*eta*eta+0.1724745300292969E4*eta*eta*eta*eta+0.153310693359375E3*eta*eta*eta-0.2435280029296875E3*eta*eta-0.18039111328125E2*eta+0.5450592041015625E1;
285  double eta6= -0.1576418005371094E4*eta*eta*eta*eta*eta*eta*eta*eta-0.46708681640625E3*eta*eta*eta*eta*eta*eta*eta+0.2361383349609375E4*eta*eta*eta*eta*eta*eta+0.67468095703125E3*eta*eta*eta*eta*eta-0.9843101806640625E3*eta*eta*eta*eta-0.26248271484375E3*eta*eta*eta+0.105248583984375E3*eta*eta+0.2338857421875E2*eta-0.40374755859375;
286  double eta7= 0.6756077165876116E3*eta*eta*eta*eta*eta*eta*eta*eta+0.3336334402901786E3*eta*eta*eta*eta*eta*eta*eta-0.908224365234375E3*eta*eta*eta*eta*eta*eta-0.43248779296875E3*eta*eta*eta*eta*eta+0.3019405517578125E3*eta*eta*eta*eta+0.13419580078125E3*eta*eta*eta-0.1731476702008929E2*eta*eta-0.6412876674107143E1*eta+0.6229248046875E-1;
287  double eta8= -0.1689019291469029E3*eta*eta*eta*eta*eta*eta*eta*eta-0.1167717041015625E3*eta*eta*eta*eta*eta*eta*eta+0.1881321899414063E3*eta*eta*eta*eta*eta*eta+0.1254214599609375E3*eta*eta*eta*eta*eta-0.4424990844726563E2*eta*eta*eta*eta-0.275332763671875E2*eta*eta*eta+0.2246343122209821E1*eta*eta+0.11647705078125E1*eta-0.7945469447544643E-2;
288  double eta9= 0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta+0.1668167201450893E2*eta*eta*eta*eta*eta*eta*eta-0.1513707275390625E2*eta*eta*eta*eta*eta*eta-0.129746337890625E2*eta*eta*eta*eta*eta+0.3136871337890625E1*eta*eta*eta*eta+0.25094970703125E1*eta*eta*eta-0.1520350864955357*eta*eta-0.1013567243303571*eta+0.5340576171875E-3;
289 
290 
291  values[0] = xi0*eta0;
292  values[1] = xi1*eta0;
293  values[2] = xi2*eta0;
294  values[3] = xi3*eta0;
295  values[4] = xi4*eta0;
296  values[5] = xi5*eta0;
297  values[6] = xi6*eta0;
298  values[7] = xi7*eta0;
299  values[8] = xi8*eta0;
300  values[9] = xi9*eta0;
301  values[10] = xi0*eta1;
302  values[11] = xi1*eta1;
303  values[12] = xi2*eta1;
304  values[13] = xi3*eta1;
305  values[14] = xi4*eta1;
306  values[15] = xi5*eta1;
307  values[16] = xi6*eta1;
308  values[17] = xi7*eta1;
309  values[18] = xi8*eta1;
310  values[19] = xi9*eta1;
311  values[20] = xi0*eta2;
312  values[21] = xi1*eta2;
313  values[22] = xi2*eta2;
314  values[23] = xi3*eta2;
315  values[24] = xi4*eta2;
316  values[25] = xi5*eta2;
317  values[26] = xi6*eta2;
318  values[27] = xi7*eta2;
319  values[28] = xi8*eta2;
320  values[29] = xi9*eta2;
321  values[30] = xi0*eta3;
322  values[31] = xi1*eta3;
323  values[32] = xi2*eta3;
324  values[33] = xi3*eta3;
325  values[34] = xi4*eta3;
326  values[35] = xi5*eta3;
327  values[36] = xi6*eta3;
328  values[37] = xi7*eta3;
329  values[38] = xi8*eta3;
330  values[39] = xi9*eta3;
331  values[40] = xi0*eta4;
332  values[41] = xi1*eta4;
333  values[42] = xi2*eta4;
334  values[43] = xi3*eta4;
335  values[44] = xi4*eta4;
336  values[45] = xi5*eta4;
337  values[46] = xi6*eta4;
338  values[47] = xi7*eta4;
339  values[48] = xi8*eta4;
340  values[49] = xi9*eta4;
341  values[50] = xi0*eta5;
342  values[51] = xi1*eta5;
343  values[52] = xi2*eta5;
344  values[53] = xi3*eta5;
345  values[54] = xi4*eta5;
346  values[55] = xi5*eta5;
347  values[56] = xi6*eta5;
348  values[57] = xi7*eta5;
349  values[58] = xi8*eta5;
350  values[59] = xi9*eta5;
351  values[60] = xi0*eta6;
352  values[61] = xi1*eta6;
353  values[62] = xi2*eta6;
354  values[63] = xi3*eta6;
355  values[64] = xi4*eta6;
356  values[65] = xi5*eta6;
357  values[66] = xi6*eta6;
358  values[67] = xi7*eta6;
359  values[68] = xi8*eta6;
360  values[69] = xi9*eta6;
361  values[70] = xi0*eta7;
362  values[71] = xi1*eta7;
363  values[72] = xi2*eta7;
364  values[73] = xi3*eta7;
365  values[74] = xi4*eta7;
366  values[75] = xi5*eta7;
367  values[76] = xi6*eta7;
368  values[77] = xi7*eta7;
369  values[78] = xi8*eta7;
370  values[79] = xi9*eta7;
371  values[80] = xi0*eta8;
372  values[81] = xi1*eta8;
373  values[82] = xi2*eta8;
374  values[83] = xi3*eta8;
375  values[84] = xi4*eta8;
376  values[85] = xi5*eta8;
377  values[86] = xi6*eta8;
378  values[87] = xi7*eta8;
379  values[88] = xi8*eta8;
380  values[89] = xi9*eta8;
381  values[90] = xi0*eta9;
382  values[91] = xi1*eta9;
383  values[92] = xi2*eta9;
384  values[93] = xi3*eta9;
385  values[94] = xi4*eta9;
386  values[95] = xi5*eta9;
387  values[96] = xi6*eta9;
388  values[97] = xi7*eta9;
389  values[98] = xi8*eta9;
390  values[99] = xi9*eta9;
391 }
392 
393 // values of the derivatives in xi-xi direction
394 static void C_Q_Q9_2D_DeriveXiXi(double xi, double eta, double *values)
395 {
396  double xi0= -0.1501350481305804E3*xi*xi*xi*xi*xi*xi*xi+0.1167717041015625E3*xi*xi*xi*xi*xi*xi+0.908224365234375E2*xi*xi*xi*xi*xi-0.648731689453125E2*xi*xi*xi*xi-0.125474853515625E2*xi*xi*xi+0.75284912109375E1*xi*xi+0.3040701729910714*xi-0.1013567243303571;
397  double xi1= 0.1351215433175223E4*xi*xi*xi*xi*xi*xi*xi-0.8174019287109375E3*xi*xi*xi*xi*xi*xi-0.1128793139648438E4*xi*xi*xi*xi*xi+0.6271072998046875E3*xi*xi*xi*xi+0.1769996337890625E3*xi*xi*xi-0.825998291015625E2*xi*xi-0.4492686244419643E1*xi+0.11647705078125E1;
398  double xi2= -0.5404861732700893E4*xi*xi*xi*xi*xi*xi*xi+0.233543408203125E4*xi*xi*xi*xi*xi*xi+0.544934619140625E4*xi*xi*xi*xi*xi-0.216243896484375E4*xi*xi*xi*xi-0.120776220703125E4*xi*xi*xi+0.40258740234375E3*xi*xi+0.3462953404017857E2*xi-0.6412876674107143E1;
399  double xi3= 0.1261134404296875E5*xi*xi*xi*xi*xi*xi*xi-0.326960771484375E4*xi*xi*xi*xi*xi*xi-0.1416830009765625E5*xi*xi*xi*xi*xi+0.337340478515625E4*xi*xi*xi*xi+0.393724072265625E4*xi*xi*xi-0.78744814453125E3*xi*xi-0.21049716796875E3*xi+0.2338857421875E2;
400  double xi4= -0.1891701606445313E5*xi*xi*xi*xi*xi*xi*xi+0.1634803857421875E4*xi*xi*xi*xi*xi*xi+0.2234231938476563E5*xi*xi*xi*xi*xi-0.1773199951171875E4*xi*xi*xi*xi-0.6898981201171875E4*xi*xi*xi+0.459932080078125E3*xi*xi+0.487056005859375E3*xi-0.18039111328125E2;
401  double xi5= 0.1891701606445313E5*xi*xi*xi*xi*xi*xi*xi+0.1634803857421875E4*xi*xi*xi*xi*xi*xi-0.2234231938476563E5*xi*xi*xi*xi*xi-0.1773199951171875E4*xi*xi*xi*xi+0.6898981201171875E4*xi*xi*xi+0.459932080078125E3*xi*xi-0.487056005859375E3*xi-0.18039111328125E2;
402  double xi6= -0.1261134404296875E5*xi*xi*xi*xi*xi*xi*xi-0.326960771484375E4*xi*xi*xi*xi*xi*xi+0.1416830009765625E5*xi*xi*xi*xi*xi+0.337340478515625E4*xi*xi*xi*xi-0.393724072265625E4*xi*xi*xi-0.78744814453125E3*xi*xi+0.21049716796875E3*xi+0.2338857421875E2;
403  double xi7= 0.5404861732700893E4*xi*xi*xi*xi*xi*xi*xi+0.233543408203125E4*xi*xi*xi*xi*xi*xi-0.544934619140625E4*xi*xi*xi*xi*xi-0.216243896484375E4*xi*xi*xi*xi+0.120776220703125E4*xi*xi*xi+0.40258740234375E3*xi*xi-0.3462953404017857E2*xi-0.6412876674107143E1;
404  double xi8= -0.1351215433175223E4*xi*xi*xi*xi*xi*xi*xi-0.8174019287109375E3*xi*xi*xi*xi*xi*xi+0.1128793139648438E4*xi*xi*xi*xi*xi+0.6271072998046875E3*xi*xi*xi*xi-0.1769996337890625E3*xi*xi*xi-0.825998291015625E2*xi*xi+0.4492686244419643E1*xi+0.11647705078125E1;
405  double xi9= 0.1501350481305804E3*xi*xi*xi*xi*xi*xi*xi+0.1167717041015625E3*xi*xi*xi*xi*xi*xi-0.908224365234375E2*xi*xi*xi*xi*xi-0.648731689453125E2*xi*xi*xi*xi+0.125474853515625E2*xi*xi*xi+0.75284912109375E1*xi*xi-0.3040701729910714*xi-0.1013567243303571;
406 
407  double eta0= -0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta+0.216243896484375E1*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta-0.627374267578125*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta+0.5067836216517857E-1*eta*eta*eta-0.5067836216517857E-1*eta*eta-0.5340576171875E-3*eta+0.5340576171875E-3;
408  double eta1= 0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.1459646301269531E2*eta*eta*eta*eta*eta*eta*eta*eta-0.2687602713448661E2*eta*eta*eta*eta*eta*eta*eta+0.2090357666015625E2*eta*eta*eta*eta*eta*eta+0.8849981689453125E1*eta*eta*eta*eta*eta-0.6883319091796875E1*eta*eta*eta*eta-0.7487810407366071*eta*eta*eta+0.58238525390625*eta*eta+0.7945469447544643E-2*eta-0.61798095703125E-2;
409  double eta2= -0.7506752406529018E2*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.4170418003627232E2*eta*eta*eta*eta*eta*eta*eta*eta+0.129746337890625E3*eta*eta*eta*eta*eta*eta*eta-0.72081298828125E2*eta*eta*eta*eta*eta*eta-0.603881103515625E2*eta*eta*eta*eta*eta+0.335489501953125E2*eta*eta*eta*eta+0.5771589006696429E1*eta*eta*eta-0.3206438337053571E1*eta*eta-0.6229248046875E-1*eta+0.3460693359375E-1;
410  double eta3= 0.1751575561523438E3*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.5838585205078125E2*eta*eta*eta*eta*eta*eta*eta*eta-0.337340478515625E3*eta*eta*eta*eta*eta*eta*eta+0.112446826171875E3*eta*eta*eta*eta*eta*eta+0.1968620361328125E3*eta*eta*eta*eta*eta-0.656206787109375E2*eta*eta*eta*eta-0.35082861328125E2*eta*eta*eta+0.11694287109375E2*eta*eta+0.40374755859375*eta-0.13458251953125;
411  double eta4= -0.2627363342285156E3*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2919292602539063E2*eta*eta*eta*eta*eta*eta*eta*eta+0.5319599853515625E3*eta*eta*eta*eta*eta*eta*eta-0.591066650390625E2*eta*eta*eta*eta*eta*eta-0.3449490600585938E3*eta*eta*eta*eta*eta+0.3832767333984375E2*eta*eta*eta*eta+0.811760009765625E2*eta*eta*eta-0.90195556640625E1*eta*eta-0.5450592041015625E1*eta+0.605621337890625;
412  double eta5= 0.2627363342285156E3*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2919292602539063E2*eta*eta*eta*eta*eta*eta*eta*eta-0.5319599853515625E3*eta*eta*eta*eta*eta*eta*eta-0.591066650390625E2*eta*eta*eta*eta*eta*eta+0.3449490600585938E3*eta*eta*eta*eta*eta+0.3832767333984375E2*eta*eta*eta*eta-0.811760009765625E2*eta*eta*eta-0.90195556640625E1*eta*eta+0.5450592041015625E1*eta+0.605621337890625;
413  double eta6= -0.1751575561523438E3*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.5838585205078125E2*eta*eta*eta*eta*eta*eta*eta*eta+0.337340478515625E3*eta*eta*eta*eta*eta*eta*eta+0.112446826171875E3*eta*eta*eta*eta*eta*eta-0.1968620361328125E3*eta*eta*eta*eta*eta-0.656206787109375E2*eta*eta*eta*eta+0.35082861328125E2*eta*eta*eta+0.11694287109375E2*eta*eta-0.40374755859375*eta-0.13458251953125;
414  double eta7= 0.7506752406529018E2*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.4170418003627232E2*eta*eta*eta*eta*eta*eta*eta*eta-0.129746337890625E3*eta*eta*eta*eta*eta*eta*eta-0.72081298828125E2*eta*eta*eta*eta*eta*eta+0.603881103515625E2*eta*eta*eta*eta*eta+0.335489501953125E2*eta*eta*eta*eta-0.5771589006696429E1*eta*eta*eta-0.3206438337053571E1*eta*eta+0.6229248046875E-1*eta+0.3460693359375E-1;
415  double eta8= -0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta*eta-0.1459646301269531E2*eta*eta*eta*eta*eta*eta*eta*eta+0.2687602713448661E2*eta*eta*eta*eta*eta*eta*eta+0.2090357666015625E2*eta*eta*eta*eta*eta*eta-0.8849981689453125E1*eta*eta*eta*eta*eta-0.6883319091796875E1*eta*eta*eta*eta+0.7487810407366071*eta*eta*eta+0.58238525390625*eta*eta-0.7945469447544643E-2*eta-0.61798095703125E-2;
416  double eta9= 0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta*eta+0.2085209001813616E1*eta*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta*eta-0.216243896484375E1*eta*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta*eta+0.627374267578125*eta*eta*eta*eta-0.5067836216517857E-1*eta*eta*eta-0.5067836216517857E-1*eta*eta+0.5340576171875E-3*eta+0.5340576171875E-3;
417 
418 
419  values[0] = xi0*eta0;
420  values[1] = xi1*eta0;
421  values[2] = xi2*eta0;
422  values[3] = xi3*eta0;
423  values[4] = xi4*eta0;
424  values[5] = xi5*eta0;
425  values[6] = xi6*eta0;
426  values[7] = xi7*eta0;
427  values[8] = xi8*eta0;
428  values[9] = xi9*eta0;
429  values[10] = xi0*eta1;
430  values[11] = xi1*eta1;
431  values[12] = xi2*eta1;
432  values[13] = xi3*eta1;
433  values[14] = xi4*eta1;
434  values[15] = xi5*eta1;
435  values[16] = xi6*eta1;
436  values[17] = xi7*eta1;
437  values[18] = xi8*eta1;
438  values[19] = xi9*eta1;
439  values[20] = xi0*eta2;
440  values[21] = xi1*eta2;
441  values[22] = xi2*eta2;
442  values[23] = xi3*eta2;
443  values[24] = xi4*eta2;
444  values[25] = xi5*eta2;
445  values[26] = xi6*eta2;
446  values[27] = xi7*eta2;
447  values[28] = xi8*eta2;
448  values[29] = xi9*eta2;
449  values[30] = xi0*eta3;
450  values[31] = xi1*eta3;
451  values[32] = xi2*eta3;
452  values[33] = xi3*eta3;
453  values[34] = xi4*eta3;
454  values[35] = xi5*eta3;
455  values[36] = xi6*eta3;
456  values[37] = xi7*eta3;
457  values[38] = xi8*eta3;
458  values[39] = xi9*eta3;
459  values[40] = xi0*eta4;
460  values[41] = xi1*eta4;
461  values[42] = xi2*eta4;
462  values[43] = xi3*eta4;
463  values[44] = xi4*eta4;
464  values[45] = xi5*eta4;
465  values[46] = xi6*eta4;
466  values[47] = xi7*eta4;
467  values[48] = xi8*eta4;
468  values[49] = xi9*eta4;
469  values[50] = xi0*eta5;
470  values[51] = xi1*eta5;
471  values[52] = xi2*eta5;
472  values[53] = xi3*eta5;
473  values[54] = xi4*eta5;
474  values[55] = xi5*eta5;
475  values[56] = xi6*eta5;
476  values[57] = xi7*eta5;
477  values[58] = xi8*eta5;
478  values[59] = xi9*eta5;
479  values[60] = xi0*eta6;
480  values[61] = xi1*eta6;
481  values[62] = xi2*eta6;
482  values[63] = xi3*eta6;
483  values[64] = xi4*eta6;
484  values[65] = xi5*eta6;
485  values[66] = xi6*eta6;
486  values[67] = xi7*eta6;
487  values[68] = xi8*eta6;
488  values[69] = xi9*eta6;
489  values[70] = xi0*eta7;
490  values[71] = xi1*eta7;
491  values[72] = xi2*eta7;
492  values[73] = xi3*eta7;
493  values[74] = xi4*eta7;
494  values[75] = xi5*eta7;
495  values[76] = xi6*eta7;
496  values[77] = xi7*eta7;
497  values[78] = xi8*eta7;
498  values[79] = xi9*eta7;
499  values[80] = xi0*eta8;
500  values[81] = xi1*eta8;
501  values[82] = xi2*eta8;
502  values[83] = xi3*eta8;
503  values[84] = xi4*eta8;
504  values[85] = xi5*eta8;
505  values[86] = xi6*eta8;
506  values[87] = xi7*eta8;
507  values[88] = xi8*eta8;
508  values[89] = xi9*eta8;
509  values[90] = xi0*eta9;
510  values[91] = xi1*eta9;
511  values[92] = xi2*eta9;
512  values[93] = xi3*eta9;
513  values[94] = xi4*eta9;
514  values[95] = xi5*eta9;
515  values[96] = xi6*eta9;
516  values[97] = xi7*eta9;
517  values[98] = xi8*eta9;
518  values[99] = xi9*eta9;
519 }
520 
521 // values of the derivatives in xi-eta direction
522 static void C_Q_Q9_2D_DeriveXiEta(double xi, double eta, double *values)
523 {
524  double xi0= -0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi+0.1668167201450893E2*xi*xi*xi*xi*xi*xi*xi+0.1513707275390625E2*xi*xi*xi*xi*xi*xi-0.129746337890625E2*xi*xi*xi*xi*xi-0.3136871337890625E1*xi*xi*xi*xi+0.25094970703125E1*xi*xi*xi+0.1520350864955357*xi*xi-0.1013567243303571*xi-0.5340576171875E-3;
525  double xi1= 0.1689019291469029E3*xi*xi*xi*xi*xi*xi*xi*xi-0.1167717041015625E3*xi*xi*xi*xi*xi*xi*xi-0.1881321899414063E3*xi*xi*xi*xi*xi*xi+0.1254214599609375E3*xi*xi*xi*xi*xi+0.4424990844726563E2*xi*xi*xi*xi-0.275332763671875E2*xi*xi*xi-0.2246343122209821E1*xi*xi+0.11647705078125E1*xi+0.7945469447544643E-2;
526  double xi2= -0.6756077165876116E3*xi*xi*xi*xi*xi*xi*xi*xi+0.3336334402901786E3*xi*xi*xi*xi*xi*xi*xi+0.908224365234375E3*xi*xi*xi*xi*xi*xi-0.43248779296875E3*xi*xi*xi*xi*xi-0.3019405517578125E3*xi*xi*xi*xi+0.13419580078125E3*xi*xi*xi+0.1731476702008929E2*xi*xi-0.6412876674107143E1*xi-0.6229248046875E-1;
527  double xi3= 0.1576418005371094E4*xi*xi*xi*xi*xi*xi*xi*xi-0.46708681640625E3*xi*xi*xi*xi*xi*xi*xi-0.2361383349609375E4*xi*xi*xi*xi*xi*xi+0.67468095703125E3*xi*xi*xi*xi*xi+0.9843101806640625E3*xi*xi*xi*xi-0.26248271484375E3*xi*xi*xi-0.105248583984375E3*xi*xi+0.2338857421875E2*xi+0.40374755859375;
528  double xi4= -0.2364627008056641E4*xi*xi*xi*xi*xi*xi*xi*xi+0.233543408203125E3*xi*xi*xi*xi*xi*xi*xi+0.3723719897460938E4*xi*xi*xi*xi*xi*xi-0.354639990234375E3*xi*xi*xi*xi*xi-0.1724745300292969E4*xi*xi*xi*xi+0.153310693359375E3*xi*xi*xi+0.2435280029296875E3*xi*xi-0.18039111328125E2*xi-0.5450592041015625E1;
529  double xi5= 0.2364627008056641E4*xi*xi*xi*xi*xi*xi*xi*xi+0.233543408203125E3*xi*xi*xi*xi*xi*xi*xi-0.3723719897460938E4*xi*xi*xi*xi*xi*xi-0.354639990234375E3*xi*xi*xi*xi*xi+0.1724745300292969E4*xi*xi*xi*xi+0.153310693359375E3*xi*xi*xi-0.2435280029296875E3*xi*xi-0.18039111328125E2*xi+0.5450592041015625E1;
530  double xi6= -0.1576418005371094E4*xi*xi*xi*xi*xi*xi*xi*xi-0.46708681640625E3*xi*xi*xi*xi*xi*xi*xi+0.2361383349609375E4*xi*xi*xi*xi*xi*xi+0.67468095703125E3*xi*xi*xi*xi*xi-0.9843101806640625E3*xi*xi*xi*xi-0.26248271484375E3*xi*xi*xi+0.105248583984375E3*xi*xi+0.2338857421875E2*xi-0.40374755859375;
531  double xi7= 0.6756077165876116E3*xi*xi*xi*xi*xi*xi*xi*xi+0.3336334402901786E3*xi*xi*xi*xi*xi*xi*xi-0.908224365234375E3*xi*xi*xi*xi*xi*xi-0.43248779296875E3*xi*xi*xi*xi*xi+0.3019405517578125E3*xi*xi*xi*xi+0.13419580078125E3*xi*xi*xi-0.1731476702008929E2*xi*xi-0.6412876674107143E1*xi+0.6229248046875E-1;
532  double xi8= -0.1689019291469029E3*xi*xi*xi*xi*xi*xi*xi*xi-0.1167717041015625E3*xi*xi*xi*xi*xi*xi*xi+0.1881321899414063E3*xi*xi*xi*xi*xi*xi+0.1254214599609375E3*xi*xi*xi*xi*xi-0.4424990844726563E2*xi*xi*xi*xi-0.275332763671875E2*xi*xi*xi+0.2246343122209821E1*xi*xi+0.11647705078125E1*xi-0.7945469447544643E-2;
533  double xi9= 0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi+0.1668167201450893E2*xi*xi*xi*xi*xi*xi*xi-0.1513707275390625E2*xi*xi*xi*xi*xi*xi-0.129746337890625E2*xi*xi*xi*xi*xi+0.3136871337890625E1*xi*xi*xi*xi+0.25094970703125E1*xi*xi*xi-0.1520350864955357*xi*xi-0.1013567243303571*xi+0.5340576171875E-3;
534 
535  double eta0= -0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta+0.1668167201450893E2*eta*eta*eta*eta*eta*eta*eta+0.1513707275390625E2*eta*eta*eta*eta*eta*eta-0.129746337890625E2*eta*eta*eta*eta*eta-0.3136871337890625E1*eta*eta*eta*eta+0.25094970703125E1*eta*eta*eta+0.1520350864955357*eta*eta-0.1013567243303571*eta-0.5340576171875E-3;
536  double eta1= 0.1689019291469029E3*eta*eta*eta*eta*eta*eta*eta*eta-0.1167717041015625E3*eta*eta*eta*eta*eta*eta*eta-0.1881321899414063E3*eta*eta*eta*eta*eta*eta+0.1254214599609375E3*eta*eta*eta*eta*eta+0.4424990844726563E2*eta*eta*eta*eta-0.275332763671875E2*eta*eta*eta-0.2246343122209821E1*eta*eta+0.11647705078125E1*eta+0.7945469447544643E-2;
537  double eta2= -0.6756077165876116E3*eta*eta*eta*eta*eta*eta*eta*eta+0.3336334402901786E3*eta*eta*eta*eta*eta*eta*eta+0.908224365234375E3*eta*eta*eta*eta*eta*eta-0.43248779296875E3*eta*eta*eta*eta*eta-0.3019405517578125E3*eta*eta*eta*eta+0.13419580078125E3*eta*eta*eta+0.1731476702008929E2*eta*eta-0.6412876674107143E1*eta-0.6229248046875E-1;
538  double eta3= 0.1576418005371094E4*eta*eta*eta*eta*eta*eta*eta*eta-0.46708681640625E3*eta*eta*eta*eta*eta*eta*eta-0.2361383349609375E4*eta*eta*eta*eta*eta*eta+0.67468095703125E3*eta*eta*eta*eta*eta+0.9843101806640625E3*eta*eta*eta*eta-0.26248271484375E3*eta*eta*eta-0.105248583984375E3*eta*eta+0.2338857421875E2*eta+0.40374755859375;
539  double eta4= -0.2364627008056641E4*eta*eta*eta*eta*eta*eta*eta*eta+0.233543408203125E3*eta*eta*eta*eta*eta*eta*eta+0.3723719897460938E4*eta*eta*eta*eta*eta*eta-0.354639990234375E3*eta*eta*eta*eta*eta-0.1724745300292969E4*eta*eta*eta*eta+0.153310693359375E3*eta*eta*eta+0.2435280029296875E3*eta*eta-0.18039111328125E2*eta-0.5450592041015625E1;
540  double eta5= 0.2364627008056641E4*eta*eta*eta*eta*eta*eta*eta*eta+0.233543408203125E3*eta*eta*eta*eta*eta*eta*eta-0.3723719897460938E4*eta*eta*eta*eta*eta*eta-0.354639990234375E3*eta*eta*eta*eta*eta+0.1724745300292969E4*eta*eta*eta*eta+0.153310693359375E3*eta*eta*eta-0.2435280029296875E3*eta*eta-0.18039111328125E2*eta+0.5450592041015625E1;
541  double eta6= -0.1576418005371094E4*eta*eta*eta*eta*eta*eta*eta*eta-0.46708681640625E3*eta*eta*eta*eta*eta*eta*eta+0.2361383349609375E4*eta*eta*eta*eta*eta*eta+0.67468095703125E3*eta*eta*eta*eta*eta-0.9843101806640625E3*eta*eta*eta*eta-0.26248271484375E3*eta*eta*eta+0.105248583984375E3*eta*eta+0.2338857421875E2*eta-0.40374755859375;
542  double eta7= 0.6756077165876116E3*eta*eta*eta*eta*eta*eta*eta*eta+0.3336334402901786E3*eta*eta*eta*eta*eta*eta*eta-0.908224365234375E3*eta*eta*eta*eta*eta*eta-0.43248779296875E3*eta*eta*eta*eta*eta+0.3019405517578125E3*eta*eta*eta*eta+0.13419580078125E3*eta*eta*eta-0.1731476702008929E2*eta*eta-0.6412876674107143E1*eta+0.6229248046875E-1;
543  double eta8= -0.1689019291469029E3*eta*eta*eta*eta*eta*eta*eta*eta-0.1167717041015625E3*eta*eta*eta*eta*eta*eta*eta+0.1881321899414063E3*eta*eta*eta*eta*eta*eta+0.1254214599609375E3*eta*eta*eta*eta*eta-0.4424990844726563E2*eta*eta*eta*eta-0.275332763671875E2*eta*eta*eta+0.2246343122209821E1*eta*eta+0.11647705078125E1*eta-0.7945469447544643E-2;
544  double eta9= 0.1876688101632254E2*eta*eta*eta*eta*eta*eta*eta*eta+0.1668167201450893E2*eta*eta*eta*eta*eta*eta*eta-0.1513707275390625E2*eta*eta*eta*eta*eta*eta-0.129746337890625E2*eta*eta*eta*eta*eta+0.3136871337890625E1*eta*eta*eta*eta+0.25094970703125E1*eta*eta*eta-0.1520350864955357*eta*eta-0.1013567243303571*eta+0.5340576171875E-3;
545 
546 
547  values[0] = xi0*eta0;
548  values[1] = xi1*eta0;
549  values[2] = xi2*eta0;
550  values[3] = xi3*eta0;
551  values[4] = xi4*eta0;
552  values[5] = xi5*eta0;
553  values[6] = xi6*eta0;
554  values[7] = xi7*eta0;
555  values[8] = xi8*eta0;
556  values[9] = xi9*eta0;
557  values[10] = xi0*eta1;
558  values[11] = xi1*eta1;
559  values[12] = xi2*eta1;
560  values[13] = xi3*eta1;
561  values[14] = xi4*eta1;
562  values[15] = xi5*eta1;
563  values[16] = xi6*eta1;
564  values[17] = xi7*eta1;
565  values[18] = xi8*eta1;
566  values[19] = xi9*eta1;
567  values[20] = xi0*eta2;
568  values[21] = xi1*eta2;
569  values[22] = xi2*eta2;
570  values[23] = xi3*eta2;
571  values[24] = xi4*eta2;
572  values[25] = xi5*eta2;
573  values[26] = xi6*eta2;
574  values[27] = xi7*eta2;
575  values[28] = xi8*eta2;
576  values[29] = xi9*eta2;
577  values[30] = xi0*eta3;
578  values[31] = xi1*eta3;
579  values[32] = xi2*eta3;
580  values[33] = xi3*eta3;
581  values[34] = xi4*eta3;
582  values[35] = xi5*eta3;
583  values[36] = xi6*eta3;
584  values[37] = xi7*eta3;
585  values[38] = xi8*eta3;
586  values[39] = xi9*eta3;
587  values[40] = xi0*eta4;
588  values[41] = xi1*eta4;
589  values[42] = xi2*eta4;
590  values[43] = xi3*eta4;
591  values[44] = xi4*eta4;
592  values[45] = xi5*eta4;
593  values[46] = xi6*eta4;
594  values[47] = xi7*eta4;
595  values[48] = xi8*eta4;
596  values[49] = xi9*eta4;
597  values[50] = xi0*eta5;
598  values[51] = xi1*eta5;
599  values[52] = xi2*eta5;
600  values[53] = xi3*eta5;
601  values[54] = xi4*eta5;
602  values[55] = xi5*eta5;
603  values[56] = xi6*eta5;
604  values[57] = xi7*eta5;
605  values[58] = xi8*eta5;
606  values[59] = xi9*eta5;
607  values[60] = xi0*eta6;
608  values[61] = xi1*eta6;
609  values[62] = xi2*eta6;
610  values[63] = xi3*eta6;
611  values[64] = xi4*eta6;
612  values[65] = xi5*eta6;
613  values[66] = xi6*eta6;
614  values[67] = xi7*eta6;
615  values[68] = xi8*eta6;
616  values[69] = xi9*eta6;
617  values[70] = xi0*eta7;
618  values[71] = xi1*eta7;
619  values[72] = xi2*eta7;
620  values[73] = xi3*eta7;
621  values[74] = xi4*eta7;
622  values[75] = xi5*eta7;
623  values[76] = xi6*eta7;
624  values[77] = xi7*eta7;
625  values[78] = xi8*eta7;
626  values[79] = xi9*eta7;
627  values[80] = xi0*eta8;
628  values[81] = xi1*eta8;
629  values[82] = xi2*eta8;
630  values[83] = xi3*eta8;
631  values[84] = xi4*eta8;
632  values[85] = xi5*eta8;
633  values[86] = xi6*eta8;
634  values[87] = xi7*eta8;
635  values[88] = xi8*eta8;
636  values[89] = xi9*eta8;
637  values[90] = xi0*eta9;
638  values[91] = xi1*eta9;
639  values[92] = xi2*eta9;
640  values[93] = xi3*eta9;
641  values[94] = xi4*eta9;
642  values[95] = xi5*eta9;
643  values[96] = xi6*eta9;
644  values[97] = xi7*eta9;
645  values[98] = xi8*eta9;
646  values[99] = xi9*eta9;
647 }
648 
649 // values of the derivatives in eta-eta direction
650 static void C_Q_Q9_2D_DeriveEtaEta(double xi, double eta, double *values)
651 {
652  double xi0= -0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi+0.216243896484375E1*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi-0.627374267578125*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi+0.5067836216517857E-1*xi*xi*xi-0.5067836216517857E-1*xi*xi-0.5340576171875E-3*xi+0.5340576171875E-3;
653  double xi1= 0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.1459646301269531E2*xi*xi*xi*xi*xi*xi*xi*xi-0.2687602713448661E2*xi*xi*xi*xi*xi*xi*xi+0.2090357666015625E2*xi*xi*xi*xi*xi*xi+0.8849981689453125E1*xi*xi*xi*xi*xi-0.6883319091796875E1*xi*xi*xi*xi-0.7487810407366071*xi*xi*xi+0.58238525390625*xi*xi+0.7945469447544643E-2*xi-0.61798095703125E-2;
654  double xi2= -0.7506752406529018E2*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.4170418003627232E2*xi*xi*xi*xi*xi*xi*xi*xi+0.129746337890625E3*xi*xi*xi*xi*xi*xi*xi-0.72081298828125E2*xi*xi*xi*xi*xi*xi-0.603881103515625E2*xi*xi*xi*xi*xi+0.335489501953125E2*xi*xi*xi*xi+0.5771589006696429E1*xi*xi*xi-0.3206438337053571E1*xi*xi-0.6229248046875E-1*xi+0.3460693359375E-1;
655  double xi3= 0.1751575561523438E3*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.5838585205078125E2*xi*xi*xi*xi*xi*xi*xi*xi-0.337340478515625E3*xi*xi*xi*xi*xi*xi*xi+0.112446826171875E3*xi*xi*xi*xi*xi*xi+0.1968620361328125E3*xi*xi*xi*xi*xi-0.656206787109375E2*xi*xi*xi*xi-0.35082861328125E2*xi*xi*xi+0.11694287109375E2*xi*xi+0.40374755859375*xi-0.13458251953125;
656  double xi4= -0.2627363342285156E3*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2919292602539063E2*xi*xi*xi*xi*xi*xi*xi*xi+0.5319599853515625E3*xi*xi*xi*xi*xi*xi*xi-0.591066650390625E2*xi*xi*xi*xi*xi*xi-0.3449490600585938E3*xi*xi*xi*xi*xi+0.3832767333984375E2*xi*xi*xi*xi+0.811760009765625E2*xi*xi*xi-0.90195556640625E1*xi*xi-0.5450592041015625E1*xi+0.605621337890625;
657  double xi5= 0.2627363342285156E3*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2919292602539063E2*xi*xi*xi*xi*xi*xi*xi*xi-0.5319599853515625E3*xi*xi*xi*xi*xi*xi*xi-0.591066650390625E2*xi*xi*xi*xi*xi*xi+0.3449490600585938E3*xi*xi*xi*xi*xi+0.3832767333984375E2*xi*xi*xi*xi-0.811760009765625E2*xi*xi*xi-0.90195556640625E1*xi*xi+0.5450592041015625E1*xi+0.605621337890625;
658  double xi6= -0.1751575561523438E3*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.5838585205078125E2*xi*xi*xi*xi*xi*xi*xi*xi+0.337340478515625E3*xi*xi*xi*xi*xi*xi*xi+0.112446826171875E3*xi*xi*xi*xi*xi*xi-0.1968620361328125E3*xi*xi*xi*xi*xi-0.656206787109375E2*xi*xi*xi*xi+0.35082861328125E2*xi*xi*xi+0.11694287109375E2*xi*xi-0.40374755859375*xi-0.13458251953125;
659  double xi7= 0.7506752406529018E2*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.4170418003627232E2*xi*xi*xi*xi*xi*xi*xi*xi-0.129746337890625E3*xi*xi*xi*xi*xi*xi*xi-0.72081298828125E2*xi*xi*xi*xi*xi*xi+0.603881103515625E2*xi*xi*xi*xi*xi+0.335489501953125E2*xi*xi*xi*xi-0.5771589006696429E1*xi*xi*xi-0.3206438337053571E1*xi*xi+0.6229248046875E-1*xi+0.3460693359375E-1;
660  double xi8= -0.1876688101632254E2*xi*xi*xi*xi*xi*xi*xi*xi*xi-0.1459646301269531E2*xi*xi*xi*xi*xi*xi*xi*xi+0.2687602713448661E2*xi*xi*xi*xi*xi*xi*xi+0.2090357666015625E2*xi*xi*xi*xi*xi*xi-0.8849981689453125E1*xi*xi*xi*xi*xi-0.6883319091796875E1*xi*xi*xi*xi+0.7487810407366071*xi*xi*xi+0.58238525390625*xi*xi-0.7945469447544643E-2*xi-0.61798095703125E-2;
661  double xi9= 0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi*xi+0.2085209001813616E1*xi*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi*xi-0.216243896484375E1*xi*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi*xi+0.627374267578125*xi*xi*xi*xi-0.5067836216517857E-1*xi*xi*xi-0.5067836216517857E-1*xi*xi+0.5340576171875E-3*xi+0.5340576171875E-3;
662 
663  double eta0= -0.1501350481305804E3*eta*eta*eta*eta*eta*eta*eta+0.1167717041015625E3*eta*eta*eta*eta*eta*eta+0.908224365234375E2*eta*eta*eta*eta*eta-0.648731689453125E2*eta*eta*eta*eta-0.125474853515625E2*eta*eta*eta+0.75284912109375E1*eta*eta+0.3040701729910714*eta-0.1013567243303571;
664  double eta1= 0.1351215433175223E4*eta*eta*eta*eta*eta*eta*eta-0.8174019287109375E3*eta*eta*eta*eta*eta*eta-0.1128793139648438E4*eta*eta*eta*eta*eta+0.6271072998046875E3*eta*eta*eta*eta+0.1769996337890625E3*eta*eta*eta-0.825998291015625E2*eta*eta-0.4492686244419643E1*eta+0.11647705078125E1;
665  double eta2= -0.5404861732700893E4*eta*eta*eta*eta*eta*eta*eta+0.233543408203125E4*eta*eta*eta*eta*eta*eta+0.544934619140625E4*eta*eta*eta*eta*eta-0.216243896484375E4*eta*eta*eta*eta-0.120776220703125E4*eta*eta*eta+0.40258740234375E3*eta*eta+0.3462953404017857E2*eta-0.6412876674107143E1;
666  double eta3= 0.1261134404296875E5*eta*eta*eta*eta*eta*eta*eta-0.326960771484375E4*eta*eta*eta*eta*eta*eta-0.1416830009765625E5*eta*eta*eta*eta*eta+0.337340478515625E4*eta*eta*eta*eta+0.393724072265625E4*eta*eta*eta-0.78744814453125E3*eta*eta-0.21049716796875E3*eta+0.2338857421875E2;
667  double eta4= -0.1891701606445313E5*eta*eta*eta*eta*eta*eta*eta+0.1634803857421875E4*eta*eta*eta*eta*eta*eta+0.2234231938476563E5*eta*eta*eta*eta*eta-0.1773199951171875E4*eta*eta*eta*eta-0.6898981201171875E4*eta*eta*eta+0.459932080078125E3*eta*eta+0.487056005859375E3*eta-0.18039111328125E2;
668  double eta5= 0.1891701606445313E5*eta*eta*eta*eta*eta*eta*eta+0.1634803857421875E4*eta*eta*eta*eta*eta*eta-0.2234231938476563E5*eta*eta*eta*eta*eta-0.1773199951171875E4*eta*eta*eta*eta+0.6898981201171875E4*eta*eta*eta+0.459932080078125E3*eta*eta-0.487056005859375E3*eta-0.18039111328125E2;
669  double eta6= -0.1261134404296875E5*eta*eta*eta*eta*eta*eta*eta-0.326960771484375E4*eta*eta*eta*eta*eta*eta+0.1416830009765625E5*eta*eta*eta*eta*eta+0.337340478515625E4*eta*eta*eta*eta-0.393724072265625E4*eta*eta*eta-0.78744814453125E3*eta*eta+0.21049716796875E3*eta+0.2338857421875E2;
670  double eta7= 0.5404861732700893E4*eta*eta*eta*eta*eta*eta*eta+0.233543408203125E4*eta*eta*eta*eta*eta*eta-0.544934619140625E4*eta*eta*eta*eta*eta-0.216243896484375E4*eta*eta*eta*eta+0.120776220703125E4*eta*eta*eta+0.40258740234375E3*eta*eta-0.3462953404017857E2*eta-0.6412876674107143E1;
671  double eta8= -0.1351215433175223E4*eta*eta*eta*eta*eta*eta*eta-0.8174019287109375E3*eta*eta*eta*eta*eta*eta+0.1128793139648438E4*eta*eta*eta*eta*eta+0.6271072998046875E3*eta*eta*eta*eta-0.1769996337890625E3*eta*eta*eta-0.825998291015625E2*eta*eta+0.4492686244419643E1*eta+0.11647705078125E1;
672  double eta9= 0.1501350481305804E3*eta*eta*eta*eta*eta*eta*eta+0.1167717041015625E3*eta*eta*eta*eta*eta*eta-0.908224365234375E2*eta*eta*eta*eta*eta-0.648731689453125E2*eta*eta*eta*eta+0.125474853515625E2*eta*eta*eta+0.75284912109375E1*eta*eta-0.3040701729910714*eta-0.1013567243303571;
673 
674 
675  values[0] = xi0*eta0;
676  values[1] = xi1*eta0;
677  values[2] = xi2*eta0;
678  values[3] = xi3*eta0;
679  values[4] = xi4*eta0;
680  values[5] = xi5*eta0;
681  values[6] = xi6*eta0;
682  values[7] = xi7*eta0;
683  values[8] = xi8*eta0;
684  values[9] = xi9*eta0;
685  values[10] = xi0*eta1;
686  values[11] = xi1*eta1;
687  values[12] = xi2*eta1;
688  values[13] = xi3*eta1;
689  values[14] = xi4*eta1;
690  values[15] = xi5*eta1;
691  values[16] = xi6*eta1;
692  values[17] = xi7*eta1;
693  values[18] = xi8*eta1;
694  values[19] = xi9*eta1;
695  values[20] = xi0*eta2;
696  values[21] = xi1*eta2;
697  values[22] = xi2*eta2;
698  values[23] = xi3*eta2;
699  values[24] = xi4*eta2;
700  values[25] = xi5*eta2;
701  values[26] = xi6*eta2;
702  values[27] = xi7*eta2;
703  values[28] = xi8*eta2;
704  values[29] = xi9*eta2;
705  values[30] = xi0*eta3;
706  values[31] = xi1*eta3;
707  values[32] = xi2*eta3;
708  values[33] = xi3*eta3;
709  values[34] = xi4*eta3;
710  values[35] = xi5*eta3;
711  values[36] = xi6*eta3;
712  values[37] = xi7*eta3;
713  values[38] = xi8*eta3;
714  values[39] = xi9*eta3;
715  values[40] = xi0*eta4;
716  values[41] = xi1*eta4;
717  values[42] = xi2*eta4;
718  values[43] = xi3*eta4;
719  values[44] = xi4*eta4;
720  values[45] = xi5*eta4;
721  values[46] = xi6*eta4;
722  values[47] = xi7*eta4;
723  values[48] = xi8*eta4;
724  values[49] = xi9*eta4;
725  values[50] = xi0*eta5;
726  values[51] = xi1*eta5;
727  values[52] = xi2*eta5;
728  values[53] = xi3*eta5;
729  values[54] = xi4*eta5;
730  values[55] = xi5*eta5;
731  values[56] = xi6*eta5;
732  values[57] = xi7*eta5;
733  values[58] = xi8*eta5;
734  values[59] = xi9*eta5;
735  values[60] = xi0*eta6;
736  values[61] = xi1*eta6;
737  values[62] = xi2*eta6;
738  values[63] = xi3*eta6;
739  values[64] = xi4*eta6;
740  values[65] = xi5*eta6;
741  values[66] = xi6*eta6;
742  values[67] = xi7*eta6;
743  values[68] = xi8*eta6;
744  values[69] = xi9*eta6;
745  values[70] = xi0*eta7;
746  values[71] = xi1*eta7;
747  values[72] = xi2*eta7;
748  values[73] = xi3*eta7;
749  values[74] = xi4*eta7;
750  values[75] = xi5*eta7;
751  values[76] = xi6*eta7;
752  values[77] = xi7*eta7;
753  values[78] = xi8*eta7;
754  values[79] = xi9*eta7;
755  values[80] = xi0*eta8;
756  values[81] = xi1*eta8;
757  values[82] = xi2*eta8;
758  values[83] = xi3*eta8;
759  values[84] = xi4*eta8;
760  values[85] = xi5*eta8;
761  values[86] = xi6*eta8;
762  values[87] = xi7*eta8;
763  values[88] = xi8*eta8;
764  values[89] = xi9*eta8;
765  values[90] = xi0*eta9;
766  values[91] = xi1*eta9;
767  values[92] = xi2*eta9;
768  values[93] = xi3*eta9;
769  values[94] = xi4*eta9;
770  values[95] = xi5*eta9;
771  values[96] = xi6*eta9;
772  values[97] = xi7*eta9;
773  values[98] = xi8*eta9;
774  values[99] = xi9*eta9;
775 }
776 
777 
778 // ***********************************************************************
779 
780 TBaseFunct2D *BF_C_Q_Q9_2D_Obj = new TBaseFunct2D
781  (100, BF_C_Q_Q9_2D, BFUnitSquare,
782  C_Q_Q9_2D_Funct, C_Q_Q9_2D_DeriveXi,
783  C_Q_Q9_2D_DeriveEta, C_Q_Q9_2D_DeriveXiXi,
784  C_Q_Q9_2D_DeriveXiEta, C_Q_Q9_2D_DeriveEtaEta, 9, 9,
785  0, NULL);
TBaseFunct2D
Definition: BaseFunct2D.h:27
Generated on Tue May 10 2016 09:07:57 for ParMooN by   doxygen 1.8.5