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
|
// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause
#include "vtkPResampleToImage.h"
#include "vtkActor.h"
#include "vtkCamera.h"
#include "vtkClipDataSet.h"
#include "vtkCompositeRenderManager.h"
#include "vtkContourFilter.h"
#include "vtkMPIController.h"
#include "vtkNew.h"
#include "vtkPExtractVOI.h"
#include "vtkPieceScalars.h"
#include "vtkPolyDataMapper.h"
#include "vtkRTAnalyticSource.h"
#include "vtkRegressionTestImage.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkSmartPointer.h"
// clang-format off
#include "vtk_diy2.h" // must include this before any diy header
#include VTK_DIY2(diy/mpi.hpp)
// clang-format on
int TestPResampleToImage(int argc, char* argv[])
{
diy::mpi::environment mpienv(argc, argv);
vtkNew<vtkMPIController> controller;
controller->Initialize(&argc, &argv, true);
diy::mpi::communicator world;
// Setup parallel rendering
vtkNew<vtkCompositeRenderManager> prm;
vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::Take(prm->MakeRenderer());
vtkSmartPointer<vtkRenderWindow> renWin =
vtkSmartPointer<vtkRenderWindow>::Take(prm->MakeRenderWindow());
renWin->AddRenderer(renderer);
renWin->DoubleBufferOn();
renWin->SetMultiSamples(0);
vtkNew<vtkRenderWindowInteractor> iren;
iren->SetRenderWindow(renWin);
prm->SetRenderWindow(renWin);
prm->SetController(controller);
// Create Pipeline
vtkNew<vtkRTAnalyticSource> wavelet;
wavelet->SetWholeExtent(0, 31, 0, 31, 0, 31);
wavelet->SetCenter(16, 16, 16);
vtkNew<vtkClipDataSet> clip;
clip->SetInputConnection(wavelet->GetOutputPort());
clip->SetValue(157);
vtkNew<vtkPResampleToImage> resample;
resample->SetUseInputBounds(true);
resample->SetSamplingDimensions(64, 64, 64);
resample->SetInputConnection(clip->GetOutputPort());
vtkNew<vtkPExtractVOI> voi;
voi->SetVOI(4, 59, 4, 59, 4, 59);
voi->SetInputConnection(resample->GetOutputPort());
vtkNew<vtkContourFilter> contour;
contour->SetValue(0, 200);
contour->ComputeNormalsOn();
contour->SetInputConnection(voi->GetOutputPort());
vtkNew<vtkPieceScalars> pieceScalars;
pieceScalars->SetInputConnection(contour->GetOutputPort());
pieceScalars->SetScalarModeToCellData();
// Execute pipeline and render
vtkNew<vtkPolyDataMapper> mapper;
mapper->SetInputConnection(pieceScalars->GetOutputPort());
mapper->SetScalarModeToUseCellFieldData();
mapper->SelectColorArray("Piece");
mapper->SetScalarRange(0, world.size() - 1);
mapper->SetPiece(world.rank());
mapper->SetNumberOfPieces(world.size());
mapper->Update();
vtkNew<vtkActor> actor;
actor->SetMapper(mapper);
renderer->AddActor(actor);
int retVal;
if (world.rank() == 0)
{
prm->ResetAllCameras();
renderer->GetActiveCamera()->Azimuth(90);
renWin->Render();
retVal = vtkRegressionTester::Test(argc, argv, renWin, 10);
if (retVal == vtkRegressionTester::DO_INTERACTOR)
{
prm->StartInteractor();
}
controller->TriggerBreakRMIs();
}
else
{
prm->StartServices();
}
world.barrier();
diy::mpi::broadcast(world, retVal, 0);
controller->Finalize(true);
return !retVal;
}
|
