mirror of
https://gitlab.kitware.com/vtk/vtk-m
synced 2024-09-20 02:55:47 +00:00
102 lines
3.3 KiB
C++
102 lines
3.3 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.
|
|
//
|
|
// Copyright 2014 National Technology & Engineering Solutions of Sandia, LLC (NTESS).
|
|
// Copyright 2014 UT-Battelle, LLC.
|
|
// Copyright 2014 Los Alamos National Security.
|
|
//
|
|
// Under the terms of Contract DE-NA0003525 with NTESS,
|
|
// the U.S. Government retains certain rights in this software.
|
|
//
|
|
// Under the terms of Contract DE-AC52-06NA25396 with Los Alamos National
|
|
// Laboratory (LANL), the U.S. Government retains certain rights in
|
|
// this software.
|
|
//============================================================================
|
|
|
|
#ifndef vtk_m_filter_MeshQuality_h
|
|
#define vtk_m_filter_MeshQuality_h
|
|
|
|
#include <vtkm/CellShape.h>
|
|
#include <vtkm/filter/FilterCell.h>
|
|
#include <vtkm/worklet/FieldStatistics.h>
|
|
#include <vtkm/worklet/MeshQuality.h>
|
|
|
|
namespace vtkm
|
|
{
|
|
namespace filter
|
|
{
|
|
|
|
//Names of the available cell metrics, for use in
|
|
//the output dataset fields
|
|
static const std::string MetricNames[] = {
|
|
"area", "aspectGamma", "aspectRatio", "condition", "diagonalRatio", "jacobian",
|
|
"minAngle", "maxAngle", "oddy", "relativeSize", "scaledJacobian", "shape",
|
|
"shear", "skew", "stretch", "taper", "volume", "warpage"
|
|
};
|
|
|
|
//Different cell metrics available to use
|
|
//This must follow the same order as the MetricNames above
|
|
enum class CellMetric
|
|
{
|
|
AREA,
|
|
ASPECT_GAMMA,
|
|
ASPECT_RATIO,
|
|
CONDITION,
|
|
DIAGONAL_RATIO,
|
|
JACOBIAN,
|
|
MIN_ANGLE,
|
|
MAX_ANGLE,
|
|
ODDY,
|
|
RELATIVE_SIZE,
|
|
SCALED_JACOBIAN,
|
|
SHAPE,
|
|
SHEAR,
|
|
SKEW,
|
|
STRETCH,
|
|
TAPER,
|
|
VOLUME,
|
|
WARPAGE,
|
|
NUMBER_OF_CELL_METRICS,
|
|
EMPTY
|
|
};
|
|
|
|
/** \brief Computes the quality of an unstructured cell-based mesh. The quality is defined in terms of the
|
|
* summary statistics (frequency, mean, variance, min, max) of metrics computed over the mesh
|
|
* cells. One of several different metrics can be specified for a given cell type, and the mesh
|
|
* can consist of one or more different cell types. The resulting mesh quality is stored as one
|
|
* or more new fields in the output dataset of this filter, with a separate field for each cell type.
|
|
* Each field contains the metric summary statistics for the cell type.
|
|
* Summary statists with all 0 values imply that the specified metric does not support the cell type.
|
|
*/
|
|
class MeshQuality : public vtkm::filter::FilterCell<MeshQuality>
|
|
{
|
|
public:
|
|
using SupportedTypes = vtkm::TypeListTagFieldVec3;
|
|
|
|
VTKM_CONT MeshQuality(CellMetric);
|
|
void SetOutputName(const std::string& s) { this->OutputName = s; };
|
|
|
|
template <typename T, typename StorageType, typename DerivedPolicy>
|
|
VTKM_CONT vtkm::cont::DataSet DoExecute(
|
|
const vtkm::cont::DataSet& input,
|
|
const vtkm::cont::ArrayHandle<vtkm::Vec<T, 3>, StorageType>& points,
|
|
const vtkm::filter::FieldMetadata& fieldMeta,
|
|
const vtkm::filter::PolicyBase<DerivedPolicy>& policy);
|
|
|
|
private:
|
|
CellMetric MyMetric;
|
|
std::string OutputName;
|
|
};
|
|
|
|
} // namespace filter
|
|
} // namespace vtkm
|
|
|
|
#include <vtkm/filter/MeshQuality.hxx>
|
|
|
|
#endif // vtk_m_filter_MeshQuality_h
|