MaxwellBC.C

Go to the documentation of this file.

#include <iostream>
#include "Simulation.H"



namespace Maxwell {

  namespace bc {


    int ComputeDirichletBC(const std::vector<size_t> &bufferIndices, double *variable, const double &value)
    {

      std::vector<size_t>::const_iterator bufferIndicesIt = bufferIndices.begin();
      while(bufferIndicesIt != bufferIndices.end()) {
        size_t bufferIndex(*bufferIndicesIt++);
        variable[bufferIndex] = value;
      }

      return(0);

    }


    int ComputeSATDirichletBC(const std::vector<size_t> &bufferIndices, double *rhsVar, double *stateVar,
                              const double &value, const double &weight)
    {

      std::vector<size_t>::const_iterator bufferIndicesIt = bufferIndices.begin();
      while(bufferIndicesIt != bufferIndices.end()) {
        size_t bufferIndex(*bufferIndicesIt++);
        rhsVar[bufferIndex] = rhsVar[bufferIndex]-weight*(stateVar[bufferIndex]-value);
      }

      return(0);

    }


  }

}