<HTML>
<HEAD>
<TITLE>Re: [Paraview] problems generating ghost cells, update on PKdTree questions</TITLE>
</HEAD>
<BODY>
<FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'>I think setting UPDATE_NUMBER_OF_GHOST_LEVELS on the output of the filter before calling Update is the right thing to do. However, in your example code I do not see where that value is getting set. I only see the update extent being set, and that is not sufficient to produce ghost cells.<BR>
<BR>
-Ken<BR>
<BR>
<BR>
On 11/19/09 4:50 PM, "Christine Corbett Moran" <<a href="corbett@physik.uzh.ch">corbett@physik.uzh.ch</a>> wrote:<BR>
<BR>
</SPAN></FONT><BLOCKQUOTE><FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'>I'm having trouble with generating ghost cells, either on my own or<BR>
with D3 (via the duplicate boundary points mode). In neither case am I<BR>
able to see a single ghost level of value 1 for any of my points on a<BR>
test data set, which has 10,000 cells, one point per cell, distributed<BR>
on 4 processors. I try both generating ghost cells at level one within<BR>
my d3 call in my reader as follows:<BR>
vtkSmartPointer<vtkDistributedDataFilter> d3 = \<BR>
vtkSmartPointer<vtkDistributedDataFilter>::New();<BR>
d3->AddInput(tipsyReadInitialOutput);<BR>
d3->UpdateInformation();<BR>
vtkStreamingDemandDrivenPipeline* exec = \<BR>
static_cast<vtkStreamingDemandDrivenPipeline*>(d3->GetExecutive()); <BR>
// adds one ghostlevel<BR>
exec->SetUpdateExtent(exec->GetOutputInformation(0), piece, numPieces, 1);<BR>
d3->Update();<BR>
// changing output to output of d3<BR>
output->ShallowCopy(d3->GetOutput());<BR>
output->GetInformation()->Set(<BR>
vtkDataObject::DATA_NUMBER_OF_GHOST_LEVELS(), 1);<BR>
And, as far as I can tell, the procedure for a filter requesting ghost<BR>
levels is:<BR>
// Requesting one level of ghost cells<BR>
vtkInformation* inInfo =<BR>
inputVector[0]->GetInformationObject(0);<BR>
inInfo->Set(<BR>
vtkStreamingDemandDrivenPipeline::UPDATE_NUMBER_OF_GHOST_LEVELS(),1);<BR>
And my reader or D3 would then have to access the inInfo act on this<BR>
request. However neither my filter nor D3 are generating any ghost<BR>
cells, as is verified by examining:<BR>
output->GetPointData()->GetArray("vtkGhostLevels")->GetTuple(<BR>
nextHaloId)[0]<BR>
Can any one see what's going on? I've read a chunk of vtk source<BR>
referencing vtkGhostLevels or using UPDATE_NUMBER_OF_GHOST_LEVELS but<BR>
haven't found the clue yet.<BR>
<BR>
An update on my D3/PKd tree question: I decided to run D3<BR>
automatically if the user checks a button in my readers and strongly<BR>
recommend they run D3 if they don't use my readers. I found that<BR>
BuildLocator() or BuildLocatorFromPoints(points) never worked as I had<BR>
hoped from the documentation (I wanted to call it once/or on all<BR>
processes on the PkD tree generated by D3 and have it build a locator<BR>
which works across all processes, i.e. can also find points belonging<BR>
to a neighbor process), so instead I build a KdTree locator local to<BR>
each process, separate from the PKdTree generated by D3, and<BR>
coordinate a search if necessary, I can still make this highly data<BR>
parallel for my applications, e.g. finding the mass of all points<BR>
within a given radius by doing local searches/consolidation and just<BR>
sending the results for a final consolidation to the root.<BR>
<BR>
Christine<BR>
_______________________________________________<BR>
Powered by www.kitware.com<BR>
<BR>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a><BR>
<BR>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView">http://paraview.org/Wiki/ParaView</a><BR>
<BR>
Follow this link to subscribe/unsubscribe:<BR>
<a href="http://www.paraview.org/mailman/listinfo/paraview">http://www.paraview.org/mailman/listinfo/paraview</a><BR>
<BR>
<BR>
</SPAN></FONT></BLOCKQUOTE><FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'><BR>
</SPAN></FONT><FONT SIZE="1"><FONT FACE="Monaco, Courier New"><SPAN STYLE='font-size:7.5pt'><BR>
**** Kenneth Moreland<BR>
*** Sandia National Laboratories<BR>
*********** <BR>
*** *** *** email: <a href="kmorel@sandia.gov">kmorel@sandia.gov</a><BR>
** *** ** phone: (505) 844-8919<BR>
*** web: <a href="http://www.cs.unm.edu/~kmorel">http://www.cs.unm.edu/~kmorel</a><BR>
</SPAN></FONT></FONT><FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'><BR>
</SPAN></FONT>
</BODY>
</HTML>