7.1
general documentation
cs_quadrature.h
Go to the documentation of this file.
1 #ifndef __CS_QUADRATURE_H__
2 #define __CS_QUADRATURE_H__
3 
4 /*============================================================================
5  * Functions to handle quadrature rules
6  *============================================================================*/
7 
8 /*
9  This file is part of Code_Saturne, a general-purpose CFD tool.
10 
11  Copyright (C) 1998-2021 EDF S.A.
12 
13  This program is free software; you can redistribute it and/or modify it under
14  the terms of the GNU General Public License as published by the Free Software
15  Foundation; either version 2 of the License, or (at your option) any later
16  version.
17 
18  This program is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
20  FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
21  details.
22 
23  You should have received a copy of the GNU General Public License along with
24  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
25  Street, Fifth Floor, Boston, MA 02110-1301, USA.
26 */
27 
28 /*----------------------------------------------------------------------------
29  * Local headers
30  *----------------------------------------------------------------------------*/
31 
32 #include <bft_error.h>
33 
34 #include "cs_base.h"
35 #include "cs_defs.h"
36 #include "cs_flag.h"
37 #include "cs_math.h"
38 #include "cs_param_types.h"
39 
40 /*----------------------------------------------------------------------------*/
41 
42 BEGIN_C_DECLS
43 
44 /*============================================================================
45  * Macro definitions
46  *============================================================================*/
47 
48 /*============================================================================
49  * Type definitions
50  *============================================================================*/
51 
52 typedef enum {
53 
54  CS_QUADRATURE_NONE,
55  CS_QUADRATURE_BARY, /* Value at the barycenter * meas */
56  CS_QUADRATURE_BARY_SUBDIV, /* Value at the barycenter * meas on a sub-mesh */
57  CS_QUADRATURE_HIGHER, /* Unique weight but several Gauss points */
58  CS_QUADRATURE_HIGHEST, /* Specific weight for each Gauss points */
59  CS_QUADRATURE_N_TYPES
60 
61 } cs_quadrature_type_t;
62 
63 /*----------------------------------------------------------------------------*/
74 /*----------------------------------------------------------------------------*/
75 
76 typedef void
77 (cs_quadrature_edge_t) (const cs_real_3_t v1,
78  const cs_real_3_t v2,
79  double len,
80  cs_real_3_t gpts[],
81  double *weights);
82 
83 /*----------------------------------------------------------------------------*/
95 /*----------------------------------------------------------------------------*/
96 
97 typedef void
98 (cs_quadrature_tria_t) (const cs_real_3_t v1,
99  const cs_real_3_t v2,
100  const cs_real_3_t v3,
101  double area,
102  cs_real_3_t gpts[],
103  double *weights);
104 
105 /*----------------------------------------------------------------------------*/
117 /*----------------------------------------------------------------------------*/
118 
119 typedef void
120 (cs_quadrature_tet_t) (const cs_real_3_t v1,
121  const cs_real_3_t v2,
122  const cs_real_3_t v3,
123  const cs_real_3_t v4,
124  double vol,
125  cs_real_3_t gpts[],
126  double weights[]);
127 
128 /*----------------------------------------------------------------------------*/
141 /*----------------------------------------------------------------------------*/
142 
143 typedef void
144 (cs_quadrature_edge_integral_t)(double tcur,
145  const cs_real_3_t v1,
146  const cs_real_3_t v2,
147  double len,
148  cs_analytic_func_t *ana,
149  void *input,
150  double results[]);
151 
152 /*----------------------------------------------------------------------------*/
166 /*----------------------------------------------------------------------------*/
167 
168 typedef void
169 (cs_quadrature_tria_integral_t)(double tcur,
170  const cs_real_3_t v1,
171  const cs_real_3_t v2,
172  const cs_real_3_t v3,
173  double area,
174  cs_analytic_func_t *ana,
175  void *input,
176  double results[]);
177 
178 /*----------------------------------------------------------------------------*/
193 /*----------------------------------------------------------------------------*/
194 
195 typedef void
196 (cs_quadrature_tetra_integral_t)(double tcur,
197  const cs_real_3_t v1,
198  const cs_real_3_t v2,
199  const cs_real_3_t v3,
200  const cs_real_3_t v4,
201  double vol,
202  cs_analytic_func_t *ana,
203  void *input,
204  double results[]);
205 
206 
207 /*============================================================================
208  * Public function prototypes
209  *============================================================================*/
210 
211 /*----------------------------------------------------------------------------*/
215 /*----------------------------------------------------------------------------*/
216 
217 void
218 cs_quadrature_setup(void);
219 
220 /*----------------------------------------------------------------------------*/
228 /*----------------------------------------------------------------------------*/
229 
230 const char *
231 cs_quadrature_get_type_name(const cs_quadrature_type_t type);
232 
233 /*----------------------------------------------------------------------------*/
244 /*----------------------------------------------------------------------------*/
245 
246 static inline void
247 cs_quadrature_edge_1pt(const cs_real_3_t v1,
248  const cs_real_3_t v2,
249  double len,
250  cs_real_3_t gpts[],
251  double *w)
252 {
253  gpts[0][0] = 0.5*(v1[0] + v2[0]);
254  gpts[0][1] = 0.5*(v1[1] + v2[1]);
255  gpts[0][2] = 0.5*(v1[2] + v2[2]);
256  w[0] = len;
257 }
258 
259 /*----------------------------------------------------------------------------*/
270 /*----------------------------------------------------------------------------*/
271 
272 void
273 cs_quadrature_edge_2pts(const cs_real_3_t v1,
274  const cs_real_3_t v2,
275  double len,
276  cs_real_3_t gpts[],
277  double *w);
278 
279 /*----------------------------------------------------------------------------*/
290 /*----------------------------------------------------------------------------*/
291 
292 void
293 cs_quadrature_edge_3pts(const cs_real_3_t v1,
294  const cs_real_3_t v2,
295  double len,
296  cs_real_3_t gpts[],
297  double w[]);
298 
299 /*----------------------------------------------------------------------------*/
311 /*----------------------------------------------------------------------------*/
312 
313 static inline void
314 cs_quadrature_tria_1pt(const cs_real_3_t v1,
315  const cs_real_3_t v2,
316  const cs_real_3_t v3,
317  double area,
318  cs_real_3_t gpts[],
319  double *w)
320 {
321  gpts[0][0] = cs_math_1ov3 * (v1[0] + v2[0] + v3[0]);
322  gpts[0][1] = cs_math_1ov3 * (v1[1] + v2[1] + v3[1]);
323  gpts[0][2] = cs_math_1ov3 * (v1[2] + v2[2] + v3[2]);
324  w[0] = area;
325 }
326 
327 /*----------------------------------------------------------------------------*/
339 /*----------------------------------------------------------------------------*/
340 
341 void
342 cs_quadrature_tria_3pts(const cs_real_3_t v1,
343  const cs_real_3_t v2,
344  const cs_real_3_t v3,
345  double area,
346  cs_real_3_t gpts[],
347  double *w);
348 
349 /*----------------------------------------------------------------------------*/
361 /*----------------------------------------------------------------------------*/
362 
363 void
364 cs_quadrature_tria_4pts(const cs_real_3_t v1,
365  const cs_real_3_t v2,
366  const cs_real_3_t v3,
367  double area,
368  cs_real_3_t gpts[],
369  double w[]);
370 
371 /*----------------------------------------------------------------------------*/
383 /*----------------------------------------------------------------------------*/
384 
385 void
386 cs_quadrature_tria_7pts(const cs_real_3_t v1,
387  const cs_real_3_t v2,
388  const cs_real_3_t v3,
389  double area,
390  cs_real_3_t gpts[],
391  double w[]);
392 
393 /*----------------------------------------------------------------------------*/
406 /*----------------------------------------------------------------------------*/
407 
408 static inline void
409 cs_quadrature_tet_1pt(const cs_real_3_t v1,
410  const cs_real_3_t v2,
411  const cs_real_3_t v3,
412  const cs_real_3_t v4,
413  double vol,
414  cs_real_3_t gpts[],
415  double weight[])
416 {
417  gpts[0][0] = 0.25 * (v1[0] + v2[0] + v3[0] + v4[0]);
418  gpts[0][1] = 0.25 * (v1[1] + v2[1] + v3[1] + v4[1]);
419  gpts[0][2] = 0.25 * (v1[2] + v2[2] + v3[2] + v4[2]);
420  weight[0] = vol;
421 }
422 
423 /*----------------------------------------------------------------------------*/
436 /*----------------------------------------------------------------------------*/
437 
438 void
439 cs_quadrature_tet_4pts(const cs_real_3_t v1,
440  const cs_real_3_t v2,
441  const cs_real_3_t v3,
442  const cs_real_3_t v4,
443  double vol,
444  cs_real_3_t gpts[],
445  double weights[]);
446 
447 /*----------------------------------------------------------------------------*/
460 /*----------------------------------------------------------------------------*/
461 
462 void
463 cs_quadrature_tet_5pts(const cs_real_3_t v1,
464  const cs_real_3_t v2,
465  const cs_real_3_t v3,
466  const cs_real_3_t v4,
467  double vol,
468  cs_real_3_t gpts[],
469  double weights[]);
470 
471 /*----------------------------------------------------------------------------*/
484 /*----------------------------------------------------------------------------*/
485 
486 void
487 cs_quadrature_tet_15pts(const cs_real_3_t v1,
488  const cs_real_3_t v2,
489  const cs_real_3_t v3,
490  const cs_real_3_t v4,
491  double vol,
492  cs_real_3_t gpts[],
493  double weights[]);
494 
495 /*----------------------------------------------------------------------------*/
509 /*----------------------------------------------------------------------------*/
510 
511 static inline void
512 cs_quadrature_edge_1pt_scal(double tcur,
513  const cs_real_3_t v1,
514  const cs_real_3_t v2,
515  double len,
516  cs_analytic_func_t *ana,
517  void *input,
518  double results[])
519 {
520  cs_real_3_t xg;
521  double feval;
522 
523  /* Copied from cs_quadrature_1pt */
524  xg[0] = .5 * (v1[0] + v2[0]);
525  xg[1] = .5 * (v1[1] + v2[1]);
526  xg[2] = .5 * (v1[2] + v2[2]);
527 
528  /* Evaluate the function at the Gauss points */
529  ana(tcur, 1, NULL, xg, false, input, &feval);
530 
531  /* Update the result with the quadrature rule */
532  *results += len * feval;
533 }
534 
535 /*----------------------------------------------------------------------------*/
549 /*----------------------------------------------------------------------------*/
550 
551 static inline void
552 cs_quadrature_edge_2pts_scal(double tcur,
553  const cs_real_3_t v1,
554  const cs_real_3_t v2,
555  double len,
556  cs_analytic_func_t *ana,
557  void *input,
558  double results[])
559 {
560  cs_real_3_t gauss_pts[2];
561  double feval[2], weights[2];
562 
563  /* Compute Gauss points and its unique weight */
564  cs_quadrature_edge_2pts(v1, v2, len, gauss_pts, weights);
565 
566  /* Evaluate the function at the Gauss points */
567  ana(tcur, 2, NULL, (const cs_real_t *)gauss_pts, false, input, feval);
568 
569  /* Update the result with the quadrature rule */
570  *results += weights[0] * feval[0] + weights[1] * feval[1];
571 }
572 
573 /*----------------------------------------------------------------------------*/
587 /*----------------------------------------------------------------------------*/
588 
589 static inline void
590 cs_quadrature_edge_3pts_scal(double tcur,
591  const cs_real_3_t v1,
592  const cs_real_3_t v2,
593  double len,
594  cs_analytic_func_t *ana,
595  void *input,
596  double results[])
597 {
598  cs_real_3_t gauss_pts[3];
599  double feval[3], weights[3];
600 
601  /* Compute Gauss points and its weights */
602  cs_quadrature_edge_3pts(v1, v2, len, gauss_pts, weights);
603 
604  /* Evaluate the function at the Gauss points */
605  ana(tcur, 3, NULL, (const cs_real_t *)gauss_pts, false, input, feval);
606 
607  /* Update the result with the quadrature rule */
608  *results += weights[0]*feval[0] + weights[1]*feval[1] + weights[2]*feval[2];
609 }
610 
611 /*----------------------------------------------------------------------------*/
625 /*----------------------------------------------------------------------------*/
626 
627 static inline void
628 cs_quadrature_edge_1pt_vect(double tcur,
629  const cs_real_3_t v1,
630  const cs_real_3_t v2,
631  double len,
632  cs_analytic_func_t *ana,
633  void *input,
634  double results[])
635 {
636  cs_real_3_t xg;
637  double feval[3];
638 
639  /* Copied from cs_quadrature_1pt */
640  xg[0] = .5 * (v1[0] + v2[0]);
641  xg[1] = .5 * (v1[1] + v2[1]);
642  xg[2] = .5 * (v1[2] + v2[2]);
643 
644  /* Evaluate the function at the Gauss points */
645  ana(tcur, 1, NULL, xg, false, input, feval);
646 
647  /* Update the result with the quadrature rule */
648  results[0] += len * feval[0];
649  results[1] += len * feval[1];
650  results[2] += len * feval[2];
651 }
652 
653 /*----------------------------------------------------------------------------*/
667 /*----------------------------------------------------------------------------*/
668 
669 static inline void
670 cs_quadrature_edge_2pts_vect(double tcur,
671  const cs_real_3_t v1,
672  const cs_real_3_t v2,
673  double len,
674  cs_analytic_func_t *ana,
675  void *input,
676  double results[])
677 {
678  cs_real_3_t gauss_pts[2];
679  double feval[6], weights[2];
680 
681  /* Compute Gauss points and its unique weight */
682  cs_quadrature_edge_2pts(v1, v2, len, gauss_pts, weights);
683 
684  /* Evaluate the function at the Gauss points */
685  ana(tcur, 2, NULL, (const cs_real_t *)gauss_pts, false, input, feval);
686 
687  /* Update the result with the quadrature rule */
688  results[0] += weights[0] * feval[0] + weights[1] * feval[3];
689  results[1] += weights[0] * feval[1] + weights[1] * feval[4];
690  results[2] += weights[0] * feval[2] + weights[1] * feval[5];
691 }
692 
693 /*----------------------------------------------------------------------------*/
707 /*----------------------------------------------------------------------------*/
708 
709 static inline void
710 cs_quadrature_edge_3pts_vect(double tcur,
711  const cs_real_3_t v1,
712  const cs_real_3_t v2,
713  double len,
714  cs_analytic_func_t *ana,
715  void *input,
716  double results[])
717 {
718  cs_real_3_t gauss_pts[3];
719  double feval[9], weights[3];
720 
721  /* Compute Gauss points and its weights */
722  cs_quadrature_edge_3pts(v1, v2, len, gauss_pts, weights);
723 
724  /* Evaluate the function at the Gauss points */
725  ana(tcur, 3, NULL, (const cs_real_t *)gauss_pts, false, input, feval);
726 
727  /* Update the result with the quadrature rule */
728  for (int p = 0; p < 3; p++) {
729  results[0] += weights[p] * feval[3*p ];
730  results[1] += weights[p] * feval[3*p+1];
731  results[2] += weights[p] * feval[3*p+2];
732  }
733 }
734 
735 /*----------------------------------------------------------------------------*/
750 /*----------------------------------------------------------------------------*/
751 
752 static inline void
753 cs_quadrature_tria_1pt_scal(double tcur,
754  const cs_real_3_t v1,
755  const cs_real_3_t v2,
756  const cs_real_3_t v3,
757  double area,
758  cs_analytic_func_t *ana,
759  void *input,
760  double results[])
761 {
762  cs_real_3_t xg;
763  double evaluation;
764 
765  /* Copied from cs_quadrature_1pt */
766  xg[0] = cs_math_1ov3 * (v1[0] + v2[0] + v3[0]);
767  xg[1] = cs_math_1ov3 * (v1[1] + v2[1] + v3[1]);
768  xg[2] = cs_math_1ov3 * (v1[2] + v2[2] + v3[2]);
769 
770  ana(tcur, 1, NULL, xg, false, input, &evaluation);
771 
772  *results += area * evaluation;
773 }
774 
775 /*----------------------------------------------------------------------------*/
790 /*----------------------------------------------------------------------------*/
791 
792 static inline void
793 cs_quadrature_tria_3pts_scal(double tcur,
794  const cs_real_3_t v1,
795  const cs_real_3_t v2,
796  const cs_real_3_t v3,
797  double area,
798  cs_analytic_func_t *ana,
799  void *input,
800  double results[])
801 {
802  cs_real_3_t gauss_pts[3];
803  double evaluation[3], weights[3];
804 
805  /* Compute Gauss points and its unique weight */
806  cs_quadrature_tria_3pts(v1, v2, v3, area, gauss_pts, weights);
807 
808  ana(tcur, 3, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
809 
810  /* Return results */
811  *results += weights[0] * evaluation[0] + weights[1] * evaluation[1] +
812  weights[2] * evaluation[2];
813 }
814 
815 /*----------------------------------------------------------------------------*/
830 /*----------------------------------------------------------------------------*/
831 
832 static inline void
833 cs_quadrature_tria_4pts_scal(double tcur,
834  const cs_real_3_t v1,
835  const cs_real_3_t v2,
836  const cs_real_3_t v3,
837  double area,
838  cs_analytic_func_t *ana,
839  void *input,
840  double results[])
841 {
842  cs_real_3_t gauss_pts[4];
843  double evaluation[4], weights[4];
844 
845  /* Compute Gauss points and its weights */
846  cs_quadrature_tria_4pts(v1, v2, v3, area, gauss_pts, weights);
847 
848  ana(tcur, 4, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
849 
850  *results += weights[0] * evaluation[0] + weights[1] * evaluation[1] +
851  weights[2] * evaluation[2] + weights[3] * evaluation[3];
852 }
853 
854 /*----------------------------------------------------------------------------*/
869 /*----------------------------------------------------------------------------*/
870 
871 static inline void
872 cs_quadrature_tria_7pts_scal(double tcur,
873  const cs_real_3_t v1,
874  const cs_real_3_t v2,
875  const cs_real_3_t v3,
876  double area,
877  cs_analytic_func_t *ana,
878  void *input,
879  double results[])
880 {
881  cs_real_3_t gauss_pts[7];
882  double evaluation[7], weights[7];
883 
884  /* Compute Gauss points and its weights */
885  cs_quadrature_tria_7pts(v1, v2, v3, area, gauss_pts, weights);
886 
887  ana(tcur, 7, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
888 
889  *results += weights[0] * evaluation[0] + weights[1] * evaluation[1] +
890  weights[2] * evaluation[2] + weights[3] * evaluation[3] +
891  weights[4] * evaluation[4] + weights[5] * evaluation[5] +
892  weights[6] * evaluation[6] ;
893 }
894 
895 /*----------------------------------------------------------------------------*/
910 /*----------------------------------------------------------------------------*/
911 
912 static inline void
913 cs_quadrature_tria_1pt_vect(double tcur,
914  const cs_real_3_t v1,
915  const cs_real_3_t v2,
916  const cs_real_3_t v3,
917  double area,
918  cs_analytic_func_t *ana,
919  void *input,
920  double results[])
921 {
922  cs_real_3_t xg;
923  double evaluation[3];
924 
925  /* Copied from cs_quadrature_1pt */
926  xg[0] = cs_math_1ov3 * (v1[0] + v2[0] + v3[0]);
927  xg[1] = cs_math_1ov3 * (v1[1] + v2[1] + v3[1]);
928  xg[2] = cs_math_1ov3 * (v1[2] + v2[2] + v3[2]);
929 
930  ana(tcur, 1, NULL, xg, false, input, evaluation);
931 
932  results[0] += area * evaluation[0];
933  results[1] += area * evaluation[1];
934  results[2] += area * evaluation[2];
935 }
936 
937 /*----------------------------------------------------------------------------*/
952 /*----------------------------------------------------------------------------*/
953 
954 static inline void
955 cs_quadrature_tria_3pts_vect(double tcur,
956  const cs_real_3_t v1,
957  const cs_real_3_t v2,
958  const cs_real_3_t v3,
959  double area,
960  cs_analytic_func_t *ana,
961  void *input,
962  double results[])
963 {
964  cs_real_3_t gauss_pts[3];
965  double evaluation[3*3], weights[3];
966 
967  /* Compute Gauss points and its unique weight */
968  cs_quadrature_tria_3pts(v1, v2, v3, area, gauss_pts, weights);
969 
970  ana(tcur, 3, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
971 
972  for (int p = 0; p < 3; p++) {
973  results[0] += weights[p] * evaluation[3*p ];
974  results[1] += weights[p] * evaluation[3*p+1];
975  results[2] += weights[p] * evaluation[3*p+2];
976  }
977 }
978 
979 /*----------------------------------------------------------------------------*/
994 /*----------------------------------------------------------------------------*/
995 
996 static inline void
997 cs_quadrature_tria_4pts_vect(double tcur,
998  const cs_real_3_t v1,
999  const cs_real_3_t v2,
1000  const cs_real_3_t v3,
1001  double area,
1002  cs_analytic_func_t *ana,
1003  void *input,
1004  double results[])
1005 {
1006  cs_real_3_t gauss_pts[4];
1007  double evaluation[3*4], weights[4];
1008 
1009  /* Compute Gauss points and its weights */
1010  cs_quadrature_tria_4pts(v1, v2, v3, area, gauss_pts, weights);
1011 
1012  ana(tcur, 4, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1013 
1014  for (int p = 0; p < 4; p++) {
1015  results[0] += weights[p] * evaluation[3*p ];
1016  results[1] += weights[p] * evaluation[3*p+1];
1017  results[2] += weights[p] * evaluation[3*p+2];
1018  }
1019 }
1020 
1021 /*----------------------------------------------------------------------------*/
1036 /*----------------------------------------------------------------------------*/
1037 
1038 static inline void
1039 cs_quadrature_tria_7pts_vect(double tcur,
1040  const cs_real_3_t v1,
1041  const cs_real_3_t v2,
1042  const cs_real_3_t v3,
1043  double area,
1044  cs_analytic_func_t *ana,
1045  void *input,
1046  double results[])
1047 {
1048  cs_real_3_t gauss_pts[7];
1049  double evaluation[3*7], weights[7];
1050 
1051  /* Compute Gauss points and its weights */
1052  cs_quadrature_tria_7pts(v1, v2, v3, area, gauss_pts, weights);
1053 
1054  ana(tcur, 7, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1055 
1056  for (int p = 0; p < 7; p++) {
1057  results[0] += weights[p] * evaluation[3*p ];
1058  results[1] += weights[p] * evaluation[3*p+1];
1059  results[2] += weights[p] * evaluation[3*p+2];
1060  }
1061 }
1062 
1063 /*----------------------------------------------------------------------------*/
1078 /*----------------------------------------------------------------------------*/
1079 
1080 static inline void
1081 cs_quadrature_tria_1pt_tens(double tcur,
1082  const cs_real_3_t v1,
1083  const cs_real_3_t v2,
1084  const cs_real_3_t v3,
1085  double area,
1086  cs_analytic_func_t *ana,
1087  void *input,
1088  double results[])
1089 {
1090  cs_real_3_t xg;
1091  double evaluation[9];
1092 
1093  /* Copied from cs_quadrature_1pt */
1094  xg[0] = cs_math_1ov3 * (v1[0] + v2[0] + v3[0]);
1095  xg[1] = cs_math_1ov3 * (v1[1] + v2[1] + v3[1]);
1096  xg[2] = cs_math_1ov3 * (v1[2] + v2[2] + v3[2]);
1097 
1098  ana(tcur, 1, NULL, xg, false, input, evaluation);
1099 
1100  for (short int ij = 0; ij < 9; ij++)
1101  results[ij] += area * evaluation[ij];
1102 }
1103 
1104 /*----------------------------------------------------------------------------*/
1119 /*----------------------------------------------------------------------------*/
1120 
1121 static inline void
1122 cs_quadrature_tria_3pts_tens(double tcur,
1123  const cs_real_3_t v1,
1124  const cs_real_3_t v2,
1125  const cs_real_3_t v3,
1126  double area,
1127  cs_analytic_func_t *ana,
1128  void *input,
1129  double results[])
1130 {
1131  cs_real_3_t gauss_pts[3];
1132  double evaluation[9*3], weights[3];
1133 
1134  /* Compute Gauss points and its unique weight */
1135  cs_quadrature_tria_3pts(v1, v2, v3, area, gauss_pts, weights);
1136 
1137  ana(tcur, 3, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1138 
1139  for (int p = 0; p < 3; p++) {
1140  const double wp = weights[p];
1141  double *eval_p = evaluation + 9*p;
1142  for (short int ij = 0; ij < 9; ij++)
1143  results[ij] += wp * eval_p[ij];
1144  }
1145 }
1146 
1147 /*----------------------------------------------------------------------------*/
1162 /*----------------------------------------------------------------------------*/
1163 
1164 static inline void
1165 cs_quadrature_tria_4pts_tens(double tcur,
1166  const cs_real_3_t v1,
1167  const cs_real_3_t v2,
1168  const cs_real_3_t v3,
1169  double area,
1170  cs_analytic_func_t *ana,
1171  void *input,
1172  double results[])
1173 {
1174  cs_real_3_t gauss_pts[4];
1175  double evaluation[9*4], weights[4];
1176 
1177  /* Compute Gauss points and its weights */
1178  cs_quadrature_tria_4pts(v1, v2, v3, area, gauss_pts, weights);
1179 
1180  ana(tcur, 4, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1181 
1182  for (int p = 0; p < 4; p++) {
1183  const double wp = weights[p];
1184  double *eval_p = evaluation + 9*p;
1185  for (short int ij = 0; ij < 9; ij++)
1186  results[ij] += wp * eval_p[ij];
1187  }
1188 }
1189 
1190 /*----------------------------------------------------------------------------*/
1205 /*----------------------------------------------------------------------------*/
1206 
1207 static inline void
1208 cs_quadrature_tria_7pts_tens(double tcur,
1209  const cs_real_3_t v1,
1210  const cs_real_3_t v2,
1211  const cs_real_3_t v3,
1212  double area,
1213  cs_analytic_func_t *ana,
1214  void *input,
1215  double results[])
1216 {
1217  cs_real_3_t gauss_pts[7];
1218  double evaluation[9*7], weights[7];
1219 
1220  /* Compute Gauss points and its weights */
1221  cs_quadrature_tria_7pts(v1, v2, v3, area, gauss_pts, weights);
1222 
1223  ana(tcur, 7, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1224 
1225  for (int p = 0; p < 7; p++) {
1226  const double wp = weights[p];
1227  double *eval_p = evaluation + 9*p;
1228  for (short int ij = 0; ij < 9; ij++)
1229  results[ij] += wp * eval_p[ij];
1230  }
1231 }
1232 
1233 /*----------------------------------------------------------------------------*/
1249 /*----------------------------------------------------------------------------*/
1250 
1251 static inline void
1252 cs_quadrature_tet_1pt_scal(double tcur,
1253  const cs_real_3_t v1,
1254  const cs_real_3_t v2,
1255  const cs_real_3_t v3,
1256  const cs_real_3_t v4,
1257  double vol,
1258  cs_analytic_func_t *ana,
1259  void *input,
1260  double results[])
1261 {
1262  cs_real_3_t xg;
1263  double evaluation;
1264 
1265  /* Copied from cs_quadrature_tet_1pt */
1266  xg[0] = 0.25 * (v1[0] + v2[0] + v3[0] + v4[0]);
1267  xg[1] = 0.25 * (v1[1] + v2[1] + v3[1] + v4[1]);
1268  xg[2] = 0.25 * (v1[2] + v2[2] + v3[2] + v4[2]);
1269 
1270  ana(tcur, 1, NULL, xg, false, input, &evaluation);
1271 
1272  *results += vol * evaluation;
1273 }
1274 
1275 /*----------------------------------------------------------------------------*/
1291 /*----------------------------------------------------------------------------*/
1292 
1293 static inline void
1294 cs_quadrature_tet_4pts_scal(double tcur,
1295  const cs_real_3_t v1,
1296  const cs_real_3_t v2,
1297  const cs_real_3_t v3,
1298  const cs_real_3_t v4,
1299  double vol,
1300  cs_analytic_func_t *ana,
1301  void *input,
1302  double results[])
1303 {
1304  cs_real_3_t gauss_pts[4];
1305  double evaluation[4], weights[4];
1306 
1307  /* Compute Gauss points and its unique weight */
1308  cs_quadrature_tet_4pts(v1, v2, v3, v4, vol, gauss_pts, weights);
1309 
1310  ana(tcur, 4, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1311 
1312  /* Return results */
1313  *results += weights[0] * evaluation[0] + weights[1] * evaluation[1] +
1314  weights[2] * evaluation[2] + weights[3] * evaluation[3];
1315 }
1316 
1317 /*----------------------------------------------------------------------------*/
1333 /*----------------------------------------------------------------------------*/
1334 
1335 static inline void
1336 cs_quadrature_tet_5pts_scal(double tcur,
1337  const cs_real_3_t v1,
1338  const cs_real_3_t v2,
1339  const cs_real_3_t v3,
1340  const cs_real_3_t v4,
1341  double vol,
1342  cs_analytic_func_t *ana,
1343  void *input,
1344  double results[])
1345 {
1346  cs_real_3_t gauss_pts[5];
1347  double evaluation[5], weights[5];
1348 
1349  /* Compute Gauss points and its weights */
1350  cs_quadrature_tet_5pts(v1, v2, v3, v4, vol, gauss_pts, weights);
1351 
1352  ana(tcur, 5, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1353 
1354  *results += weights[0] * evaluation[0] + weights[1] * evaluation[1] +
1355  weights[2] * evaluation[2] + weights[3] * evaluation[3] +
1356  weights[4] * evaluation[4];
1357 }
1358 
1359 /*----------------------------------------------------------------------------*/
1375 /*----------------------------------------------------------------------------*/
1376 
1377 static inline void
1378 cs_quadrature_tet_1pt_vect(double tcur,
1379  const cs_real_3_t v1,
1380  const cs_real_3_t v2,
1381  const cs_real_3_t v3,
1382  const cs_real_3_t v4,
1383  double vol,
1384  cs_analytic_func_t *ana,
1385  void *input,
1386  double results[])
1387 {
1388  cs_real_3_t xg;
1389  double evaluation[3];
1390 
1391  /* Copied from cs_quadrature_tet_1pt */
1392  xg[0] = 0.25 * (v1[0] + v2[0] + v3[0] + v4[0]);
1393  xg[1] = 0.25 * (v1[1] + v2[1] + v3[1] + v4[1]);
1394  xg[2] = 0.25 * (v1[2] + v2[2] + v3[2] + v4[2]);
1395 
1396  ana(tcur, 1, NULL, xg, false, input, evaluation);
1397 
1398  results[0] += vol * evaluation[0];
1399  results[1] += vol * evaluation[1];
1400  results[2] += vol * evaluation[2];
1401 }
1402 
1403 /*----------------------------------------------------------------------------*/
1419 /*----------------------------------------------------------------------------*/
1420 
1421 static inline void
1422 cs_quadrature_tet_4pts_vect(double tcur,
1423  const cs_real_3_t v1,
1424  const cs_real_3_t v2,
1425  const cs_real_3_t v3,
1426  const cs_real_3_t v4,
1427  double vol,
1428  cs_analytic_func_t *ana,
1429  void *input,
1430  double results[])
1431 {
1432  cs_real_3_t gauss_pts[4];
1433  double evaluation[3*4], weights[4];
1434 
1435  /* Compute Gauss points and its unique weight */
1436  cs_quadrature_tet_4pts(v1, v2, v3, v4, vol, gauss_pts, weights);
1437 
1438  ana(tcur, 4, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1439 
1440  for (int p = 0; p < 4; p++) {
1441  results[0] += weights[p] * evaluation[3*p ];
1442  results[1] += weights[p] * evaluation[3*p+1];
1443  results[2] += weights[p] * evaluation[3*p+2];
1444  }
1445 }
1446 
1447 /*----------------------------------------------------------------------------*/
1463 /*----------------------------------------------------------------------------*/
1464 
1465 static inline void
1466 cs_quadrature_tet_5pts_vect(double tcur,
1467  const cs_real_3_t v1,
1468  const cs_real_3_t v2,
1469  const cs_real_3_t v3,
1470  const cs_real_3_t v4,
1471  double vol,
1472  cs_analytic_func_t *ana,
1473  void *input,
1474  double results[])
1475 {
1476  cs_real_3_t gauss_pts[5];
1477  double evaluation[3*5], weights[5];
1478 
1479  /* Compute Gauss points and its weights */
1480  cs_quadrature_tet_5pts(v1, v2, v3, v4, vol, gauss_pts, weights);
1481 
1482  ana(tcur, 5, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1483 
1484  for (int p = 0; p < 5; p++) {
1485  results[0] += weights[p] * evaluation[3*p ];
1486  results[1] += weights[p] * evaluation[3*p+1];
1487  results[2] += weights[p] * evaluation[3*p+2];
1488  }
1489 }
1490 
1491 /*----------------------------------------------------------------------------*/
1507 /*----------------------------------------------------------------------------*/
1508 
1509 static inline void
1510 cs_quadrature_tet_1pt_tens(double tcur,
1511  const cs_real_3_t v1,
1512  const cs_real_3_t v2,
1513  const cs_real_3_t v3,
1514  const cs_real_3_t v4,
1515  double vol,
1516  cs_analytic_func_t *ana,
1517  void *input,
1518  double results[])
1519 {
1520  cs_real_3_t xg;
1521  double evaluation[9];
1522 
1523  /* Copied from cs_quadrature_tet_1pt */
1524  xg[0] = 0.25 * (v1[0] + v2[0] + v3[0] + v4[0]);
1525  xg[1] = 0.25 * (v1[1] + v2[1] + v3[1] + v4[1]);
1526  xg[2] = 0.25 * (v1[2] + v2[2] + v3[2] + v4[2]);
1527 
1528  ana(tcur, 1, NULL, xg, false, input, evaluation);
1529 
1530  for (short int ij = 0; ij < 9; ij++)
1531  results[ij] += vol * evaluation[ij];
1532 }
1533 
1534 /*----------------------------------------------------------------------------*/
1550 /*----------------------------------------------------------------------------*/
1551 
1552 static inline void
1553 cs_quadrature_tet_4pts_tens(double tcur,
1554  const cs_real_3_t v1,
1555  const cs_real_3_t v2,
1556  const cs_real_3_t v3,
1557  const cs_real_3_t v4,
1558  double vol,
1559  cs_analytic_func_t *ana,
1560  void *input,
1561  double results[])
1562 {
1563  cs_real_3_t gauss_pts[4];
1564  double evaluation[9*4], weights[4];
1565 
1566  /* Compute Gauss points and its unique weight */
1567  cs_quadrature_tet_4pts(v1, v2, v3, v4, vol, gauss_pts, weights);
1568 
1569  ana(tcur, 4, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1570 
1571  for (int p = 0; p < 4; p++) {
1572  const double wp = weights[p];
1573  double *eval_p = evaluation + 9*p;
1574  for (short int ij = 0; ij < 9; ij++)
1575  results[ij] += wp * eval_p[ij];
1576  }
1577 }
1578 
1579 /*----------------------------------------------------------------------------*/
1595 /*----------------------------------------------------------------------------*/
1596 
1597 static inline void
1598 cs_quadrature_tet_5pts_tens(double tcur,
1599  const cs_real_3_t v1,
1600  const cs_real_3_t v2,
1601  const cs_real_3_t v3,
1602  const cs_real_3_t v4,
1603  double vol,
1604  cs_analytic_func_t *ana,
1605  void *input,
1606  double results[])
1607 {
1608  cs_real_3_t gauss_pts[5];
1609  double evaluation[9*5], weights[5];
1610 
1611  /* Compute Gauss points and its weights */
1612  cs_quadrature_tet_5pts(v1, v2, v3, v4, vol, gauss_pts, weights);
1613 
1614  ana(tcur, 5, NULL, (const cs_real_t *)gauss_pts, false, input, evaluation);
1615 
1616  for (int p = 0; p < 5; p++) {
1617  const double wp = weights[p];
1618  double *eval_p = evaluation + 9*p;
1619  for (short int ij = 0; ij < 9; ij++)
1620  results[ij] += wp * eval_p[ij];
1621  }
1622 }
1623 
1624 /*----------------------------------------------------------------------------*/
1635 /*----------------------------------------------------------------------------*/
1636 
1637 static inline cs_quadrature_edge_integral_t *
1638 cs_quadrature_get_edge_integral(int dim,
1639  cs_quadrature_type_t qtype)
1640 {
1641  switch (dim) {
1642 
1643  case 1: /* Scalar-valued integral */
1644 
1645  switch (qtype) {
1646 
1647  case CS_QUADRATURE_BARY:
1648  case CS_QUADRATURE_BARY_SUBDIV:
1649  return cs_quadrature_edge_1pt_scal;
1650  case CS_QUADRATURE_HIGHER:
1651  return cs_quadrature_edge_2pts_scal;
1652  case CS_QUADRATURE_HIGHEST:
1653  return cs_quadrature_edge_3pts_scal;
1654 
1655  default:
1656  bft_error(__FILE__, __LINE__, 0,
1657  " %s: Invalid quadrature type\n", __func__);
1658  }
1659  break;
1660 
1661  case 3: /* Vector-valued case */
1662 
1663  switch (qtype) {
1664 
1665  case CS_QUADRATURE_BARY:
1666  case CS_QUADRATURE_BARY_SUBDIV:
1667  return cs_quadrature_edge_1pt_vect;
1668  case CS_QUADRATURE_HIGHER:
1669  return cs_quadrature_edge_2pts_vect;
1670  case CS_QUADRATURE_HIGHEST:
1671  return cs_quadrature_edge_3pts_vect;
1672 
1673  default:
1674  bft_error(__FILE__, __LINE__, 0,
1675  " %s: Invalid quadrature type\n", __func__);
1676  }
1677  break;
1678 
1679  default:
1680  bft_error(__FILE__, __LINE__, 0,
1681  " %s: Invalid dimension value %d. Only 1 and 3 are valid.\n",
1682  __func__, dim);
1683 
1684  } /* switch on dim */
1685 
1686  return NULL; /* Should not go to this stage */
1687 }
1688 
1689 /*----------------------------------------------------------------------------*/
1699 /*----------------------------------------------------------------------------*/
1700 
1701 static inline cs_quadrature_tria_integral_t *
1702 cs_quadrature_get_tria_integral(int dim,
1703  cs_quadrature_type_t qtype)
1704 {
1705  switch (dim) {
1706 
1707  case 1: /* Scalar-valued integral */
1708 
1709  switch (qtype) {
1710 
1711  case CS_QUADRATURE_BARY:
1712  case CS_QUADRATURE_BARY_SUBDIV:
1713  return cs_quadrature_tria_1pt_scal;
1714  case CS_QUADRATURE_HIGHER:
1715  return cs_quadrature_tria_4pts_scal;
1716  case CS_QUADRATURE_HIGHEST:
1717  return cs_quadrature_tria_7pts_scal;
1718 
1719  default:
1720  bft_error(__FILE__, __LINE__, 0,
1721  " %s: Invalid quadrature type\n", __func__);
1722  }
1723  break;
1724 
1725  case 3: /* Vector-valued case */
1726 
1727  switch (qtype) {
1728 
1729  case CS_QUADRATURE_BARY:
1730  case CS_QUADRATURE_BARY_SUBDIV:
1731  return cs_quadrature_tria_1pt_vect;
1732  case CS_QUADRATURE_HIGHER:
1733  return cs_quadrature_tria_4pts_vect;
1734  case CS_QUADRATURE_HIGHEST:
1735  return cs_quadrature_tria_7pts_vect;
1736 
1737  default:
1738  bft_error(__FILE__, __LINE__, 0,
1739  " %s: Invalid quadrature type\n", __func__);
1740  }
1741  break;
1742 
1743  case 9: /* Tensor-valued case */
1744 
1745  switch (qtype) {
1746 
1747  case CS_QUADRATURE_BARY:
1748  case CS_QUADRATURE_BARY_SUBDIV:
1749  return cs_quadrature_tria_1pt_tens;
1750  case CS_QUADRATURE_HIGHER:
1751  return cs_quadrature_tria_4pts_tens;
1752  case CS_QUADRATURE_HIGHEST:
1753  return cs_quadrature_tria_7pts_tens;
1754 
1755  default:
1756  bft_error(__FILE__, __LINE__, 0,
1757  " %s: Invalid quadrature type\n", __func__);
1758  }
1759  break;
1760 
1761  default:
1762  bft_error(__FILE__, __LINE__, 0,
1763  " %s: Invalid dimension value %d. Only 1, 3 and 9 are valid.\n",
1764  __func__, dim);
1765 
1766  } /* switch on dim */
1767 
1768  return NULL; /* Should not go to this stage */
1769 }
1770 
1771 /*----------------------------------------------------------------------------*/
1781 /*----------------------------------------------------------------------------*/
1782 
1783 static inline cs_quadrature_tetra_integral_t *
1784 cs_quadrature_get_tetra_integral(int dim,
1785  cs_quadrature_type_t qtype)
1786 {
1787  switch (dim) {
1788 
1789  case 1: /* Scalar-valued case */
1790 
1791  switch (qtype) {
1792 
1793  case CS_QUADRATURE_BARY:
1794  case CS_QUADRATURE_BARY_SUBDIV:
1795  return cs_quadrature_tet_1pt_scal;
1796  case CS_QUADRATURE_HIGHER:
1797  return cs_quadrature_tet_4pts_scal;
1798  case CS_QUADRATURE_HIGHEST:
1799  return cs_quadrature_tet_5pts_scal;
1800 
1801  default:
1802  bft_error(__FILE__, __LINE__, 0,
1803  " %s: Invalid quadrature type\n", __func__);
1804  }
1805  break;
1806 
1807  case 3: /* Vector-valued case */
1808 
1809  switch (qtype) {
1810 
1811  case CS_QUADRATURE_BARY:
1812  case CS_QUADRATURE_BARY_SUBDIV:
1813  return cs_quadrature_tet_1pt_vect;
1814  case CS_QUADRATURE_HIGHER:
1815  return cs_quadrature_tet_4pts_vect;
1816  case CS_QUADRATURE_HIGHEST:
1817  return cs_quadrature_tet_5pts_vect;
1818 
1819  default:
1820  bft_error(__FILE__, __LINE__, 0,
1821  " %s: Invalid quadrature type\n", __func__);
1822  }
1823  break;
1824 
1825  case 9: /* Tensor-valued case */
1826 
1827  switch (qtype) {
1828 
1829  case CS_QUADRATURE_BARY:
1830  case CS_QUADRATURE_BARY_SUBDIV:
1831  return cs_quadrature_tet_1pt_tens;
1832  case CS_QUADRATURE_HIGHER:
1833  return cs_quadrature_tet_4pts_tens;
1834  case CS_QUADRATURE_HIGHEST:
1835  return cs_quadrature_tet_5pts_tens;
1836 
1837  default:
1838  bft_error(__FILE__, __LINE__, 0,
1839  " %s: Invalid quadrature type\n", __func__);
1840  }
1841  break;
1842 
1843  default:
1844  bft_error(__FILE__, __LINE__, 0,
1845  " %s: Invalid dimension value %d. Only 1, 3 and 9 are valid.\n",
1846  __func__, dim);
1847 
1848  } /* Switch on dim */
1849 
1850  /* Avoid no return warning */
1851  return NULL;
1852 }
1853 
1854 /*----------------------------------------------------------------------------*/
1866 /*----------------------------------------------------------------------------*/
1867 
1868 cs_eflag_t
1869 cs_quadrature_get_flag(const cs_quadrature_type_t qtype,
1870  const cs_flag_t loc);
1871 
1872 /*----------------------------------------------------------------------------*/
1873 
1874 END_C_DECLS
1875 
1876 #endif /* __CS_QUADRATURE_H__ */
cs_quadrature_tria_integral_t
void() cs_quadrature_tria_integral_t(double tcur, const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, double area, cs_analytic_func_t *ana, void *input, double results[])
Compute the integral over a triangle based on a specified quadrature rule and add it to results...
Definition: cs_quadrature.h:169
cs_quadrature_tet_4pts
void cs_quadrature_tet_4pts(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, const cs_real_3_t v4, double vol, cs_real_3_t gpts[], double weights[])
Compute the quadrature in a tetrehedra. Exact for 2nd order polynomials (order 3).
Definition: cs_quadrature.c:341
cs_quadrature_tet_5pts
void cs_quadrature_tet_5pts(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, const cs_real_3_t v4, double vol, cs_real_3_t gpts[], double weights[])
Compute the quadrature in a tetrehedra. Exact for 3rd order polynomials (order 4).
Definition: cs_quadrature.c:381
cs_base.h
cs_quadrature_tetra_integral_t
void() cs_quadrature_tetra_integral_t(double tcur, const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, const cs_real_3_t v4, double vol, cs_analytic_func_t *ana, void *input, double results[])
Compute the integral over a tetrahedron based on a specified quadrature rule and add it to results...
Definition: cs_quadrature.h:196
cs_param_types.h
cs_defs.h
BEGIN_C_DECLS
#define BEGIN_C_DECLS
Definition: cs_defs.h:510
CS_QUADRATURE_BARY_SUBDIV
Definition: cs_quadrature.h:56
p
Definition: cs_field_pointer.h:67
cs_quadrature_edge_3pts
void cs_quadrature_edge_3pts(const cs_real_3_t v1, const cs_real_3_t v2, double len, cs_real_3_t gpts[], double w[])
Compute quadrature points for an edge from v1 -> v2 (3 points) Exact for polynomial function up to or...
Definition: cs_quadrature.c:198
cs_quadrature_edge_2pts
void cs_quadrature_edge_2pts(const cs_real_3_t v1, const cs_real_3_t v2, double len, cs_real_3_t gpts[], double *w)
Compute quadrature points for an edge from v1 -> v2 (2 points) Exact for polynomial function up to or...
CS_QUADRATURE_N_TYPES
Definition: cs_quadrature.h:59
bft_error
void bft_error(const char *const file_name, const int line_num, const int sys_error_code, const char *const format,...)
Calls the error handler (set by bft_error_handler_set() or default).
Definition: bft_error.c:193
cs_quadrature_get_flag
cs_eflag_t cs_quadrature_get_flag(const cs_quadrature_type_t qtype, const cs_flag_t loc)
Get the flags adapted to the given quadrature type qtype and the location on which the quadrature sho...
Definition: cs_quadrature.c:488
cs_real_t
double cs_real_t
Floating-point value.
Definition: cs_defs.h:322
CS_QUADRATURE_HIGHEST
Definition: cs_quadrature.h:58
cs_quadrature_tria_3pts
void cs_quadrature_tria_3pts(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, double area, cs_real_3_t gpts[], double *w)
Compute quadrature points for a triangle (3 points) Exact for polynomial function up to order 2...
cs_quadrature_get_type_name
const char * cs_quadrature_get_type_name(const cs_quadrature_type_t type)
Return th name associated to a type of quadrature.
Definition: cs_quadrature.c:149
cs_math_1ov3
const cs_real_t cs_math_1ov3
cs_quadrature_type_t
cs_quadrature_type_t
Definition: cs_quadrature.h:52
cs_quadrature_tria_4pts
void cs_quadrature_tria_4pts(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, double area, cs_real_3_t gpts[], double w[])
Compute quadrature points for a triangle (4 points) Exact for polynomial function up to order 3...
Definition: cs_quadrature.c:265
cs_analytic_func_t
void() cs_analytic_func_t(cs_real_t time, cs_lnum_t n_elts, const cs_lnum_t *elt_ids, const cs_real_t *coords, bool dense_output, void *input, cs_real_t *retval)
Generic function pointer for an evaluation relying on an analytic function.
Definition: cs_param_types.h:100
cs_flag.h
bft_error.h
cs_quadrature_setup
void cs_quadrature_setup(void)
Compute constant weights for all quadratures used.
Definition: cs_quadrature.c:107
cs_quadrature_edge_integral_t
void() cs_quadrature_edge_integral_t(double tcur, const cs_real_3_t v1, const cs_real_3_t v2, double len, cs_analytic_func_t *ana, void *input, double results[])
Compute the integral over an edge based on a specified quadrature rule and add it to results...
Definition: cs_quadrature.h:144
cs_quadrature_edge_t
void() cs_quadrature_edge_t(const cs_real_3_t v1, const cs_real_3_t v2, double len, cs_real_3_t gpts[], double *weights)
Generic function pointer to compute the quadrature points for an edge from v1 -> v2.
Definition: cs_quadrature.h:77
cs_real_3_t
cs_real_t cs_real_3_t[3]
vector of 3 floating-point values
Definition: cs_defs.h:335
cs_quadrature_tet_t
void() cs_quadrature_tet_t(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, const cs_real_3_t v4, double vol, cs_real_3_t gpts[], double weights[])
Generic function to compute the quadrature points in a tetrehedra.
Definition: cs_quadrature.h:120
cs_quadrature_tria_7pts
void cs_quadrature_tria_7pts(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, double area, cs_real_3_t gpts[], double w[])
Compute quadrature points for a triangle (7 points) Exact for polynomial function up to order 5...
Definition: cs_quadrature.c:300
cs_math.h
END_C_DECLS
#define END_C_DECLS
Definition: cs_defs.h:511
cs_flag_t
unsigned short int cs_flag_t
Definition: cs_defs.h:324
cs_quadrature_tet_15pts
void cs_quadrature_tet_15pts(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, const cs_real_3_t v4, double vol, cs_real_3_t gpts[], double weights[])
Compute the quadrature in a tetrehedra. Exact for 5th order polynomials (order 6).
Definition: cs_quadrature.c:425
cs_eflag_t
unsigned int cs_eflag_t
Definition: cs_flag.h:170
cs_quadrature_tria_t
void() cs_quadrature_tria_t(const cs_real_3_t v1, const cs_real_3_t v2, const cs_real_3_t v3, double area, cs_real_3_t gpts[], double *weights)
Generic functoin pointer to compute the quadrature points for a triangle.
Definition: cs_quadrature.h:98
CS_QUADRATURE_HIGHER
Definition: cs_quadrature.h:57
CS_QUADRATURE_NONE
Definition: cs_quadrature.h:54
CS_QUADRATURE_BARY
Definition: cs_quadrature.h:55