mirror of
https://gitlab.kitware.com/vtk/vtk-m
synced 2024-09-19 10:35:42 +00:00
54 lines
1.7 KiB
C++
54 lines
1.7 KiB
C++
//============================================================================
|
|
// Copyright (c) Kitware, Inc.
|
|
// All rights reserved.
|
|
// See LICENSE.txt 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 <vtkm/cont/Initialize.h>
|
|
|
|
#include <vtkm/io/VTKDataSetReader.h>
|
|
#include <vtkm/io/VTKDataSetWriter.h>
|
|
|
|
#include <vtkm/filter/contour/ClipWithField.h>
|
|
|
|
#include <cstdlib>
|
|
#include <iostream>
|
|
|
|
int main(int argc, char* argv[])
|
|
{
|
|
vtkm::cont::Initialize(argc, argv);
|
|
|
|
if ((argc < 4) || (argc > 5))
|
|
{
|
|
std::cerr << "Usage: " << argv[0] << " in_data.vtk field_name clip_value [out_data.vtk]\n\n";
|
|
std::cerr << "For example, you could use the ucd3d.vtk that comes with the VTK-m source:\n\n";
|
|
std::cerr << " " << argv[0]
|
|
<< " <path-to-vtkm-source>/data/data/unstructured/ucd3d.vtk v 0.3\n";
|
|
return 1;
|
|
}
|
|
std::string infilename = argv[1];
|
|
std::string infield = argv[2];
|
|
double fieldValue = std::atof(argv[3]);
|
|
std::string outfilename = "out_data.vtk";
|
|
if (argc == 5)
|
|
{
|
|
outfilename = argv[4];
|
|
}
|
|
|
|
vtkm::io::VTKDataSetReader reader(infilename);
|
|
vtkm::cont::DataSet input = reader.ReadDataSet();
|
|
|
|
vtkm::filter::contour::ClipWithField clipFilter;
|
|
clipFilter.SetActiveField(infield);
|
|
clipFilter.SetClipValue(fieldValue);
|
|
vtkm::cont::DataSet output = clipFilter.Execute(input);
|
|
|
|
vtkm::io::VTKDataSetWriter writer(outfilename);
|
|
writer.WriteDataSet(output);
|
|
|
|
return 0;
|
|
}
|