Trop_Volume {TML}

Estimate the volume of a tropical polytope

Description

This function uses tropical HAR with a uniform target distribution to estimate the volume of a tropical polytope

Usage

Trop_Volume(B, P, x0, S, I, R)

Arguments

B

matrix of points defining a minimum enclosing ball for a polytope P; rows are the points

P

matrix of points defining a tropical polytope; rows are the points

x0

initial point used for the HAR sampler

S

number of points to sample from the minimum enclosing ball

I

number of iterations for the HAR sampler

R

radius of the minimum enclosing tropical ball

Value

list containing ratio of points falling in P; volume of the tropical ball; volume estimate of P

Author(s)

David Barnhill [e-mail]

References

Barnhill, David, Ruriko Yoshida and Keiji Miura (2023). Maximum Inscribed and Minimum Enclosing Tropical Balls of Tropical Polytopes and Applications to Volume Estimation and Uniform Sampling.

Examples

P <-matrix(c(0,0,0,0,3,1,0,2,5),3,3,TRUE)
BR<-min_enc_ball(P)
B<-trop_bal.vert(BR[[1]],BR[[2]])
x0<-c(0,1.5,.4)
S<-200
I<-50
R<-BR[[2]]
Trop_Volume(B,P,x0,S,I,R)