#include "cinstancereduite.h"
#include "cinstance.h"

using std::vector;
vector<double> CInstanceReduite::m_DEM ;
vector<double> CInstanceReduite::m_DEM_min ;
vector<double> CInstanceReduite::m_EDEM ;
vector<double> CInstanceReduite::m_EDEM_min ;

vector<std::vector<double> > CInstanceReduite::m_EMAX ;
vector<std::vector<double> > CInstanceReduite::m_PMAX ;
int CInstanceReduite::m_week_duration = 0 ;
// Ajout Daniel
vector<double> CInstanceReduite::m_ThProdMax ;
vector<double> CInstanceReduite::m_ThProdMin ;
vector<double> CInstanceReduite::m_ThCout ;

void CInstanceReduite::init()
{
    // supposition : D(t) est constant
    m_week_duration = T()*::D()/H();

#ifdef RED_MOYENNE
    calculeMoyenneDemandeParSemaine();
    calculeThProdParSemaine();
    calculeThCoutParSemaine();
    cacluleEMAXParSemaine();
#endif
}

void CInstanceReduite::calculeMoyenneDemandeParSemaine()
{
    static timestep_t week_length = (timestep_t)( T()/H() ) ;
    m_DEM.resize( H(), 0 ) ;
    m_DEM_min.resize( H(), INFINI);
    m_EDEM.resize( H(), 0 ) ;
    m_EDEM_min.resize( H(), INFINI );
    for (week_t w=(week_t)0; w<H(); ++w)
    {
        for (scenario_t s=(scenario_t)0; s<S(); s=(scenario_t) (s+1))
        {
            double demande_scenario, energy_scenario = 0 ;
            for (timestep_t t = (timestep_t) (w * week_length);
                 t<(timestep_t)(w*week_length+week_length);
                 t=(timestep_t) (t+1))
            {
                demande_scenario += ::DEM(t,s) ;
                energy_scenario += ::EDEM(t,s) ;
            }
            m_DEM[w]     += demande_scenario ;
            m_DEM_min[w] = std::min(m_DEM_min[w], demande_scenario);
            m_EDEM[w]     += energy_scenario ;
            m_EDEM_min[w] = std::min(m_EDEM_min[w], energy_scenario);
        }
        m_DEM[w]  /= (S()*week_length) ;
        m_EDEM[w] /= S() ;
    }
}


/**
 * calcule EMAX pour toutes les centrales de type 2 et pour toutes les semaines
 */
void CInstanceReduite::cacluleEMAXParSemaine()
{
    static timestep_t week_length = (timestep_t)( T()/H() ) ;
    vector<double> proto( H(), 0);
    m_EMAX.resize( I(), proto );
    m_PMAX.resize( I(), proto );

    for (powerplants2_t i=(powerplants2_t)0; i<I(); i=(powerplants2_t)(i+1))
    {
        for (week_t w=(week_t)0; w<H(); w=(week_t) (w+1))
        {
            for (timestep_t t = (timestep_t) (w * week_length);
                 t<(timestep_t)(w*week_length+week_length);
                 t=(timestep_t) (t+1))
            {
                EMAX(i,w) += ::EMAX(i,t);
                PMAX(i,w) += ::PMAX(i,t);
            }
            PMAX(i,w) /= week_length ;
        }
    }
}


/**
 * calcule ThProdMin et ThProdMax pour toutes les centrales de type 1 et pour toutes les semaines
 */
void CInstanceReduite::calculeThProdParSemaine()
{
    static timestep_t week_length = (timestep_t)( T()/H() ) ;
    m_ThProdMax.resize( H(), 0 );
    m_ThProdMin.resize( H(), 0 );

    for (week_t w=(week_t)0; w<H(); w=(week_t) (w+1))
    {
        double som_max = 0;
        double som_min = 0;
        for (powerplants1_t j=(powerplants1_t)0; j<J(); j=(powerplants1_t)(j+1))
        {
            for (timestep_t t = (timestep_t) (w * week_length);
                 t<(timestep_t)(w*week_length+week_length);
                 t=(timestep_t) (t+1))
            {
                for (scenario_t s = (scenario_t) 0; s < S();s=(scenario_t)(s+1))
                {
                    som_max += ::EMAX(j,t,s);
                    som_min += ::EMIN(j,t,s);
                }
            }
        }
        som_max/=S();som_min/=S();
        m_ThProdMax[w] = som_max;
        m_ThProdMin[w] = som_min;
    }

    #ifdef TEST_PMAX
    for (week_t w=(week_t)0; w<H(); w=(week_t) (w+1))
    {
        printf( "%f ",::ThProdMax(w));
    }
    printf( "\n");
    #endif

}

void CInstanceReduite::calculeThCoutParSemaine()
{
    static timestep_t week_length = (timestep_t)( T()/H() ) ;
    m_ThCout.resize( H(), 0 );

    for (week_t w=(week_t)0; w<H(); w=(week_t) (w+1))
    {
        double num = 0;
        for (powerplants1_t j=(powerplants1_t)0; j<J(); j=(powerplants1_t)(j+1))
        {
            for (timestep_t t = (timestep_t) (w * week_length);
                 t<(timestep_t)(w*week_length+week_length);
                 t=(timestep_t) (t+1))
            {
                for (scenario_t s = (scenario_t) 0; s < S();s=(scenario_t)(s+1))
                {
                    num += ::EMAX(j,t,s)*C(j,t,s);
                }
            }
        }
        num /= S()*m_ThProdMax[w];
        m_ThCout[w] = num;
        std::cout << num << " " << std::endl ;
    }//*/

}