2019-08-08 22:20:11 +00:00
|
|
|
//=============================================================================
|
|
|
|
//
|
|
|
|
// 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.
|
|
|
|
//
|
|
|
|
//=============================================================================
|
|
|
|
|
2019-08-05 20:15:10 +00:00
|
|
|
#include <cmath>
|
|
|
|
#include <string>
|
|
|
|
#include <vector>
|
|
|
|
|
|
|
|
#include <vtkm/cont/ArrayCopy.h>
|
|
|
|
#include <vtkm/cont/DataSet.h>
|
|
|
|
#include <vtkm/cont/DataSetBuilderUniform.h>
|
2019-08-09 17:34:01 +00:00
|
|
|
#include <vtkm/cont/Initialize.h>
|
|
|
|
|
2020-04-30 12:12:04 +00:00
|
|
|
#include <vtkm/io/VTKDataSetReader.h>
|
|
|
|
#include <vtkm/io/VTKDataSetWriter.h>
|
2019-08-05 20:15:10 +00:00
|
|
|
|
|
|
|
#include <vtkm/filter/LagrangianStructures.h>
|
|
|
|
|
|
|
|
int main(int argc, char** argv)
|
|
|
|
{
|
2019-08-09 21:03:26 +00:00
|
|
|
vtkm::cont::Initialize(argc, argv);
|
2019-08-09 17:34:01 +00:00
|
|
|
|
2019-08-05 20:15:10 +00:00
|
|
|
if (argc < 3)
|
|
|
|
{
|
|
|
|
std::cout << "Usage : flte <input dataset> <vector field name>" << std::endl;
|
|
|
|
}
|
|
|
|
std::string datasetName(argv[1]);
|
|
|
|
std::string variableName(argv[2]);
|
|
|
|
|
|
|
|
std::cout << "Reading input dataset" << std::endl;
|
|
|
|
vtkm::cont::DataSet input;
|
2020-04-30 12:12:04 +00:00
|
|
|
vtkm::io::VTKDataSetReader reader(datasetName);
|
2019-08-05 20:15:10 +00:00
|
|
|
input = reader.ReadDataSet();
|
|
|
|
std::cout << "Read input dataset" << std::endl;
|
|
|
|
|
|
|
|
vtkm::filter::LagrangianStructures lcsFilter;
|
2019-08-13 16:19:32 +00:00
|
|
|
lcsFilter.SetStepSize(0.025f);
|
2019-08-05 20:15:10 +00:00
|
|
|
lcsFilter.SetNumberOfSteps(500);
|
2019-08-13 16:19:32 +00:00
|
|
|
lcsFilter.SetAdvectionTime(0.025f * 500);
|
|
|
|
lcsFilter.SetOutputFieldName("gradient");
|
2019-08-05 20:15:10 +00:00
|
|
|
lcsFilter.SetActiveField(variableName);
|
|
|
|
|
|
|
|
vtkm::cont::DataSet output = lcsFilter.Execute(input);
|
2020-04-30 12:12:04 +00:00
|
|
|
vtkm::io::VTKDataSetWriter writer("out.vtk");
|
2019-08-05 20:15:10 +00:00
|
|
|
writer.WriteDataSet(output);
|
|
|
|
std::cout << "Written output dataset" << std::endl;
|
|
|
|
|
|
|
|
return 0;
|
|
|
|
}
|