CoarseIfaceInfo.h

Go to the documentation of this file.

/***************************************************************************
 *  ThunderEgg, a library for solvers on adaptively refined block-structured
 *  Cartesian grids.
 *
 *  Copyright (c) 2020-2021 Scott Aiton
 *
 *  This program 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.
 *
 *  This program 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 this program.  If not, see <https://www.gnu.org/licenses/>.
 ***************************************************************************/
 
#ifndef THUNDEREGG_SCHUR_COARSEIFACEINFO_H
#define THUNDEREGG_SCHUR_COARSEIFACEINFO_H
 
#include <ThunderEgg/PatchInfo.h>
#include <ThunderEgg/Schur/IfaceInfo.h>
namespace ThunderEgg {
namespace Schur {
template<int D>
class CoarseIfaceInfo : public IfaceInfo<D>
{
private:
  static int GetId(const PatchInfo<D>& pinfo, Side<D> s)
  {
    return (int)(pinfo.id * Side<D>::number_of + s.getIndex());
  }
 
public:
  Orthant<D - 1> orth_on_coarse;
  int coarse_rank;
  int coarse_id;
  int coarse_col_local_index = -1;
  int coarse_global_index = -1;
  CoarseIfaceInfo(const PatchInfo<D>& pinfo, Side<D> s)
    : IfaceInfo<D>(pinfo.rank, GetId(pinfo, s))
  {
    // fine and coarse interfaces always belong to their patches
    auto nbr_info = pinfo.getCoarseNbrInfo(s);
    orth_on_coarse = nbr_info.orth_on_coarse;
    coarse_rank = nbr_info.rank;
    coarse_id = nbr_info.id * Side<D>::number_of + s.opposite().getIndex();
  }
};
} // namespace Schur
} // namespace ThunderEgg
#endif