olb::ConcreteBlockRefinementO< T, DESCRIPTOR, PLATFORM, OPERATOR > Struct Template Reference

#include <operator.h>

Collaboration diagram for olb::ConcreteBlockRefinementO< T, DESCRIPTOR, PLATFORM, OPERATOR >:

Public Member Functions
void	apply (ConcreteBlockRefinementContextD< T, DESCRIPTOR, PLATFORM > &data)

Detailed Description

template<typename T, typename DESCRIPTOR, Platform PLATFORM, typename OPERATOR>
requires (isPlatformCPU(PLATFORM))
struct olb::ConcreteBlockRefinementO< T, DESCRIPTOR, PLATFORM, OPERATOR >

Definition at line 35 of file operator.h.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , Platform PLATFORM, typename OPERATOR >

void olb::ConcreteBlockRefinementO< T, DESCRIPTOR, PLATFORM, OPERATOR >::apply ( ConcreteBlockRefinementContextD< T, DESCRIPTOR, PLATFORM > & data )

inline

Definition at line 36 of file operator.h.

                                                                           {
    const auto& cellIdCoarse = data.getConcreteData().template getField<fields::refinement::CELL_ID_COARSE>();
    const auto& cellIdFine = data.getConcreteData().template getField<fields::refinement::CELL_ID_FINE>();
    auto parameters = data.getConcreteData().template getData<OperatorParameters<OPERATOR>>().parameters;
 
           if constexpr (OPERATOR::scope == refinement::OperatorScope::PerCoarseCell) {
      #ifdef PARALLEL_MODE_OMP
      #pragma omp parallel for schedule(static)
      #endif
      for (CellID i=0; i < cellIdFine.getSize(); ++i) {
        if (cellIdCoarse[0][i] != 0) {
          cpu::Cell cCell{data.getCoarseLattice(), cellIdCoarse[0][i]};
          cpu::Cell fCell{data.getFineLattice(), cellIdFine[0][i]};
          cpu::Row  cData{data.getConcreteData(), i};
          OPERATOR().apply(cCell, fCell, cData, parameters);
        }
      }
    } else if constexpr (OPERATOR::scope == refinement::OperatorScope::PerFineCell) {
      if constexpr (OPERATOR::data::template contains<fields::refinement::CONTEXT_NEIGHBORS>()) {
        #ifdef PARALLEL_MODE_OMP
        #pragma omp parallel for schedule(static)
        #endif
        for (CellID i=0; i < cellIdFine.getSize(); ++i) {
          cpu::Cell fCell{data.getFineLattice(), cellIdFine[0][i]};
          refinement::CoarseCell cCell{cpu::Cell{data.getCoarseLattice(), cellIdCoarse[0][i]},
                                       data.getFineLattice().getLatticeR(cellIdFine[0][i])};
          refinement::ContextDataWithNeighbors cData{cpu::Row{data.getConcreteData(), i}};
          OPERATOR().apply(cCell, fCell, cData, parameters);
        }
      } else {
        #ifdef PARALLEL_MODE_OMP
        #pragma omp parallel for schedule(static)
        #endif
        for (CellID i=0; i < cellIdFine.getSize(); ++i) {
          cpu::Cell fCell{data.getFineLattice(), cellIdFine[0][i]};
          refinement::CoarseCell cCell{cpu::Cell{data.getCoarseLattice(), cellIdCoarse[0][i]},
                                       data.getFineLattice().getLatticeR(cellIdFine[0][i])};
          refinement::ContextData cData{cpu::Row{data.getConcreteData(), i}};
          OPERATOR().apply(cCell, fCell, cData, parameters);
        }
      }
    }
  }

References olb::ConcreteBlockRefinementContextD< T, DESCRIPTOR, PLATFORM >::getCoarseLattice(), olb::ConcreteBlockRefinementContextD< T, DESCRIPTOR, PLATFORM >::getConcreteData(), olb::ConcreteBlockRefinementContextD< T, DESCRIPTOR, PLATFORM >::getFineLattice(), OPERATOR, olb::refinement::PerCoarseCell, and olb::refinement::PerFineCell.

Here is the call graph for this function:

---

The documentation for this struct was generated from the following file:

• src/refinement/platform/cpu/operator.h