CD_DataOpsImplem.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_DATAOPSIMPLEM_H
#define CD_DATAOPSIMPLEM_H
 
// Chombo includes
#include <CH_Timer.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
template <typename T>
void
DataOps::setValue(LevelData<MFInterfaceFAB<T>>& a_lhs, const T& a_value)
{
  CH_TIME("DataOps::setValue");
 
  const DisjointBoxLayout& dbl = a_lhs.disjointBoxLayout();
  const DataIterator&      dit = dbl.dataIterator();
 
  const int nbox = dit.size();
 
#pragma omp parallel for schedule(runtime)
  for (int mybox = 0; mybox < nbox; mybox++) {
    const DataIndex& din = dit[mybox];
 
    a_lhs[din].setVal(a_value);
  }
}
 
template <typename T>
int
DataOps::sgn(const T a_value)
{
  return ((a_value > 0.0) - (a_value < 0.0));
}
 
template <typename T>
void
DataOps::copy(EBAMRData<T>& a_dst, const EBAMRData<T>& a_src, const Interval a_dstInterv, const Interval a_srcInterv)
{
  CH_TIME("DataOps::copy(EBAMRData)");
 
  CH_assert(a_dst.size() == a_src.size());
  CH_assert(a_dstInterv.size() == a_srcInterv.size());
 
  for (int lvl = 0; lvl < a_dst.size(); lvl++) {
    CH_assert(a_dst[lvl]->nComp() > a_dstInterv.end());
    CH_assert(a_src[lvl]->nComp() > a_srcInterv.end());
 
    a_src[lvl]->localCopyTo(a_srcInterv, *a_dst[lvl], a_dstInterv);
  }
}
 
#include <CD_NamespaceFooter.H>
 
#endif