face_geom.hpp

Go to the documentation of this file.

// core — FaceGeom: the static (sub)face geometry CSR of the collocated AMR projection.
//
// Factored out of flow.hpp into its own header so the device face-geometry assembler
// (facegeom_assembly.hpp) and the flow driver (flow.hpp) can both name the type without a
// circular include (flow.hpp consumes the assembler, the assembler produces a FaceGeom).
//
// Per (sub)face (matching AmrPoisson::forEachFaceFull, 2:1 sub-faces): neighbour, axis, dir, α·area
// (divergence weight), raw area (advection flux), face-normal distance, openness α (gradient gate),
// and the two SOU upstream-probe leaves. Per cell: 1/V and a fluid flag. Plus, since ROADMAP
// B5, the seam-reconstruction descriptors and records of docs/amr_cf_convective.md §5.2 (built
// by seam_recon.hpp, read by the two advective kernels only; empty where there is no 2:1 seam).
//
// Requires a Kokkos build (View).
#ifndef PECLET_AMR_FACE_GEOM_HPP
#define PECLET_AMR_FACE_GEOM_HPP
 
#include "peclet/amr/common.hpp"
#include "peclet/core/common/types.hpp"
#include "peclet/core/common/view.hpp"
 
namespace peclet::amr {
 
struct FaceGeom {
  View<Index> start;       
  View<Index> nbr;         
  View<int> axis;          
  View<int> dir;           
  View<double> alphaArea;  
  View<double> rawArea;    
  View<double> dist;       
  View<double> alpha;      
  View<Index> upupI;  
  View<Index> upupJ;  
  View<double> invVol;  
  View<char> fluid;     
  Index n = 0;
  // ---- seam reconstruction of the ADVECTED value (docs/amr_cf_convective.md §5.2) -------------
  // Empty unless the seam builder ran (seam_recon.hpp; setSolid). `seam` is the per-slot
  // descriptor id, −1 = plain slot; the descriptor arrays are indexed by it and the records by a
  // descriptor's four record ids. Read ONLY by the two advective kernels.
  View<Index> seam;      
  View<Index> sampI;     
  View<Index> sampJ;     
  View<Index> uuRecI;    
  View<Index> uuRecJ;    
  View<double> d1I;      
  View<double> d1J;      
  View<Index> recStart;  
  View<double> recDist;  
  View<Index> recCell;   
  View<double> recW;     
};
 
struct SeamView {
  View<const Index> seam, sampI, sampJ, uuRecI, uuRecJ, recStart, recCell;
  View<const double> d1I, d1J, recDist, recW;
  bool on = false;
 
  template <class FldFn>
  KOKKOS_INLINE_FUNCTION double gather(Index r, const FldFn& fld) const {
    double v = 0.0;
    for (Index e = recStart(r); e < recStart(r + 1); ++e)
      v += recW(e) * fld(recCell(e));
    return v;
  }
};
 
inline SeamView makeSeamView(const FaceGeom& g, bool on) {
  SeamView sv;
  sv.on = on && g.seam.extent(0) > 0;
  sv.seam = g.seam;
  sv.sampI = g.sampI;
  sv.sampJ = g.sampJ;
  sv.uuRecI = g.uuRecI;
  sv.uuRecJ = g.uuRecJ;
  sv.recStart = g.recStart;
  sv.recCell = g.recCell;
  sv.d1I = g.d1I;
  sv.d1J = g.d1J;
  sv.recDist = g.recDist;
  sv.recW = g.recW;
  return sv;
}
 
}  // namespace peclet::amr
 
#endif  // PECLET_AMR_FACE_GEOM_HPP