Commit ee294536 authored by Pierre Aubert's avatar Pierre Aubert
Browse files

Remove debug stuff for intrinsics Gray Scott

parent 7f9259bd
......@@ -62,9 +62,8 @@ void grayscott_propagation(float * outMatVecU, float * outMatVecV, const float *
long stencilIndexCol(0l);
for(long l(firstColStencil); l < lastColStencil; ++l){
// float deltaSquare(matDeltaSquare[stencilIndexRow*nbStencilCol + stencilIndexCol]);
PRegVecf vecDeltaSquare(vecOne);
// PRegVecf vecDeltaSquare(plib_load_ps(matBroadcastDeltaSquare +
// (stencilIndexRow*nbStencilCol + stencilIndexCol)*PLIB_VECTOR_SIZE_FLOAT));
PRegVecf vecDeltaSquare(plib_load_ps(matBroadcastDeltaSquare +
(stencilIndexRow*nbStencilCol + stencilIndexCol)*PLIB_VECTOR_SIZE_FLOAT));
PRegVecf vecKLU(plib_load_ps(matVecVecU + (k*nbVecCol + l)*PLIB_VECTOR_SIZE_FLOAT));
PRegVecf vecKLV(plib_load_ps(matVecVecV + (k*nbVecCol + l)*PLIB_VECTOR_SIZE_FLOAT));
......
......@@ -93,13 +93,6 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte
float * ptrVecMatStencil = vecMatDeltaSquare.getData();
size_t nbVecRow(tmpVecInV.getFullNbRow()), nbVecCol(tmpVecInV.getNbCol());
MatrixHdf5 fullVecMat;
fullVecMat.setAllDim(nbVecCol, nbVecRow);
fullVecMat.resize(nbImage);
//TODO : create an other MatrixHdf5 to store the vectorial matrices and see what's going on
PTensor<float> tmpScalOutV(AllocMode::ALIGNED);
ProgressTime progress(nbImage);
......@@ -120,12 +113,9 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte
}
if(tmpVecV1 != tmpVecOutV.getData()){
tmpScalOutV.fromVecToScalNeigbhour(tmpVecOutV);
fullVecMat.setRow(i, tmpVecOutV.getData());
}else{
tmpScalOutV.fromVecToScalNeigbhour(tmpVecInV); //The pointers were swaped
fullVecMat.setRow(i, tmpVecInV.getData());
}
fullMat.setRow(i, tmpScalOutV.getData());
// fullMat.setRow(i, tmpVecV1);
......@@ -135,7 +125,6 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte
std::cerr << "Done" << std::endl;
//Let's save the output file
fullMat.write(outputFile);
fullVecMat.write("./output_vec.h5");
return true;
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment