polyMesh::findCell slow in parallel, but meshSearch::findCell sometimes fails

pbachant · November 13, 2015, 02:04

I am working on an fvOption source that moves around the mesh (actuator line), so I often need to lookup what cell an element is inside to read local velocity and cell volume. The code works great in serial using polyMesh::findCell, but this slows down a lot in parallel.

I tried replacing these calls with meshSearch::findCell(position, 0), and things get much faster, but sometimes cells are not found on any processors when run in parallel, which shouldn't be possible, right?

Any ideas on how I might either fix up meshSearch::findCell or speed up polyMesh::findCell?

hamedhiv · August 5, 2017, 08:57

mesh search class does not support parallel proccessing.
you should use findCell function.

localCellId [Pstream::myProcNo()]= mesh().findCell(location);
reduce(localCellId, sumOp<List<label> >());

November 13, 2015, 02:04	polyMesh::findCell slow in parallel, but meshSearch::findCell sometimes fails	#1
pbachant Senior Member Pete Bachant Join Date: Jun 2012 Location: Boston, MA Posts: 173 Rep Power: 14	I am working on an fvOption source that moves around the mesh (actuator line), so I often need to lookup what cell an element is inside to read local velocity and cell volume. The code works great in serial using polyMesh::findCell, but this slows down a lot in parallel. I tried replacing these calls with meshSearch::findCell(position, 0), and things get much faster, but sometimes cells are not found on any processors when run in parallel, which shouldn't be possible, right? Any ideas on how I might either fix up meshSearch::findCell or speed up polyMesh::findCell? __________________ Home \| Twitter \| GitHub

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Explicitly filtered LES	saeedi	Main CFD Forum	16	October 14, 2015 12:58
simpleFoam parallel	AndrewMortimer	OpenFOAM Running, Solving & CFD	12	August 7, 2015 19:45
[swak4Foam] fails in parallel with -otherTime?	Phicau	OpenFOAM Community Contributions	3	June 26, 2013 14:00
rhoCentralFoam solver with Slip BCs fails in Parallel Only	JLight	OpenFOAM Running, Solving & CFD	2	October 11, 2012 22:08
Grid Check Fails in Parallel Processing Mode	askance	Main CFD Forum	0	October 20, 2010 11:11

August 5, 2017, 08:57		#2
hamedhiv New Member hamed eshraghi Join Date: Jan 2016 Posts: 3 Rep Power: 10	mesh search class does not support parallel proccessing. you should use findCell function. localCellId [Pstream::myProcNo()]= mesh().findCell(location); reduce(localCellId, sumOp<List<label> >());