freegs.quadrature

Contents

freegs.quadrature#

Quadrature rules for averaging over polygons

Note: integration weights set so that sum of weights is 1 and gives the average of a function over the polygon, not the integral.

License#

Copyright 2019 Ben Dudson, University of York. Email: benjamin.dudson@york.ac.uk

This file is part of FreeGS.

FreeGS is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

FreeGS is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public License along with FreeGS. If not, see <http://www.gnu.org/licenses/>.

Functions

average(func, quad)

Average func(r,z) using given quadrature

polygon_quad(polygon[, n])

Given a polygon in the form [(r1,z1), (r2,z2), (r3,z3), ...] calculates a set of quadrature points and weights, by splitting the polygon into triangles.

triangle_quad(triangle[, n])

Given a triangle in the form [(r1,z1), (r2,z2), (r3,z3)] returns a list of evaluation points [(r,z,weight),...]