Function burke

Source
pub fn burke(psi: &Psi) -> Result<(Weights, f64)>
Expand description

Applies Burke’s Interior Point Method (IPM) to solve a convex optimization problem.

The objective function to maximize is: f(x) = Σ(log(Σ(ψ_ij * x_j))) for i = 1 to n_sub

subject to: 1. x_j ≥ 0 for all j = 1 to n_point, 2. Σ(x_j) = 1,

where ψ is an n_sub×n_point matrix with non-negative entries and x is a probability vector.

§Arguments

  • psi - A reference to a Psi structure containing the input matrix.

§Returns

On success, returns a tuple (weights, obj) where:

  • Weights contains the optimized weights (probabilities) for each support point.
  • obj is the value of the objective function at the solution.

§Errors

This function returns an error if any step in the optimization (e.g. Cholesky factorization) fails.