From f1e9a090389a41fe09340284fe09515122e56f45 Mon Sep 17 00:00:00 2001
From: michael <michael@localhost.localdomain>
Date: Fri, 5 Nov 2021 23:40:52 +0100
Subject: [PATCH] Added function getElementComponent in class Field

---
 CDMATH/mesh/inc/Field.hxx |  2 ++
 CDMATH/mesh/src/Field.cxx | 49 +++++++++++++++++++++++++++++++++++++++
 2 files changed, 51 insertions(+)

diff --git a/CDMATH/mesh/inc/Field.hxx b/CDMATH/mesh/inc/Field.hxx
index bc4ff85..542a688 100755
--- a/CDMATH/mesh/inc/Field.hxx
+++ b/CDMATH/mesh/inc/Field.hxx
@@ -239,6 +239,8 @@ class Field
 
     EntityType getTypeOfField ( void ) const ;
 
+	// returns the x, y or z component of the element (node cell or face) with number i
+    double getElementComponent(int i, int comp) const;
     /**
      * return the MEDCouplingField pointer
      * return _field
diff --git a/CDMATH/mesh/src/Field.cxx b/CDMATH/mesh/src/Field.cxx
index b7a4e94..6887ea0 100755
--- a/CDMATH/mesh/src/Field.cxx
+++ b/CDMATH/mesh/src/Field.cxx
@@ -741,6 +741,55 @@ Field::getTypeOfField ( void ) const
 	return _typeField;
 }
 
+double 
+Field::getElementComponent(int i, int comp) const
+{
+	switch( _typeField )
+	{
+		case CELLS:
+			switch( comp )
+			{
+				case 0:
+					return _mesh.getCell(i).x();
+				case 1:
+					return _mesh.getCell(i).y();
+				case 2:
+					return _mesh.getCell(i).z();
+				default:
+					cout<<"Wrong component number "<< comp <<" , dimension is "<< _mesh.getSpaceDimension() << ", field values are on CELLS" <<endl;
+					throw CdmathException("Field::getElementComponent : Wrong component number");
+			}
+		case NODES:
+			switch( comp )
+			{
+				case 0:
+					return _mesh.getNode(i).x();
+				case 1:
+					return _mesh.getNode(i).y();
+				case 2:
+					return _mesh.getNode(i).z();
+				default:
+					cout<<"Wrong component number "<< comp <<" , dimension is "<< _mesh.getSpaceDimension() << ", field values are on NODES" <<endl;
+					throw CdmathException("Field::getElementComponent : Wrong component number");
+			}
+		case FACES:
+			switch( comp )
+			{
+				case 0:
+					return _mesh.getFace(i).x();
+				case 1:
+					return _mesh.getFace(i).y();
+				case 2:
+					return _mesh.getFace(i).z();
+				default:
+					cout<<"Wrong component number "<< comp <<" , dimension is "<< _mesh.getSpaceDimension() << ", field values are on FACES"<<endl;
+					throw CdmathException("Field::getElementComponent : Wrong component number");
+			}
+		default:
+			throw CdmathException("Accepted field supports are CELLS, NODES and FACES");
+	}
+	
+}
 //----------------------------------------------------------------------
 void
 Field::setName ( const string fieldName )
-- 
2.39.2