slim2d.hpp

#pragma once
 
#include <Eigen/Core>
 
#include <array>
#include <vector>
 
namespace wmtk::components::tracking {
 
struct SlimElement
{
    std::array<int, 3> v;
    Eigen::Matrix2d rest_inv = Eigen::Matrix2d::Identity();
    double weight = 0;
};
 
SlimElement slim_element(
    const std::array<int, 3>& v,
    const Eigen::Vector3d& p0,
    const Eigen::Vector3d& p1,
    const Eigen::Vector3d& p2,
    double weight);
 
SlimElement slim_element(
    const std::array<int, 3>& v,
    const Eigen::Vector2d& q0,
    const Eigen::Vector2d& q1,
    const Eigen::Vector2d& q2,
    double weight);
 
double slim_energy(const std::vector<Eigen::Vector2d>& x, const std::vector<SlimElement>& elements);
 
double slim(
    std::vector<Eigen::Vector2d>& x,
    const std::vector<char>& fixed,
    const std::vector<SlimElement>& elements,
    int iterations);
 
} // namespace wmtk::components::tracking