From: Anders P. <an...@op...> - 2008-11-23 20:24:54
|
On 21 nov 2008, at 12.40, Anders Peterson wrote: > On 21 nov 2008, at 10.07, Oddvar Grønning wrote: > >> I have a dynamic optimization problem to solve that mostly is >> overdetermined, so the optimization problem has a set of solutions >> instead of one global solution. I guess this is the reason of the >> numerical difficulties of my systems. > > Yes, I think the problem is that you have several (near) redundant > constraints. Don't think so anymore... I've long felt that the QP solvers needed a pre-solver that analysed and fixed the model/problem for better numerical performance. I think that the Matlab solver has one, and that's the difference. By experimenting a bit I managed to improve the quality of the solution to your latest model/problem: {{0.0}, {0.0175}, {-0.0175}, {1.4639}, {5.0}, {4.8768}, {4.4581}, {-6.7724}, {0.2257}} This solution is VERY close to the what you got from Matlab! Until ojAlgo actually has such a pre-solver you have to think about how you construct your matrices. /Anders |