BoundaryConstraint.cpp

#include "BoundaryConstraint.h"
#include "DartConstraint.h"
#include "Mesh.h"

bool BoundaryConstraint::is_uniform(Mesh const &m) const {
    DartConstraint const &dc = m.dart_constraint(DartConstraints[0]);
    double_t delta_0{dc.delta()};

    for(size_t i : DartConstraints) {
        double delta_i{m.dart_constraint(i).delta()};
        if (std::abs(delta_0 - delta_i) > FLT_EPSILON) {
            return false;
        }
    }

    return true;
}

double BoundaryConstraint::smallest_parametric_edge_length(Mesh &m) const {
    double_t delta_min{DBL_MAX};

    for (size_t i : DartConstraints) {
        DartConstraint const &dc = m.dart_constraint(i);
        delta_min = std::min(delta_min, std::abs(dc.S1 - dc.S0));
    }

    return delta_min;
}

double BoundaryConstraint::queue_uniform(Mesh &m, double delta_min) const {
    double delta_max{delta_min};

    for(size_t i : DartConstraints) {
        DartConstraint const &dc = m.dart_constraint(i);
        double_t delta = std::abs(dc.S1 - dc.S0);
        if (delta > delta_min) {
            m.add_to_queue(std::make_unique<MidpointQueuer>(dc.dart()));

            delta_max = std::max(delta_max, delta);
        }
    }

    return delta_max;
}

void BoundaryConstraint::add_to_queue(Mesh &m, size_t dart) const {
    if (UniformDiscretization) {
        queue_uniform(m, smallest_parametric_edge_length(m) / 2.0);
    } else {
        m.add_to_queue(std::make_unique<MidpointQueuer>(dart));
    }
}


void BoundaryConstraint::make_uniform(Mesh &m) { make_uniform(m, smallest_parametric_edge_length(m)); };

void BoundaryConstraint::make_uniform(Mesh &m, double delta_min) {
    delta_min *= (1 + FLT_EPSILON);
    double_t delta_max{DBL_MAX};
    while (delta_max > delta_min) {
        delta_max = queue_uniform(m, delta_min);
        m.insert_from_queue();
    }
}

void BoundaryConstraint::refine(Mesh &m, double_t tol) {
    double_t err_max{DBL_MAX};

    while (err_max > tol) {
        err_max = 0.0;
        for (size_t i : DartConstraints) {
            DartConstraint const &dc = m.dart_constraint(i);

            double_t len_c = ConstraintCurve->length(dc.S0, dc.S1);
            double_t len_e = m.edge(dc.dart()).length(m);
            double_t err_i = std::abs(len_c - len_e) / len_c;

            if (err_i > tol) {
                m.add_to_queue(std::make_unique<MidpointQueuer>(dc.dart()));

                err_max = std::max(err_max, err_i);
            }
        }

        m.insert_from_queue();
    }
}

void BoundaryConstraint::sort_dart_constraints(Mesh const &m) {
    std::sort(DartConstraints.begin(), DartConstraints.end(),
              [&](size_t i, size_t j) { return m.dart_constraint(i).S0 < m.dart_constraint(j).S0; }
    );
}

std::vector<size_t> BoundaryConstraint::nodes(Mesh const &m) {
    sort_dart_constraints(m);
    std::vector<size_t> n;
    n.reserve(DartConstraints.size());
    for (size_t i : DartConstraints) {
        n.push_back(m.dart_constraint(i).base(m));
    }

    return n;
}