search_swept_face.h

Go to the documentation of this file.

/*
This file is part of the Ristra portage project.
Please see the license file at the root of this repository, or at:
    https://github.com/laristra/portage/blob/master/LICENSE
*/

#ifndef PORTAGE_SEARCH_SEARCH_SWEPT_FACE_H_
#define PORTAGE_SEARCH_SEARCH_SWEPT_FACE_H_

#include <vector>
#include <algorithm>

// portage includes
#include "portage/support/portage.h"

namespace Portage {

template <int D, Entity_kind on_what,
          typename SourceMeshType, typename TargetMeshType>
class SearchSweptFace {
public:
  SearchSweptFace() = delete;

  // Constructor with Meshes
  SearchSweptFace(const SourceMeshType & source_mesh,
                  const TargetMeshType & target_mesh)
    : sourceMesh_(source_mesh), targetMesh_(target_mesh)  {
#ifdef DEBUG
    int numSourceCells = sourceMesh_.num_owned_cells() +
      sourceMesh_.num_ghost_cells();
    int numTargetCells = targetMesh_.num_owned_cells() +
      targetMesh_.num_ghost_cells();
    assert(numSourceCells == numTargetCells);
#endif    
  }  // SearchSweptFace::SearchSweptFace

  SearchSweptFace(const SearchSweptFace&) = delete;

  SearchSweptFace & operator = (const SearchSweptFace &) = delete;

  ~SearchSweptFace() = default;

  std::vector<int> operator() (const int entityId) const {
    std::vector<int> candidates;
    std::cerr << "Swept face search is not implemented for a generic entity kind" << std::endl;
    return candidates;
  }

private:
  const SourceMeshType &sourceMesh_;
  const TargetMeshType &targetMesh_;
};  // class SearchSweptFace

template <int D, typename SourceMeshType, typename TargetMeshType>
class SearchSweptFace<D, Entity_kind::NODE, SourceMeshType, TargetMeshType> {
public:
  SearchSweptFace() = delete;

  // Constructor with Meshes
  SearchSweptFace(const SourceMeshType & source_mesh,
                  const TargetMeshType & target_mesh)
    : sourceMesh_(source_mesh), targetMesh_(target_mesh)  {
#ifdef DEBUG
    int numSourceCells = sourceMesh_.num_owned_cells() +
      sourceMesh_.num_ghost_cells();
    int numTargetCells = targetMesh_.num_owned_cells() +
      targetMesh_.num_ghost_cells();
    assert(numSourceCells == numTargetCells);
#endif    
  }  // SearchSweptFace::SearchSweptFace

  SearchSweptFace & operator = (const SearchSweptFace &) = delete;

  ~SearchSweptFace() = default;

  std::vector<int> operator() (const int entityId) const {
    std::vector<int> candidates;
    std::cerr << "Swept face search is not implemented for nodal control volumes" << std::endl;
    return candidates;
  }

private:
  const SourceMeshType &sourceMesh_;
  const TargetMeshType &targetMesh_;
};  // class SearchSweptFace

template <int D, typename SourceMeshType, typename TargetMeshType>
class SearchSweptFace<D, Entity_kind::CELL, SourceMeshType, TargetMeshType> {
public:
  SearchSweptFace() = delete;

  // Constructor with Meshes
  SearchSweptFace(const SourceMeshType & source_mesh,
                  const TargetMeshType & target_mesh)
    : sourceMesh_(source_mesh), targetMesh_(target_mesh)  {
#ifdef DEBUG
    int numSourceCells = sourceMesh_.num_owned_cells() +
      sourceMesh_.num_ghost_cells();
    int numTargetCells = targetMesh_.num_owned_cells() +
      targetMesh_.num_ghost_cells();
    assert(numSourceCells == numTargetCells);
#endif    
  }  // SearchSweptFace::SearchSweptFace

  SearchSweptFace & operator = (const SearchSweptFace &) = delete;

  ~SearchSweptFace() = default;

  std::vector<int> operator() (const int entityId) const {
    std::vector<int> candidates;
    sourceMesh_.cell_get_face_adj_cells(entityId, Entity_type::ALL, &candidates);
#ifdef DEBUG
    std::vector<int> alt_candidates;
    targetMesh_.cell_get_face_adj_cells(entityId, Entity_type::ALL, &alt_candidates);
    assert(candidates.size() == alt_candidates.size());
    int const num_candidates = candidates.size();
    for (int icc = 0; icc < num_candidates; icc++)
      assert(candidates[icc] == alt_candidates[icc]);
#endif  

    return candidates;
  }

private:
  const SourceMeshType &sourceMesh_;
  const TargetMeshType &targetMesh_;
};  // class SearchSweptFace

}  // namespace Portage

#endif  // PORTAGE_SEARCH_SEARCH_SWEPT_FACE_H_