MeshMath.hpp

Go to the documentation of this file.

#pragma once
 
#include "vc/core/types/ITKMesh.hpp"
 
namespace volcart::meshmath
{
 
template <
    typename T,
    std::enable_if_t<std::is_floating_point<T>::value, bool> = true>
T TriangleArea(T a, T b, T c)
{
    // Sort the side lengths so that a >= b >= c
    auto nc = std::min(a, std::min(b, c));
    auto na = std::max(a, std::max(b, c));
    auto nb = a + b + c - na - nc;
 
    // Calculate the area
    auto p = (na + (nb + nc)) * (nc - (na - nb)) * (nc + (na - nb)) *
             (na + (nb - nc));
 
    return 0.25 * std::sqrt(p);
}
 
double SurfaceArea(const ITKMesh::Pointer& mesh);
 
}  // namespace volcart::meshmath