/*
 Egy eroforrást tartalmazó gyártócella ütemezése.
 Cél: az átlagos készletszint minimalizálása.
 Megoldás: SPT ütemezési szabály alkalmazása

*/

#include<stdio.h>
#include<stdlib.h>

//---------------------------------------------------------------------------
typedef struct {
                int id;  //azonosito
                long ProcT; //muveleti ido
                long StartT; //inditasi idopont
                long EndT;   //befejezesi idopont
               }T_JOB;

void Simulation( T_JOB* job, int NJ, int* s, long t0);
void Evalute( T_JOB* job, int NJ, double* obj_f );
void print_obj_f( double* obj_f );



int main(int argc, char* argv[])
{
  int NJ; //munkak szama, number of jobs
  T_JOB* job;  //munkak dinamikusan letrehozott tombje
  int* s;      //utemterv
  int i;       //munka indexe
  double obj_f[3];  //3 celfuggveny

  printf("\n SPT utemezesi szabaly demo");
  printf("\n Munkak szam = ");
  scanf("%d", &NJ);  //ellenorzes nincs

  //mem. foglalas
  job = (T_JOB*) calloc( NJ, sizeof( T_JOB ) );  //strukturatomb
  s = (int*) calloc( NJ, sizeof( int ) );  //int tomb

  //termelesinformatikai rendszer
  //veletlenszeru adatokkal helyettesitjuk
  for ( i=0; i<NJ; i++ )
    {
      job[i].id = i;  //decimalis azonosito
      job[i].ProcT = 1 + rand() % 100; //veletlenszam
      //printf("\n ProcT=%ld", job[i].ProcT);
      s[i] = i;  //ad-hoc sorrend
    }

  //ad-hoc sorrend
  Simulation( job, NJ, s, 0);
  Evalute( job, NJ, obj_f );
  print_obj_f( obj_f );


  /*
  rendezes
  szimulacio
  ertekeles
  eredmeny megjelenitese
  */


  //mem. felszabaditas
  free( job );
  free( s );

  getch();  //biil.
  return 0;
}

void Simulation( T_JOB* job, int NJ, int* s, long t0)
{
  //vegrahajtas szimulacioja
  //s mutatja a sorrendet
  //az eroforras egyidejuleg csak egyetlen munkan dolgozhat
  //az atallitas ideje bele van integralva a muveleti idobbe

  int pos;
  for( pos=0; pos<NJ; pos++ )
     {
       if ( pos == 0 ) //kezdo muvelet
            job[ s[pos] ].StartT = t0; //ref. idopontban indul
       else
            job[ s[pos] ].StartT = job[ s[pos-1] ].EndT;

       job[s[pos]].EndT=job[s[pos]].StartT+job[s[pos]].ProcT;

       //printf("\n %d \t %ld \t %ld \t %ld", pos, job[ s[pos] ].StartT, job[ s[pos] ].ProcT, job[ s[pos] ].EndT);
     }
}

void Evalute( T_JOB* job, int NJ, double* obj_f )
{
   double Na;    //atlagos keszletszin WIP work in process
   double Cmax;  //befejezesi idopont, completion time
   double Csum = 0;  //munkak befejezesi idopontjainak osszege, sum of job completion time

   int i;
   for ( i=0; i<NJ; i++ )
     {
       //BUG =
       Csum += job[i].EndT;

       if ( i==0 )
          Cmax = job[i].EndT;
       else
            if ( Cmax < job[i].EndT )
                 Cmax = job[i].EndT;
     }

     Na = Csum / Cmax;

   obj_f[0] = Cmax;
   obj_f[1] = Csum;
   obj_f[2] = Na;
}

void print_obj_f( double* obj_f )
{
   printf("\n Celfuggveny ertekek:");
   printf("\n Cmax = %lf ", obj_f[0]);
   printf("\n Csum = %lf ", obj_f[1]);
   printf("\n WIP = %lf ", obj_f[2]);

}