v0.15.0
Loading...
Searching...
No Matches
Functions
tetrahedron_nco_rule.c File Reference
#include "tetrahedron_nco_rule.h"

Go to the source code of this file.

Functions

double r8mat_det_4d (double a[])
 
void reference_to_physical_t4 (double t[], int n, double ref[], double phy[])
 
int tetrahedron_nco_degree (int rule)
 
int tetrahedron_nco_order_num (int rule)
 
void tetrahedron_nco_rule (int rule, int order_num, double xyz[], double w[])
 
int tetrahedron_nco_rule_num ()
 
int * tetrahedron_nco_suborder (int rule, int suborder_num)
 
int tetrahedron_nco_suborder_num (int rule)
 
void tetrahedron_nco_subrule (int rule, int suborder_num, double suborder_xyz[], double suborder_w[])
 
void tetrahedron_nco_subrule_01 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
void tetrahedron_nco_subrule_02 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
void tetrahedron_nco_subrule_03 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
void tetrahedron_nco_subrule_04 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
void tetrahedron_nco_subrule_05 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
void tetrahedron_nco_subrule_06 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
void tetrahedron_nco_subrule_07 (int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
 
double tetrahedron_volume (double tetra[3 *4])
 

Function Documentation

◆ r8mat_det_4d()

double r8mat_det_4d ( double  a[])

Definition at line 6 of file tetrahedron_nco_rule.c.

37{
38 double det;
39
40 det =
41 a[0+0*4] * (
42 a[1+1*4] * ( a[2+2*4] * a[3+3*4] - a[2+3*4] * a[3+2*4] )
43 - a[1+2*4] * ( a[2+1*4] * a[3+3*4] - a[2+3*4] * a[3+1*4] )
44 + a[1+3*4] * ( a[2+1*4] * a[3+2*4] - a[2+2*4] * a[3+1*4] ) )
45 - a[0+1*4] * (
46 a[1+0*4] * ( a[2+2*4] * a[3+3*4] - a[2+3*4] * a[3+2*4] )
47 - a[1+2*4] * ( a[2+0*4] * a[3+3*4] - a[2+3*4] * a[3+0*4] )
48 + a[1+3*4] * ( a[2+0*4] * a[3+2*4] - a[2+2*4] * a[3+0*4] ) )
49 + a[0+2*4] * (
50 a[1+0*4] * ( a[2+1*4] * a[3+3*4] - a[2+3*4] * a[3+1*4] )
51 - a[1+1*4] * ( a[2+0*4] * a[3+3*4] - a[2+3*4] * a[3+0*4] )
52 + a[1+3*4] * ( a[2+0*4] * a[3+1*4] - a[2+1*4] * a[3+0*4] ) )
53 - a[0+3*4] * (
54 a[1+0*4] * ( a[2+1*4] * a[3+2*4] - a[2+2*4] * a[3+1*4] )
55 - a[1+1*4] * ( a[2+0*4] * a[3+2*4] - a[2+2*4] * a[3+0*4] )
56 + a[1+2*4] * ( a[2+0*4] * a[3+1*4] - a[2+1*4] * a[3+0*4] ) );
57
58 return det;
59}
a
constexpr double a
Definition approx_sphere.cpp:30

◆ reference_to_physical_t4()

void reference_to_physical_t4 ( double  t[],
int  n,
double  ref[],
double  phy[] 
)

Definition at line 62 of file tetrahedron_nco_rule.c.

108{
109 int i;
110 int j;
111
112 for ( i = 0; i < 3; i++ )
113 {
114 for ( j = 0; j < n; j++ )
115 {
116 phy[i+j*3] = t[i+0*3] * ( 1.0 - ref[0+j*3] - ref[1+j*3] - ref[2+j*3] )
117 + t[i+1*3] * + ref[0+j*3]
118 + t[i+2*3] * + ref[1+j*3]
119 + t[i+3*3] * + ref[2+j*3];
120 }
121 }
122
123 return;
124}
i
FTensor::Index< 'i', SPACE_DIM > i
Definition hcurl_divergence_operator_2d.cpp:27
n
const double n
refractive index of diffusive medium
Definition initial_diffusion.cpp:38
j
FTensor::Index< 'j', 3 > j
Definition matrix_function.cpp:19
t
constexpr double t
plate stiffness
Definition plate.cpp:58

◆ tetrahedron_nco_degree()

int tetrahedron_nco_degree ( int  rule)

Definition at line 127 of file tetrahedron_nco_rule.c.

161{
162 int degree;
163
164 if ( 1 <= rule && rule <= 7 )
165 {
166 degree = rule - 1;
167 }
168 else
169 {
170 degree = -1;
171 fprintf ( stderr, "\n" );
172 fprintf ( stderr, "TETRAHEDRON_NCO_DEGREE - Fatal error!\n" );
173 fprintf ( stderr, " Illegal RULE = %d\n", rule );
174 exit ( 1 );
175 }
176
177 return degree;
178}

◆ tetrahedron_nco_order_num()

int tetrahedron_nco_order_num ( int  rule)

Definition at line 181 of file tetrahedron_nco_rule.c.

215{
216 int order;
217 int order_num;
218 int *suborder;
219 int suborder_num;
220
221 suborder_num = tetrahedron_nco_suborder_num ( rule );
222
223 suborder = tetrahedron_nco_suborder ( rule, suborder_num );
224
225 order_num = 0;
226 for ( order = 0; order < suborder_num; order++ )
227 {
228 order_num = order_num + suborder[order];
229 }
230
231 free ( suborder );
232
233 return order_num;
234}
order
constexpr int order
Definition dg_projection.cpp:19
tetrahedron_nco_suborder_num
int tetrahedron_nco_suborder_num(int rule)
Definition tetrahedron_nco_rule.c:811
tetrahedron_nco_suborder
int * tetrahedron_nco_suborder(int rule, int suborder_num)
Definition tetrahedron_nco_rule.c:710

◆ tetrahedron_nco_rule()

void tetrahedron_nco_rule ( int  rule,
int  order_num,
double  xyz[],
double  w[] 
)

Definition at line 237 of file tetrahedron_nco_rule.c.

274{
275 int o;
276 int s;
277 int *suborder;
278 int suborder_num;
279 double *suborder_w;
280 double *suborder_xyz;
281/*
282 Get the suborder information.
283*/
284 suborder_num = tetrahedron_nco_suborder_num ( rule );
285
286 suborder_xyz = ( double * ) malloc ( 4 * suborder_num * sizeof ( double ) );
287 suborder_w = ( double * ) malloc ( suborder_num * sizeof ( double ) );
288
289 suborder = tetrahedron_nco_suborder ( rule, suborder_num );
290
291 tetrahedron_nco_subrule ( rule, suborder_num, suborder_xyz, suborder_w );
292/*
293 Expand the suborder information to a full order rule.
294*/
295 o = 0;
296
297 for ( s = 0; s < suborder_num; s++ )
298 {
299 if ( suborder[s] == 1 )
300 {
301 xyz[0+o*3] = suborder_xyz[0+s*4];
302 xyz[1+o*3] = suborder_xyz[1+s*4];
303 xyz[2+o*3] = suborder_xyz[2+s*4];
304 w[o] = suborder_w[s];
305 o = o + 1;
306 }
307/*
308 Fourfold symmetry on (A,A,A,B)
309
310 123 AAA
311 124 AAB
312 142 ABA
313 412 BAA
314*/
315 else if ( suborder[s] == 4 )
316 {
317 xyz[0+o*3] = suborder_xyz[0+s*4];
318 xyz[1+o*3] = suborder_xyz[1+s*4];
319 xyz[2+o*3] = suborder_xyz[2+s*4];
320 w[o] = suborder_w[s];
321 o = o + 1;
322
323 xyz[0+o*3] = suborder_xyz[0+s*4];
324 xyz[1+o*3] = suborder_xyz[1+s*4];
325 xyz[2+o*3] = suborder_xyz[3+s*4];
326 w[o] = suborder_w[s];
327 o = o + 1;
328
329 xyz[0+o*3] = suborder_xyz[0+s*4];
330 xyz[1+o*3] = suborder_xyz[3+s*4];
331 xyz[2+o*3] = suborder_xyz[1+s*4];
332 w[o] = suborder_w[s];
333 o = o + 1;
334
335 xyz[0+o*3] = suborder_xyz[3+s*4];
336 xyz[1+o*3] = suborder_xyz[0+s*4];
337 xyz[2+o*3] = suborder_xyz[1+s*4];
338 w[o] = suborder_w[s];
339 o = o + 1;
340 }
341/*
342 Sixfold symmetry on (A,A,B,B):
343
344 123 (A,A,B)
345 132 (A,B,A),
346 134 (A,B,B)
347 312 (B,A,A)
348 314 (B,A,B)
349 341 (B,B,A)
350*/
351 else if ( suborder[s] == 6 )
352 {
353 xyz[0+o*3] = suborder_xyz[0+s*4];
354 xyz[1+o*3] = suborder_xyz[1+s*4];
355 xyz[2+o*3] = suborder_xyz[2+s*4];
356 w[o] = suborder_w[s];
357 o = o + 1;
358
359 xyz[0+o*3] = suborder_xyz[0+s*4];
360 xyz[1+o*3] = suborder_xyz[2+s*4];
361 xyz[2+o*3] = suborder_xyz[1+s*4];
362 w[o] = suborder_w[s];
363 o = o + 1;
364
365 xyz[0+o*3] = suborder_xyz[0+s*4];
366 xyz[1+o*3] = suborder_xyz[2+s*4];
367 xyz[2+o*3] = suborder_xyz[3+s*4];
368 w[o] = suborder_w[s];
369 o = o + 1;
370
371 xyz[0+o*3] = suborder_xyz[2+s*4];
372 xyz[1+o*3] = suborder_xyz[0+s*4];
373 xyz[2+o*3] = suborder_xyz[1+s*4];
374 w[o] = suborder_w[s];
375 o = o + 1;
376
377 xyz[0+o*3] = suborder_xyz[2+s*4];
378 xyz[1+o*3] = suborder_xyz[0+s*4];
379 xyz[2+o*3] = suborder_xyz[3+s*4];
380 w[o] = suborder_w[s];
381 o = o + 1;
382
383 xyz[0+o*3] = suborder_xyz[2+s*4];
384 xyz[1+o*3] = suborder_xyz[3+s*4];
385 xyz[2+o*3] = suborder_xyz[0+s*4];
386 w[o] = suborder_w[s];
387 o = o + 1;
388 }
389/*
390 Twelvefold symmetry on (A,A,B,C):
391
392 123 (A,A,B)
393 124 (A,A,C)
394 132 (A,B,A)
395 134 (A,B,C)
396 142 (A,C,A)
397 143 (A,C,B)
398 312 (B,A,A)
399 314 (B,A,C)
400 341 (B,C,A)
401 412 (C,A,A)
402 413 (C,A,B)
403 431 (C,B,A)
404*/
405 else if ( suborder[s] == 12 )
406 {
407 xyz[0+o*3] = suborder_xyz[0+s*4];
408 xyz[1+o*3] = suborder_xyz[1+s*4];
409 xyz[2+o*3] = suborder_xyz[2+s*4];
410 w[o] = suborder_w[s];
411 o = o + 1;
412
413 xyz[0+o*3] = suborder_xyz[0+s*4];
414 xyz[1+o*3] = suborder_xyz[1+s*4];
415 xyz[2+o*3] = suborder_xyz[3+s*4];
416 w[o] = suborder_w[s];
417 o = o + 1;
418
419 xyz[0+o*3] = suborder_xyz[0+s*4];
420 xyz[1+o*3] = suborder_xyz[2+s*4];
421 xyz[2+o*3] = suborder_xyz[1+s*4];
422 w[o] = suborder_w[s];
423 o = o + 1;
424
425 xyz[0+o*3] = suborder_xyz[0+s*4];
426 xyz[1+o*3] = suborder_xyz[2+s*4];
427 xyz[2+o*3] = suborder_xyz[3+s*4];
428 w[o] = suborder_w[s];
429 o = o + 1;
430
431 xyz[0+o*3] = suborder_xyz[0+s*4];
432 xyz[1+o*3] = suborder_xyz[3+s*4];
433 xyz[2+o*3] = suborder_xyz[1+s*4];
434 w[o] = suborder_w[s];
435 o = o + 1;
436
437 xyz[0+o*3] = suborder_xyz[0+s*4];
438 xyz[1+o*3] = suborder_xyz[3+s*4];
439 xyz[2+o*3] = suborder_xyz[2+s*4];
440 w[o] = suborder_w[s];
441 o = o + 1;
442
443 xyz[0+o*3] = suborder_xyz[2+s*4];
444 xyz[1+o*3] = suborder_xyz[0+s*4];
445 xyz[2+o*3] = suborder_xyz[1+s*4];
446 w[o] = suborder_w[s];
447 o = o + 1;
448
449 xyz[0+o*3] = suborder_xyz[2+s*4];
450 xyz[1+o*3] = suborder_xyz[0+s*4];
451 xyz[2+o*3] = suborder_xyz[3+s*4];
452 w[o] = suborder_w[s];
453 o = o + 1;
454
455 xyz[0+o*3] = suborder_xyz[2+s*4];
456 xyz[1+o*3] = suborder_xyz[3+s*4];
457 xyz[2+o*3] = suborder_xyz[1+s*4];
458 w[o] = suborder_w[s];
459 o = o + 1;
460
461 xyz[0+o*3] = suborder_xyz[3+s*4];
462 xyz[1+o*3] = suborder_xyz[0+s*4];
463 xyz[2+o*3] = suborder_xyz[1+s*4];
464 w[o] = suborder_w[s];
465 o = o + 1;
466
467 xyz[0+o*3] = suborder_xyz[3+s*4];
468 xyz[1+o*3] = suborder_xyz[0+s*4];
469 xyz[2+o*3] = suborder_xyz[2+s*4];
470 w[o] = suborder_w[s];
471 o = o + 1;
472
473 xyz[0+o*3] = suborder_xyz[3+s*4];
474 xyz[1+o*3] = suborder_xyz[2+s*4];
475 xyz[2+o*3] = suborder_xyz[0+s*4];
476 w[o] = suborder_w[s];
477 o = o + 1;
478 }
479/*
480 24 fold symmetry on (A,B,C,D):
481
482 123 (A,B,C)
483 124 (A,B,D)
484 132 (A,C,B)
485 134 (A,C,D)
486 142 (A,D,B)
487 143 (A,D,C)
488 213 (B,A,C)
489 214 (B,A,D)
490 231 (B,C,A)
491 234 (B,C,D)
492 241 (B,D,A)
493 243 (B,D,C)
494 312 (C,A,B)
495 314 (C,A,D)
496 321 (C,B,A)
497 324 (C,B,D)
498 341 (C,D,A)
499 342 (C,D,B)
500 412 (D,A,B)
501 413 (D,A,C)
502 421 (D,B,A)
503 423 (D,B,C)
504 431 (D,C,A)
505 432 (D,C,B)
506*/
507 else if ( suborder[s] == 24 )
508 {
509 xyz[0+o*3] = suborder_xyz[0+s*4];
510 xyz[1+o*3] = suborder_xyz[1+s*4];
511 xyz[2+o*3] = suborder_xyz[2+s*4];
512 w[o] = suborder_w[s];
513 o = o + 1;
514
515 xyz[0+o*3] = suborder_xyz[0+s*4];
516 xyz[1+o*3] = suborder_xyz[1+s*4];
517 xyz[2+o*3] = suborder_xyz[3+s*4];
518 w[o] = suborder_w[s];
519 o = o + 1;
520
521 xyz[0+o*3] = suborder_xyz[0+s*4];
522 xyz[1+o*3] = suborder_xyz[2+s*4];
523 xyz[2+o*3] = suborder_xyz[1+s*4];
524 w[o] = suborder_w[s];
525 o = o + 1;
526
527 xyz[0+o*3] = suborder_xyz[0+s*4];
528 xyz[1+o*3] = suborder_xyz[2+s*4];
529 xyz[2+o*3] = suborder_xyz[3+s*4];
530 w[o] = suborder_w[s];
531 o = o + 1;
532
533 xyz[0+o*3] = suborder_xyz[0+s*4];
534 xyz[1+o*3] = suborder_xyz[3+s*4];
535 xyz[2+o*3] = suborder_xyz[1+s*4];
536 w[o] = suborder_w[s];
537 o = o + 1;
538
539 xyz[0+o*3] = suborder_xyz[0+s*4];
540 xyz[1+o*3] = suborder_xyz[3+s*4];
541 xyz[2+o*3] = suborder_xyz[2+s*4];
542 w[o] = suborder_w[s];
543 o = o + 1;
544
545 xyz[0+o*3] = suborder_xyz[1+s*4];
546 xyz[1+o*3] = suborder_xyz[0+s*4];
547 xyz[2+o*3] = suborder_xyz[3+s*4];
548 w[o] = suborder_w[s];
549 o = o + 1;
550
551 xyz[0+o*3] = suborder_xyz[1+s*4];
552 xyz[1+o*3] = suborder_xyz[0+s*4];
553 xyz[2+o*3] = suborder_xyz[4+s*4];
554 w[o] = suborder_w[s];
555 o = o + 1;
556
557 xyz[0+o*3] = suborder_xyz[1+s*4];
558 xyz[1+o*3] = suborder_xyz[2+s*4];
559 xyz[2+o*3] = suborder_xyz[0+s*4];
560 w[o] = suborder_w[s];
561 o = o + 1;
562
563 xyz[0+o*3] = suborder_xyz[1+s*4];
564 xyz[1+o*3] = suborder_xyz[2+s*4];
565 xyz[2+o*3] = suborder_xyz[3+s*4];
566 w[o] = suborder_w[s];
567 o = o + 1;
568
569 xyz[0+o*3] = suborder_xyz[1+s*4];
570 xyz[1+o*3] = suborder_xyz[3+s*4];
571 xyz[2+o*3] = suborder_xyz[0+s*4];
572 w[o] = suborder_w[s];
573 o = o + 1;
574
575 xyz[0+o*3] = suborder_xyz[1+s*4];
576 xyz[1+o*3] = suborder_xyz[3+s*4];
577 xyz[2+o*3] = suborder_xyz[2+s*4];
578 w[o] = suborder_w[s];
579 o = o + 1;
580
581 xyz[0+o*3] = suborder_xyz[2+s*4];
582 xyz[1+o*3] = suborder_xyz[0+s*4];
583 xyz[2+o*3] = suborder_xyz[1+s*4];
584 w[o] = suborder_w[s];
585 o = o + 1;
586
587 xyz[0+o*3] = suborder_xyz[2+s*4];
588 xyz[1+o*3] = suborder_xyz[0+s*4];
589 xyz[2+o*3] = suborder_xyz[3+s*4];
590 w[o] = suborder_w[s];
591 o = o + 1;
592
593 xyz[0+o*3] = suborder_xyz[2+s*4];
594 xyz[1+o*3] = suborder_xyz[1+s*4];
595 xyz[2+o*3] = suborder_xyz[0+s*4];
596 w[o] = suborder_w[s];
597 o = o + 1;
598
599 xyz[0+o*3] = suborder_xyz[2+s*4];
600 xyz[1+o*3] = suborder_xyz[1+s*4];
601 xyz[2+o*3] = suborder_xyz[3+s*4];
602 w[o] = suborder_w[s];
603 o = o + 1;
604
605 xyz[0+o*3] = suborder_xyz[2+s*4];
606 xyz[1+o*3] = suborder_xyz[3+s*4];
607 xyz[2+o*3] = suborder_xyz[0+s*4];
608 w[o] = suborder_w[s];
609 o = o + 1;
610
611 xyz[0+o*3] = suborder_xyz[2+s*4];
612 xyz[1+o*3] = suborder_xyz[3+s*4];
613 xyz[2+o*3] = suborder_xyz[1+s*4];
614 w[o] = suborder_w[s];
615 o = o + 1;
616
617 xyz[0+o*3] = suborder_xyz[3+s*4];
618 xyz[1+o*3] = suborder_xyz[0+s*4];
619 xyz[2+o*3] = suborder_xyz[1+s*4];
620 w[o] = suborder_w[s];
621 o = o + 1;
622
623 xyz[0+o*3] = suborder_xyz[3+s*4];
624 xyz[1+o*3] = suborder_xyz[0+s*4];
625 xyz[2+o*3] = suborder_xyz[2+s*4];
626 w[o] = suborder_w[s];
627 o = o + 1;
628
629 xyz[0+o*3] = suborder_xyz[3+s*4];
630 xyz[1+o*3] = suborder_xyz[1+s*4];
631 xyz[2+o*3] = suborder_xyz[0+s*4];
632 w[o] = suborder_w[s];
633 o = o + 1;
634
635 xyz[0+o*3] = suborder_xyz[3+s*4];
636 xyz[1+o*3] = suborder_xyz[1+s*4];
637 xyz[2+o*3] = suborder_xyz[2+s*4];
638 w[o] = suborder_w[s];
639 o = o + 1;
640
641 xyz[0+o*3] = suborder_xyz[3+s*4];
642 xyz[1+o*3] = suborder_xyz[2+s*4];
643 xyz[2+o*3] = suborder_xyz[0+s*4];
644 w[o] = suborder_w[s];
645 o = o + 1;
646
647 xyz[0+o*3] = suborder_xyz[3+s*4];
648 xyz[1+o*3] = suborder_xyz[2+s*4];
649 xyz[2+o*3] = suborder_xyz[1+s*4];
650 w[o] = suborder_w[s];
651 o = o + 1;
652 }
653 else
654 {
655 fprintf ( stderr, "\n" );
656 fprintf ( stderr, "TETRAHEDRON_NCO_RULE - Fatal error!\n" );
657 fprintf ( stderr, " Illegal SUBORDER(%d) = %d\n", s, suborder[s] );
658 exit ( 1 );
659 }
660 }
661
662 free ( suborder );
663 free ( suborder_xyz );
664 free ( suborder_w );
665
666 return;
667}
sdf_wavy_2d.w
int w
Definition sdf_wavy_2d.py:6
tetrahedron_nco_subrule
void tetrahedron_nco_subrule(int rule, int suborder_num, double suborder_xyz[], double suborder_w[])
Definition tetrahedron_nco_rule.c:889

◆ tetrahedron_nco_rule_num()

int tetrahedron_nco_rule_num ( )

Definition at line 670 of file tetrahedron_nco_rule.c.

701{
702 int rule_num;
703
704 rule_num = 7;
705
706 return rule_num;
707}

◆ tetrahedron_nco_suborder()

int * tetrahedron_nco_suborder ( int  rule,
int  suborder_num 
)

Definition at line 710 of file tetrahedron_nco_rule.c.

746{
747 int *suborder;
748
749 suborder = ( int * ) malloc ( suborder_num * sizeof ( int ) );
750
751 if ( rule == 1 )
752 {
753 suborder[0] = 1;
754 }
755 else if ( rule == 2 )
756 {
757 suborder[0] = 4;
758 }
759 else if ( rule == 3 )
760 {
761 suborder[0] = 4;
762 suborder[1] = 6;
763 }
764 else if ( rule == 4 )
765 {
766 suborder[0] = 4;
767 suborder[1] = 12;
768 suborder[2] = 4;
769 }
770 else if ( rule == 5 )
771 {
772 suborder[0] = 4;
773 suborder[1] = 12;
774 suborder[2] = 6;
775 suborder[3] = 12;
776 suborder[4] = 1;
777 }
778 else if ( rule == 6 )
779 {
780 suborder[0] = 4;
781 suborder[1] = 12;
782 suborder[2] = 12;
783 suborder[3] = 12;
784 suborder[4] = 12;
785 suborder[5] = 4;
786 }
787 else if ( rule == 7 )
788 {
789 suborder[0] = 4;
790 suborder[1] = 12;
791 suborder[2] = 12;
792 suborder[3] = 12;
793 suborder[4] = 6;
794 suborder[5] = 24;
795 suborder[6] = 4;
796 suborder[7] = 4;
797 suborder[8] = 6;
798 }
799 else
800 {
801 fprintf ( stderr, "\n" );
802 fprintf ( stderr, "TETRAHEDRON_NCO_SUBORDER - Fatal error!\n" );
803 fprintf ( stderr, " Illegal RULE = %d\n", rule );
804 exit ( 1 );
805 }
806
807 return suborder;
808}

◆ tetrahedron_nco_suborder_num()

int tetrahedron_nco_suborder_num ( int  rule)

Definition at line 811 of file tetrahedron_nco_rule.c.

845{
846 int suborder_num;
847
848 if ( rule == 1 )
849 {
850 suborder_num = 1;
851 }
852 else if ( rule == 2 )
853 {
854 suborder_num = 1;
855 }
856 else if ( rule == 3 )
857 {
858 suborder_num = 2;
859 }
860 else if ( rule == 4 )
861 {
862 suborder_num = 3;
863 }
864 else if ( rule == 5 )
865 {
866 suborder_num = 5;
867 }
868 else if ( rule == 6 )
869 {
870 suborder_num = 6;
871 }
872 else if ( rule == 7 )
873 {
874 suborder_num = 9;
875 }
876 else
877 {
878 suborder_num = -1;
879 fprintf ( stderr, "\n" );
880 fprintf ( stderr, "TETRAHEDRON_NCO_SUBORDER_NUM - Fatal error!\n" );
881 fprintf ( stderr, " Illegal RULE = %d\n", rule );
882 exit ( 1 );
883 }
884
885 return suborder_num;
886}

◆ tetrahedron_nco_subrule()

void tetrahedron_nco_subrule ( int  rule,
int  suborder_num,
double  suborder_xyz[],
double  suborder_w[] 
)

Definition at line 889 of file tetrahedron_nco_rule.c.

928{
929 int i;
930 int s;
931 int suborder_w_d;
932 int *suborder_w_n;
933 int suborder_xyz_d;
934 int *suborder_xyz_n;
935
936 suborder_xyz_n = ( int * ) malloc ( 4 * suborder_num * sizeof ( int ) );
937 suborder_w_n = ( int * ) malloc ( suborder_num * sizeof ( int ) );
938
939 if ( rule == 1 )
940 {
941 tetrahedron_nco_subrule_01 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
942 suborder_w_n, &suborder_w_d );
943 }
944 else if ( rule == 2 )
945 {
946 tetrahedron_nco_subrule_02 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
947 suborder_w_n, &suborder_w_d );
948 }
949 else if ( rule == 3 )
950 {
951 tetrahedron_nco_subrule_03 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
952 suborder_w_n, &suborder_w_d );
953 }
954 else if ( rule == 4 )
955 {
956 tetrahedron_nco_subrule_04 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
957 suborder_w_n, &suborder_w_d );
958 }
959 else if ( rule == 5 )
960 {
961 tetrahedron_nco_subrule_05 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
962 suborder_w_n, &suborder_w_d );
963 }
964 else if ( rule == 6 )
965 {
966 tetrahedron_nco_subrule_06 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
967 suborder_w_n, &suborder_w_d );
968 }
969 else if ( rule == 7 )
970 {
971 tetrahedron_nco_subrule_07 ( suborder_num, suborder_xyz_n, &suborder_xyz_d,
972 suborder_w_n, &suborder_w_d );
973 }
974 else
975 {
976 fprintf ( stderr, "\n" );
977 fprintf ( stderr, "TETRAHEDRON_NCO_SUBRULE - Fatal error!\n" );
978 fprintf ( stderr, " Illegal RULE = %d\n", rule );
979 exit ( 1 );
980 }
981
982 for ( s = 0; s < suborder_num; s++ )
983 {
984 for ( i = 0; i < 4; i++ )
985 {
986 suborder_xyz[i+s*4] =
987 ( double ) ( 1 + suborder_xyz_n[i+s*4] )
988 / ( double ) ( 4 + suborder_xyz_d );
989 }
990 }
991 for ( s = 0; s < suborder_num; s++ )
992 {
993 suborder_w[s] = ( double ) suborder_w_n[s] / ( double ) suborder_w_d;
994 }
995
996 free ( suborder_w_n );
997 free ( suborder_xyz_n );
998
999 return;
1000}
tetrahedron_nco_subrule_04
void tetrahedron_nco_subrule_04(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1223
tetrahedron_nco_subrule_03
void tetrahedron_nco_subrule_03(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1149
tetrahedron_nco_subrule_05
void tetrahedron_nco_subrule_05(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1298
tetrahedron_nco_subrule_06
void tetrahedron_nco_subrule_06(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1375
tetrahedron_nco_subrule_01
void tetrahedron_nco_subrule_01(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1003
tetrahedron_nco_subrule_07
void tetrahedron_nco_subrule_07(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1453
tetrahedron_nco_subrule_02
void tetrahedron_nco_subrule_02(int suborder_num, int suborder_xyz_n[], int *suborder_xyz_d, int suborder_w_n[], int *suborder_w_d)
Definition tetrahedron_nco_rule.c:1076

◆ tetrahedron_nco_subrule_01()

void tetrahedron_nco_subrule_01 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1003 of file tetrahedron_nco_rule.c.

1047{
1048 int i;
1049 int s;
1050 int suborder_xyz_n_01[4*1] = {
1051 0, 0, 0, 0
1052 };
1053 int suborder_xyz_d_01 = 0;
1054 int suborder_w_n_01[1] = { 1 };
1055 int suborder_w_d_01 = 1;
1056
1057 for ( s = 0; s < suborder_num; s++ )
1058 {
1059 for ( i = 0; i < 4; i++ )
1060 {
1061 suborder_xyz_n[i+s*4] = suborder_xyz_n_01[i+s*4];
1062 }
1063 }
1064 *suborder_xyz_d = suborder_xyz_d_01;
1065
1066 for ( s = 0; s < suborder_num; s++ )
1067 {
1068 suborder_w_n[s] = suborder_w_n_01[s];
1069 }
1070 *suborder_w_d = suborder_w_d_01;
1071
1072 return;
1073}

◆ tetrahedron_nco_subrule_02()

void tetrahedron_nco_subrule_02 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1076 of file tetrahedron_nco_rule.c.

1120{
1121 int i;
1122 int s;
1123 int suborder_xyz_n_02[4*1] = {
1124 0, 0, 0, 1
1125 };
1126 int suborder_xyz_d_02 = 1;
1127 int suborder_w_n_02[1] = { 1 };
1128 int suborder_w_d_02 = 4;
1129
1130 for ( s = 0; s < suborder_num; s++ )
1131 {
1132 for ( i = 0; i < 4; i++ )
1133 {
1134 suborder_xyz_n[i+s*4] = suborder_xyz_n_02[i+s*4];
1135 }
1136 }
1137 *suborder_xyz_d = suborder_xyz_d_02;
1138
1139 for ( s = 0; s < suborder_num; s++ )
1140 {
1141 suborder_w_n[s] = suborder_w_n_02[s];
1142 }
1143 *suborder_w_d = suborder_w_d_02;
1144
1145 return;
1146}

◆ tetrahedron_nco_subrule_03()

void tetrahedron_nco_subrule_03 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1149 of file tetrahedron_nco_rule.c.

1193{
1194 int i;
1195 int s;
1196 int suborder_xyz_n_03[4*2] = {
1197 0, 0, 0, 2,
1198 1, 1, 0, 0
1199 };
1200 int suborder_xyz_d_03 = 2;
1201 int suborder_w_n_03[2] = { 11, -4 };
1202 int suborder_w_d_03 = 20;
1203
1204 for ( s = 0; s < suborder_num; s++ )
1205 {
1206 for ( i = 0; i < 4; i++ )
1207 {
1208 suborder_xyz_n[i+s*4] = suborder_xyz_n_03[i+s*4];
1209 }
1210 }
1211 *suborder_xyz_d = suborder_xyz_d_03;
1212
1213 for ( s = 0; s < suborder_num; s++ )
1214 {
1215 suborder_w_n[s] = suborder_w_n_03[s];
1216 }
1217 *suborder_w_d = suborder_w_d_03;
1218
1219 return;
1220}

◆ tetrahedron_nco_subrule_04()

void tetrahedron_nco_subrule_04 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1223 of file tetrahedron_nco_rule.c.

1267{
1268 int i;
1269 int s;
1270 int suborder_xyz_n_04[4*3] = {
1271 0, 0, 0, 3,
1272 0, 0, 1, 2,
1273 1, 1, 1, 0
1274 };
1275 int suborder_xyz_d_04 = 3;
1276 int suborder_w_n_04[3] = { 20, 13, -29 };
1277 int suborder_w_d_04 = 120;
1278
1279 for ( s = 0; s < suborder_num; s++ )
1280 {
1281 for ( i = 0; i < 4; i++ )
1282 {
1283 suborder_xyz_n[i+s*4] = suborder_xyz_n_04[i+s*4];
1284 }
1285 }
1286 *suborder_xyz_d = suborder_xyz_d_04;
1287
1288 for ( s = 0; s < suborder_num; s++ )
1289 {
1290 suborder_w_n[s] = suborder_w_n_04[s];
1291 }
1292 *suborder_w_d = suborder_w_d_04;
1293
1294 return;
1295}

◆ tetrahedron_nco_subrule_05()

void tetrahedron_nco_subrule_05 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1298 of file tetrahedron_nco_rule.c.

1342{
1343 int i;
1344 int s;
1345 int suborder_xyz_n_05[4*5] = {
1346 0, 0, 0, 4,
1347 0, 0, 3, 1,
1348 2, 2, 0, 0,
1349 1, 1, 0, 2,
1350 1, 1, 1, 1
1351 };
1352 int suborder_xyz_d_05 = 4;
1353 int suborder_w_n_05[5] = { 79, -68, 142, -12, 2 };
1354 int suborder_w_d_05 = 210;
1355
1356 for ( s = 0; s < suborder_num; s++ )
1357 {
1358 for ( i = 0; i < 4; i++ )
1359 {
1360 suborder_xyz_n[i+s*4] = suborder_xyz_n_05[i+s*4];
1361 }
1362 }
1363 *suborder_xyz_d = suborder_xyz_d_05;
1364
1365 for ( s = 0; s < suborder_num; s++ )
1366 {
1367 suborder_w_n[s] = suborder_w_n_05[s];
1368 }
1369 *suborder_w_d = suborder_w_d_05;
1370
1371 return;
1372}

◆ tetrahedron_nco_subrule_06()

void tetrahedron_nco_subrule_06 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1375 of file tetrahedron_nco_rule.c.

1419{
1420 int i;
1421 int s;
1422 int suborder_xyz_n_06[4*6] = {
1423 0, 0, 0, 5,
1424 0, 0, 4, 1,
1425 0, 0, 3, 2,
1426 1, 1, 0, 3,
1427 2, 2, 1, 0,
1428 1, 1, 1, 2
1429 };
1430 int suborder_xyz_d_06 = 5;
1431 int suborder_w_n_06[6] = { 277, 97, 223, -713, 505, -53 };
1432 int suborder_w_d_06 = 2240;
1433
1434 for ( s = 0; s < suborder_num; s++ )
1435 {
1436 for ( i = 0; i < 4; i++ )
1437 {
1438 suborder_xyz_n[i+s*4] = suborder_xyz_n_06[i+s*4];
1439 }
1440 }
1441 *suborder_xyz_d = suborder_xyz_d_06;
1442
1443 for ( s = 0; s < suborder_num; s++ )
1444 {
1445 suborder_w_n[s] = suborder_w_n_06[s];
1446 }
1447 *suborder_w_d = suborder_w_d_06;
1448
1449 return;
1450}

◆ tetrahedron_nco_subrule_07()

void tetrahedron_nco_subrule_07 ( int  suborder_num,
int  suborder_xyz_n[],
int *  suborder_xyz_d,
int  suborder_w_n[],
int *  suborder_w_d 
)

Definition at line 1453 of file tetrahedron_nco_rule.c.

1497{
1498 int i;
1499 int s;
1500 int suborder_xyz_n_07[4*9] = {
1501 0, 0, 0, 6,
1502 0, 0, 5, 1,
1503 0, 0, 4, 2,
1504 1, 1, 0, 4,
1505 3, 3, 0, 0,
1506 3, 2, 1, 0,
1507 1, 1, 1, 3,
1508 2, 2, 2, 0,
1509 2, 2, 1, 1
1510 };
1511 int suborder_xyz_d_07 = 6;
1512 int suborder_w_n_07[9] = { 430, -587, 1327, 187, -1298, -398, 22, 1537, -38 };
1513 int suborder_w_d_07 = 1512;
1514
1515 for ( s = 0; s < suborder_num; s++ )
1516 {
1517 for ( i = 0; i < 4; i++ )
1518 {
1519 suborder_xyz_n[i+s*4] = suborder_xyz_n_07[i+s*4];
1520 }
1521 }
1522 *suborder_xyz_d = suborder_xyz_d_07;
1523
1524 for ( s = 0; s < suborder_num; s++ )
1525 {
1526 suborder_w_n[s] = suborder_w_n_07[s];
1527 }
1528 *suborder_w_d = suborder_w_d_07;
1529
1530 return;
1531}

◆ tetrahedron_volume()

double tetrahedron_volume ( double  tetra[3 *4])

Definition at line 1534 of file tetrahedron_nco_rule.c.

1560{
1561 double a[4*4];
1562 int i;
1563 int j;
1564 double volume;
1565
1566 for ( i = 0; i < 3; i++ )
1567 {
1568 for ( j = 0; j < 4; j++ )
1569 {
1570 a[i+j*4] = tetra[i+j*3];
1571 }
1572 }
1573
1574 i = 3;
1575 for ( j = 0; j < 4; j++ )
1576 {
1577 a[i+j*4] = 1.0;
1578 }
1579
1580 volume = fabs ( r8mat_det_4d ( a ) ) / 6.0;
1581
1582 return volume;
1583}
r8mat_det_4d
double r8mat_det_4d(double a[])
Definition tetrahedron_nco_rule.c:6
Generated by Doxygen 1.9.8 and hosted at University of Glasgow