Contours and Quadrature Rules

Contours

We plan to support three kinds of contours: ellipse, circle, and rectangle (not support yet).

We use the composite trapezoidal rule on ellipses and circles. We can get the quadrature points by using get_quadpts. For example,

# Create a ellipse
elp = Cim.ellipse([150, 0], 148, 148)
# Get the quadrature points associated to the trapezoidal rule
pts = get_quadpts(ctr, 50)