Skip to main content

Contour

Struct Contour 

Source
pub struct Contour<F: ComplexScalar> { /* private fields */ }
Expand description

Ordered integration contour.

A Contour stores a sequence of ContourSegments. The integrator evaluates each segment independently and sums the resulting contributions.

The current concrete contour type is designed for complex-valued contours and supports heterogeneous built-in pieces such as line segments and circular arcs.

Implementations§

Source§

impl<F> Contour<F>

Source

pub fn indent( self, pole: F::Complex, radius: F, side: IndentSide, tolerance: F, ) -> Self

Source§

impl<F: ComplexScalar> Contour<F>

Source

pub fn from_pieces(pieces: Vec<ContourSegment<F>>) -> Self

Creates a contour from explicit contour pieces.

§Panics

Panics if pieces is empty.

Source

pub fn pieces(&self) -> &[ContourSegment<F>]

Source

pub fn into_pieces(self) -> Vec<ContourSegment<F>>

Source

pub fn reverse(self) -> Self

Source

pub fn close(self) -> Self

Source

pub fn with_principal_value( self, pole: F::Complex, radius: F, side: IndentSide, tolerance: F, ) -> Self

Source

pub fn indent_many( self, singularities: impl IntoIterator<Item = (F::Complex, F, IndentSide)>, tolerance: F, ) -> Self

Source§

impl<F> Contour<F>
where F: ComplexScalar,

Source

pub fn piecewise_linear(points: Vec<F::Complex>) -> Self

Creates a piecewise-linear contour through the supplied points.

Consecutive points are joined by line segments.

§Panics

Panics if fewer than two points are supplied.

Source

pub fn upper_half_disk_centered(center: F::Complex, radius: F) -> Self
where F: FromPrimitive,

Counter-clockwise upper half-disk contour.

Path:

center - R  →  center + R
center + R  →  center - R   through the upper half-plane
Source

pub fn lower_half_disk_centered(center: F::Complex, radius: F) -> Self
where F: FromPrimitive,

Clockwise lower half-disk contour.

Path:

center - R  →  center + R
center + R  →  center - R   through the lower half-plane
Source

pub fn real_axis(radius: F) -> Self
where F: FromPrimitive,

Constructs a contour following the real axis from -radius to radius.

This is the canonical finite approximation to the real line used when numerically evaluating improper integrals over (-∞, ∞).

The contour consists of a single straight line segment.

§Orientation

The contour is traversed from left to right.

§Notes

The infinite real axis is recovered in the limit radius → ∞.

Source

pub fn real_axis_offset(radius: F, imaginary_offset: F) -> Self
where F: FromPrimitive,

Constructs a straight contour parallel to the real axis.

The contour runs from

-radius + i·offset

to

radius + i·offset.
§Orientation

The contour is traversed from left to right.

§Applications

Offset contours occur frequently in physics and applied mathematics, including:

  • causal Green’s functions (i0⁺ prescriptions),
  • Laplace and Fourier inversion,
  • contour deformation to avoid poles,
  • regularisation of principal-value integrals.
Source

pub fn upper_half_disk(radius: F) -> Self
where F: FromPrimitive,

Constructs a counter-clockwise upper half-disk.

The contour consists of

  1. a straight line along the real axis from -radius to radius,
  2. a circular arc returning through the upper half-plane.
       ●
    .-' '-.
  .'       '.
-R-----------R
§Orientation

The resulting contour is positively oriented (counter-clockwise).

§Applications

This contour is commonly used with:

  • the residue theorem,
  • Jordan’s lemma,
  • Fourier transform evaluation,
  • contour integration in wave propagation.
§Notes

The infinite upper-half-plane contour is recovered by taking radius → ∞.

Source

pub fn lower_half_disk(radius: F) -> Self
where F: FromPrimitive,

Constructs a clockwise lower half-disk.

The contour consists of

  1. a straight line along the real axis from -radius to radius,
  2. a circular arc returning through the lower half-plane.
-R-----------R
  '.       .'
    '-._.-'
§Orientation

The resulting contour is negatively oriented (clockwise).

§Applications

Useful when applying the residue theorem to integrands that decay in the lower half-plane, such as Fourier integrals with negative arguments.

§Notes

The infinite lower-half-plane contour is recovered by taking radius → ∞.

Source

pub fn upper_half_disk_offset(radius: F, imaginary_offset: F) -> Self
where F: FromPrimitive,

Constructs an upper half-disk translated vertically.

The contour is identical to upper_half_disk except that it is centred at

i · imaginary_offset.

Consequently, the straight segment lies along

Im(z) = imaginary_offset.
§Applications

Shifted contours are useful when implementing

  • i0⁺ prescriptions,
  • contour regularisation,
  • displaced Bromwich contours,
  • Green’s function calculations.
Source

pub fn lower_half_disk_offset(radius: F, imaginary_offset: F) -> Self
where F: FromPrimitive,

Constructs a lower half-disk translated vertically.

The contour is identical to lower_half_disk except that it is centred at

i · imaginary_offset.

The straight segment therefore lies on

Im(z) = imaginary_offset.

This contour is frequently used when closing contours below the real axis while avoiding nearby singularities.

Trait Implementations§

Source§

impl<F: Clone + ComplexScalar> Clone for Contour<F>

Source§

fn clone(&self) -> Contour<F>

Returns a duplicate of the value. Read more
1.0.0 (const: unstable) · Source§

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source. Read more
Source§

impl<F: Debug + ComplexScalar> Debug for Contour<F>

Source§

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter. Read more

Auto Trait Implementations§

§

impl<F> Freeze for Contour<F>

§

impl<F> RefUnwindSafe for Contour<F>

§

impl<F> Send for Contour<F>
where Vec<ContourSegment<F>>: Send,

§

impl<F> Sync for Contour<F>
where Vec<ContourSegment<F>>: Sync,

§

impl<F> Unpin for Contour<F>

§

impl<F> UnsafeUnpin for Contour<F>

§

impl<F> UnwindSafe for Contour<F>

Blanket Implementations§

Source§

impl<T> Any for T
where T: 'static + ?Sized,

Source§

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more
Source§

impl<T> Borrow<T> for T
where T: ?Sized,

Source§

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more
Source§

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source§

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more
Source§

impl<T> CloneToUninit for T
where T: Clone,

Source§

unsafe fn clone_to_uninit(&self, dest: *mut u8)

🔬This is a nightly-only experimental API. (clone_to_uninit)
Performs copy-assignment from self to dest. Read more
Source§

impl<T> From<T> for T

Source§

fn from(t: T) -> T

Returns the argument unchanged.

Source§

impl<Proc> GenerateBuilder for Proc

Source§

fn build_for<P>(self, problem: P) -> Builder<Infallible<Proc>, P, Uninitialised>
where Proc: Procedure<P>, <Proc as Procedure<P>>::State: UserState,

Source§

impl<Proc> GenerateBuilderFallible for Proc

Source§

fn build_for<P>(self, problem: P) -> Builder<Proc, P, Uninitialised>
where Proc: FallibleProcedure<P>, <Proc as FallibleProcedure<P>>::State: UserState,

Source§

impl<T> Instrument for T

Source§

fn instrument(self, span: Span) -> Instrumented<Self>

Instruments this type with the provided Span, returning an Instrumented wrapper. Read more
Source§

fn in_current_span(self) -> Instrumented<Self>

Instruments this type with the current Span, returning an Instrumented wrapper. Read more
Source§

impl<T, U> Into<U> for T
where U: From<T>,

Source§

fn into(self) -> U

Calls U::from(self).

That is, this conversion is whatever the implementation of From<T> for U chooses to do.

Source§

impl<T> Same for T

Source§

type Output = T

Should always be Self
Source§

impl<SS, SP> SupersetOf<SS> for SP
where SS: SubsetOf<SP>,

Source§

fn to_subset(&self) -> Option<SS>

The inverse inclusion map: attempts to construct self from the equivalent element of its superset. Read more
Source§

fn is_in_subset(&self) -> bool

Checks if self is actually part of its subset T (and can be converted to it).
Source§

fn to_subset_unchecked(&self) -> SS

Use with care! Same as self.to_subset but without any property checks. Always succeeds.
Source§

fn from_subset(element: &SS) -> SP

The inclusion map: converts self to the equivalent element of its superset.
Source§

impl<T> ToOwned for T
where T: Clone,

Source§

type Owned = T

The resulting type after obtaining ownership.
Source§

fn to_owned(&self) -> T

Creates owned data from borrowed data, usually by cloning. Read more
Source§

fn clone_into(&self, target: &mut T)

Uses borrowed data to replace owned data, usually by cloning. Read more
Source§

impl<T, U> TryFrom<U> for T
where U: Into<T>,

Source§

type Error = !

The type returned in the event of a conversion error.
Source§

fn try_from(value: U) -> Result<T, !>

Performs the conversion.
Source§

impl<T, U> TryInto<U> for T
where U: TryFrom<T>,

Source§

type Error = <U as TryFrom<T>>::Error

The type returned in the event of a conversion error.
Source§

fn try_into(self) -> Result<U, <U as TryFrom<T>>::Error>

Performs the conversion.
Source§

impl<T> WithSubscriber for T

Source§

fn with_subscriber<S>(self, subscriber: S) -> WithDispatch<Self>
where S: Into<Dispatch>,

Attaches the provided Subscriber to this type, returning a WithDispatch wrapper. Read more
Source§

fn with_current_subscriber(self) -> WithDispatch<Self>

Attaches the current default Subscriber to this type, returning a WithDispatch wrapper. Read more