Update Engine prototype

Glen Whitney 2024-02-08 22:09:15 +00:00
parent e854cc7d16
commit c49cc775b4

@ -31,7 +31,7 @@ Some of these are obvious, but we may as well write them down. Note we may want
## Overall solver approaches
* Find the Gröbner basis of the constraints. This definitely tells us a lot about the possible configurations. For example, the basis will be 1 if and only if there are no solutions (correct?). And we should be able to tell if the system is rigid, i.e., has only finitely many solutions (correct?). In this case, can we find all of the algebraic number solutions? And hopefully in the general case, we still stand to find some rational points. We're of course also fine with finding say quadratic or other low-degree algebraic numbers that satisfy the constraints. (I'm no algebraic geometer; can a 1-dimensional variety have no rational points? I assume for any fixed degree, it can have only finitely many algebraic points of that degree over Q?
* Find the Gröbner basis of the constraints. This definitely tells us a lot about the possible configurations. For example, the basis will be 1 if and only if there are no solutions (correct?). And we should be able to tell if the system is rigid, i.e., has only finitely many solutions (correct?). In this case, can we find all of the algebraic number solutions? And hopefully in the general case, we still stand to find some rational points. We're of course also fine with finding say quadratic or other low-degree algebraic numbers that satisfy the constraints. (I'm no algebraic geometer; can a 1-dimensional variety have no rational points? I assume for any fixed degree, it can have only finitely many algebraic points of that degree over Q?)
* Use Newton's method to minimize the sum of the squares of the constraints. If we do find an approximate zero this way, is there a way to find a nearby algebraic zero? Sadly, if we don't find a zero this way, we may have just gotten stuck at a local minimum, so that's no proof the constraints are inconsistent.
* Use homotopy continuation (see https://www.juliahomotopycontinuation.org/) to get numerical information about the solution variety (and generate numerical approximations to points on it, correct?) If I understand correctly, an advantage of this method is that we get an almost-sure proof that the system is inconsistent (if it is), correct? And become almost sure whether it is rigid (has only isolated solutions), correct? As with Newton's method, we still have the question of whether having approximate numerical solutions helps us find exact algebraic ones. Also, if I understand correctly, this approach comes with a natural way to update solutions as one element is dragged, is that right?
* Since generally speaking each constraint is quadratic, there _might_ be some useful information here: https://math.stackexchange.com/questions/2119007/how-do-you-solve-a-system-of-quadratic-equations