#pragma once
#include
namespace internal
{
// local: cylinder face x^2+y^2-1=0
struct cylinder_face_t final : subface {
static constexpr uint64_t max_degree = 2;
static constexpr equation_system_type eq_sys_type = equation_system_type::implicit;
std::function fetch_sdf_evaluator() const override;
std::function fetch_sdf_grad_evaluator() const override;
// u: planar angle from x-axis to z-axis
// v: depth/height from xz-plane to y-axis
std::function fetch_point_by_param_evaluator() const override;
std::function fetch_curve_constraint_evaluator(parameter_u_t constraint_var_type,
double u) const override;
std::function fetch_curve_constraint_evaluator(parameter_v_t constraint_var_type,
double v) const override;
std::function fetch_solver_evaluator() const override;
};
struct cylinder_paired_model_matrix {
internal::paired_model_matrix *data{};
};
} // namespace internal
namespace detail
{
template <>
struct hasher {
size_t operator()(const internal::cylinder_paired_model_matrix &block) const
{
Eigen::Matrix character_rows = block.data->world_to_local.matrix().topRows<2>();
return XXH3_64bits(character_rows.data(), sizeof(decltype(character_rows)));
}
};
template <>
struct default_elem_ctor {
internal::cylinder_face_t operator()(const internal::cylinder_paired_model_matrix &k) const
{
internal::cylinder_face_t res{};
res.model_matrices = const_cast(k.data);
return res;
}
};
} // namespace detail