Beispiel #1
0
bool Line1D::is_coord_in_element(const CoordsT& coord, const NodesT& nodes)
{
  MappedCoordsT mapped_coord = mapped_coordinate(coord,nodes);
  if( (mapped_coord[KSI] >= -0.5) &&
      (mapped_coord[KSI] <= 0.5) )
  {
    return true;
  }
  else
  {
    return false;
  }
}
Beispiel #2
0
bool Line1D::is_coord_in_element(const CoordsT& coord, const NodesT& nodes)
{
  static const Real tolerance = 1e-12;
  MappedCoordsT mapped_coord = mapped_coordinate(coord,nodes);
  if( (mapped_coord[KSI] >= -0.5 - tolerance) &&
      (mapped_coord[KSI] <= 0.5 + tolerance) )
  {
    return true;
  }
  else
  {
    return false;
  }
}
Beispiel #3
0
bool Tetra3D::is_coord_in_element(const CoordsT& coord, const NodesT& nodes)
{
  MappedCoordsT mapped_coord = mapped_coordinate(coord,nodes);
  const Real tolerance = 10*math::Consts::eps();
  if((mapped_coord[KSI] >= -tolerance) &&
     (mapped_coord[ETA] >= -tolerance) &&
     (mapped_coord[ZTA] >= -tolerance) &&
     (mapped_coord.sum() <= 1.+tolerance))
  {
    return true;
  }
  else
  {
    return false;
  }
}