induced_subcomplex.h

Go to the documentation of this file.

/*
 *  Copyright (C) 2014 Josh Bialkowski (jbialk@mit.edu)
 *
 *  This file is part of mpblocks.
 *
 *  mpblocks is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  mpblocks is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with mpblocks.  If not, see <http://www.gnu.org/licenses/>.
 */
#ifndef EDELSBRUNNER96_INDUCED_SUBCOMPLEX_H_
#define EDELSBRUNNER96_INDUCED_SUBCOMPLEX_H_

#include <array>
#include <set>

namespace edelsbrunner96 {


template <typename Traits>
class InducedSubcomplex {
 public:
  typedef typename Traits::Scalar     Scalar;
  typedef typename Traits::Storage    Storage;
  typedef typename Traits::PointRef   PointRef;
  typedef typename Traits::SimplexRef SimplexRef;
  typedef typename Traits::Simplex    Simplex;

  struct BreadthFirstNode {
    SimplexRef s_ref;
    PointRef   v_missing;

    BreadthFirstNode(SimplexRef s_ref, PointRef v_missing)
        : s_ref(s_ref),
          v_missing(v_missing) {
    }
  };

  std::array<PointRef, Traits::NDim + 2> V;

  std::array<SimplexRef, Traits::NDim + 2> S;

  std::array<PointRef, Traits::NDim> f;

  std::array<PointRef, 2> e;

  std::array<SimplexRef, 2>  s;

  void Init(Storage& storage, SimplexRef s0_ref, SimplexRef s1_ref);

  void Init(Storage& storage, const std::array<SimplexRef, 2>& s);

  bool IsLocallyRegular(Storage& storage);

  void Build(Storage& storage);


  bool IsFlippable(Storage& storage);

  uint8_t FindIndexOfPeakVertex(Simplex& s);

  void Flip(Storage& storage);
};

}  // namespace edelsbrunner96

#endif  // EDELSBRUNNER96_INDUCED_SUBCOMPLEX_H_