ShapeTet4-impl.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
namespace NumLib
{
template <class T_X, class T_N>
void ShapeTet4::computeShapeFunction(const T_X& r, T_N& N)
{
    N[0] = 1. - r[0] - r[1] - r[2];
    N[1] = r[0];
    N[2] = r[1];
    N[3] = r[2];
}
 
template <class T_X, class T_N>
void ShapeTet4::computeGradShapeFunction(const T_X& /*r*/, T_N& dN)
{
    // dr
    dN[0] = -1.0;
    dN[1] = 1.0;
    dN[2] = 0.0;
    dN[3] = 0.0;
 
    // ds
    dN[4] = -1.0;
    dN[5] = 0.0;
    dN[6] = 1.0;
    dN[7] = 0.0;
 
    // dt
    dN[8] = -1.0;
    dN[9] = 0.0;
    dN[10] = 0.0;
    dN[11] = 1.0;
}
 
}  // namespace NumLib