Comments (6)
- we can only solve polynomial system, so the square root is not allowed
- Your system has the same zeros if you just omit the sqrt :)
from homotopycontinuation.jl.
Hi,
Thanks for the reply @saschatimme .
Regarding 1) - I thought in this issue: #498 it was stated that there was support for sqrt
on main?
No this is unfortunately not so straightforward. We internally have a fairly complex system for the evaluation of systems where the input are taylor polynomials. We need this for the predictor step to compute higher order derivatives of
x(t)
. See Chapter 13 of Evaluating Derivatives from Griewank and Walther.On the main branch we have already support for
sqrt
but not for general real powers. For this we need to first implement the taylor polynomial expansions ofexp
andlog
since you can rewritex^a
asexp(a * log(x))
and reuse theexp
andlog
primitives. Here are the recursive relationships forexp
andlog
(p. 306)If somebody is interested in giving this a shot I can also provide more guidance.
and re 2) - Yes true, although I was just providing a minimal example to show the error. My actual system is more complex..
from homotopycontinuation.jl.
Regarding 1) - I thought in this issue: #498 it was stated that there was support for sqrt on main?
Mhm maybe there is a bug then. I will look more into this over the weekend. As for whether you can use sqrt or not: Our math only works when sqrt
is used for parameters. So no variables inside sqrt
allowed.
from homotopycontinuation.jl.
Suppose I can rearrange my system so that the square roots are on one side of the equation, is it valid in general to square both sides and solve the new system using homotopy continuation methods? I really need to solve a system where the variables sit inside the square root and this is the only way I can think of, but I'm not sure if its valid...
from homotopycontinuation.jl.
I think this should be a valid way to rewrite the system.
A more general way of rewriting a system with a sqrt is the following pattern:
Starting from
f(x) + sqrt(g(x)) = 0
add one additional variable y
and write the system
f(x) + a = 0
a^2 = g(x)
Your approach as well as mine both increases the number of solution but this is unavoidable.
I wrote some time ago a paper where the true problem also had a sqrt in it's formulation. Maybe the general techniques used there can be useful for you https://arxiv.org/pdf/2009.13408.pdf
Otherwise, if you can share more about your problem (here or via email) with me or @PBrdng, then we can maybe give more concrete advice. We are both always happy to learn about applications of HC :)
from homotopycontinuation.jl.
Sure, would be great to discuss with you. I will send you an email to [email protected] ?
from homotopycontinuation.jl.
Related Issues (20)
- Monodromy, how to choose u? HOT 3
- Type inconsistency/failure to evaluate expressions HOT 1
- solve errors for real target subspace HOT 7
- :invalid_startvalue on simple system HOT 2
- Converting from Nemo.fmpq_mpoly to HomotopyContinuation.Expression HOT 3
- cache_min_length failing with empty iterator HOT 1
- Seed for mondromy not working as intended
- Symbolic linear system solving not supported by ModelKit HOT 2
- Issues dealing with BigFloat coefficients HOT 6
- find_start_pair failing on some linearly-parametrized square systems HOT 2
- `OP_POW` not defined HOT 5
- Add checkpoints for monodromy_solve
- HC fails to solve sanity check HOT 7
- Strange behaviour of numerical irreducible decomposition HOT 4
- Can't build a system with Unicode symbols in variable names HOT 2
- `Threads.resize_nthreads!` is deprecated HOT 1
- Non-real solution incorrectly certified as real
- Start System Complaints: Incorrect Rank Calculation HOT 2
- What is the best way to approximate trigonometric functions in HC HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from homotopycontinuation.jl.