power_rule_test

power_rule_test, a C code which calls power_rule() to create a power rule, that is, a multidimensional quadrature rule formed as a repeated product of a single one-dimensional quadrature rule.

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Related Data and Programs:

power_rule, a C code which constructs a power rule, that is, a product quadrature rule from identical 1D factor rules.

Source Code:

power_rule_test.sh, runs all the tests.
power_rule_test.txt, the output file.

As an example, we start with a 1D Clenshaw Curtis rule of order 3:

cc3_x.txt, the abscissas of a 1D Clenshaw Curtis rule of order 3.
cc3_w.txt, the weights of a 1D Clenshaw Curtis rule of order 3.
cc3_r.txt, defines the range of the region.

We give the command

power_rule cc3 2

to create a Clenshaw Curtis product rule of order 9 for a 2D region. The files which are created include:

power_x.txt, the abscissas of a 2D Clenshaw Curtis power rule of order 9.
power_w.txt, the weights of a 2D Clenshaw Curtis power rule of order 9.
power_r.txt, defines the range of the power region.

Last revised on 27 July 2019.