diff --git a/src/Makefile.am b/src/Makefile.am
index 7aa6946c8f77151e07a9d31229b205831decaff0..d631e1b8b63ff3e9792375ef401153d2beb7d2e8 100755
--- a/src/Makefile.am
+++ b/src/Makefile.am
@@ -2,7 +2,7 @@
 LINKLIBS = $(OPENMP_LDFLAGS) $(LIBHILBERT_LIB) $(PETSC_LIB) $(SUITESPARSE_LIBS) $(LAPACK_LIBS) $(BLAS_LIBS)  $(METIS_LIB) $(PARMETIS_LIB) $(DEFAULT_LIB) $(PTHREAD_LIBS) $(OPT_LIBS) $(BOOST_LDFLAGS) $(BOOST_PROGRAM_OPTIONS_LIB) $(BOOST_IOSTREAMS_LIB) $(HDF5_LDFLAGS)  $(HDF5_LIBS) $(LIBQUADMATH) $(OPENMP_LDFLAGS) $(LIBIFCORE)
 
 noinst_PROGRAMS = numerics
-numerics_SOURCES = main.cpp ../../openfpm_vcluster/src/VCluster.cpp ../../openfpm_devices/src/memory/HeapMemory.cpp ../../openfpm_devices/src/memory/PtrMemory.cpp ../../openfpm_devices/src/Memleak_check.cpp
+numerics_SOURCES = main.cpp ../../openfpm_vcluster/src/VCluster/VCluster.cpp ../../openfpm_devices/src/memory/HeapMemory.cpp ../../openfpm_devices/src/memory/PtrMemory.cpp ../../openfpm_devices/src/Memleak_check.cpp
 numerics_CXXFLAGS = $(OPENMP_CFLAGS) $(LIBHILBERT_INCLUDE) $(AM_CXXFLAGS) $(HDF5_CPPFLAGS) $(INCLUDES_PATH) $(BOOST_CPPFLAGS) $(SUITESPARSE_INCLUDE) $(METIS_INCLUDE) $(PARMETIS_INCLUDE) $(EIGEN_INCLUDE) $(PETSC_INCLUDE) -Wno-deprecated-declarations -Wno-unused-local-typedefs
 numerics_CFLAGS = $(CUDA_CFLAGS)
 numerics_LDADD = $(LINKLIBS) -lparmetis -lmetis
diff --git a/src/Matrix/SparseMatrix_Eigen.hpp b/src/Matrix/SparseMatrix_Eigen.hpp
index 9a5fc42ab86e39e2f04828b7fe0b1cdd2bbc6a7e..ce65579714f51f156b5b47114ee0d38b298b6df1 100644
--- a/src/Matrix/SparseMatrix_Eigen.hpp
+++ b/src/Matrix/SparseMatrix_Eigen.hpp
@@ -10,7 +10,7 @@
 
 #include "Vector/map_vector.hpp"
 #include <boost/mpl/int.hpp>
-#include "VCluster.hpp"
+#include "VCluster/VCluster.hpp"
 
 #define EIGEN_TRIPLET 1
 
diff --git a/src/Operators/Vector/vector_dist_operators_apply_kernel.hpp b/src/Operators/Vector/vector_dist_operators_apply_kernel.hpp
index 59a51645eaadfeabc78a919b78a128c87f46adb4..447cfcede8044095c1103360be1e785041544d54 100644
--- a/src/Operators/Vector/vector_dist_operators_apply_kernel.hpp
+++ b/src/Operators/Vector/vector_dist_operators_apply_kernel.hpp
@@ -256,7 +256,7 @@ class vector_dist_expression_op<exp1,vector_type,VECT_APPLYKER_IN>
 	//! The vector that contain the particles
 	const vector_orig & vd;
 
-	//! Get the return type of apply the kernel to a particle
+	//! Get the return type of applying the kernel to a particle
 	typedef typename apply_kernel_rtype<decltype(o1.value(vect_dist_key_dx(0)))>::rtype rtype;
 
 public:
diff --git a/src/Operators/Vector/vector_dist_operators_functions.hpp b/src/Operators/Vector/vector_dist_operators_functions.hpp
index 800ab69fa3af90da898370334bf8122016cc4d30..d1634e83b33eafe2d7b2bc9dc60055757eb731a7 100644
--- a/src/Operators/Vector/vector_dist_operators_functions.hpp
+++ b/src/Operators/Vector/vector_dist_operators_functions.hpp
@@ -207,15 +207,17 @@ CREATE_VDIST_ARG2_FUNC(pmul,pmul,VECT_PMUL)
 
 ////////// Special function reduce /////////////////////////
 
-/*! \brief Expression that create a reduction
+
+/*! \brief expression that encapsulate a vector reduction expression
  *
- * \tparam exp1 expression
- * \tparam vector_type type of the vector
+ * \tparam exp1 expression 1
+ * \tparam vector_type type of vector on which the expression is acting
  *
  */
 template <typename exp1, typename vector_type>
 class vector_dist_expression_op<exp1,vector_type,VECT_SUM_REDUCE>
 {
+
 	//! expression on which apply the reduction
 	const exp1 o1;
 
@@ -230,13 +232,13 @@ class vector_dist_expression_op<exp1,vector_type,VECT_SUM_REDUCE>
 
 public:
 
-	//! Constructor from expression and the vector of particles
+	//! constructor from an epxression exp1 and a vector vd
 	vector_dist_expression_op(const exp1 & o1, const vector_type & vd)
 	:o1(o1), vd(vd)
 	{}
 
 	//! sum reduction require initialization where we calculate the reduction
-	//! this produce a cache for the calculated value
+	// this produce a cache for the calculated value
 	inline void init() const
 	{
 		o1.init();
@@ -269,7 +271,7 @@ public:
 	}
 };
 
-
+//! Reduce function (it generate an expression)
 template<typename exp1, typename exp2_, unsigned int op1, typename vector_type>
 inline vector_dist_expression_op<vector_dist_expression_op<exp1,exp2_,op1>,vector_type,VECT_SUM_REDUCE>
 rsum(const vector_dist_expression_op<exp1,exp2_,op1> & va, const vector_type & vd)
@@ -279,6 +281,7 @@ rsum(const vector_dist_expression_op<exp1,exp2_,op1> & va, const vector_type & v
 	return exp_sum;
 }
 
+//! Reduce function (It generate an expression)
 template<unsigned int prp1, typename v1, typename vector_type>
 inline vector_dist_expression_op<vector_dist_expression<prp1,v1>,vector_type,VECT_SUM_REDUCE>
 rsum(const vector_dist_expression<prp1,v1> & va, const vector_type & vd)
diff --git a/src/Operators/Vector/vector_dist_operators_unit_tests.hpp b/src/Operators/Vector/vector_dist_operators_unit_tests.hpp
index 2bd0d4aac0ab1c8a4a68ff4ec74864045969995a..f1f5e5ef42beb9933e173fa75b9fd8354729cf9c 100644
--- a/src/Operators/Vector/vector_dist_operators_unit_tests.hpp
+++ b/src/Operators/Vector/vector_dist_operators_unit_tests.hpp
@@ -652,10 +652,10 @@ typedef vector_dist<3,float,aggregate<float,float,float,VectorS<3,float>,VectorS
 //! Exponential kernel
 struct exp_kernel
 {
-	//! variance
+	//! variance of the exponential kernel
 	float var;
 
-	//! Set the variance of the exponential kernel
+	//! Exponential kernel giving variance
 	exp_kernel(float var)
 	:var(var)
 	{}
diff --git a/src/unit_test_init_cleanup.hpp b/src/unit_test_init_cleanup.hpp
index 5128804068f705666fba9e5a79b0cb258bd161cb..7cf89f01b1dac63dbc8753048851ff822e90f0fa 100644
--- a/src/unit_test_init_cleanup.hpp
+++ b/src/unit_test_init_cleanup.hpp
@@ -8,7 +8,7 @@
 #ifndef UNIT_TEST_INIT_CLEANUP_HPP_
 #define UNIT_TEST_INIT_CLEANUP_HPP_
 
-#include "VCluster.hpp"
+#include "VCluster/VCluster.hpp"
 
 struct ut_start {
     ut_start()   {