prism_jaskowiec_rule
prism_jaskowiec_rule,
a MATLAB code which
returns symmetric quadrature rules, with exactness up to total degree 20,
over the interior of a prism with triangular base,
by Jan Jaskowiec, Natarajan Sukumar.
Licensing:
The computer code and data files made available on this
web page are distributed under
the MIT license
Languages:
prism_jaskowiec_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.
Related Data and Programs:
prism_jaskowiec_rule_test
prism_witherden_rule,
a MATLAB code which
returns a symmetric Witherden quadrature rule for a prism with
triangular base, with exactness up to total degree 10.
Reference:
-
Jan Jaskowiec, Natarajan Sukumar,
High order cubature rules for tetrahedra and pyramids,
International Journal of Numerical Methods in Engineering,
Volume 121, Number 11, pages 2418-2436, 15 June 2020.
-
Jan Jaskowiec, Natarajan Sukumar,
High order symmetric cubature rules for tetrahedra and pyramids,
International Journal of Numerical Methods in Engineering,
Volume 122, Number 1, pages 148-171, 24 August 2020.
Source Code:
-
comp_next.m,
returns the next composition of an integer.
-
monomial_value.m,
evaluates a multidimensional monomial.
-
prism_jaskowiec_rule.m,
returns a quadrature rule of given precision for the unit
prism.
-
prism_unit_monomial_integral.m,
returns the exact integral of a given monomial over the unit prism.
-
prism_unit_volume.m,
computes the volume of a unit prism.
-
rule_order.m,
returns the number of quadrature points used in a rule
of given degree.
-
rule00.m,
returns the rule of degree 0.
-
rule01.m,
returns the rule of degree 1.
-
rule02.m,
returns the rule of degree 2.
-
rule03.m,
returns the rule of degree 3.
-
rule04.m,
returns the rule of degree 4.
-
rule05.m,
returns the rule of degree 5.
-
rule06.m,
returns the rule of degree 6.
-
rule07.m,
returns the rule of degree 7.
-
rule08.m,
returns the rule of degree 8.
-
rule09.m,
returns the rule of degree 9.
-
rule10.m,
returns the rule of degree 10.
-
rule11.m,
returns the rule of degree 11.
-
rule12.m,
returns the rule of degree 12.
-
rule13.m,
returns the rule of degree 13.
-
rule14.m,
returns the rule of degree 14.
-
rule15.m,
returns the rule of degree 15.
-
rule16.m,
returns the rule of degree 16.
-
rule17.m,
returns the rule of degree 17.
-
rule18.m,
returns the rule of degree 18.
-
rule19.m,
returns the rule of degree 19.
-
rule20.m,
returns the rule of degree 20.
-
prism_jaskowiec_rule_to_matlab.m,
a MATLAB script which reads a data file describing one of
the rules (with 128 digits) and writes a function
to return points and weights.
Last revised on 02 April 2023.