96 vector<Array<OneD, NekDouble> > xyz(4);
97 vector<NodeSharedPtr> ns =
m_el->GetVertexList();
98 for (
int i = 0; i < 4; i++)
100 Array<OneD, NekDouble> x(3);
116 J(0, 0) = -0.5 * b1 * xyz[0][0] + 0.5 * b1 * xyz[1][0] +
117 0.5 * b2 * xyz[2][0] - 0.5 * b2 * xyz[3][0];
118 J(1, 0) = -0.5 * b1 * xyz[0][1] + 0.5 * b1 * xyz[1][1] +
119 0.5 * b2 * xyz[2][1] - 0.5 * b2 * xyz[3][1];
121 J(0, 1) = -0.5 * a1 * xyz[0][0] - 0.5 * a2 * xyz[1][0] +
122 0.5 * a2 * xyz[2][0] + 0.5 * a1 * xyz[3][0];
123 J(1, 1) = -0.5 * a1 * xyz[0][1] - 0.5 * a2 * xyz[1][1] +
124 0.5 * a2 * xyz[2][1] + 0.5 * a1 * xyz[3][1];
128 Array<OneD, NekDouble> r(10, 0.0);
130 r[9] = 1.0 / (J(0, 0) * J(1, 1) - J(0, 1) * J(1, 0));
153 J(0, 0) = -0.5 * b1 * xyz[0][0] + 0.5 * b1 * xyz[1][0] +
154 0.5 * b2 * xyz[2][0] - 0.5 * b2 * xyz[3][0];
155 J(1, 0) = -0.5 * b1 * xyz[0][1] + 0.5 * b1 * xyz[1][1] +
156 0.5 * b2 * xyz[2][1] - 0.5 * b2 * xyz[3][1];
158 J(0, 1) = -0.5 * a1 * xyz[0][0] - 0.5 * a2 * xyz[1][0] +
159 0.5 * a2 * xyz[2][0] + 0.5 * a1 * xyz[3][0];
160 J(1, 1) = -0.5 * a1 * xyz[0][1] - 0.5 * a2 * xyz[1][1] +
161 0.5 * a2 * xyz[2][1] + 0.5 * a1 * xyz[3][1];
165 Array<OneD, NekDouble> r(10, 0.0);
167 r[9] = 1.0 / (J(0, 0) * J(1, 1) - J(0, 1) * J(1, 0));
199 Array<OneD, NekDouble> r(10, 0.0);
201 r[9] = 1.0 / (J(0, 0) * J(1, 1) - J(0, 1) * J(1, 0));
239 Array<OneD, NekDouble> r(10, 0.0);
241 r[9] = 1.0 / (J(0, 0) * (J(1, 1) * J(2, 2) - J(2, 1) * J(1, 2)) -
242 J(0, 1) * (J(1, 0) * J(2, 2) - J(2, 0) * J(1, 2)) +
243 J(0, 2) * (J(1, 0) * J(2, 1) - J(2, 0) * J(1, 1)));
263 Array<OneD, NekDouble> u1, v1, u2, v2, w1, w2;
267 vector<Array<OneD, NekDouble> > xyz(6);
268 vector<NodeSharedPtr> ns =
m_el->GetVertexList();
269 for (
int i = 0; i < 6; i++)
271 Array<OneD, NekDouble> x(3);
288 J(0, 0) = -0.5 * b1 * xyz[0][0] + 0.5 * b1 * xyz[1][0] +
289 0.5 * b2 * xyz[2][0] - 0.5 * b2 * xyz[3][0];
290 J(1, 0) = -0.5 * b1 * xyz[0][1] + 0.5 * b1 * xyz[1][1] +
291 0.5 * b2 * xyz[2][1] - 0.5 * b2 * xyz[3][1];
292 J(2, 0) = -0.5 * b1 * xyz[0][2] + 0.5 * b1 * xyz[1][2] +
293 0.5 * b2 * xyz[2][2] - 0.5 * b2 * xyz[3][2];
295 J(0, 1) = 0.5 * d * xyz[0][0] - 0.5 * a2 * xyz[1][0] +
296 0.5 * a2 * xyz[2][0] - 0.5 * d * xyz[3][0] -
297 0.5 * c2 * xyz[4][0] + 0.5 * c2 * xyz[5][0];
298 J(1, 1) = 0.5 * d * xyz[0][1] - 0.5 * a2 * xyz[1][1] +
299 0.5 * a2 * xyz[2][1] - 0.5 * d * xyz[3][1] -
300 0.5 * c2 * xyz[4][1] + 0.5 * c2 * xyz[5][1];
301 J(2, 1) = 0.5 * d * xyz[0][2] - 0.5 * a2 * xyz[1][2] +
302 0.5 * a2 * xyz[2][2] - 0.5 * d * xyz[3][2] -
303 0.5 * c2 * xyz[4][2] + 0.5 * c2 * xyz[5][2];
305 J(0, 2) = -0.5 * b1 * xyz[0][0] - 0.5 * b2 * xyz[3][0] +
306 0.5 * b1 * xyz[4][0] + 0.5 * b2 * xyz[5][0];
307 J(1, 2) = -0.5 * b1 * xyz[0][1] - 0.5 * b2 * xyz[3][1] +
308 0.5 * b1 * xyz[4][1] + 0.5 * b2 * xyz[5][1];
309 J(2, 2) = -0.5 * b1 * xyz[0][2] - 0.5 * b2 * xyz[3][2] +
310 0.5 * b1 * xyz[4][2] + 0.5 * b2 * xyz[5][2];
314 Array<OneD, NekDouble> r(10, 0.0);
316 r[9] = 1.0 / (J(0, 0) * (J(1, 1) * J(2, 2) - J(2, 1) * J(1, 2)) -
317 J(0, 1) * (J(1, 0) * J(2, 2) - J(2, 0) * J(1, 2)) +
318 J(0, 2) * (J(1, 0) * J(2, 1) - J(2, 0) * J(1, 1)));
341 J(0, 0) = -0.5 * b1 * xyz[0][0] + 0.5 * b1 * xyz[1][0] +
342 0.5 * b2 * xyz[2][0] - 0.5 * b2 * xyz[3][0];
343 J(1, 0) = -0.5 * b1 * xyz[0][1] + 0.5 * b1 * xyz[1][1] +
344 0.5 * b2 * xyz[2][1] - 0.5 * b2 * xyz[3][1];
345 J(2, 0) = -0.5 * b1 * xyz[0][2] + 0.5 * b1 * xyz[1][2] +
346 0.5 * b2 * xyz[2][2] - 0.5 * b2 * xyz[3][2];
348 J(0, 1) = 0.5 * d * xyz[0][0] - 0.5 * a2 * xyz[1][0] +
349 0.5 * a2 * xyz[2][0] - 0.5 * d * xyz[3][0] -
350 0.5 * c2 * xyz[4][0] + 0.5 * c2 * xyz[5][0];
351 J(1, 1) = 0.5 * d * xyz[0][1] - 0.5 * a2 * xyz[1][1] +
352 0.5 * a2 * xyz[2][1] - 0.5 * d * xyz[3][1] -
353 0.5 * c2 * xyz[4][1] + 0.5 * c2 * xyz[5][1];
354 J(2, 1) = 0.5 * d * xyz[0][2] - 0.5 * a2 * xyz[1][2] +
355 0.5 * a2 * xyz[2][2] - 0.5 * d * xyz[3][2] -
356 0.5 * c2 * xyz[4][2] + 0.5 * c2 * xyz[5][2];
358 J(0, 2) = -0.5 * b1 * xyz[0][0] - 0.5 * b2 * xyz[3][0] +
359 0.5 * b1 * xyz[4][0] + 0.5 * b2 * xyz[5][0];
360 J(1, 2) = -0.5 * b1 * xyz[0][1] - 0.5 * b2 * xyz[3][1] +
361 0.5 * b1 * xyz[4][1] + 0.5 * b2 * xyz[5][1];
362 J(2, 2) = -0.5 * b1 * xyz[0][2] - 0.5 * b2 * xyz[3][2] +
363 0.5 * b1 * xyz[4][2] + 0.5 * b2 * xyz[5][2];
367 Array<OneD, NekDouble> r(10, 0.0);
369 r[9] = 1.0 / (J(0, 0) * (J(1, 1) * J(2, 2) - J(2, 1) * J(1, 2)) -
370 J(0, 1) * (J(1, 0) * J(2, 2) - J(2, 0) * J(1, 2)) +
371 J(0, 2) * (J(1, 0) * J(2, 1) - J(2, 0) * J(1, 1)));
397 vector<Array<OneD, NekDouble> > xyz(8);
398 vector<NodeSharedPtr> ns =
m_el->GetVertexList();
399 for (
int i = 0; i < 8; i++)
401 Array<OneD, NekDouble> x(3);
419 J(0, 0) = -0.5 * b1 * c1 * xyz[0][0] + 0.5 * b1 * c1 * xyz[1][0] +
420 0.5 * b2 * c1 * xyz[2][0] - 0.5 * b2 * c1 * xyz[3][0] -
421 0.5 * b1 * c2 * xyz[5][0] + 0.5 * b1 * c2 * xyz[5][0] +
422 0.5 * b2 * c2 * xyz[6][0] - 0.5 * b2 * c2 * xyz[7][0];
423 J(1, 0) = -0.5 * b1 * c1 * xyz[0][1] + 0.5 * b1 * c1 * xyz[1][1] +
424 0.5 * b2 * c1 * xyz[2][1] - 0.5 * b2 * c1 * xyz[3][1] -
425 0.5 * b1 * c2 * xyz[5][1] + 0.5 * b1 * c2 * xyz[5][1] +
426 0.5 * b2 * c2 * xyz[6][1] - 0.5 * b2 * c2 * xyz[7][1];
427 J(2, 0) = -0.5 * b1 * c1 * xyz[0][2] + 0.5 * b1 * c1 * xyz[1][2] +
428 0.5 * b2 * c1 * xyz[2][2] - 0.5 * b2 * c1 * xyz[3][2] -
429 0.5 * b1 * c2 * xyz[5][2] + 0.5 * b1 * c2 * xyz[5][2] +
430 0.5 * b2 * c2 * xyz[6][2] - 0.5 * b2 * c2 * xyz[7][2];
432 J(0, 1) = -0.5 * a1 * c1 * xyz[0][0] - 0.5 * a2 * c1 * xyz[1][0] +
433 0.5 * a2 * c1 * xyz[2][0] + 0.5 * a1 * c1 * xyz[3][0] -
434 0.5 * a1 * c2 * xyz[5][0] - 0.5 * a2 * c2 * xyz[5][0] +
435 0.5 * a2 * c2 * xyz[6][0] + 0.5 * a1 * c2 * xyz[7][0];
436 J(1, 1) = -0.5 * a1 * c1 * xyz[0][1] - 0.5 * a2 * c1 * xyz[1][1] +
437 0.5 * a2 * c1 * xyz[2][1] + 0.5 * a1 * c1 * xyz[3][1] -
438 0.5 * a1 * c2 * xyz[5][1] - 0.5 * a2 * c2 * xyz[5][1] +
439 0.5 * a2 * c2 * xyz[6][1] + 0.5 * a1 * c2 * xyz[7][1];
440 J(2, 1) = -0.5 * a1 * c1 * xyz[0][2] - 0.5 * a2 * c1 * xyz[1][2] +
441 0.5 * a2 * c1 * xyz[2][2] + 0.5 * a1 * c1 * xyz[3][2] -
442 0.5 * a1 * c2 * xyz[5][2] - 0.5 * a2 * c2 * xyz[5][2] +
443 0.5 * a2 * c2 * xyz[6][2] + 0.5 * a1 * c2 * xyz[7][2];
445 J(0, 0) = -0.5 * b1 * a1 * xyz[0][0] - 0.5 * b1 * a2 * xyz[1][0] -
446 0.5 * b2 * a2 * xyz[2][0] - 0.5 * b2 * a1 * xyz[3][0] +
447 0.5 * b1 * a1 * xyz[5][0] + 0.5 * b1 * a2 * xyz[5][0] +
448 0.5 * b2 * a2 * xyz[6][0] + 0.5 * b2 * a1 * xyz[7][0];
449 J(1, 0) = -0.5 * b1 * a1 * xyz[0][1] - 0.5 * b1 * a2 * xyz[1][1] -
450 0.5 * b2 * a2 * xyz[2][1] - 0.5 * b2 * a1 * xyz[3][1] +
451 0.5 * b1 * a1 * xyz[5][1] + 0.5 * b1 * a2 * xyz[5][1] +
452 0.5 * b2 * a2 * xyz[6][1] + 0.5 * b2 * a1 * xyz[7][1];
453 J(2, 0) = -0.5 * b1 * a1 * xyz[0][2] - 0.5 * b1 * a2 * xyz[1][2] -
454 0.5 * b2 * a2 * xyz[2][2] - 0.5 * b2 * a1 * xyz[3][2] +
455 0.5 * b1 * a1 * xyz[5][2] + 0.5 * b1 * a2 * xyz[5][2] +
456 0.5 * b2 * a2 * xyz[6][2] + 0.5 * b2 * a1 * xyz[7][2];
460 Array<OneD, NekDouble> r(10, 0.0);
462 r[9] = 1.0 / (J(0, 0) * (J(1, 1) * J(2, 2) - J(2, 1) * J(1, 2)) -
463 J(0, 1) * (J(1, 0) * J(2, 2) - J(2, 0) * J(1, 2)) +
464 J(0, 2) * (J(1, 0) * J(2, 1) - J(2, 0) * J(1, 1)));
489 J(0, 0) = -0.5 * b1 * c1 * xyz[0][0] + 0.5 * b1 * c1 * xyz[1][0] +
490 0.5 * b2 * c1 * xyz[2][0] - 0.5 * b2 * c1 * xyz[3][0] -
491 0.5 * b1 * c2 * xyz[5][0] + 0.5 * b1 * c2 * xyz[5][0] +
492 0.5 * b2 * c2 * xyz[6][0] - 0.5 * b2 * c2 * xyz[7][0];
493 J(1, 0) = -0.5 * b1 * c1 * xyz[0][1] + 0.5 * b1 * c1 * xyz[1][1] +
494 0.5 * b2 * c1 * xyz[2][1] - 0.5 * b2 * c1 * xyz[3][1] -
495 0.5 * b1 * c2 * xyz[5][1] + 0.5 * b1 * c2 * xyz[5][1] +
496 0.5 * b2 * c2 * xyz[6][1] - 0.5 * b2 * c2 * xyz[7][1];
497 J(2, 0) = -0.5 * b1 * c1 * xyz[0][2] + 0.5 * b1 * c1 * xyz[1][2] +
498 0.5 * b2 * c1 * xyz[2][2] - 0.5 * b2 * c1 * xyz[3][2] -
499 0.5 * b1 * c2 * xyz[5][2] + 0.5 * b1 * c2 * xyz[5][2] +
500 0.5 * b2 * c2 * xyz[6][2] - 0.5 * b2 * c2 * xyz[7][2];
502 J(0, 1) = -0.5 * a1 * c1 * xyz[0][0] - 0.5 * a2 * c1 * xyz[1][0] +
503 0.5 * a2 * c1 * xyz[2][0] + 0.5 * a1 * c1 * xyz[3][0] -
504 0.5 * a1 * c2 * xyz[5][0] - 0.5 * a2 * c2 * xyz[5][0] +
505 0.5 * a2 * c2 * xyz[6][0] + 0.5 * a1 * c2 * xyz[7][0];
506 J(1, 1) = -0.5 * a1 * c1 * xyz[0][1] - 0.5 * a2 * c1 * xyz[1][1] +
507 0.5 * a2 * c1 * xyz[2][1] + 0.5 * a1 * c1 * xyz[3][1] -
508 0.5 * a1 * c2 * xyz[5][1] - 0.5 * a2 * c2 * xyz[5][1] +
509 0.5 * a2 * c2 * xyz[6][1] + 0.5 * a1 * c2 * xyz[7][1];
510 J(2, 1) = -0.5 * a1 * c1 * xyz[0][2] - 0.5 * a2 * c1 * xyz[1][2] +
511 0.5 * a2 * c1 * xyz[2][2] + 0.5 * a1 * c1 * xyz[3][2] -
512 0.5 * a1 * c2 * xyz[5][2] - 0.5 * a2 * c2 * xyz[5][2] +
513 0.5 * a2 * c2 * xyz[6][2] + 0.5 * a1 * c2 * xyz[7][2];
515 J(0, 0) = -0.5 * b1 * a1 * xyz[0][0] - 0.5 * b1 * a2 * xyz[1][0] -
516 0.5 * b2 * a2 * xyz[2][0] - 0.5 * b2 * a1 * xyz[3][0] +
517 0.5 * b1 * a1 * xyz[5][0] + 0.5 * b1 * a2 * xyz[5][0] +
518 0.5 * b2 * a2 * xyz[6][0] + 0.5 * b2 * a1 * xyz[7][0];
519 J(1, 0) = -0.5 * b1 * a1 * xyz[0][1] - 0.5 * b1 * a2 * xyz[1][1] -
520 0.5 * b2 * a2 * xyz[2][1] - 0.5 * b2 * a1 * xyz[3][1] +
521 0.5 * b1 * a1 * xyz[5][1] + 0.5 * b1 * a2 * xyz[5][1] +
522 0.5 * b2 * a2 * xyz[6][1] + 0.5 * b2 * a1 * xyz[7][1];
523 J(2, 0) = -0.5 * b1 * a1 * xyz[0][2] - 0.5 * b1 * a2 * xyz[1][2] -
524 0.5 * b2 * a2 * xyz[2][2] - 0.5 * b2 * a1 * xyz[3][2] +
525 0.5 * b1 * a1 * xyz[5][2] + 0.5 * b1 * a2 * xyz[5][2] +
526 0.5 * b2 * a2 * xyz[6][2] + 0.5 * b2 * a1 * xyz[7][2];
530 Array<OneD, NekDouble> r(10, 0.0);
532 r[9] = 1.0 / (J(0, 0) * (J(1, 1) * J(2, 2) - J(2, 1) * J(1, 2)) -
533 J(0, 1) * (J(1, 0) * J(2, 2) - J(2, 0) * J(1, 2)) +
534 J(0, 2) * (J(1, 0) * J(2, 1) - J(2, 0) * J(1, 1)));
#define ASSERTL0(condition, msg)
std::vector< std::vector< NekDouble * > > nodes
PointsManagerT & PointsManager(void)
NekMatrix< NekDouble, StandardMatrixTag > DNekMat
std::vector< Array< OneD, NekDouble > > mapsStd
std::vector< Array< OneD, NekDouble > > maps
DerivUtilSharedPtr m_derivUtil
3D electrostatically spaced points on a Prism