hexahedron_witherden_rule, a Fortran90 code which returns a Witherden quadrature rule, with exactness up to total degree 11, over the interior of a hexahedron.
The unit hexahedron has vertices (0,0,0), (1,0,0), (1,1,0), (0,1,0), (0,0,1), (1,0,1), (1,1,1), (0,1,1).
The computer code and data files made available on this web page are distributed under the MIT license
hexahedron_witherden_rule is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.
hexahedron_witherden_rule_test