#ifndef _FERESULT_H
#define _FERESULT_H

#include "vtkDoubleArray.h"
#include "vtkPointData.h"
#include "vtkPoints.h"
#include "vtkUnstructuredGrid.h"
#include "vtkPoints.h"
#include "vtkTetra.h"
#include "vtkHexahedron.h"
#include "vtkCellArray.h"

#include "stdio.h"
#include <iostream>
#include <fstream>
#include <string> 
#include <vector>
#include <list>
#include <map>

using namespace std;

struct tetraCell {
		long id;
		long node1;
		long node2;
		long node3;
		long node4;		
};

struct resultValue {
	long elementNr;
	long value;
};

struct hexaCell {
		long id;
		long node1;
		long node2;
		long node3;
		long node4;	
		long node5;
		long node6;
		long node7;
		long node8;
};



class FEresult{
private:
	//NR typedef für die maps
	typedef std::map<long, long> lmap;

public:
	
	FEresult(std::ifstream &); //NR Referenz auf ifstream
	
	ifstream &istr; //NR Referenz auf ifstream
	vtkPoints* nodeCoordinates;	
	vtkUnstructuredGrid* uGrid;
	vtkDoubleArray* scalars;
	vtkLookupTable* lut;
	hexaCell* hexa;
	tetraCell* tetra;

	vector<hexaCell*>  hexVector;
	vector<tetraCell*> tetraVector;
	vector<int> numbers; 

	//NR Deklaration der maps
	lmap pointIndexLUT;
	lmap tetraIndexLUT;
	lmap hexaIndexLUT;

	int i, stop;
	long nodeNumber, cellNumber, numNodes,  elementType;
	long numCells;
	long numStressValues, n, j;
	float d;

	int cellNr;
	int cellContent;
	double value;

	int *cells;
	
	void readInputFile();

	void buildGeometry(int elementType);
	void buildHexa();
	void buildTetra();
	void buildLookUpTable();
	void buildScalars();
	void buildTetraScalars();
	void buildHexaScalars();

	void setNumNodes(int numberOfNodes);
	void setNumCells(int numberOfCells);
	void setElementType(int element);
	void setCellNumber(int num);

	int  getCellNumber();
	long getElementType();
	long getNumCells();
	long getNumNodes();	
};

#endif _FERESULT_H