load_data.cpp

#include <iostream>
#include <iomanip>
#include <cmath>
#include <fstream>
#include <sstream>
#include <string>
#include <vector>
#include <unordered_map>
#include "data.h"

using namespace std;

typedef double number;

enum Section { NONE, NODE, ELEMENTS, BC };

void load_data(int& isDataLoaded, GlobalData& data, vector<Node>& nodes, vector<Element>& elements, Grid& grid, vector<int>& borders) //wczytywanie danych z pliku
{
	//Wczytywanie pliku:
	ifstream file("../../../Test1_4_4.txt");
	//ifstream file("../../../Test2_4_4_MixGrid.txt");
	//ifstream file("../../../Test3_31_31_kwadrat.txt");
	//ifstream file("../../../Test1_4_4_HG_bez_i_z_BC_wektor_PG.txt"); 
	//ifstream file("../../../Test2_4_4_MixGrid HG_bez_i_z_BC_wektor_PG.txt");

	if (!file.is_open())
	{
		cout << "Error: Could not open the file or it does not exist!" << endl;
		isDataLoaded = 1;
		return;
	}

	nodes.clear();
	elements.clear();
	borders.clear();

	string line;
	Section curr_section = NONE;

	while (getline(file, line))
	{
		if (line.empty())
		{
			continue;
		}

		//Sekcje
		if (line.find("*Node") != string::npos)
		{
			curr_section = NODE;
			continue;
		}
		else if (line.find("*Element") != string::npos)
		{
			curr_section = ELEMENTS;
			continue;
		}
		else if (line.find("*BC") != string::npos)
		{
			curr_section = BC;
			continue;
		}

		if (curr_section == NODE)
		{
			istringstream iss(line);
			int index;
			number x, y;
			char comma;
			if (iss >> index >> comma >> x >> comma >> y)
			{
				nodes.push_back({ x, y, 0.0 });
			}
		}
		else if (curr_section == ELEMENTS)
		{
			istringstream iss(line);
			int index;
			char comma;
			Element e;
			if (iss >> index >> comma >> e.ID[0] >> comma >> e.ID[1] >> comma >> e.ID[2] >> comma >> e.ID[3])
			{
				elements.push_back(e);
			}
		}
		else if (curr_section == BC)
		{
			istringstream iss(line);
			int nodeID;

			while (iss >> nodeID) //bo w lini s? boarders
			{
				borders.push_back(nodeID); //miedzy przecinkami

				if (nodeID - 1 < nodes.size())
				{
					nodes[nodeID - 1].BC = 1.0;
				}

				if (iss.peek() == ',') iss.ignore();
			}
		}
		//Pozostale parametry
		else if (curr_section == NONE)
		{
			istringstream iss(line);
			string parameterName;

			if (line.find("SimulationTime") != string::npos)
			{
				if (iss >> parameterName >> data.SimulationTime)
					continue;
			}
			if (line.find("SimulationStepTime") != string::npos)
			{
				if (iss >> parameterName >> data.TimeStep)
					continue;
			}
			if (line.find("Conductivity") != string::npos)
			{
				if (iss >> parameterName >> data.Conductivity)
					continue;
			}
			if (line.find("Alfa") != string::npos)
			{
				if (iss >> parameterName >> data.Alfa)
					continue;
			}
			if (line.find("Tot") != string::npos)
			{
				if (iss >> parameterName >> data.Tot)
					continue;
			}
			if (line.find("InitialTemp") != string::npos)
			{
				if (iss >> parameterName >> data.InitialTemp)
					continue;
			}
			if (line.find("Density") != string::npos)
			{
				if (iss >> parameterName >> data.Density)
					continue;
			}
			if (line.find("SpecificHeat") != string::npos)
			{
				if (iss >> parameterName >> data.SpecificHeat)
					continue;
			}
			if (line.find("Nodes number") != string::npos)
			{
				string parameterName2;
				if (iss >> parameterName >> parameterName2 >> data.nN)
					continue;
			}
			if (line.find("Elements number") != string::npos)
			{
				string parameterName2;
				if (iss >> parameterName >> parameterName2 >> data.nE)
					continue;
			}
		}

		//grid.nN = nodes;
		//grid.nE = elements;
	}

	file.close();

	cout << "Wczytano " << nodes.size() << " wezlow i " << elements.size() << " elementow!" << endl;

	isDataLoaded = 0;

	return;
}

void printData(int& isDataLoaded, GlobalData& data, vector<Node>& nodes, vector<Element>& elements, Grid& grid, vector<int>& borders) //wypisywanie danych
{
	if (isDataLoaded == 1)
	{
		return;
	}

	cout << fixed << setprecision(2);
	cout << "Simulation Time: " << data.SimulationTime << endl;
	cout << "Time Step: " << data.TimeStep << endl;
	cout << "Conductivity: " << data.Conductivity << endl;
	cout << "Alfa: " << data.Alfa << endl;
	cout << "Tot: " << data.Tot << endl;
	cout << "Initial Temp: " << data.InitialTemp << endl;
	cout << "Density: " << data.Density << endl;
	cout << "Specific Heat: " << data.SpecificHeat << endl;
	cout << "Number of Nodes: " << data.nN << endl;
	cout << "Number of Elements: " << data.nE << endl;
	cout << "\nNodes:" << endl;
	for (size_t i = 0; i < nodes.size(); ++i)
	{
		cout << "Node " << i + 1 << ": (" << nodes[i].x << ", " << nodes[i].y << ")" << endl;
	}
	cout << "\nElements:" << endl;
	for (size_t i = 0; i < elements.size(); ++i)
	{
		cout << "Element " << i + 1 << ": IDs(" << elements[i].ID[0] << ", " << elements[i].ID[1] << ", "
			<< elements[i].ID[2] << ", " << elements[i].ID[3] << ")" << endl;
	}
	cout << "\nBoundary Conditions (BC) Node IDs:" << endl;
	for (size_t i = 0; i < borders.size(); ++i)
	{
		cout << borders[i] << " ";
	}
	cout << endl;

	isDataLoaded = 0;
	return;
}