mirror of
https://gitlab.kitware.com/vtk/vtk-m
synced 2024-09-19 10:35:42 +00:00
478f765db2
Previously, point fields compressed by ZFP were attached as point fields on the output. However, using them as a point field would cause problems. So, instead attache them as `WholeDataSet` fields. Also fixed a problem where the 1D decompressor created an output of the wrong size.
48 lines
1.8 KiB
C++
48 lines
1.8 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/filter/zfp/ZFPCompressor3D.h>
|
|
#include <vtkm/filter/zfp/worklet/ZFPCompressor.h>
|
|
|
|
namespace vtkm
|
|
{
|
|
namespace filter
|
|
{
|
|
namespace zfp
|
|
{
|
|
//-----------------------------------------------------------------------------
|
|
VTKM_CONT vtkm::cont::DataSet ZFPCompressor3D::DoExecute(const vtkm::cont::DataSet& input)
|
|
{
|
|
vtkm::cont::CellSetStructured<3> cellSet;
|
|
input.GetCellSet().AsCellSet(cellSet);
|
|
vtkm::Id3 pointDimensions = cellSet.GetPointDimensions();
|
|
|
|
vtkm::cont::ArrayHandle<vtkm::Int64> compressed;
|
|
|
|
vtkm::worklet::ZFPCompressor compressor;
|
|
using SupportedTypes = vtkm::List<vtkm::Int32, vtkm::Float32, vtkm::Float64>;
|
|
this->GetFieldFromDataSet(input)
|
|
.GetData()
|
|
.CastAndCallForTypesWithFloatFallback<SupportedTypes, VTKM_DEFAULT_STORAGE_LIST>(
|
|
[&](const auto& concrete) {
|
|
compressed = compressor.Compress(concrete, rate, pointDimensions);
|
|
});
|
|
|
|
// Note: the compressed array is set as a WholeDataSet field. It is really associated with
|
|
// the points, but the size does not match and problems will occur if the user attempts to
|
|
// use it as a point data set. The decompressor will place the data back as a point field.
|
|
// (This might cause issues if cell fields are ever supported.)
|
|
return this->CreateResultField(
|
|
input, "compressed", vtkm::cont::Field::Association::WholeDataSet, compressed);
|
|
}
|
|
} // namespace zfp
|
|
} // namespace filter
|
|
} // namespace vtkm
|