cogfaces.cxx

有限元学习研究用源代码(老外的),供科研人员参考

#include "cogfaces.hxx"

void CogeometryForBoundary::setDelta(cogFloat delta)
{
  Cogeometry::setDelta(delta);
  old->setDelta(delta);
}

cogIndex	CogeometryForBoundary::Line (cogFlag1& f, const cogLine& s) const
{
  wzIndex rc = Cogeometry::Line(f,s);
  if(rc != cogRCFaceFound) return rc;
  BoundaryCondition(f);
  return cogRCFaceFound;
}

cogIndex 	CogeometryForBoundary::BoundaryCondition(cogFlag1& flag) const
{
  face->Point(flag.p1);
  flag.p1.segment() = cogFace(flag.p1.segment().index());
  return cogRCConditionFound;
}

void CogeometryForFaces::setDelta(cogFloat delta)
{
  Cogeometry::setDelta(delta);
  old->setDelta(delta);
}

cogIndex CogeometryForFaces::Line (cogFlag1& flag, const cogLine& s) const
{
  cogIndex rc = old->Line(flag,s);
  face->Point(flag.p1);
  flag.p1.segment() = cogFace(flag.p1.segment().index());
  return rc;
}

void CogeometryWithBoundary::setDelta(cogFloat delta)
{
  Cogeometry::setDelta(delta);
  old->setDelta(delta);
}

cogIndex	CogeometryWithBoundary::Line (cogFlag1& f, const cogLine& s) const
{
  wzIndex rc = old->Line(f,s);
  if(rc != cogRCFaceFound) return rc;
  BoundaryCondition(f);
  return rc;
}

cogIndex 	CogeometryWithBoundary::BoundaryCondition(cogFlag1& flag) const
{
  flag.p1.segment() = f(flag.p1);
  return cogRCConditionFound;
}

void CogeometryWithFaces::setDelta(cogFloat delta)
{
  Cogeometry::setDelta(delta);
  old->setDelta(delta);
}

cogIndex CogeometryWithFaces::Line (cogFlag1& flag, const cogLine& s) const
{
  cogIndex rc = old->Line(flag,s);
  if(rc) flag.p1.segment()=f(flag.p1);
  return rc;
}