Add todo for debugging

cmake/phoenix_base_project.cmake

@@ -37,9 +37,9 @@ function(phoenix_base_project programName programVersion programDescritpion prog
 			set(CMAKE_CXX_FLAGS "--std=c++2a -Wall -Werror -g -O0  -fprofile-arcs -ftest-coverage" PARENT_SCOPE)
 		else()
 			if(GPU_MODE)
-				set(CMAKE_CXX_FLAGS "--std=c++17 -Wall -Werror -g -O2" PARENT_SCOPE)
+				set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} --std=c++17 -Wall -Werror -g -O2" PARENT_SCOPE)
 			else()
-				set(CMAKE_CXX_FLAGS "--std=c++2a -Wall -Werror -g -O2" PARENT_SCOPE)
+				set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} --std=c++2a -Wall -Werror -g -O2" PARENT_SCOPE)
 			endif()
 		endif()
 		

src/cpp20/Sequential/GrayScott/Program/main_intrinsics.cpp

@@ -59,6 +59,8 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte
 	fullMat.setAllDim(nbCol, nbRow);
 	fullMat.resize(nbImage);
 	
+	//TODO : create an other MatrixHdf5 to store the vectorial matrices and see what's going on
+	
 	PTensor<float> tmpInU, tmpInV, tmpOutU, tmpOutV;
 	float *tmpU1 = NULL, *tmpU2 = NULL, *tmpV1 = NULL, *tmpV2 = NULL;
 	allocate_temporary(tmpU1, tmpU2, tmpV1, tmpV2, tmpInU, tmpInV, tmpOutU, tmpOutV, nbRow, nbCol);
@@ -102,8 +104,8 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte
 					ptrVecMatStencil, nbStencilRow, nbStencilCol,
 					diffudionRateU, diffusionRateV, feedRate, killRate, dt);
 			//Let's update the dupplicated values
-			reshuffle_averageDupplicateVecNeighbour(tmpU2, nbVecRow, nbVecCol);
-			reshuffle_averageDupplicateVecNeighbour(tmpV2, nbVecRow, nbVecCol);
+			reshuffle_updateDupplicateVecNeighbour(tmpU2, nbVecRow, nbVecCol, PLIB_VECTOR_SIZE_FLOAT);
+			reshuffle_updateDupplicateVecNeighbour(tmpV2, nbVecRow, nbVecCol, PLIB_VECTOR_SIZE_FLOAT);
 			
 			///Let's swap the pointer
 			swapValue(tmpU1, tmpU2);