1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
/*=========================================================================
Program: Visualization Toolkit
Module: $RCSfile: TestExodusTime.cxx,v $
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRegressionTestImage.h"
#include "vtkCamera.h"
#include "vtkExodusReader.h"
#include "vtkTemporalShiftScale.h"
#include "vtkCompositeDataPipeline.h"
#include "vtkActor.h"
#include "vtkMultiGroupDataGeometryFilter.h"
#include "vtkSmartPointer.h"
#include "vtkTemporalDataSet.h"
#include "vtkThreshold.h"
#include "vtkTemporalInterpolator.h"
#include "vtkPolyDataMapper.h"
//-------------------------------------------------------------------------
int TestExodusTime(int argc, char *argv[])
{
// we have to use a compsite pipeline
vtkCompositeDataPipeline* prototype = vtkCompositeDataPipeline::New();
vtkAlgorithm::SetDefaultExecutivePrototype(prototype);
prototype->Delete();
// create the reader
vtkSmartPointer<vtkExodusReader> reader =
vtkSmartPointer<vtkExodusReader>::New();
reader->SetFileName("C:/can.ex2");
// shift and scale the time range to that it run from -0.5 to 0.5
vtkSmartPointer<vtkTemporalShiftScale> tempss =
vtkSmartPointer<vtkTemporalShiftScale>::New();
tempss->SetScale(232.5);
tempss->SetInputConnection(reader->GetOutputPort());
// interpolate if needed
vtkSmartPointer<vtkTemporalInterpolator> interp =
vtkSmartPointer<vtkTemporalInterpolator>::New();
interp->SetInputConnection(tempss->GetOutputPort());
vtkSmartPointer<vtkThreshold> contour =
vtkSmartPointer<vtkThreshold>::New();
contour->SetInputConnection(interp->GetOutputPort());
contour->ThresholdByUpper(0.5);
vtkSmartPointer<vtkMultiGroupDataGeometryFilter> geom =
vtkSmartPointer<vtkMultiGroupDataGeometryFilter>::New();
geom->SetInputConnection(contour->GetOutputPort());
// map them
vtkSmartPointer<vtkPolyDataMapper> mapper =
vtkSmartPointer<vtkPolyDataMapper>::New();
mapper->SetInputConnection(geom->GetOutputPort());
mapper->ScalarVisibilityOn();
mapper->SetScalarModeToUseCellFieldData();
mapper->SelectColorArray("BlockId");
mapper->SetScalarRange(0,3);
vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
actor->SetMapper(mapper);
vtkSmartPointer<vtkRenderer> renderer =
vtkSmartPointer<vtkRenderer>::New();
vtkSmartPointer<vtkRenderWindow> renWin =
vtkSmartPointer<vtkRenderWindow>::New();
vtkSmartPointer<vtkRenderWindowInteractor> iren =
vtkSmartPointer<vtkRenderWindowInteractor>::New();
renderer->AddActor( actor );
renderer->SetBackground(0.5, 0.5, 0.5);
renWin->AddRenderer( renderer );
renWin->SetSize( 300, 300 );
iren->SetRenderWindow( renWin );
renWin->Render();
renderer->GetActiveCamera()->Elevation(-120);
// ask for some specific data points
vtkStreamingDemandDrivenPipeline *sdd =
vtkStreamingDemandDrivenPipeline::SafeDownCast(geom->GetExecutive());
double times[1];
times[0] = 0;
int i;
for (i = 0; i < 100; ++i)
{
times[0] = i/100.0;
sdd->SetUpdateTimeSteps(0, times, 1);
mapper->Modified();
renderer->ResetCamera();
//renderer->SetBackground(0.5*(i%2), 0.5, 0.5*((i/2)%2));
renWin->Render();
}
int retVal = vtkRegressionTestImage( renWin );
if ( retVal == vtkRegressionTester::DO_INTERACTOR)
{
iren->Start();
}
vtkAlgorithm::SetDefaultExecutivePrototype(0);
return !retVal;
}
|
