diff --git a/DEMSystems/sphereDEMSystem/sphereDEMSystem.H b/DEMSystems/sphereDEMSystem/sphereDEMSystem.H
index 723311ee..7e6b7521 100644
--- a/DEMSystems/sphereDEMSystem/sphereDEMSystem.H
+++ b/DEMSystems/sphereDEMSystem/sphereDEMSystem.H
@@ -56,8 +56,14 @@ protected:
 
 	uniquePtr<interaction> 			interaction_ = nullptr;
 
-	
 
+	int32 									numDomains_;
+
+	VectorDual<box> 				domains_;
+
+	Vector<int32Vector_H> 	parIndexInDomain_;
+
+	
 	auto& Control()
 	{
 		return Control_();
@@ -102,6 +108,11 @@ public:
 			Control_->g().z()};
 	}
 
+	bool inline usingDoulle()const
+	{
+		return pFlow::usingDouble__;
+	}
+
 };
 
 
diff --git a/src/phasicFlow/types/basicTypes/builtinTypes.H b/src/phasicFlow/types/basicTypes/builtinTypes.H
index 3f0347dd..6896231f 100755
--- a/src/phasicFlow/types/basicTypes/builtinTypes.H
+++ b/src/phasicFlow/types/basicTypes/builtinTypes.H
@@ -31,34 +31,36 @@ namespace pFlow
 #ifdef pFlow_Build_Double 
   #define useDouble 1
 inline const char* floatingPointType__ = "double";
+inline const bool usingDouble__ = true;
 #else
   #define useDouble 0
 inline const char* floatingPointType__ = "float";
+inline const bool usingDouble__ = false;
 #endif
 
 
 // scalars
 #if useDouble
-	using real 	= double;
+  using real  = double;
 #else
-	using real 	= float;
+  using real  = float;
 #endif
 
-using int8 		= signed char;
+using int8    = signed char;
 
-using int16 	=	short int;
+using int16   = short int;
 
-using int32 	=	int;
+using int32   = int;
 
-using int64 	=	long long int;
+using int64   = long long int;
 
 using uint16  = unsigned short int  ;
 
-using uint32 	=	unsigned int;
+using uint32  = unsigned int;
 
 using label   = std::size_t;
 
-using word 		= std::string;
+using word    = std::string;
 
 inline 
 auto floatingPointDescription()