Commit 9a8544e3 authored by Pierre Aubert's avatar Pierre Aubert
Browse files

Seems to be bugs in the intrinsics version of the Gray Scott reaction

parent 82250241
...@@ -13,8 +13,8 @@ ...@@ -13,8 +13,8 @@
/** @param nbElement : number of elements of the tables /** @param nbElement : number of elements of the tables
*/ */
void evaluateGrayScott(size_t nbElement){ void evaluateGrayScott(size_t nbElement){
size_t nbRow(nbElement*PLIB_VECTOR_SIZE_FLOAT/2lu); size_t nbRow(nbElement*PLIB_VECTOR_SIZE_FLOAT);
size_t nbCol(nbElement*PLIB_VECTOR_SIZE_FLOAT); size_t nbCol(nbElement*PLIB_VECTOR_SIZE_FLOAT*2lu);
nbElement = nbRow*nbCol; nbElement = nbRow*nbCol;
......
...@@ -15,8 +15,8 @@ ...@@ -15,8 +15,8 @@
/** @param nbElement : number of elements of the tables /** @param nbElement : number of elements of the tables
*/ */
void evaluateGrayScott(size_t nbElement){ void evaluateGrayScott(size_t nbElement){
size_t nbRow(nbElement*PLIB_VECTOR_SIZE_FLOAT/2lu); size_t nbRow(nbElement*PLIB_VECTOR_SIZE_FLOAT);
size_t nbCol(nbElement*PLIB_VECTOR_SIZE_FLOAT); size_t nbCol(nbElement*PLIB_VECTOR_SIZE_FLOAT*2lu);
nbElement = nbRow*nbCol; nbElement = nbRow*nbCol;
......
...@@ -14,8 +14,8 @@ ...@@ -14,8 +14,8 @@
/** @param nbElement : number of elements of the tables /** @param nbElement : number of elements of the tables
*/ */
void evaluateGrayScott(size_t nbElement){ void evaluateGrayScott(size_t nbElement){
size_t nbRow(nbElement*PLIB_VECTOR_SIZE_FLOAT/2lu); size_t nbRow(nbElement*PLIB_VECTOR_SIZE_FLOAT);
size_t nbCol(nbElement*PLIB_VECTOR_SIZE_FLOAT); size_t nbCol(nbElement*PLIB_VECTOR_SIZE_FLOAT*2lu);
nbElement = nbRow*nbCol; nbElement = nbRow*nbCol;
......
...@@ -115,10 +115,10 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte ...@@ -115,10 +115,10 @@ bool simulateImage(size_t nbRow, size_t nbCol, size_t nbImage, size_t nbExtraSte
}else{ }else{
tmpScalOutV.fromVecToScalNeigbhour(tmpVecInV); //The pointers were swaped tmpScalOutV.fromVecToScalNeigbhour(tmpVecInV); //The pointers were swaped
} }
// fullMat.setRow(i, tmpScalOutV.getData()); fullMat.setRow(i, tmpScalOutV.getData());
// fullMat.setRow(i, tmpV1); // fullMat.setRow(i, tmpV1);
fullMat.setRow(i, tmpV2); // fullMat.setRow(i, tmpV2);
} }
progress.finish(); progress.finish();
std::cerr << "Done" << std::endl; std::cerr << "Done" << std::endl;
......
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