Optimizer: no way to check for infeasibility!

[gchabert]

When the optimizer returns, how the user can be informed that infeasibility occurred (no feasible point at all)?

[bneveu]

It is reported in Optimizer::report
there are 2 possibilities

1/ the problem is infeasible
it is reported " infeasible problem"
2/ no feasible point has been found, but the problem may be feasible
cout << " no feasible point found " << endl;

maybe the optimize function itself should return a status object and not void.

[gchabert]

Yes, that was my point. It would be nice to have a return code (or some other status function) to check infeasibility programatically.
Gilles

[gchabert]

OK c'est fait par contre ce serait bien que tu y jettes un coup d'oeil.
En particulier, on déclarait un problème "infeasible" dans report() si uplo_of_epsboxes =+oo et loup==+oo. Je pense que la condition (uplo_of_epsboxes =+oo) suffit et cela me permet ainsi de prendre en compte le cas où une borne initiale est donnée par l'utilisateur et que celle-ci s'avère en deça du minimum global réel (cette borne impose du coup un loup < +oo). Voir les commentaires de optimize() dans le header Optimizer.h.
Merci de valider...

[bneveu]

On 18/09/2014 10:47, Gilles Chabert wrote:

OK c'est fait par contre ce serait bien que tu y jettes un coup d'oeil.
En particulier, on déclarait un problème "infeasible" dans report() si
uplo_of_epsboxes =+oo et loup==+oo. Je pense que la condition
(uplo_of_epsboxes =+oo) suffit et cela me permet ainsi de prendre en
compte le cas où une borne initiale est donnée par l'utilisateur et
que celle-ci s'avère en deça du minimum global réel (cette borne
impose du coup un loup < +oo). Voir les commentaires de optimize()
dans le header Optimizer.h.
Merci de valider...

—
Reply to this email directly or view it on GitHub
#8 (comment).

le report() ne marche plus si tu enlèves la condition loup==+oo

il trouve infeasible tout problème qui finit bien.
si loup différent de +oo sert à une borne donnée par l'utilisateur et
n'indique plus qu'aucune solution n'a été trouvée ,
il faut un autre indicateur pour cela (en particulier savoir si
loup_point est une solution);

../build/examples/optimizer04
../benchs/benchs-optim/coconutbenchmark-library1/ex2_1_9.bch acidhc4 xn
smearsumrel 1.e-8 1.e-8 1000 1

avec trace =1
...
simplex1] loup update -0.375000007218 loup point (0 ; 0 ; 0 ;
0.249993311176 ; 0.250017558483 ; 0.249993939414 ; 0.24999520083 ; 0 ; 0
; 0)
ymax=-0.375000017218 uplo= -0.37500008873
[simplex1] loup update -0.375000007218 loup point (0 ; 0 ; 0 ;
0.249993311176 ; 0.250017558483 ; 0.249993939414 ; 0.24999520083 ; 0 ; 0
; 0)
ymax=-0.375000017218 uplo= -0.375000084636
infeasible problem
cpu time used 4.82119800001s.
number of cells 1822
Automatic cleanup of gaol... done

A +,

Bertrand.

[gchabert]

Ah oui c'est à cause du upperbounding j'imagine... on peut du coup trouver des points faisables et n'avoir aucune boîte epsilonesques -> uplo_of_epsboxes=+oo.
C'est ça?

[bneveu]

On 18/09/2014 14:55, Gilles Chabert wrote:

Ah oui c'est à cause du upperbounding j'imagine... on peut du coup
trouver des points faisables et n'avoir aucune boîte epsilonesques ->
uplo_of_epsboxes=+oo.
C'est ça?

—
Reply to this email directly or view it on GitHub
#8 (comment).

oui, c'est que la condition d'arrêt sur l'objectif s'est déclenchée : on
a trouvé un point faisable et loup-uplo < eps_obj
avant d'avoir des boîtes epsilonesques , ce qui arrive la plupart du
temps (heureusement) .

Bertrand.

[gchabert]

J'ai rajouté du coup la condition (loup==initial_loup) qui généralise simplement l'ancienne (loup==+oo).
En effet, on accepte un nouveau loup-point que s'il diminue strictement le loup (donc il ne peut pas y avoir de loup-point avec loup=initial_loup).
A priori ça devrait être bon pour toi, du coup je ferme le ticket.
Gilles