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

feat: support rational functions in HomotopyContinuationProblem #3151

Merged

Conversation

AayushSabharwal
Copy link
Member

Requires SciML/SciMLBase.jl#829

Checklist

  • Appropriate tests were added
  • Any code changes were done in a way that does not break public API
  • All documentation related to code changes were updated
  • The new code follows the
    contributor guidelines, in particular the SciML Style Guide and
    COLPRAC.
  • Any new documentation only uses public API

Additional context

Add any other context about the problem here.

jac = true, eval_expression, eval_module)

denominator = MTK.build_explicit_observed_function(sys, denoms)
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

don't you just want the numerator?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The denominator is used to rule out roots that cause a 0/0. A trivial example would be 0 ~ (x - y) / (x - 3), 0 ~ y - 3. If we take the numerators and solve, we get x = y = 3, but that causes 0 / 0 in the first polynomial and hence the system is actually infeasible.

@AayushSabharwal AayushSabharwal force-pushed the as/homotopy-rational-poly branch from 92da576 to 77ceea5 Compare October 26, 2024 08:28
@AayushSabharwal AayushSabharwal force-pushed the as/homotopy-rational-poly branch from 77ceea5 to 7aae63d Compare October 26, 2024 15:06
@ChrisRackauckas ChrisRackauckas merged commit 1266976 into SciML:master Oct 26, 2024
35 of 39 checks passed
@AayushSabharwal AayushSabharwal deleted the as/homotopy-rational-poly branch October 26, 2024 18:06
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants