Marching.hpp

Go to the documentation of this file.

#pragma once
 
#include <wmtk/Mesh.hpp>
 
namespace wmtk::components {
 
class Marching
{
public:
    Marching(
        attribute::MeshAttributeHandle& pos_handle,
        std::map<PrimitiveType, attribute::MeshAttributeHandle>& label_handles,
        const std::vector<int64_t>& input_values,
        const int64_t output_value);
 
    void process();
 
    void add_filter(const attribute::MeshAttributeHandle& label, const int64_t value);
 
    void add_pass_through(const attribute::MeshAttributeHandle& pass_through);
 
    void add_pass_through(const std::vector<attribute::MeshAttributeHandle>& pass_through);
 
    void add_isovalue(const attribute::MeshAttributeHandle& scalar_field, const double isovalue);
 
    void set_isovalue_linesearch_iterations(const int64_t linesearch_iterations);
 
    void invert_filter();
 
private:
    Mesh& m_mesh;
 
    attribute::MeshAttributeHandle m_pos_handle;
    attribute::MeshAttributeHandle m_vertex_tag_handle;
    std::optional<attribute::MeshAttributeHandle> m_edge_tag_handle;
    std::optional<attribute::MeshAttributeHandle> m_face_tag_handle;
    std::vector<int64_t> m_input_values;
    int64_t m_output_value;
 
    std::vector<attribute::MeshAttributeHandle> m_filter_labels;
    std::vector<int64_t> m_filter_values;
 
    std::vector<attribute::MeshAttributeHandle> m_pass_through_attributes;
 
    std::optional<attribute::MeshAttributeHandle> m_scalar_field;
    double m_isovalue = std::numeric_limits<double>::lowest();
 
    int64_t m_linesearch_iterations = 10;
 
    bool m_invert_filter = false;
};
 
} // namespace wmtk::components