Imported Upstream version 0.0.34

1 files changed, 88 insertions, 0 deletions

diff --git a/CSXCAD/src/CSPropDiscMaterial.h b/CSXCAD/src/CSPropDiscMaterial.h
new file mode 100644
index 0000000..5317ac7
--- /dev/null
+++ b/CSXCAD/src/CSPropDiscMaterial.h
@@ -0,0 +1,88 @@
+/*
+*	Copyright (C) 2008-2012 Thorsten Liebig (Thorsten.Liebig@gmx.de)
+*
+*	This program is free software: you can redistribute it and/or modify
+*	it under the terms of the GNU Lesser General Public License as published
+*	by the Free Software Foundation, either version 3 of the License, or
+*	(at your option) any later version.
+*
+*	This program is distributed in the hope that it will be useful,
+*	but WITHOUT ANY WARRANTY; without even the implied warranty of
+*	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+*	GNU Lesser General Public License for more details.
+*
+*	You should have received a copy of the GNU Lesser General Public License
+*	along with this program.  If not, see <http://www.gnu.org/licenses/>.
+*/
+
+#pragma once
+
+#include "CSProperties.h"
+#include "CSPropMaterial.h"
+
+typedef unsigned char uint8;
+
+class vtkPolyData;
+
+//! Continuous Structure Discrete Material Property
+/*!
+  This Property reads a discrete material distribution from a file. (currently only HDF5)
+  */
+class CSXCAD_EXPORT CSPropDiscMaterial : public CSPropMaterial
+{
+public:
+	CSPropDiscMaterial(ParameterSet* paraSet);
+	CSPropDiscMaterial(CSProperties* prop);
+	CSPropDiscMaterial(unsigned int ID, ParameterSet* paraSet);
+	virtual ~CSPropDiscMaterial();
+
+	virtual const std::string GetTypeXMLString() const {return std::string("Discrete-Material");}
+
+	virtual double GetEpsilonWeighted(int ny, const double* coords);
+	virtual double GetMueWeighted(int ny, const double* coords);
+	virtual double GetKappaWeighted(int ny, const double* coords);
+	virtual double GetSigmaWeighted(int ny, const double* coords);
+
+	virtual double GetDensityWeighted(const double* coords);
+
+	//! Set true if database index 0 is used as background material (default), or false if CSPropMaterial should be used as index 0
+	virtual void SetUseDataBaseForBackground(bool val) {m_DB_Background=val;}
+
+	CSTransform* GetTransform() {return m_Transform;}
+
+	double GetScale() {return m_Scale;}
+
+	virtual void Init();
+
+	virtual bool Write2XML(TiXmlNode& root, bool parameterised=true, bool sparse=false);
+	virtual bool ReadFromXML(TiXmlNode &root);
+
+	bool ReadHDF5(std::string filename);
+
+	virtual void ShowPropertyStatus(std::ostream& stream);
+
+	//! Create a vtkPolyData surface that separates the discrete material from background material
+	virtual vtkPolyData* CreatePolyDataModel() const;
+
+protected:
+	unsigned int GetWeightingPos(const double* coords);
+	int GetDBPos(const double* coords);
+
+	int m_FileType;
+	std::string m_Filename;
+	unsigned int m_Size[3];
+	unsigned int m_DB_size;
+	uint8* m_Disc_Ind;
+	float *m_mesh[3];
+	float *m_Disc_epsR;
+	float *m_Disc_kappa;
+	float *m_Disc_mueR;
+	float *m_Disc_sigma;
+	float *m_Disc_Density;
+	double m_Scale;
+	bool m_DB_Background;
+	CSTransform* m_Transform;
+
+	void* ReadDataSet(std::string filename, std::string d_name, int type_id, int &rank, unsigned int &size, bool debug=false);
+};
+