#include <iostream>
#include "main.h"
#include "cube_array.h"
#include "heat_cube.h"

using namespace std;

void Set_pipe_temp(Cube_array& big_cube,int input_temp[]);
void List_pipe(Cube_array& big_cube);

int main()
{
    cerr << "Heat simulation" << endl;

    Cube_array big_cube;  //make a 60x60x200 cube
    big_cube.SetRatio(0,0,0);//air/air
    big_cube.SetRatio(0,1,883);//air/rock 8.83%
    big_cube.SetRatio(1,0,883);//rock/air
    big_cube.SetRatio(1,1,27);//rock/rock 0.27% down=*0,66=0.18% up=*1.33=0.36%


    for (int i=0;i<15;i++) big_cube.Element(30,30,i).SetType(air);//set the central pipe to air at 1.5m

    for (int i=0;i<20;i++) big_cube.Element(30,20,i).SetType(air);//set the no.2 pipe to air at 2m
    for (int i=0;i<20;i++) big_cube.Element(30,40,i).SetType(air);//set the no.2 pipe to air at 2m
    for (int i=0;i<20;i++) big_cube.Element(20,30,i).SetType(air);//set the no.2 pipe to air at 2m
    for (int i=0;i<20;i++) big_cube.Element(40,30,i).SetType(air);//set the no.2 pipe to air at 2m

    for (int i=0;i<25;i++) big_cube.Element(10,30,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(50,30,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(30,50,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(30,10,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(16,16,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(16,44,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(44,16,i).SetType(air);//set the no.3 pipe to air at 2.5m
    for (int i=0;i<25;i++) big_cube.Element(44,44,i).SetType(air);//set the no.3 pipe to air at 2.5m

    int input_temp[13];
    Set_pipe_temp(big_cube,input_temp);
    for (int i=0;i<13;i++){
        //cout<<"pipe no."<<i<<" input temp="<<input_temp[i]<<endl;
        input_temp[i]=0;
    }


    for (int j=1;j<61;j++){
        cout<<"simulation no."<<j<<": "<<endl;

        big_cube.Flow();//Flow the energy in the system
        big_cube.Update(); //Update the temperatures
        Set_pipe_temp(big_cube,input_temp);
        for (int i=0;i<13;i++){
            cout<<"pipe no."<<i<<" input p="<<(input_temp[i]*2000/60/j)<<endl;
        }
        cout<<"rock temperature(30,25,10)="<<big_cube.Element(30,25,10).GetTemperature()<<endl;
    }



    big_cube.Flow();//Flow the energy in the system

    big_cube.Update(); //Update the temperatures

    List_pipe(big_cube);

  //big_cube.List();//list the inner cube temperatures




    //big_cube.List();//list the inner cube temperatures
    //big_cube.List();//list the inner cube temperatures



    return 0;
}

void Set_pipe_temp(Cube_array& big_cube,int input_temp[])
{

for (int i=0;i<15;i++) input_temp[0]+=big_cube.Element(30,30,i).SetTemperature((int)(500-(i*13)));//set the no.1 pipe 500 to 305celsius

for (int i=0;i<20;i++) input_temp[1]+=big_cube.Element(30,20,i).SetTemperature((int)(290-(i*6.5)));//set the no.2 pipe 290 to 160celsius
for (int i=0;i<20;i++) input_temp[2]+=big_cube.Element(30,40,i).SetTemperature((int)(290-(i*6.5)));//set the no.2 pipe 290 to 160celsius
for (int i=0;i<20;i++) input_temp[3]+=big_cube.Element(20,30,i).SetTemperature((int)(290-(i*6.5)));//set the no.2 pipe 290 to 160celsius
for (int i=0;i<20;i++) input_temp[4]+=big_cube.Element(40,30,i).SetTemperature((int)(290-(i*6.5)));//set the no.2 pipe 290 to 160celsius

for (int i=0;i<25;i++) input_temp[5]+=big_cube.Element(10,30,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[6]+=big_cube.Element(50,30,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[7]+=big_cube.Element(30,50,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[8]+=big_cube.Element(30,10,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[9]+=big_cube.Element(16,16,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[10]+=big_cube.Element(16,44,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[11]+=big_cube.Element(44,16,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius
for (int i=0;i<25;i++) input_temp[12]+=big_cube.Element(44,44,i).SetTemperature((int)(150-(i*3.6)));//set the no.2 pipe 150 to 60celsius

}

void List_pipe(Cube_array& big_cube)
{

cout<<"no.1 pipe(central)"<<endl;
for (int i=0;i<15;i++) cout <<"@-"<<i*10<<"cm="<<big_cube.Element(30,30,i).GetTemperature()<<" C"<<endl;//get the no.1 pipe temperature
cout<<"no.2 pipe"<<endl;
for (int i=0;i<20;i++) cout <<"@-"<<i*10<<"cm="<<big_cube.Element(20,30,i).GetTemperature()<<" C"<<endl;//get the no.2 pipe temperature
cout<<"no.3 pipe"<<endl;
for (int i=0;i<25;i++) cout <<"@-"<<i*10<<"cm="<<big_cube.Element(10,30,i).GetTemperature()<<" C"<<endl;//get the no.3 pipe temperature


}
/**calculations for HeatFlowRatio

air-rock:

octave:18> deltaT=200
deltaT = 200
octave:19> d=0.125
d = 0.12500
octave:20> A=d*pi / 4 * 0.1
A = 0.0098175
octave:21> kappa = 30
kappa = 30
octave:22> Patadas=deltaT * kappa *A
Patadas = 58.905
octave:23> Patadas * 86400
ans =  5.0894e+06
octave:24> Enapi=Patadas * 86400
Enapi =  5.0894e+06
octave:25> Enapi / 2000
ans = 2544.7
octave:26> Eperc = Patadas*60
Eperc = 3534.3
octave:27> terfogatifajho=2000
terfogatifajho = 2000
octave:28> Tvalt = Eperc / terfogatifajho
Tvalt = 1.7671
octave:29>

rock-rock:

octave:1> hovezetesiegyutthato=0.9
hovezetesiegyutthato = 0.90000
octave:2> deltaT=100
deltaT = 100
octave:3> TeljesitmenyAtfolyas=deltaT * hovezetesiegyutthato * 0.1
TeljesitmenyAtfolyas = 9
octave:4> Eperc=TeljesitmenyAtfolyas*60
Eperc = 540
octave:5>  terfogatifajho=2000
terfogatifajho = 2000
octave:6> Tvaltozas=Eperc /  terfogatifajho
Tvaltozas = 0.27000
octave:7>

*/