oomph::TElementShape< 3, 2 > Class Reference

Return local coordinates of node j. More...

#include <Telements.h>

Public Member Functions
void	local_coordinate_of_node (const unsigned &j, Vector< double > &s) const
void	shape (const Vector< double > &s, Shape &psi) const
	Shape function for specific TElement<3,2> More...
void	dshape_local (const Vector< double > &s, Shape &psi, DShape &dpsids) const
	Derivatives of shape functions for specific TElement<3,2> More...
void	d2shape_local (const Vector< double > &s, Shape &psi, DShape &dpsids, DShape &d2psids) const

Detailed Description

Return local coordinates of node j.

//////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////

Member Function Documentation

d2shape_local()

void oomph::TElementShape< 3, 2 >::d2shape_local ( const Vector< double > &  s, Shape &  psi, DShape &  dpsids, DShape &  d2psids  ) const

inline

Second derivatives of shape functions for specific TElement<3,2>: d2psids(i,0) = \( \partial^2 \psi_j / \partial s_0^2 \) d2psids(i,1) = \( \partial^2 \psi_j / \partial s_1^2 \) d2psids(i,2) = \( \partial^2 \psi_j / \partial s_2^2 \) d2psids(i,3) = \( \partial^2 \psi_j / \partial s_0 \partial s_1 \) d2psids(i,4) = \( \partial^2 \psi_j / \partial s_0 \partial s_2 \) d2psids(i,5) = \( \partial^2 \psi_j / \partial s_1 \partial s_2 \)

    {
      // ALH: Don't know why object qualifier is needed
      this->dshape_local(s, psi, dpsids);
 
      for (unsigned i = 0; i < 4; i++)
      {
        d2psids(i, 0) = 0.0;
        d2psids(i, 1) = 0.0;
        d2psids(i, 2) = 0.0;
        d2psids(i, 3) = 0.0;
        d2psids(i, 4) = 0.0;
        d2psids(i, 5) = 0.0;
      }
    }

References i, and s.

dshape_local()

void oomph::TElementShape< 3, 2 >::dshape_local ( const Vector< double > &  s, Shape &  psi, DShape &  dpsids  ) const

inline

Derivatives of shape functions for specific TElement<3,2>

    {
      // ALH: Don't know why object qualifier is needed
      this->shape(s, psi);
 
      // Derivatives
      dpsids(0, 0) = 1.0;
      dpsids(0, 1) = 0.0;
      dpsids(0, 2) = 0.0;
 
      dpsids(1, 0) = 0.0;
      dpsids(1, 1) = 1.0;
      dpsids(1, 2) = 0.0;
 
      dpsids(2, 0) = 0.0;
      dpsids(2, 1) = 0.0;
      dpsids(2, 2) = 1.0;
 
      dpsids(3, 0) = -1.0;
      dpsids(3, 1) = -1.0;
      dpsids(3, 2) = -1.0;
    }

References oomph::OneDimLagrange::shape().

local_coordinate_of_node()

void oomph::TElementShape< 3, 2 >::local_coordinate_of_node ( const unsigned &  j, Vector< double > &  s  ) const

inline

    {
      s.resize(3);
 
      switch (j)
      {
        case 0:
          s[0] = 1.0;
          s[1] = 0.0;
          s[2] = 0.0;
          break;
 
        case 1:
          s[0] = 0.0;
          s[1] = 1.0;
          s[2] = 0.0;
          break;
 
        case 2:
          s[0] = 0.0;
          s[1] = 0.0;
          s[2] = 1.0;
          break;
 
        case 3:
          s[0] = 0.0;
          s[1] = 0.0;
          s[2] = 0.0;
          break;
 
        default:
          std::ostringstream error_message;
          error_message
            << "Element only has four nodes; called with node number " << j
            << std::endl;
 
          throw OomphLibError(error_message.str(),
                              OOMPH_CURRENT_FUNCTION,
                              OOMPH_EXCEPTION_LOCATION);
      }
    }

References j, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and s.

shape()

void oomph::TElementShape< 3, 2 >::shape ( const Vector< double > &  s, Shape &  psi  ) const

inline

Shape function for specific TElement<3,2>

    {
      psi[0] = s[0];
      psi[1] = s[1];
      psi[2] = s[2];
      psi[3] = 1.0 - s[0] - s[1] - s[2];
    }

References s.

---

The documentation for this class was generated from the following file:

• Telements.h