MeshTests.C File Reference

Go to the source code of this file.

Functions
void	triRodPositionsRef (std::vector< Point > &positions, Real nrings, Real pitch, Point center)
bool	relativeEq (double x, double y)
bool	pointLess (const Point &a, const Point &b)
	TEST (MeshTests, triRodCoordinates)

Variables
const double	tol = 1e-8
const double	eps = 1e-9

Function Documentation

pointLess()

bool pointLess	(	const Point &	a,
		const Point &	b
	)

Definition at line 57 of file MeshTests.C.

Referenced by TEST().

{
  for (int i = 0; i < LIBMESH_DIM; i++)
  {
    if (relativeEq(a(i), b(i)))
      continue;
    return a(i) < b(i);
  }
  return false;
}

relativeEq()

bool relativeEq	(	double	x,
		double	y
	)

Definition at line 44 of file MeshTests.C.

Referenced by pointLess().

{
  if (std::abs(x) < eps && std::abs(y) < eps)
    return true;

  double diff = std::abs(x - y);
  x = std::abs(x);
  y = std::abs(y);
  double largest = (y > x) ? y : x;
  return diff <= largest * tol;
}

TEST()

TEST	(	MeshTests	,
		triRodCoordinates
	)

Definition at line 68 of file MeshTests.C.

{

  int nrings = 4;
  Real pitch = 1;
  Point center(0, 0, 0);

  std::vector<Point> positions;
  triRodPositionsRef(positions, nrings, pitch, center);
  std::sort(positions.begin(), positions.end(), pointLess);

  std::vector<Point> positions2;
  TriSubChannelMesh::rodPositions(positions2, nrings, pitch, center);
  std::sort(positions2.begin(), positions2.end(), pointLess);

  ASSERT_EQ(positions.size(), positions2.size());

  std::stringstream msg;
  for (size_t i = 0; i < positions.size(); i++)
    if (!positions[i].absolute_fuzzy_equals(positions2[i]))
      msg << "point " << i + 1 << " differs: " << positions[i] << " != " << positions2[i] << "\n";

  if (msg.str().size() > 0)
    FAIL() << msg.str();
}

triRodPositionsRef()

void triRodPositionsRef	(	std::vector< Point > &	positions,
		Real	nrings,
		Real	pitch,
		Point	center
	)

Definition at line 20 of file MeshTests.C.

Referenced by TEST().

{
  positions.clear();

  const Real start_r = (nrings - 1) * pitch;
  const Real theta0 = 2 * libMesh::pi / 3;
  const Real startx = start_r * std::cos(theta0);
  const Real starty = start_r * std::sin(theta0);
  for (int i = 0; i < nrings; i++)
  {
    int n_rods_in_row = nrings + i;
    const Real y = starty - i * pitch * std::sin(theta0);
    const Real x_row = startx + i * pitch * std::cos(theta0);
    for (int j = 0; j < n_rods_in_row; j++)
    {
      const Real x = x_row + j * pitch;
      positions.emplace_back(center(0) + x, center(1) + y);
      if (i < nrings - 1)
        positions.emplace_back(center(0) + x, center(1) - y);
    }
  }
}

Variable Documentation

eps

const double eps = 1e-9

Definition at line 17 of file MeshTests.C.

Referenced by relativeEq().

tol

const double tol = 1e-8

Definition at line 16 of file MeshTests.C.

Referenced by relativeEq().