Example #1
0
void DC2DT6::BoundaryRHS(UserData<real_t>& ud)
{
   _ln->setLocal(-1.0);
   sRHS(1) += OFELI_THIRD*_ln->getDet()*ud.SurfaceForce((*_theSide)(1)->getCoord(),_theSide->getCode(1),_time);
   _ln->setLocal(0.0);
   sRHS(2) += OFELI_THIRD*_ln->getDet()*ud.SurfaceForce((*_theSide)(2)->getCoord(),_theSide->getCode(1),_time);
   _ln->setLocal(1.0);
   sRHS(3) += 4*OFELI_THIRD*_ln->getDet()*ud.SurfaceForce((*_theSide)(3)->getCoord(),_theSide->getCode(1),_time);
}