#ifndef _TEMPARRAY_H_
#define _TEMPARRAY_H_

#include "Cubehole.h"
#include "gl.h"
#include "Triangle.h"
#include <list>
#include <math.h>

class Temparray
{
  public:
    Temparray(float initialtemp, int x0, int y0, int z0);
    Temparray(){}
    ~Temparray(){}
    void calcTemp();
    float& temperatureold(int x, int y, int z, int line, int pos) {
      return temparrayold[x*sy*sz*6*4 + y*sz*6*4 + z*6*4 + line*4 + pos];
    }
    float& temperaturenew(int x, int y, int z, int line, int pos) {
      return temparraynew[x*sy*sz*6*4 + y*sz*6*4 + z*6*4 + line*4 + pos];
    }
    std::list<Triangle> getTriangles();
    void coloring();
    
    void setprobetemp(float t){probetemp = t;}
    void setinput(float input0){input = input0;}

  private:
    float *temparrayold, *temparraynew;
    Cubehole *cubearray;
    int sx, sy, sz;
    float averagetemp, earthtemp, probetemp, input;

    Cubehole& cubehole(int x, int y, int z, int line) {
      return cubearray[x*sy*sz*6 + y*sz*6 + z*6 + line];
    }
    void mergetemperature(){
      for(int i=0; i<sx; i++){
        for(int j=0; j<sy; j++){
          for(int k=0; k<sz; k++){
            for(int l=0; l<6; l++){
              for(int m=0; m<4; m++){temperatureold(i, j, k, l, m)=temperaturenew(i, j, k, l, m);}}}}}}
    
    void tempInit(float temp0, int x0, int y0, int z0){
      for(int i=0; i<x0; i++){
        for(int j=0; j<y0; j++){
          for(int k=0; k<z0; k++){
            for(int l=0; l<6; l++){
              for(int m=0; m<4; m++){
                temperaturenew(i, j, k, l, m) = temp0;
                temperatureold(i, j, k, l, m) = temp0;
              }
            }
          }
        }
      }
      earthtemp=8;
    }
};

#endif /* _TEMPARRAY_H_ */