#include <PhysicsRemesher.hpp>

Inheritance diagram for polyfem::mesh::PhysicsTetRemesher:

Collaboration diagram for polyfem::mesh::PhysicsTetRemesher:

Public Types
using	Tuple = typename Super::Tuple

Public Types inherited from polyfem::mesh::PhysicsRemesher< wmtk::TetMesh >
using	Tuple = typename Super::Tuple

using	Operations = typename Super::Operations

using	VectorNd = typename Super::VectorNd

Public Types inherited from polyfem::mesh::WildRemesher< WMTKMesh >
using	VectorNd = Eigen::Matrix< double, DIM, 1 >

using	Tuple = typename WMTKMesh::Tuple

Public Types inherited from polyfem::mesh::Remesher
template<typename T >
using	EdgeMap = std::unordered_map< std::array< size_t, 2 >, T, polyfem::utils::HashUnorderedArray< size_t, 2 >, polyfem::utils::EqualUnorderedArray< size_t, 2 > >
	Map from a (sorted) edge to an integer (ID)

template<typename T >
using	FaceMap = std::unordered_map< std::array< size_t, 3 >, T, polyfem::utils::HashUnorderedArray< size_t, 3 >, polyfem::utils::EqualUnorderedArray< size_t, 3 > >
	Map from a (sorted) edge to an integer (ID)

template<typename T >
using	TetMap = std::unordered_map< std::array< size_t, 4 >, T, polyfem::utils::HashUnorderedArray< size_t, 4 >, polyfem::utils::EqualUnorderedArray< size_t, 4 > >

template<typename T >
using	BoundaryMap = std::variant< EdgeMap< T >, FaceMap< T > >

Public Member Functions
	PhysicsTetRemesher (const State &state, const Eigen::MatrixXd &obstacle_displacements, const Eigen::MatrixXd &obstacle_vals, const double current_time, const double starting_energy)

Public Member Functions inherited from polyfem::mesh::PhysicsRemesher< wmtk::TetMesh >
	PhysicsRemesher (const State &state, const Eigen::MatrixXd &obstacle_displacements, const Eigen::MatrixXd &obstacle_vals, const double current_time, const double starting_energy)

virtual	~PhysicsRemesher ()

void	split_edges () override

bool	split_edge_before (const Tuple &t) override

bool	split_edge_after (const Tuple &t) override

void	collapse_edges () override

bool	collapse_edge_before (const Tuple &t) override

bool	collapse_edge_after (const Tuple &t) override

bool	smooth_before (const Tuple &t) override

bool	smooth_after (const Tuple &t) override

Public Member Functions inherited from polyfem::mesh::WildRemesher< WMTKMesh >
	WildRemesher (const State &state, const Eigen::MatrixXd &obstacle_displacements, const Eigen::MatrixXd &obstacle_vals, const double current_time, const double starting_energy)
	Construct a new WildRemesher object.

virtual	~WildRemesher ()=default

virtual void	init (const Eigen::MatrixXd &rest_positions, const Eigen::MatrixXd &positions, const Eigen::MatrixXi &elements, const Eigen::MatrixXd &projection_quantities, const BoundaryMap< int > &boundary_to_id, const std::vector< int > &body_ids, const EdgeMap< double > &elastic_energy, const EdgeMap< double > &contact_energy) override
	Initialize the mesh.

bool	execute () override
	Execute the remeshing.

virtual void	smooth_vertices ()

virtual void	swap_edges ()

int	dim () const override
	Dimension of the mesh.

Eigen::MatrixXd	rest_positions () const override
	Exports rest positions of the stored mesh.

Eigen::MatrixXd	displacements () const override
	Exports positions of the stored mesh.

Eigen::MatrixXd	positions () const override
	Exports displacements of the stored mesh.

Eigen::MatrixXi	edges () const override
	Exports edges of the stored mesh.

Eigen::MatrixXi	elements () const override
	Exports elements of the stored mesh.

Eigen::MatrixXi	boundary_edges () const override
	Exports boundary edges of the stored mesh.

Eigen::MatrixXi	boundary_faces () const override
	Exports boundary faces of the stored mesh.

Eigen::MatrixXd	projection_quantities () const override
	Exports projected quantities of the stored mesh.

BoundaryMap< int >	boundary_ids () const override
	Exports boundary ids of the stored mesh.

std::vector< int >	body_ids () const override
	Exports body ids of the stored mesh.

std::vector< int >	boundary_nodes (const Eigen::VectorXi &vertex_to_basis) const override
	Get the boundary nodes of the stored mesh.

int	n_quantities () const override
	Number of projection quantities (not including the position)

std::vector< Tuple >	get_facets () const
	Get a vector of all facets (edges or triangles)

std::vector< Tuple >	get_elements () const
	Get a vector of all elements (triangles or tetrahedra)

void	set_rest_positions (const Eigen::MatrixXd &positions) override
	Set rest positions of the stored mesh.

void	set_positions (const Eigen::MatrixXd &positions) override
	Set deformed positions of the stored mesh.

void	set_projection_quantities (const Eigen::MatrixXd &projection_quantities) override
	Set projected quantities of the stored mesh.

void	set_fixed (const std::vector< bool > &fixed) override
	Set if a vertex is fixed.

void	set_boundary_ids (const BoundaryMap< int > &boundary_to_id) override
	Set the boundary IDs of all edges.

void	set_body_ids (const std::vector< int > &body_ids) override
	Set the body IDs of all elements.

double	rest_edge_length (const Tuple &e) const
	Compute the length of an edge.

double	deformed_edge_length (const Tuple &e) const

VectorNd	rest_edge_center (const Tuple &e) const
	Compute the center of the edge.

VectorNd	deformed_edge_center (const Tuple &e) const

Eigen::VectorXd	edge_adjacent_element_volumes (const Tuple &e) const

double	element_volume (const Tuple &e) const
	Compute the volume (area) of an tetrahedron (triangle) element.

bool	is_boundary_vertex (const Tuple &v) const
	Is the given vertex tuple on the boundary of the mesh?

bool	is_body_boundary_vertex (const Tuple &v) const
	Is the given vertex tuple on the boundary of a body?

bool	is_boundary_edge (const Tuple &e) const
	Is the given edge tuple on the boundary of the mesh?

bool	is_body_boundary_edge (const Tuple &e) const
	Is the given edge tuple on the boundary of a body?

bool	is_boundary_facet (const Tuple &t) const
	Is the given tuple on the boundary of the mesh?

bool	is_boundary_op () const
	Is the currently cached operation a boundary operation?

std::vector< Tuple >	boundary_facets (std::vector< int > *boundary_ids=nullptr) const
	Get the boundary facets of the mesh.

std::array< Tuple, DIM >	facet_vertices (const Tuple &t) const
	Get the vertex tuples of a facet.

std::array< size_t, DIM >	facet_vids (const Tuple &t) const
	Get the vertex ids of a facet.

std::array< Tuple, VERTICES_PER_ELEMENT >	element_vertices (const Tuple &t) const
	Get the vertex tuples of an element.

std::array< size_t, VERTICES_PER_ELEMENT >	element_vids (const Tuple &t) const
	Get the vertex ids of an element.

void	element_aabb (const Tuple &t, polyfem::VectorNd &el_min, polyfem::VectorNd &el_max) const
	Get a AABB for an element.

std::array< size_t, VERTICES_PER_ELEMENT >	orient_preserve_element_reorder (const std::array< size_t, VERTICES_PER_ELEMENT > &conn, const size_t v0) const
	Reorder the element vertices so that the first vertex is v0.

std::vector< Tuple >	get_one_ring_elements_for_vertex (const Tuple &t) const
	Get the one ring of elements around a vertex.

std::vector< Tuple >	get_one_ring_boundary_edges_for_vertex (const Tuple &v) const

std::array< Tuple, 2 >	get_boundary_faces_for_edge (const Tuple &e) const

std::vector< Tuple >	get_one_ring_boundary_faces_for_vertex (const Tuple &v) const

std::vector< Tuple >	get_edges_for_elements (const std::vector< Tuple > &elements) const

size_t	facet_id (const Tuple &t) const
	Get the id of a facet (edge for triangle, triangle for tetrahedra)

size_t	element_id (const Tuple &t) const
	Get the id of an element (triangle or tetrahedra)

Tuple	tuple_from_element (size_t elem_id) const
	Get a tuple of an element.

Tuple	tuple_from_facet (size_t elem_id, int local_facet_id) const
	Get a tuple of an element with a local facet.

std::vector< Tuple >	get_incident_elements_for_edge (const Tuple &t) const
	Get the incident elements for an edge.

void	extend_local_patch (std::vector< Tuple > &patch) const
	Extend the local patch by including neighboring elements.

Tuple	opposite_vertex_on_face (const Tuple &e) const
	Get the opposite vertex on a face.

CollapseEdgeTo	collapse_boundary_edge_to (const Tuple &e) const
	Determine where to collapse an edge to.

void	write_edge_ranks_mesh (const EdgeMap< typename EdgeAttributes::EnergyRank > &elastic_ranks, const EdgeMap< typename EdgeAttributes::EnergyRank > &contact_ranks) const

EdgeAttributes &	edge_attr (const size_t e_id)
	Get a reference to an edge's attributes.

const EdgeAttributes &	edge_attr (const size_t e_id) const
	Get a const reference to an edge's attributes.

Public Member Functions inherited from polyfem::mesh::Remesher
	Remesher (const State &state, const Eigen::MatrixXd &obstacle_displacements, const Eigen::MatrixXd &obstacle_quantities, const double current_time, const double starting_energy)
	Construct a new Remesher object.

virtual	~Remesher ()=default

virtual bool	is_volume () const
	Is the mesh a volumetric mesh.

const Obstacle &	obstacle () const
	Get a reference to the collision obstacles.

const Eigen::MatrixXd &	obstacle_displacements () const
	Get a reference to the collision obstacles' displacements.

const Eigen::MatrixXd &	obstacle_quantities () const
	Get a reference to the collision obstacles' extra quantities.

void	write_mesh (const std::string &path) const
	Writes a VTU mesh file.

void	init_assembler (const std::vector< int > &body_ids) const
	Create an assembler object.

Private Types
using	Super = PhysicsRemesher< wmtk::TetMesh >

Additional Inherited Members
Static Public Member Functions inherited from polyfem::mesh::Remesher
static Eigen::MatrixXd	combine_time_integrator_quantities (const std::shared_ptr< time_integrator::ImplicitTimeIntegrator > &time_integrator)
	Combine the quantities of a time integrator into a single matrix (one column per quantity)

static void	split_time_integrator_quantities (const Eigen::MatrixXd &quantities, const int dim, Eigen::MatrixXd &x_prevs, Eigen::MatrixXd &v_prevs, Eigen::MatrixXd &a_prevs)
	Split the quantities of a time integrator into separate vectors.

static int	build_bases (const Mesh &mesh, const std::string &assembler_formulation, std::vector< polyfem::basis::ElementBases > &bases, std::vector< LocalBoundary > &local_boundary, Eigen::VectorXi &vertex_to_basis)
	Build bases for a given mesh (V, F)

static void	log_timings ()

Public Attributes inherited from polyfem::mesh::PhysicsRemesher< wmtk::TetMesh >
wmtk::AttributeCollection< VertexAttributes >	vertex_attrs

Public Attributes inherited from polyfem::mesh::WildRemesher< WMTKMesh >
wmtk::AttributeCollection< VertexAttributes >	vertex_attrs

wmtk::AttributeCollection< BoundaryAttributes >	boundary_attrs

wmtk::AttributeCollection< ElementAttributes >	element_attrs

Public Attributes inherited from polyfem::mesh::Remesher
const State &	state
	Reference to the simulation state.

int	max_op_attempts = 1

Static Public Attributes inherited from polyfem::mesh::WildRemesher< WMTKMesh >
static constexpr int	DIM

static constexpr int	VERTICES_PER_ELEMENT

static constexpr int	EDGES_PER_ELEMENT

static constexpr int	FACETS_PER_ELEMENT

static constexpr wmtk::ExecutionPolicy	EXECUTION_POLICY = wmtk::ExecutionPolicy::kSeq
	Current execuation policy (sequencial or parallel)

Static Public Attributes inherited from polyfem::mesh::Remesher
static std::unordered_map< std::string, utils::Timing >	timings
	Timings for the remeshing operations.

static double	total_time = 0

static size_t	num_solves = 0

static size_t	total_ndofs = 0

Protected Types inherited from polyfem::mesh::WildRemesher< WMTKMesh >
using	Operations = std::vector< std::pair< std::string, Tuple > >

Protected Member Functions inherited from polyfem::mesh::PhysicsRemesher< wmtk::TetMesh >
Operations	renew_neighbor_tuples (const std::string &op, const std::vector< Tuple > &tris) const override
	Renew the neighbor tuples of an operation.

bool	local_relaxation (const Tuple &t, const double acceptance_tolerance)
	Relax a local n-ring around a vertex.

bool	local_relaxation (const VectorNd &center, const double acceptance_tolerance)
	Relax a local n-ring around a vertex.

bool	local_relaxation (const std::vector< Tuple > &local_mesh_tuples, const double acceptance_tolerance)
	Relax a local mesh.

std::vector< Tuple >	local_mesh_tuples (const VectorNd &center) const
	Get the local n-ring around a vertex.

std::vector< Tuple >	local_mesh_tuples (const Tuple &v) const
	Get the local n-ring around a vertex.

double	local_mesh_energy (const VectorNd &local_mesh_center) const
	Compute the energy of a local n-ring around a vertex.

double	local_energy_before () const
	Get the energy of the local n-ring around a vertex.

double	edge_elastic_energy (const Tuple &e) const
	Compute the average elastic energy of the faces containing an edge.

void	write_priority_queue_mesh (const std::string &path, const Tuple &e) const
	Write a visualization mesh of the priority queue.

Protected Member Functions inherited from polyfem::mesh::WildRemesher< WMTKMesh >
void	init_attributes_and_connectivity (const size_t num_vertices, const Eigen::MatrixXi &elements) override
	Create an internal mesh representation and associate attributes.

virtual bool	swap_edge_before (const Tuple &t) override

virtual bool	swap_edge_after (const Tuple &t) override

bool	invariants (const std::vector< Tuple > &new_tris) override
	Check if invariants are satisfied.

bool	is_rest_inverted (const Tuple &loc) const
	Check if a triangle's rest shape is inverted.

bool	is_inverted (const Tuple &loc) const
	Check if a triangle's rest and deformed shapes are inverted.

void	cache_split_edge (const Tuple &e)
	Cache the split edge operation.

void	cache_collapse_edge (const Tuple &e, const CollapseEdgeTo collapse_to)
	Cache the edge collapse operation.

void	cache_swap_edge (const Tuple &e)
	Cache the edge swap operation.

void	map_edge_split_edge_attributes (const Tuple &t)

void	map_edge_split_boundary_attributes (const Tuple &t)

void	map_edge_split_element_attributes (const Tuple &t)

void	map_edge_collapse_vertex_attributes (const Tuple &t)

void	map_edge_collapse_boundary_attributes (const Tuple &t)

void	map_edge_collapse_edge_attributes (const Tuple &t)

void	map_edge_swap_edge_attributes (const Tuple &t)

void	map_edge_swap_element_attributes (const Tuple &t)

Protected Member Functions inherited from polyfem::mesh::Remesher
void	project_quantities ()
	Update the mesh positions and other projection quantities.

void	cache_before ()
	Cache quantities before applying an operation.

Protected Attributes inherited from polyfem::mesh::PhysicsRemesher< wmtk::TetMesh >
wmtk::ExecutePass< WildRemesher, EXECUTION_POLICY >	executor

Protected Attributes inherited from polyfem::mesh::WildRemesher< WMTKMesh >
wmtk::ExecutePass< WildRemesher, EXECUTION_POLICY >	executor

int	m_n_quantities

double	total_volume

std::conditional< std::is_same< WMTKMesh, wmtk::TriMesh >::value, std::shared_ptr< TriOperationCache >, std::shared_ptr< TetOperationCache > >::type	op_cache

Protected Attributes inherited from polyfem::mesh::Remesher
GlobalProjectionCache	global_projection_cache

const json	args
	Copy of remesh args.

const Eigen::MatrixXd	m_obstacle_displacements
	Collision obstacles' displacements.

Eigen::MatrixXd	m_obstacle_quantities
	Collision obstacles' extra quantities.

const double	current_time
	Current time.

const double	starting_energy
	Starting energy.

Detailed Description

Definition at line 140 of file PhysicsRemesher.hpp.

Member Typedef Documentation

◆ Super

using polyfem::mesh::PhysicsTetRemesher::Super = PhysicsRemesher<wmtk::TetMesh>

private

Definition at line 143 of file PhysicsRemesher.hpp.

◆ Tuple

using polyfem::mesh::PhysicsTetRemesher::Tuple = typename Super::Tuple

Definition at line 146 of file PhysicsRemesher.hpp.

Constructor & Destructor Documentation

◆ PhysicsTetRemesher()

polyfem::mesh::PhysicsTetRemesher::PhysicsTetRemesher	(	const State &	state,
		const Eigen::MatrixXd &	obstacle_displacements,
		const Eigen::MatrixXd &	obstacle_vals,
		const double	current_time,
		const double	starting_energy
	)

inline

Definition at line 148 of file PhysicsRemesher.hpp.

The documentation for this class was generated from the following file:

src/polyfem/mesh/remesh/PhysicsRemesher.hpp

Public Types

Public Member Functions