Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Implement boundary oracle for integral collocation method #121

Open
papachristoumarios opened this issue Oct 21, 2020 · 2 comments
Open

Comments

@papachristoumarios
Copy link
Collaborator

papachristoumarios commented Oct 21, 2020

The Integral Collocation method of this paper is implemented in volesti under
include/ode_solvers/integral_collocation.hpp.

The current state of the solver does not include a boundary oracle (and thus it does not currently support truncation).

The goal of this feature request is to implement a boundary oracle for the method.
More specifically the boundary oracle method has to solve the intersection of the trajectory
(computed via Lagrange interpolation at the Chebyshev nodes) which is defined by its
Chebyshev transform (via using the Chebyshev transform class provided in boost) with
the boundary of the H-polytope.

To solve it effectively, I propose converting the Chebyshev transform to a complex polynomial
of twice the degree such that the real part of its roots represent solutions of the equation.

So, all in all, we need to

  1. Compute the polynomial of twice the degree (already implemented)
  2. Compute the integral of this polynomial (a[k] x^k -> a[k] / (k + 1) x^{k + 1})
  3. Project on each of the facet normals and solve the corresponding polynomial equation numerically (e.g. with MPSolve)
  4. Keep the smallest positive of the real parts of the complex roots
@JRS296
Copy link

JRS296 commented Sep 30, 2023

@vissarion I'd like to know the status of this issue, has it been solved yet? and if not, if I can help.

@vissarion
Copy link
Member

Hi @JRS296 thanks for your interest I do not know if this issue is the right to get started.
I propose to choose one issue from https://github.com/GeomScale/volesti/issues?q=is%3Aissue+is%3Aopen+label%3A%22good+first+issue%22

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants