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

19 January 2020 09:01:50 AM CUBE_ARBQ_RULE_TEST C version Test the CUBE_ARBQ_RULE library. TEST01 Symmetric quadrature rule for a cube. Polynomial exactness degree DEGREE = 8 Number of nodes N = 42 J W X Y Z 0 0.115589 -0.447599 -0.468591 -0.977231 1 0.153073 0.119097 -0.179652 0.981507 2 0.0674205 -0.449016 0.887443 0.972053 3 0.0677107 0.865177 0.766938 -0.904582 4 0.164431 0.650543 -0.332237 -0.918244 5 0.210364 0.258921 0.411771 -0.884195 6 0.0561431 0.965699 0.267704 0.955702 7 0.185095 0.22607 -0.680921 0.722792 8 0.0752179 -0.944716 0.342837 0.894121 9 0.09141 -0.317586 0.924106 -0.899999 10 0.0368837 -0.522911 -0.990062 0.975545 11 0.105804 0.784316 -0.801611 0.865794 12 0.162443 -0.818579 0.404996 -0.87032 13 0.199256 0.282586 0.755697 -0.545666 14 0.189884 -0.760786 -0.470757 0.825128 15 0.307838 0.586623 0.136088 -0.307695 16 0.0266738 0.954863 -0.957035 -0.886927 17 0.0278913 -0.81325 -0.999929 -0.923041 18 0.202898 -0.158502 -0.84171 0.53139 19 0.238161 0.477338 0.731255 0.809511 20 0.234729 0.128405 -0.841039 -0.743166 21 0.407501 0.129189 0.426209 0.174203 22 0.366318 -0.553219 -0.2959 0.287481 23 0.411763 -0.322837 0.316518 0.68946 24 0.101878 0.681425 0.950025 -0.435345 25 0.238176 0.755778 -0.66554 -0.415336 26 0.183519 -0.900002 -0.504205 -0.622699 27 0.493286 -0.257584 -0.044014 -0.586826 28 0.12998 0.963183 0.134652 -0.601469 29 0.391024 0.69083 -0.0884841 0.557151 30 0.487986 0.182106 -0.43688 -0.0204827 31 0.170258 -0.822192 0.827084 0.491794 32 0.23396 0.858852 0.642174 0.136923 33 0.263984 -0.897234 0.134998 0.00775611 34 0.380818 -0.535893 0.693645 -0.285754 35 0.271894 -0.552999 -0.851057 -0.233014 36 0.172555 -0.00200143 0.96316 0.266956 37 0.0307028 0.937854 0.998486 0.674101 38 0.0812185 -0.971502 -0.858832 0.406603 39 0.0425638 -0.982908 0.958597 -0.556594 40 0.105808 0.995898 -0.589014 0.251071 41 0.115889 0.595292 -0.999374 0.134929 Sum 8 Volume 8 TEST02 Get a quadrature rule for the symmetric cube. Then write it to a file. Polynomial exactness degree DEGREE = 8 Quadrature rule written to file 'cube08.txt' TEST03 Get a quadrature rule for the symmetric cube. Set up GNUPLOT graphics input. Polynomial exactness degree DEGREE = 8 Created vertex file 'cube08_vertices.txt' Created node file 'cube08_nodes.txt' Created command file 'cube08_commands.txt' TEST04 Get a quadrature rule for the symmetric cube. Test its accuracy. Polynomial exactness degree DEGREE = 8 RMS error = 2.25641e-17 CUBE_ARBQ_RULE_TEST Normal end of execution. 19 January 2020 09:01:50 AM