<html>
<body>
I called the Update function on filter1 and I get the following error. I
have included the complete code below. <br>
================================================================<br>
vtkStreamingDemandDrivenPipeline (0x9f128f0): Input port 0 of algorithm
vtkAppendPolyData(0x9f112e0) has 0 connections but is not optional.<br>
================================================================<br>
Thanks. <br>
Sunil<br><br>
<i>#include "vtkContourFilter.h"<br>
#include "vtkHierarchicalDataSet.h"<br>
#include "vtkHierarchicalDataSetGeometryFilter.h"<br>
#include "vtkPolyData.h"<br>
#include "vtkXMLMultiGroupDataWriter.h"<br>
#include "vtkXMLHierarchicalDataReader.h"<br>
#include "vtkPointData.h"<br>
#include "vtkHierarchicalDataSetAlgorithm.h"<br>
#include "vtkPolyDataNormals.h"<br>
#include <iostream><br>
using namespace std;<br>
int main(int argc, char* argv[])<br>
{<br>
//read in HDS<br>
vtkXMLHierarchicalDataReader* reader =
vtkXMLHierarchicalDataReader::New();<br>
reader->SetFileName( "mb.vtu" );<br>
reader->Update();<br>
<br>
vtkMultiGroupDataSet* mgd;<br>
mgd = (reader->GetOutput(0));<br>
cout<<"Getting number of groups
:"<<endl;<br>
cout<<"Number of groups
:"<<mgd->GetNumberOfGroups()<<endl;<br>
int check = 0;<br>
for ( int grp=0;grp<mgd->GetNumberOfGroups();grp++
)<br>
{<br>
for ( int
ds=0;ds<mgd->GetNumberOfDataSets(grp);ds++ )<br>
{<br>
//activate scalar
"Density"<br>
check =
dynamic_cast<vtkDataSet*>(mgd->GetDataSet(grp,ds))-><br>
GetPointData()->SetActiveScalars( "Density" );<br>
cout<<"CHECK
is :"<<check<<endl;<br>
if ( check != -1 )<br>
cout<<"Group :"<<grp<<" Dataset
:"<<ds<<" set Density"<<endl;<br>
}<br>
}<br>
<br>
<br>
////////////setup pipelines<br>
vtkHierarchicalDataSetAlgorithm* alg = <br>
vtkHierarchicalDataSetAlgorithm::New();<br>
alg->SetInput( 0, mgd );<br>
alg->Update();<br>
<br>
<br>
vtkHierarchicalDataSetGeometryFilter* filter1 = <br>
vtkHierarchicalDataSetGeometryFilter::New( );<br>
filter1->SetInputConnection( 0, alg->GetOutputPort(0)
);<br>
filter1->Update();<br>
<br>
vtkContourFilter *contour = vtkContourFilter::New();<br>
contour->SetInputConnection( 0,
filter1->GetOutputPort(0) );<br>
contour->SetValue( 0, 0.2 );<br>
contour->UseScalarTreeOff();<br>
<br>
vtkPolyDataNormals *normals =
vtkPolyDataNormals::New();<br>
normals->SetInputConnection( 0,
contour->GetOutputPort(0) );<br>
<br>
vtkHierarchicalDataSetGeometryFilter* filter = <br>
vtkHierarchicalDataSetGeometryFilter::New( );<br>
filter->SetInput( normals->GetOutput(0) );<br><br>
int numPolys =
filter->GetOutput()->GetNumberOfPolys();<br>
std::cout << " The number
of polys is "<< numPolys << std::endl;<br><br>
//clean up<br>
reader->Delete();<br>
return 0;<br>
}<br><br>
<br><br>
</i>At 09:09 AM 6/7/2007, Marc Cotran wrote:<br>
<blockquote type=cite class=cite cite="">Hi Sunil,<br><br>
I'm not sure if this will help you as I've no experience with that
specific filter, but you could try calling Update() on it before your
std::cout call. VTK pipelines are largely demand driven, so if you don't
render the result of the pipeline, for example, the filter does not
actually execute. Calling Update() forces the object to run.<br><br>
Marc<br>
</blockquote></body>
</html>