MinEdgeLengthInvariant.cpp

Go to the documentation of this file.

 
#include "MinEdgeLengthInvariant.hpp"
#include <wmtk/Mesh.hpp>
 
namespace wmtk::invariants {
MinEdgeLengthInvariant::MinEdgeLengthInvariant(
    const Mesh& m,
    const TypedAttributeHandle<double>& coordinate,
    double threshold_squared)
    : Invariant(m, true, false, false)
    , m_coordinate_handle(coordinate)
    , m_threshold_squared(threshold_squared)
{}
bool MinEdgeLengthInvariant::before(const simplex::Simplex& t) const
{
    const attribute::Accessor<double> accessor = mesh().create_const_accessor(m_coordinate_handle);
 
    auto p0 = accessor.const_vector_attribute(t.tuple());
    auto p1 =
        accessor.const_vector_attribute(mesh().switch_tuple(t.tuple(), PrimitiveType::Vertex));
    const double l_squared = (p1 - p0).squaredNorm();
    return l_squared > m_threshold_squared;
}
} // namespace wmtk::invariants