SingleElement Namespace Reference

Helper namespace for the single element test case. More...

Functions
void	face_normal (const unsigned &face, const Vector< double > &x, Vector< double > &n)

Detailed Description

Helper namespace for the single element test case.

Function Documentation

face_normal()

void SingleElement::face_normal	(	const unsigned &	face,
		const Vector< double > &	x,
		Vector< double > &	n
	)

 {
  double N[3] = {0.0,0.0,0.0};
  switch(face)
   {
   case 0:
    N[0] = -1.0;
    N[1] = 0.0;
    N[2] = 0.0;
    break;
 
   case 1:
    N[0] = 0.0;
    N[1] = -1.0;
    N[2] = 0.0;
    break;
 
   case 2:
    N[0] = 0.0;
    N[1] = 0.0;
    N[2] = -1.0;
    break;
 
   case 3:
    Vector<double> s(3);
    //Calculate the local "bulk" (unrotated) coordinate on the face
    for(unsigned i=0;i<3;i++)
     {
      s[0] = -0.5 + 0.5*sqrt(4.0*x[0] - 3);
      s[1] = -0.5 + 0.5*sqrt(4.0*x[1] - 3);
      s[2] = -0.5 + 0.5*sqrt(4.0*x[2] - 3);
     }
      
    //Hence compute the outer normal
    N[0] = (1.0 + 2.0*s[1])*(1.0 + 2.0*s[2]);
    N[1] = (1.0 + 2.0*s[0])*(1.0 + 2.0*s[2]);
    N[2] = (1.0 + 2.0*s[0])*(1.0 + 2.0*s[1]);
    break;
   }
 
  double length = sqrt(N[0]*N[0] + N[1]*N[1] + N[2]*N[2]);
  n[0] = N[0]/length;
  n[1] = N[1]/length;
  n[2] = N[2]/length;
 }

References i, n, N, s, sqrt(), and plotDoE::x.

Referenced by TFaceTestProblem< ELEMENT >::TFaceTestProblem().