fe_test.h File Reference

Go to the source code of this file.

Classes
class	FETest< order, family, elem_type >

Functions
Number	linear_test (const Point &p, const Parameters &, const std::string &, const std::string &)
Gradient	linear_test_grad (const Point &, const Parameters &, const std::string &, const std::string &)
Number	rational_test (const Point &p, const Parameters &, const std::string &, const std::string &)
Gradient	rational_test_grad (const Point &p, const Parameters &, const std::string &, const std::string &)

Variables
static const Real	rational_w = 0.75

Function Documentation

linear_test()

Number linear_test ( const Point &  p, const Parameters &  , const std::string &  , const std::string &    )

inline

Definition at line 28 of file fe_test.h.

{
  const Real & x = p(0);
  const Real & y = (LIBMESH_DIM > 1) ? p(1) : 0;
  const Real & z = (LIBMESH_DIM > 2) ? p(2) : 0;
 
  return x + 0.25*y + 0.0625*z;
}

References libMesh::Real.

Referenced by FETest< order, family, elem_type >::setUp().

linear_test_grad()

Gradient linear_test_grad ( const Point &  , const Parameters &  , const std::string &  , const std::string &    )

inline

Definition at line 41 of file fe_test.h.

{
  Gradient grad = 1;
  if (LIBMESH_DIM > 1)
    grad(1) = 0.25;
  if (LIBMESH_DIM > 2)
    grad(2) = 0.0625;
 
  return grad;
}

Referenced by FETest< order, family, elem_type >::setUp().

rational_test()

Number rational_test ( const Point &  p, const Parameters &  , const std::string &  , const std::string &    )

inline

Definition at line 62 of file fe_test.h.

{
  const Real & x = p(0);
  const Real & y = (LIBMESH_DIM > 1) ? p(1) : 0;
  const Real & z = (LIBMESH_DIM > 2) ? p(2) : 0;
 
  const Real denom = ((1-x)*(1-x)+x*x+2*rational_w*x*(1-x))*
                     ((1-y)*(1-y)+y*y+2*rational_w*y*(1-y))*
                     ((1-z)*(1-z)+z*z+2*rational_w*z*(1-z));
 
  return (x + 0.25*y + 0.0625*z)/denom;
}

References rational_w, and libMesh::Real.

Referenced by FETest< order, family, elem_type >::setUp(), and FETest< order, family, elem_type >::testU().

rational_test_grad()

Gradient rational_test_grad ( const Point &  p, const Parameters &  , const std::string &  , const std::string &    )

inline

Definition at line 79 of file fe_test.h.

{
  const Real & x = p(0);
  const Real & y = (LIBMESH_DIM > 1) ? p(1) : 0;
  const Real & z = (LIBMESH_DIM > 2) ? p(2) : 0;
 
  const Real xpoly = (1-x)*(1-x)+x*x+2*rational_w*x*(1-x);
  const Real xderiv = -2*(1-x)+2*x+2*rational_w*(1-2*x);
  const Real ypoly = (1-y)*(1-y)+y*y+2*rational_w*y*(1-y);
  const Real yderiv = -2*(1-y)+2*y+2*rational_w*(1-2*y);
  const Real zpoly = (1-z)*(1-z)+z*z+2*rational_w*z*(1-z);
  const Real zderiv = -2*(1-z)+2*z+2*rational_w*(1-2*z);
 
  const Real denom = xpoly * ypoly * zpoly;
 
  const Real numer = (x + 0.25*y + 0.0625*z);
 
  Gradient grad_n = 1, grad_d = xderiv * ypoly * zpoly;
  if (LIBMESH_DIM > 1)
    {
      grad_n(1) = 0.25;
      grad_d(1) = xpoly * yderiv * zpoly;
    }
  if (LIBMESH_DIM > 2)
    {
      grad_n(2) = 0.0625;
      grad_d(2) = xpoly * ypoly * zderiv;
    }
 
  Gradient grad = (grad_n - numer * grad_d / denom) / denom;
 
  return grad;
}

References rational_w, and libMesh::Real.

Referenced by FETest< order, family, elem_type >::setUp(), FETest< order, family, elem_type >::testGradU(), and FETest< order, family, elem_type >::testGradUComp().

Variable Documentation

rational_w

const Real rational_w = 0.75

static

Definition at line 59 of file fe_test.h.

Referenced by rational_test(), rational_test_grad(), and FETest< order, family, elem_type >::setUp().