/*
Obliczanie naprężeń i kierunków głównych w węzłach siatki - kostka 3D z dyskretyzacją siatką regularną, analiza skalowalności - pthreads 

Miłosz Czarnota i Bzowski Krzysztof
gr 1. IS/IMiIP
 

Ogólny sens:
Klasa punkt reprezentuje jeden węzeł. Celem jest stworzenie dużej ilości punktów (umieszczonych w tablicy) a nastepnie policzenie dla nich naprężeń głównych.
Tensor reprezentowany jest przez tablicę 3x3. Tworzymy więc tablicę wskaźników składającą się z 3 komórek, a potem w każdej tworzymy kolejną tablicę trój-elementową - dostajemy macież 2-wymiarową 3x3. Zapełniamy ją losowymi liczbami z zachowaniem jednak reguły, że musi być to tensor symetryczny (jak we wzorze z pdfa). Liczymy współczynniki (zmienne I1, I2, I3) równania III rzędu - stosując wzory z PDFa. Tworzymy równanie 3-rzędu. Za pomocą wzorów Cardana obliczamy pierwiastki tego równania, które są zarazem naprężeniami głównymi w punkcie. 

Czas operacji mierzymy za pomocą dedykowanej biblioteki stworzonej przez dr Krzysztofa Banasia. 
*/


#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <pthread.h>
#include <math.h>

//------------------------------
#include "pomiar_czasu.c"
#include "punkt.h"

unsigned int rozmiar;
int debug;
void *thread_func( void *vptr_args );
punkt *cube;
pthread_mutex_t mut = PTHREAD_MUTEX_INITIALIZER;
//------------------------------
using namespace std;

main (){

do {
cout<<"Debug Mode [0/1]: ";
cin >> debug;
}
while(debug !=0 && debug !=1);

cout<<"Ilosc wezlow: ";

cin>>rozmiar;
punkt *cube = new punkt[rozmiar];

cout<<"$$ OBLICZENIA SEKWENCYJNE $$"<<endl;
inicjuj_czas();
for(int i=0; i<rozmiar; i++){
	if(debug) cube[i].show();
	cube[i].oblicz_wspolczynniki();
	if(debug) cube[i].show_rownanie();
	cube[i].oblicz_pierwiastki();
	if(debug) cube[i].show();
	}


drukuj_czas();

delete [] cube;


cout<<"$$ OBLICZENIA ROWNOLEGLE $$"<<endl;
cube = new punkt[rozmiar];
punkt :: delta = 0;

inicjuj_czas();
pthread_t watek;

pthread_create( &watek, NULL, &thread_func, NULL );

for(int i=0; i<rozmiar/2; i++){

	if(debug) cube[i].show();
	cube[i].oblicz_wspolczynniki();
	if(debug) cube[i].show_rownanie();
	pthread_mutex_lock (&mut);
	cube[i].oblicz_pierwiastki();
	pthread_mutex_unlock (&mut);
	if(debug) cube[i].show();

	}

pthread_join(watek, NULL );

drukuj_czas();
delete [] cube;
}



void *thread_func( void *vptr_args ){
cout<<"Watek 2"<<endl;
for(int j=(rozmiar/2)+1; j<rozmiar; j++){

	if(debug) cube[j].show();
	cube[j].oblicz_wspolczynniki();
	if(debug) cube[j].show_rownanie();
	
	cube[j].oblicz_pierwiastki();
	
	if(debug) cube[j].show();
	}
}