Fix milp bug #2489
Fix milp bug #2489
Conversation
|
| static_cast<uint>(std::ceil(thermalClusterProduction | ||
| / currentCluster->nominalCapacityWithSpinning))); | ||
| static_cast<uint>( | ||
| std::ceil(std::round(thermalClusterProduction |
There was a problem hiding this comment.
le ceil est-il nécessaire après l'arrondi?
There was a problem hiding this comment.
Oui parce qu'en faisant le round sur des valeurs multipliées par 1000000 on obtient bien un nombre entier mais après on divise par 1000000 donc la valeur n'est pas entière. Le round est là pour arrondir à la 7ème décimale et non à un nombre entier
C'est juste un nombre arbitrairement grand. On pourrait tout à fait choisir un nombre entre 1e4 et 1e8. |
|
Ce problème d'arrondi des solutions du solveur se pose à d'autres endroits dans la code (en amont des heuristiques thermiques par exemple), à chaque fois qu'on fait un ceil sur des variables de sortie du solveur, il faut d'abord arrondir pour éviter que 1.0000000001 ne devienne 2. |
Problème constaté par Amine sur le MILP avec Xpress. L'erreur est du à un mauvais ceiling sur des valeurs égales à n+10^-10 qui deviennent n+1 avec le ceil au lieu de n.