C 代码 返回正交规则， 精确度高达15度， 在 3D 四面体的内部.rar

21 January 2020 11:07:17 AM TETRAHEDRON_ARBQ_RULE_TEST C version Test the TETRAHEDRON_ARBQ_RULE library. TEST01 Quadrature rule for a tetrahedron. Polynomial exactness degree DEGREE = 8 Number of nodes N = 44 J W X Y Z 0 0.00164537 -0.0161177 0.00079986 1.14075 1 0.0142023 -0.00714883 0.172391 0.78607 2 0.00931167 0.458992 -0.174472 -0.399143 3 0.00954964 -0.452099 -0.263294 0.353609 4 0.0166053 0.110956 -0.0629572 0.81894 5 0.0147233 0.182522 0.167838 0.476496 6 0.0169601 -0.150749 -0.0782694 0.762897 7 0.015584 -0.287732 -0.123469 -0.388006 8 0.0212175 0.368075 -0.216988 0.399755 9 0.0120674 0.478816 -0.0462944 0.114466 10 0.0256962 -0.347878 -0.129563 0.250917 11 0.0304461 -0.0354829 -0.226811 0.417462 12 0.0265737 0.259165 0.133007 -0.274778 13 0.0227767 -0.00130504 0.491941 0.303193 14 0.0303973 -0.03494 0.3002 -0.34425 15 0.0344598 -0.166046 0.140774 0.397054 16 0.0391902 0.101009 -0.225926 -0.307089 17 0.0368287 -0.358747 -0.10832 -0.19899 18 0.0237014 -0.347307 -0.37423 0.0152626 19 0.0498396 -0.0863574 -0.147806 -0.000154235 20 0.0308556 0.290117 -0.344951 0.0412103 21 0.0526017 0.104209 0.0288367 0.324891 22 0.0503024 -0.0641939 0.255813 -0.069579 23 0.0255667 -0.201699 0.440437 -0.0334099 24 0.0378055 0.184617 0.399857 -0.0130353 25 0.0488183 0.378889 -0.0995659 -0.0839032 26 0.021443 -0.481872 -0.0150975 -0.045014 27 0.0138781 0.682172 -0.395652 -0.0745152 28 0.00650296 0.397864 0.225662 -0.399394 29 0.0160044 -0.697589 -0.389965 -0.125421 30 0.0163947 0.00274296 0.8135 -0.141582 31 0.0212691 0.0164569 -0.445104 -0.185824 32 0.022155 -0.461835 -0.453565 -0.324551 33 0.0189945 0.157509 0.665015 -0.331712 34 0.01937 0.539918 -0.467462 -0.316816 35 0.0168752 -0.140391 0.701265 -0.337215 36 0.00709101 0.044534 -0.510724 -0.395001 37 0.0206706 -0.423436 0.226332 -0.333644 38 0.0129502 -0.721418 -0.272716 -0.351939 39 0.0145084 0.727812 -0.286163 -0.322617 40 0.0111868 0.490189 0.185631 -0.278004 41 0.00240124 -0.875889 -0.556637 -0.372058 42 0.00196018 0.00763238 1.05606 -0.397275 43 0.00142736 0.931743 -0.566425 -0.382107 Sum 0.942809 Volume 0.942809 TEST02 Get a quadrature rule for the tetrahedron. Then write it to a file. Polynomial exactness degree DEGREE = 8 Quadrature rule written to file 'tetrahedron08.txt' TEST03 Get a quadrature rule for the tetrahedron. Set up GNUPLOT graphics input. Polynomial exactness degree DEGREE = 8 Created vertex file 'tetrahedron08_vertices.txt' Created node file 'tetrahedron08_nodes.txt' Created command file 'tetrahedron08_commands.txt' TEST04 Get a quadrature rule for the tetrahedron. Test its accuracy. Polynomial exactness degree DEGREE = 8 RMS error = 1.74408e-17 TETRAHEDRON_ARBQ_RULE_TEST Normal end of execution. 21 January 2020 11:07:17 AM