← Back to branch summary

~ubuntu-branches/ubuntu/raring/ffc/raring

« back to all changes in this revision

Viewing changes to test/regression/references/r_auto/VectorPoisson.h

  • Committer: Package Import Robot
  • Author(s): Johannes Ring
  • Date: 2011-10-26 17:52:20 UTC
  • mfrom: (1.1.10)
  • Revision ID: package-import@ubuntu.com-20111026175220-ope1dzqv4jn2b8pq
Tags: 1.0-beta2-1
* New upstream release. This release includes some performance
  improvements for evaluating basis functions. It also adds support
  for Bessel functions and error functions.
* debian/control: Bump version numbers for python-ufc, python-fiat,
  python-instant, and python-ufl.

Show diffs side-by-side

added added

removed removed

Lines of Context:
test/regression/references/r_auto/VectorPoisson.h
1
 
// This code conforms with the UFC specification version 2.0.0
2
 
// and was automatically generated by FFC version 1.0-beta.
 
1
// This code conforms with the UFC specification version 2.0.2
 
2
// and was automatically generated by FFC version 1.0-beta+.
3
3
// 
4
4
// This code was generated with the following parameters:
5
5
// 
130
130
      double basisvalues[3] = {0.0, 0.0, 0.0};
131
131
      
132
132
      // Declare helper variables.
133
 
      unsigned int rr = 0;
134
 
      unsigned int ss = 0;
135
133
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
136
134
      
137
135
      // Compute basisvalues.
138
136
      basisvalues[0] = 1.0;
139
137
      basisvalues[1] = tmp0;
140
 
      for (unsigned int r = 0; r < 1; r++)
141
 
      {
142
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
143
 
        ss = r*(r + 1)/2;
144
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
145
 
      }// end loop over 'r'
146
 
      for (unsigned int r = 0; r < 2; r++)
147
 
      {
148
 
        for (unsigned int s = 0; s < 2 - r; s++)
149
 
        {
150
 
          rr = (r + s)*(r + s + 1)/2 + s;
151
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
152
 
        }// end loop over 's'
153
 
      }// end loop over 'r'
 
138
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
139
      basisvalues[0] *= std::sqrt(0.5);
 
140
      basisvalues[2] *= std::sqrt(1.0);
 
141
      basisvalues[1] *= std::sqrt(3.0);
154
142
      
155
143
      // Table(s) of coefficients.
156
144
      static const double coefficients0[3] = \
170
158
      double basisvalues[3] = {0.0, 0.0, 0.0};
171
159
      
172
160
      // Declare helper variables.
173
 
      unsigned int rr = 0;
174
 
      unsigned int ss = 0;
175
161
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
176
162
      
177
163
      // Compute basisvalues.
178
164
      basisvalues[0] = 1.0;
179
165
      basisvalues[1] = tmp0;
180
 
      for (unsigned int r = 0; r < 1; r++)
181
 
      {
182
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
183
 
        ss = r*(r + 1)/2;
184
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
185
 
      }// end loop over 'r'
186
 
      for (unsigned int r = 0; r < 2; r++)
187
 
      {
188
 
        for (unsigned int s = 0; s < 2 - r; s++)
189
 
        {
190
 
          rr = (r + s)*(r + s + 1)/2 + s;
191
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
192
 
        }// end loop over 's'
193
 
      }// end loop over 'r'
 
166
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
167
      basisvalues[0] *= std::sqrt(0.5);
 
168
      basisvalues[2] *= std::sqrt(1.0);
 
169
      basisvalues[1] *= std::sqrt(3.0);
194
170
      
195
171
      // Table(s) of coefficients.
196
172
      static const double coefficients0[3] = \
210
186
      double basisvalues[3] = {0.0, 0.0, 0.0};
211
187
      
212
188
      // Declare helper variables.
213
 
      unsigned int rr = 0;
214
 
      unsigned int ss = 0;
215
189
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
216
190
      
217
191
      // Compute basisvalues.
218
192
      basisvalues[0] = 1.0;
219
193
      basisvalues[1] = tmp0;
220
 
      for (unsigned int r = 0; r < 1; r++)
221
 
      {
222
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
223
 
        ss = r*(r + 1)/2;
224
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
225
 
      }// end loop over 'r'
226
 
      for (unsigned int r = 0; r < 2; r++)
227
 
      {
228
 
        for (unsigned int s = 0; s < 2 - r; s++)
229
 
        {
230
 
          rr = (r + s)*(r + s + 1)/2 + s;
231
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
232
 
        }// end loop over 's'
233
 
      }// end loop over 'r'
 
194
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
195
      basisvalues[0] *= std::sqrt(0.5);
 
196
      basisvalues[2] *= std::sqrt(1.0);
 
197
      basisvalues[1] *= std::sqrt(3.0);
234
198
      
235
199
      // Table(s) of coefficients.
236
200
      static const double coefficients0[3] = \
369
333
      double basisvalues[3] = {0.0, 0.0, 0.0};
370
334
      
371
335
      // Declare helper variables.
372
 
      unsigned int rr = 0;
373
 
      unsigned int ss = 0;
374
336
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
375
337
      
376
338
      // Compute basisvalues.
377
339
      basisvalues[0] = 1.0;
378
340
      basisvalues[1] = tmp0;
379
 
      for (unsigned int r = 0; r < 1; r++)
380
 
      {
381
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
382
 
        ss = r*(r + 1)/2;
383
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
384
 
      }// end loop over 'r'
385
 
      for (unsigned int r = 0; r < 2; r++)
386
 
      {
387
 
        for (unsigned int s = 0; s < 2 - r; s++)
388
 
        {
389
 
          rr = (r + s)*(r + s + 1)/2 + s;
390
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
391
 
        }// end loop over 's'
392
 
      }// end loop over 'r'
 
341
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
342
      basisvalues[0] *= std::sqrt(0.5);
 
343
      basisvalues[2] *= std::sqrt(1.0);
 
344
      basisvalues[1] *= std::sqrt(3.0);
393
345
      
394
346
      // Table(s) of coefficients.
395
347
      static const double coefficients0[3] = \
527
479
      double basisvalues[3] = {0.0, 0.0, 0.0};
528
480
      
529
481
      // Declare helper variables.
530
 
      unsigned int rr = 0;
531
 
      unsigned int ss = 0;
532
482
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
533
483
      
534
484
      // Compute basisvalues.
535
485
      basisvalues[0] = 1.0;
536
486
      basisvalues[1] = tmp0;
537
 
      for (unsigned int r = 0; r < 1; r++)
538
 
      {
539
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
540
 
        ss = r*(r + 1)/2;
541
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
542
 
      }// end loop over 'r'
543
 
      for (unsigned int r = 0; r < 2; r++)
544
 
      {
545
 
        for (unsigned int s = 0; s < 2 - r; s++)
546
 
        {
547
 
          rr = (r + s)*(r + s + 1)/2 + s;
548
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
549
 
        }// end loop over 's'
550
 
      }// end loop over 'r'
 
487
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
488
      basisvalues[0] *= std::sqrt(0.5);
 
489
      basisvalues[2] *= std::sqrt(1.0);
 
490
      basisvalues[1] *= std::sqrt(3.0);
551
491
      
552
492
      // Table(s) of coefficients.
553
493
      static const double coefficients0[3] = \
685
625
      double basisvalues[3] = {0.0, 0.0, 0.0};
686
626
      
687
627
      // Declare helper variables.
688
 
      unsigned int rr = 0;
689
 
      unsigned int ss = 0;
690
628
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
691
629
      
692
630
      // Compute basisvalues.
693
631
      basisvalues[0] = 1.0;
694
632
      basisvalues[1] = tmp0;
695
 
      for (unsigned int r = 0; r < 1; r++)
696
 
      {
697
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
698
 
        ss = r*(r + 1)/2;
699
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
700
 
      }// end loop over 'r'
701
 
      for (unsigned int r = 0; r < 2; r++)
702
 
      {
703
 
        for (unsigned int s = 0; s < 2 - r; s++)
704
 
        {
705
 
          rr = (r + s)*(r + s + 1)/2 + s;
706
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
707
 
        }// end loop over 's'
708
 
      }// end loop over 'r'
 
633
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
634
      basisvalues[0] *= std::sqrt(0.5);
 
635
      basisvalues[2] *= std::sqrt(1.0);
 
636
      basisvalues[1] *= std::sqrt(3.0);
709
637
      
710
638
      // Table(s) of coefficients.
711
639
      static const double coefficients0[3] = \
1097
1025
      double basisvalues[3] = {0.0, 0.0, 0.0};
1098
1026
      
1099
1027
      // Declare helper variables.
1100
 
      unsigned int rr = 0;
1101
 
      unsigned int ss = 0;
1102
1028
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1103
1029
      
1104
1030
      // Compute basisvalues.
1105
1031
      basisvalues[0] = 1.0;
1106
1032
      basisvalues[1] = tmp0;
1107
 
      for (unsigned int r = 0; r < 1; r++)
1108
 
      {
1109
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1110
 
        ss = r*(r + 1)/2;
1111
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1112
 
      }// end loop over 'r'
1113
 
      for (unsigned int r = 0; r < 2; r++)
1114
 
      {
1115
 
        for (unsigned int s = 0; s < 2 - r; s++)
1116
 
        {
1117
 
          rr = (r + s)*(r + s + 1)/2 + s;
1118
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1119
 
        }// end loop over 's'
1120
 
      }// end loop over 'r'
 
1033
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1034
      basisvalues[0] *= std::sqrt(0.5);
 
1035
      basisvalues[2] *= std::sqrt(1.0);
 
1036
      basisvalues[1] *= std::sqrt(3.0);
1121
1037
      
1122
1038
      // Table(s) of coefficients.
1123
1039
      static const double coefficients0[3] = \
1137
1053
      double basisvalues[3] = {0.0, 0.0, 0.0};
1138
1054
      
1139
1055
      // Declare helper variables.
1140
 
      unsigned int rr = 0;
1141
 
      unsigned int ss = 0;
1142
1056
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1143
1057
      
1144
1058
      // Compute basisvalues.
1145
1059
      basisvalues[0] = 1.0;
1146
1060
      basisvalues[1] = tmp0;
1147
 
      for (unsigned int r = 0; r < 1; r++)
1148
 
      {
1149
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1150
 
        ss = r*(r + 1)/2;
1151
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1152
 
      }// end loop over 'r'
1153
 
      for (unsigned int r = 0; r < 2; r++)
1154
 
      {
1155
 
        for (unsigned int s = 0; s < 2 - r; s++)
1156
 
        {
1157
 
          rr = (r + s)*(r + s + 1)/2 + s;
1158
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1159
 
        }// end loop over 's'
1160
 
      }// end loop over 'r'
 
1061
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1062
      basisvalues[0] *= std::sqrt(0.5);
 
1063
      basisvalues[2] *= std::sqrt(1.0);
 
1064
      basisvalues[1] *= std::sqrt(3.0);
1161
1065
      
1162
1066
      // Table(s) of coefficients.
1163
1067
      static const double coefficients0[3] = \
1177
1081
      double basisvalues[3] = {0.0, 0.0, 0.0};
1178
1082
      
1179
1083
      // Declare helper variables.
1180
 
      unsigned int rr = 0;
1181
 
      unsigned int ss = 0;
1182
1084
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1183
1085
      
1184
1086
      // Compute basisvalues.
1185
1087
      basisvalues[0] = 1.0;
1186
1088
      basisvalues[1] = tmp0;
1187
 
      for (unsigned int r = 0; r < 1; r++)
1188
 
      {
1189
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1190
 
        ss = r*(r + 1)/2;
1191
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1192
 
      }// end loop over 'r'
1193
 
      for (unsigned int r = 0; r < 2; r++)
1194
 
      {
1195
 
        for (unsigned int s = 0; s < 2 - r; s++)
1196
 
        {
1197
 
          rr = (r + s)*(r + s + 1)/2 + s;
1198
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1199
 
        }// end loop over 's'
1200
 
      }// end loop over 'r'
 
1089
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1090
      basisvalues[0] *= std::sqrt(0.5);
 
1091
      basisvalues[2] *= std::sqrt(1.0);
 
1092
      basisvalues[1] *= std::sqrt(3.0);
1201
1093
      
1202
1094
      // Table(s) of coefficients.
1203
1095
      static const double coefficients0[3] = \
1217
1109
      double basisvalues[3] = {0.0, 0.0, 0.0};
1218
1110
      
1219
1111
      // Declare helper variables.
1220
 
      unsigned int rr = 0;
1221
 
      unsigned int ss = 0;
1222
1112
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1223
1113
      
1224
1114
      // Compute basisvalues.
1225
1115
      basisvalues[0] = 1.0;
1226
1116
      basisvalues[1] = tmp0;
1227
 
      for (unsigned int r = 0; r < 1; r++)
1228
 
      {
1229
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1230
 
        ss = r*(r + 1)/2;
1231
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1232
 
      }// end loop over 'r'
1233
 
      for (unsigned int r = 0; r < 2; r++)
1234
 
      {
1235
 
        for (unsigned int s = 0; s < 2 - r; s++)
1236
 
        {
1237
 
          rr = (r + s)*(r + s + 1)/2 + s;
1238
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1239
 
        }// end loop over 's'
1240
 
      }// end loop over 'r'
 
1117
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1118
      basisvalues[0] *= std::sqrt(0.5);
 
1119
      basisvalues[2] *= std::sqrt(1.0);
 
1120
      basisvalues[1] *= std::sqrt(3.0);
1241
1121
      
1242
1122
      // Table(s) of coefficients.
1243
1123
      static const double coefficients0[3] = \
1257
1137
      double basisvalues[3] = {0.0, 0.0, 0.0};
1258
1138
      
1259
1139
      // Declare helper variables.
1260
 
      unsigned int rr = 0;
1261
 
      unsigned int ss = 0;
1262
1140
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1263
1141
      
1264
1142
      // Compute basisvalues.
1265
1143
      basisvalues[0] = 1.0;
1266
1144
      basisvalues[1] = tmp0;
1267
 
      for (unsigned int r = 0; r < 1; r++)
1268
 
      {
1269
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1270
 
        ss = r*(r + 1)/2;
1271
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1272
 
      }// end loop over 'r'
1273
 
      for (unsigned int r = 0; r < 2; r++)
1274
 
      {
1275
 
        for (unsigned int s = 0; s < 2 - r; s++)
1276
 
        {
1277
 
          rr = (r + s)*(r + s + 1)/2 + s;
1278
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1279
 
        }// end loop over 's'
1280
 
      }// end loop over 'r'
 
1145
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1146
      basisvalues[0] *= std::sqrt(0.5);
 
1147
      basisvalues[2] *= std::sqrt(1.0);
 
1148
      basisvalues[1] *= std::sqrt(3.0);
1281
1149
      
1282
1150
      // Table(s) of coefficients.
1283
1151
      static const double coefficients0[3] = \
1297
1165
      double basisvalues[3] = {0.0, 0.0, 0.0};
1298
1166
      
1299
1167
      // Declare helper variables.
1300
 
      unsigned int rr = 0;
1301
 
      unsigned int ss = 0;
1302
1168
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1303
1169
      
1304
1170
      // Compute basisvalues.
1305
1171
      basisvalues[0] = 1.0;
1306
1172
      basisvalues[1] = tmp0;
1307
 
      for (unsigned int r = 0; r < 1; r++)
1308
 
      {
1309
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1310
 
        ss = r*(r + 1)/2;
1311
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1312
 
      }// end loop over 'r'
1313
 
      for (unsigned int r = 0; r < 2; r++)
1314
 
      {
1315
 
        for (unsigned int s = 0; s < 2 - r; s++)
1316
 
        {
1317
 
          rr = (r + s)*(r + s + 1)/2 + s;
1318
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1319
 
        }// end loop over 's'
1320
 
      }// end loop over 'r'
 
1173
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1174
      basisvalues[0] *= std::sqrt(0.5);
 
1175
      basisvalues[2] *= std::sqrt(1.0);
 
1176
      basisvalues[1] *= std::sqrt(3.0);
1321
1177
      
1322
1178
      // Table(s) of coefficients.
1323
1179
      static const double coefficients0[3] = \
1459
1315
      double basisvalues[3] = {0.0, 0.0, 0.0};
1460
1316
      
1461
1317
      // Declare helper variables.
1462
 
      unsigned int rr = 0;
1463
 
      unsigned int ss = 0;
1464
1318
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1465
1319
      
1466
1320
      // Compute basisvalues.
1467
1321
      basisvalues[0] = 1.0;
1468
1322
      basisvalues[1] = tmp0;
1469
 
      for (unsigned int r = 0; r < 1; r++)
1470
 
      {
1471
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1472
 
        ss = r*(r + 1)/2;
1473
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1474
 
      }// end loop over 'r'
1475
 
      for (unsigned int r = 0; r < 2; r++)
1476
 
      {
1477
 
        for (unsigned int s = 0; s < 2 - r; s++)
1478
 
        {
1479
 
          rr = (r + s)*(r + s + 1)/2 + s;
1480
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1481
 
        }// end loop over 's'
1482
 
      }// end loop over 'r'
 
1323
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1324
      basisvalues[0] *= std::sqrt(0.5);
 
1325
      basisvalues[2] *= std::sqrt(1.0);
 
1326
      basisvalues[1] *= std::sqrt(3.0);
1483
1327
      
1484
1328
      // Table(s) of coefficients.
1485
1329
      static const double coefficients0[3] = \
1617
1461
      double basisvalues[3] = {0.0, 0.0, 0.0};
1618
1462
      
1619
1463
      // Declare helper variables.
1620
 
      unsigned int rr = 0;
1621
 
      unsigned int ss = 0;
1622
1464
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1623
1465
      
1624
1466
      // Compute basisvalues.
1625
1467
      basisvalues[0] = 1.0;
1626
1468
      basisvalues[1] = tmp0;
1627
 
      for (unsigned int r = 0; r < 1; r++)
1628
 
      {
1629
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1630
 
        ss = r*(r + 1)/2;
1631
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1632
 
      }// end loop over 'r'
1633
 
      for (unsigned int r = 0; r < 2; r++)
1634
 
      {
1635
 
        for (unsigned int s = 0; s < 2 - r; s++)
1636
 
        {
1637
 
          rr = (r + s)*(r + s + 1)/2 + s;
1638
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1639
 
        }// end loop over 's'
1640
 
      }// end loop over 'r'
 
1469
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1470
      basisvalues[0] *= std::sqrt(0.5);
 
1471
      basisvalues[2] *= std::sqrt(1.0);
 
1472
      basisvalues[1] *= std::sqrt(3.0);
1641
1473
      
1642
1474
      // Table(s) of coefficients.
1643
1475
      static const double coefficients0[3] = \
1775
1607
      double basisvalues[3] = {0.0, 0.0, 0.0};
1776
1608
      
1777
1609
      // Declare helper variables.
1778
 
      unsigned int rr = 0;
1779
 
      unsigned int ss = 0;
1780
1610
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1781
1611
      
1782
1612
      // Compute basisvalues.
1783
1613
      basisvalues[0] = 1.0;
1784
1614
      basisvalues[1] = tmp0;
1785
 
      for (unsigned int r = 0; r < 1; r++)
1786
 
      {
1787
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1788
 
        ss = r*(r + 1)/2;
1789
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1790
 
      }// end loop over 'r'
1791
 
      for (unsigned int r = 0; r < 2; r++)
1792
 
      {
1793
 
        for (unsigned int s = 0; s < 2 - r; s++)
1794
 
        {
1795
 
          rr = (r + s)*(r + s + 1)/2 + s;
1796
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1797
 
        }// end loop over 's'
1798
 
      }// end loop over 'r'
 
1615
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1616
      basisvalues[0] *= std::sqrt(0.5);
 
1617
      basisvalues[2] *= std::sqrt(1.0);
 
1618
      basisvalues[1] *= std::sqrt(3.0);
1799
1619
      
1800
1620
      // Table(s) of coefficients.
1801
1621
      static const double coefficients0[3] = \
1933
1753
      double basisvalues[3] = {0.0, 0.0, 0.0};
1934
1754
      
1935
1755
      // Declare helper variables.
1936
 
      unsigned int rr = 0;
1937
 
      unsigned int ss = 0;
1938
1756
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1939
1757
      
1940
1758
      // Compute basisvalues.
1941
1759
      basisvalues[0] = 1.0;
1942
1760
      basisvalues[1] = tmp0;
1943
 
      for (unsigned int r = 0; r < 1; r++)
1944
 
      {
1945
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1946
 
        ss = r*(r + 1)/2;
1947
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1948
 
      }// end loop over 'r'
1949
 
      for (unsigned int r = 0; r < 2; r++)
1950
 
      {
1951
 
        for (unsigned int s = 0; s < 2 - r; s++)
1952
 
        {
1953
 
          rr = (r + s)*(r + s + 1)/2 + s;
1954
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1955
 
        }// end loop over 's'
1956
 
      }// end loop over 'r'
 
1761
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1762
      basisvalues[0] *= std::sqrt(0.5);
 
1763
      basisvalues[2] *= std::sqrt(1.0);
 
1764
      basisvalues[1] *= std::sqrt(3.0);
1957
1765
      
1958
1766
      // Table(s) of coefficients.
1959
1767
      static const double coefficients0[3] = \
2091
1899
      double basisvalues[3] = {0.0, 0.0, 0.0};
2092
1900
      
2093
1901
      // Declare helper variables.
2094
 
      unsigned int rr = 0;
2095
 
      unsigned int ss = 0;
2096
1902
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
2097
1903
      
2098
1904
      // Compute basisvalues.
2099
1905
      basisvalues[0] = 1.0;
2100
1906
      basisvalues[1] = tmp0;
2101
 
      for (unsigned int r = 0; r < 1; r++)
2102
 
      {
2103
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
2104
 
        ss = r*(r + 1)/2;
2105
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
2106
 
      }// end loop over 'r'
2107
 
      for (unsigned int r = 0; r < 2; r++)
2108
 
      {
2109
 
        for (unsigned int s = 0; s < 2 - r; s++)
2110
 
        {
2111
 
          rr = (r + s)*(r + s + 1)/2 + s;
2112
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
2113
 
        }// end loop over 's'
2114
 
      }// end loop over 'r'
 
1907
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1908
      basisvalues[0] *= std::sqrt(0.5);
 
1909
      basisvalues[2] *= std::sqrt(1.0);
 
1910
      basisvalues[1] *= std::sqrt(3.0);
2115
1911
      
2116
1912
      // Table(s) of coefficients.
2117
1913
      static const double coefficients0[3] = \
2249
2045
      double basisvalues[3] = {0.0, 0.0, 0.0};
2250
2046
      
2251
2047
      // Declare helper variables.
2252
 
      unsigned int rr = 0;
2253
 
      unsigned int ss = 0;
2254
2048
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
2255
2049
      
2256
2050
      // Compute basisvalues.
2257
2051
      basisvalues[0] = 1.0;
2258
2052
      basisvalues[1] = tmp0;
2259
 
      for (unsigned int r = 0; r < 1; r++)
2260
 
      {
2261
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
2262
 
        ss = r*(r + 1)/2;
2263
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
2264
 
      }// end loop over 'r'
2265
 
      for (unsigned int r = 0; r < 2; r++)
2266
 
      {
2267
 
        for (unsigned int s = 0; s < 2 - r; s++)
2268
 
        {
2269
 
          rr = (r + s)*(r + s + 1)/2 + s;
2270
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
2271
 
        }// end loop over 's'
2272
 
      }// end loop over 'r'
 
2053
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
2054
      basisvalues[0] *= std::sqrt(0.5);
 
2055
      basisvalues[2] *= std::sqrt(1.0);
 
2056
      basisvalues[1] *= std::sqrt(3.0);
2273
2057
      
2274
2058
      // Table(s) of coefficients.
2275
2059
      static const double coefficients0[3] = \