convex.bounds: Function to find the boundaries of a convex set

In HI: Simulation from Distributions Supported by Nested Hyperplanes

Function to find the boundaries of a convex set

Description

Finds the boundaries of a bounded convex set along a specified straight line, using a bisection approach. It is mainly intended for use within arms.

Usage

convex.bounds(x, dir, indFunc, ..., tol=1e-07)

Arguments

x	A point within the set
dir	A vector specifying a direction
indFunc	Indicator function of the set
...	Parameters passed to indFunc
tol	Tolerance

Details

Uses a bisection algorithm along a line having parametric representation x + t * dir.

Value

A two dimensional vector ans. The boundaries of the set are x + ans[1] * dir and x + ans[2] * dir.

Author(s)

Giovanni Petris [e-mail]

Examples

## boundaries of a unit circle
convex.bounds(c(0,0), c(1,1), indFunc=function(x) crossprod(x)<1)

HI documentation built on April 30, 2022, 1:06 a.m.