Function triangle_angle_based_quality

src/geode/geometry/quality.cpp:280–311

    template < index_t dimension >
    double triangle_angle_based_quality( const Triangle< dimension >& triangle )
    {
        try
        {
            const auto& vertices = triangle.vertices();
            std::array< double, 3 > sinus;
            for( const auto v : LRange{ 3 } )
            {
                const auto point = vertices[v].get();
                const auto point_prev = vertices[( v + 2 ) % 3].get();
                const auto point_next = vertices[( v + 1 ) % 3].get();
                const auto angle =
                    compute_angle( point, point_prev, point_next );
                DEBUG( angle );
                sinus[v] = std::sin( angle );
                DEBUG( sinus[v] );
            }
            const auto denominator = sinus[0] + sinus[1] + sinus[2];
            if( denominator == 0 )
            {
                return 0;
            }
            const auto quality =
                4 * sinus[0] * sinus[1] * sinus[2] / ( denominator );
            DEBUG( quality );
            return quality;
        }
        catch( ... )
        {
            return 0;
        }
    }

    template opengeode_geometry_api double triangle_angle_based_quality< 2 >(
        const Triangle2D& );