path: root/src/helper/geom-pointwise.cpp

/**
 * \file
 * \brief Pointwise a class to manage a vector of satellites per piecewise curve
 */ /*
    * Authors:
    * 2015 Jabier Arraiza Cenoz<jabier.arraiza@marker.es>
    *
    * This code is in public domain
    */

#include <helper/geom-pointwise.h>

namespace Geom {

/**
 * @brief Pointwise a class to manage a vector of satellites per piecewise curve
 *
 * For the moment is a per curve satellite holder not per node. This is ok for
 * much cases but not a real node satellite on open paths
 * To implement this we can:
 * add extra satellite in open paths, and take notice of current open paths
 * or put extra satellites on back for each open subpath
 *
 * Also maybe the vector of satellites become a vector of
 * optional satellites, and remove the active variable in satellites.
 *
 */
Pointwise::Pointwise(Piecewise<D2<SBasis> > pwd2,
                     std::vector<Satellite> satellites)
    : _pwd2(pwd2), _satellites(satellites), _path_info(pwd2)
{
    setStart();
}
;

Pointwise::~Pointwise() {}
;

Piecewise<D2<SBasis> > Pointwise::getPwd2() const
{
    return _pwd2;
}

void Pointwise::setPwd2(Piecewise<D2<SBasis> > pwd2_in)
{
    _pwd2 = pwd2_in;
    _path_info.setPwd2(_pwd2);
}

std::vector<Satellite> Pointwise::getSatellites() const
{
    return _satellites;
}

void Pointwise::setSatellites(std::vector<Satellite> sats)
{
    _satellites = sats;
    setStart();
}

/** Update the start satellite on ope/closed paths.
 */
void Pointwise::setStart()
{
    std::vector<std::pair<size_t, bool> > path_info = _path_info.data;
    for (size_t i = 0; i < path_info.size(); i++) {
        size_t firstNode = _path_info.first(path_info[i].first);
        size_t lastNode = _path_info.last(path_info[i].first);
        if (!_path_info.closed(lastNode)) {
            _satellites[firstNode].hidden = true;
            _satellites[firstNode].active = false;
        } else {
            _satellites[firstNode].active = true;
            _satellites[firstNode].hidden = _satellites[firstNode + 1].hidden;
        }
    }
}

/** Fired when a path is modified.
 */
void Pointwise::recalculateForNewPwd2(Piecewise<D2<SBasis> > A, Geom::PathVector B)
{
    if (_pwd2.size() > A.size()) {
        pwd2Sustract(A);
    } else if (_pwd2.size() < A.size()) {
        pwd2Append(A);
    } else {
        insertDegenerateSatellites(A,B);
    }
}

/** Some nodes/subpaths are removed.
 */
void Pointwise::pwd2Sustract(Piecewise<D2<SBasis> > A)
{
    size_t counter = 0;
    std::vector<Satellite> sats;
    Piecewise<D2<SBasis> > pwd2 = _pwd2;
    setPwd2(A);
    for (size_t i = 0; i < _satellites.size(); i++) {
        if (_path_info.last(i - counter) < i - counter ||
                !are_near(pwd2[i].at0(), A[i - counter].at0())) {
            counter++;
        } else {
            sats.push_back(_satellites[i - counter]);
        }
    }
    setSatellites(sats);
}

/** Append nodes/subpaths to current pointwise
 */
void Pointwise::pwd2Append(Piecewise<D2<SBasis> > A)
{
    size_t counter = 0;
    std::vector<Satellite> sats;
    bool reversed = false;
    bool reorder = false;
    for (size_t i = 0; i < A.size(); i++) {
        size_t first = _path_info.first(i - counter);
        size_t last = _path_info.last(i - counter);
        //Check for subpath closed. If a subpath is closed, is not reversed or moved
        //to back
        _path_info.setPwd2(A);
        size_t new_subpath_index = _path_info.subPathIndex(i);
        _path_info.setPwd2(_pwd2);
        bool subpath_is_changed = false;
        if (_pwd2.size() <= i - counter) {
            subpath_is_changed = false;
        } else {
            subpath_is_changed = new_subpath_index != _path_info.subPathIndex(i - counter);
        }
        if (!reorder && first == i - counter && !are_near(_pwd2[i - counter].at0(), A[i].at0()) && !subpath_is_changed) {
            //Send the modified subpath to back
            subpathToBack(_path_info.subPathIndex(first));
            reorder = true;
            i--;
            continue;
        }
        if (!reversed && first == i - counter && !are_near(_pwd2[i - counter].at0(), A[i].at0()) && !subpath_is_changed) {
            subpathReverse(first, last);
            reversed = true;
        }

        if (_pwd2.size() <= i - counter || !are_near(_pwd2[i - counter].at0(), A[i].at0())){
            counter++;
            bool active = true;
            bool hidden = false;
            bool is_time = _satellites[0].isTime;
            bool mirror_knots = _satellites[0].hasMirror;
            double amount = 0.0;
            double degrees = 0.0;
            int steps = 0;
            Satellite sat(_satellites[0].satelliteType, is_time, active, mirror_knots,
                          hidden, amount, degrees, steps);
            sats.push_back(sat);
        } else {
            sats.push_back(_satellites[i - counter]);
        }
    }
    setPwd2(A);
    setSatellites(sats);
}

void Pointwise::subpathToBack(size_t subpath)
{
    std::vector<Geom::Path> path_in =
        path_from_piecewise(remove_short_cuts(_pwd2, 0.1), 0.001);
    size_t subpath_counter = 0;
    size_t counter = 0;
    std::vector<Geom::Path> tmp_path;
    Geom::Path to_back;
    for (PathVector::const_iterator path_it = path_in.begin();
            path_it != path_in.end(); ++path_it) {
        if (path_it->empty()) {
            continue;
        }
        Geom::Path::const_iterator curve_it1 = path_it->begin();
        Geom::Path::const_iterator curve_endit = path_it->end_default();
        const Curve &closingline = path_it->back_closed();
        if (are_near(closingline.initialPoint(), closingline.finalPoint())) {
            curve_endit = path_it->end_open();
        }
        while (curve_it1 != curve_endit) {
            if (subpath_counter == subpath) {
                _satellites.push_back(_satellites[counter]);
                _satellites.erase(_satellites.begin() + counter);
            } else {
                counter++;
            }
            ++curve_it1;
        }
        if (subpath_counter == subpath) {
            to_back = *path_it;
        } else {
            tmp_path.push_back(*path_it);
        }
        subpath_counter++;
    }
    tmp_path.push_back(to_back);
    setPwd2(remove_short_cuts(paths_to_pw(tmp_path), 0.01));
}

void Pointwise::subpathReverse(size_t start, size_t end)
{
    start++;
    for (size_t i = end; i >= start; i--) {
        _satellites.push_back(_satellites[i]);
        _satellites.erase(_satellites.begin() + i);
    }
    std::vector<Geom::Path> path_in =
        path_from_piecewise(remove_short_cuts(_pwd2, 0.1), 0.001);
    size_t counter = 0;
    size_t subpath_counter = 0;
    size_t subpath = _path_info.subPathIndex(start);
    std::vector<Geom::Path> tmp_path;
    Geom::Path rev;
    for (PathVector::const_iterator path_it = path_in.begin();
            path_it != path_in.end(); ++path_it) {
        if (path_it->empty()) {
            continue;
        }
        counter++;
        if (subpath_counter == subpath) {
            tmp_path.push_back(path_it->reverse());
        } else {
            tmp_path.push_back(*path_it);
        }
        subpath_counter++;
    }
    setPwd2(remove_short_cuts(paths_to_pw(tmp_path), 0.01));
}


/** Fired when a path is modified duplicating a node. Piecewise ignore degenerated curves.
 */
void Pointwise::insertDegenerateSatellites(Piecewise<D2<SBasis> > A, Geom::PathVector B)
{
    size_t size_A = A.size();
    _path_info.setPathVector(B);
    size_t size_B = _path_info.numberCurves();
    size_t satellite_gap = size_B - size_A;
    if (satellite_gap == 0){
        return;
    }
    size_t counter = 0;
    size_t counter_added = 0;
    for (PathVector::const_iterator path_it = B.begin();
            path_it != B.end(); ++path_it) {
        if (path_it->empty()) {
            continue;
        }
        Geom::Path::const_iterator curve_it1 = path_it->begin();
        Geom::Path::const_iterator curve_endit = path_it->end_default();
        if (path_it->closed()) {
            const Curve &closingline = path_it->back_closed();
            if (are_near(closingline.initialPoint(), closingline.finalPoint())) {
                curve_endit = path_it->end_open();
            }
        }
        while (curve_it1 != curve_endit) {
            if ((*curve_it1).isDegenerate() && counter_added < satellite_gap){
                counter_added++;
                bool active = true;
                bool hidden = false;
                bool is_time = _satellites[0].isTime;
                bool mirror_knots = _satellites[0].hasMirror;
                double amount = 0.0;
                double degrees = 0.0;
                int steps = 0;
                Satellite sat(_satellites[0].satelliteType, is_time, active, mirror_knots,
                              hidden, amount, degrees, steps);
                _satellites.insert(_satellites.begin() + counter + 1 ,sat);
            }
            ++curve_it1;
            counter++;
        }
    }

    _path_info.setPwd2(A);
    setPwd2(A);
}

} // namespace Geom
/*
  Local Variables:
  mode:c++
  c-file-style:"stroustrup"
  c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +))
  indent-tabs-mode:nil
  fill-column:99
  End:
*/
// vim:
// filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:fileencoding=utf-8:textwidth=99
// :