CD_MFHelmholtzElectrostaticEBBCImplem.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_MFHelmholtzElectrostaticEBBCImplem_H
#define CD_MFHelmholtzElectrostaticEBBCImplem_H
 
// Std includes
#include <limits>
 
// Our includes
#include <CD_MFHelmholtzElectrostaticEBBC.H>
#include <CD_NamespaceHeader.H>
 
Real
MFHelmholtzElectrostaticEBBC::getElectrodePotential(const RealVect& a_pos) const
{
  CH_TIME("MFHelmholtzElectrostaticEBBC::getElectrodePotential(RealVect)");
 
  // Find closest electrode
  int  closestElectrode = 0;
  Real minDist          = std::numeric_limits<Real>::infinity();
 
  const std::vector<std::pair<Electrode, ElectrostaticEbBc::BcFunction>>& electrodeBCs = m_electrostaticBCs.getBcs();
 
  for (int i = 0; i < electrodeBCs.size(); i++) {
    const RefCountedPtr<BaseIF>& impFunc = electrodeBCs[i].first.getImplicitFunction();
 
    const Real curDist = std::abs(impFunc->value(a_pos));
 
    if (curDist < minDist) {
      closestElectrode = i;
      minDist          = curDist;
    }
  }
 
  // Return potential of closest electrode. Again, following the (perhaps odd) convention that the FieldSolver is "time dependent" but the operator
  // factory is not, we've passed the time in by reference to the functions in m_electrostaticBCs. Calling those functions simply ignore the time-argument
  // and uses FieldSolver::m_dt instead. So, we can just use a dummy dt here.
  constexpr Real dummyDt = 0.0;
 
  return electrodeBCs[closestElectrode].second(a_pos, dummyDt);
}
 
#include <CD_NamespaceFooter.H>
 
#endif