#ifndef DA_H_INCLUDED
#define DA_H_INCLUDED

#include <iostream>
#include <vector>
#include <string>

#include <indices.h>
#include "cplan.h"
#include "dualcalc.h"


class DualAscent : public CDualCalc {
    public :

    // Matrice A de taille 260 + 51 X nombre de colonnes
        std::vector< std::vector<kWh> * > A ;
        std::vector<std::vector<kWh> > recharge; //recharge[col][k] = quantité refuelée à la k-ème campagne, sur la col-ième colonne
    kWh a(int i, int j) const { return (*A[i])[j];}

    // vecteur de cout des colonnes
    std::vector<euros> cost ;

    // fuel remaining
    std::vector<kWh> reste; // reste[col] = fuel restant en fin de scénar

    // vecteur de demande de production + affectation des colonnes
    std::vector<kWh> b ;

    // vecteur de production MAX remplie par toutes les centrales nucléaires si elles fonctionnaient toutes à PMAX + affectation colonne
    std::vector<kWh> b_EMAX ;

    // mêmes informations pour la colonne themique
    std::vector<kWh> *      therm_A ;
    std::vector<kWh>        therm_recharge ;
    euros              therm_cost ;
    kWh                therm_reste ;

    protected:
        void ajoutcolonne( std::vector<kWh> * colonne, euros cou, std::vector<CPlan::CArc>& arcs, kWh resteCol );

   public:
       DualAscent(std::vector<kWh>  dem, powerplants2_t nb_centrales);

        virtual void ajoutcolonne(CPlan* cplan);

        /// retrait d'un plan de production nucléaire
        virtual void retraitColonne(const CPlan* cplan);


        /// Fonction à appeler pour le calcul des Ui Vi
        std::vector<euros_par_kWh> calculeDuals()  const;

            /// récupération des coûts réduits pour les classes utilisant cplex)
        virtual std::vector<euros> getReducedCosts() const { return std::vector<euros>(m_plans.size());}

    private:

        mutable std::vector<euros_par_kWh>   m_u;

        mutable std::vector <euros_par_kWh>  m_lambda;

        mutable euros                   m_UB;


        std::vector<kWh> calculePoids() const;

        std::vector<euros_par_kWh> DAcalc (const std::vector<euros_par_kWh> & lambda, std::vector<euros> & L, const std::vector<kWh> & q) const;

        std::vector<kWh> sousgrad(std::vector<coeff_t> & x, const std::vector<euros_par_kWh> & u,const std::vector<euros_par_kWh> & lambda,
                                  const std::vector<euros> & L, const std::vector<kWh> & q) const;

        // Recherche de borne sup upper bound
        bool testUB () const ;

        euros recherche_UB () const ;
        euros recherche_UB2 (const std::vector<kWh> &q) const ;

        /**
         * Supprime la colonne correspondant à la centrale thermique.
         * À n'appeler qu'une unique fois, juste après la génération de colonne,
         * et juste avant la préparation des données pour CPlex
         *
         * Appelé par la classe CPlex
         */
    public:
        void vireColonneThermique() ;
        void remetColonneThermique();


        //temp : méthode de debug
        void testL(const std::vector<double>& L, const std::string& mes) const
        {
        	for ( size_t i=0 ; i < L.size() ; i++ )
        	{
        		if ( (L[i]==INFINI) || L[i]!=L[i] )
        		{
                    std::cerr<<mes<<std::endl;
                    std::cerr<<"L["<<i<<"] = "<<L[i]<<std::endl;
        			throw std::string("exit(-1)");
        		}
        	}
        }

};

#endif // DA_H_INCLUDED