Commit 88f1aa6d authored by BRAMAS Berenger's avatar BRAMAS Berenger

do not move direct result on the file result

parent 0dad7e35
......@@ -168,17 +168,17 @@ int main(int argc, char ** argv){
octreeData << "Leaf width " << widthAtLeafLevel << std::endl ;
octreeData << "Box Corner "<< boxCorner << std::endl<< std::endl ;
do{
auto * const FRestrict cell = octreeIterator.getCurrentCell();
FTreeCoordinate coordinate = cell->getCoordinate() ;
FPoint leafCenter(FReal(coordinate.getX()) * widthAtLeafLevel + widthAtLeafLevelDiv2 + boxCorner.getX(),
FReal(coordinate.getY()) * widthAtLeafLevel + widthAtLeafLevelDiv2 + boxCorner.getX(),
FReal(coordinate.getZ()) * widthAtLeafLevel + widthAtLeafLevelDiv2 + boxCorner.getX());
octreeData << "Leaf " << cell->getMortonIndex() << std::endl
<< " Center "<< coordinate << std::endl
<< " Center "<< leafCenter
<< std::endl ;
} while(octreeIterator.moveRight());
do{
auto * const FRestrict cell = octreeIterator.getCurrentCell();
FTreeCoordinate coordinate = cell->getCoordinate() ;
FPoint leafCenter(FReal(coordinate.getX()) * widthAtLeafLevel + widthAtLeafLevelDiv2 + boxCorner.getX(),
FReal(coordinate.getY()) * widthAtLeafLevel + widthAtLeafLevelDiv2 + boxCorner.getX(),
FReal(coordinate.getZ()) * widthAtLeafLevel + widthAtLeafLevelDiv2 + boxCorner.getX());
octreeData << "Leaf " << cell->getMortonIndex() << std::endl
<< " Center "<< coordinate << std::endl
<< " Center "<< leafCenter
<< std::endl ;
} while(octreeIterator.moveRight());
//
FIOVtk vtkfile ;
vtkfile.writeOctree("octreeFile.vtk","Octree ", tree) ;
......@@ -300,7 +300,7 @@ int main(int argc, char ** argv){
dpotential += part.potential * part.physicalValue;
particlesDirect[idxTarget] = part;
// particlesDirect[idxTarget] = part;
}
dpotential *= 0.5*scaleEnergy ;
......
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