13 template <
class T_X,
class T_N>
16 const double r = x[0];
17 const double s = x[1];
18 const double t = x[2];
20 N[0] = 0.125 * (1 -
r) * (1 - s) * (1 - t);
21 N[1] = 0.125 * (1 +
r) * (1 - s) * (1 - t);
22 N[2] = 0.125 * (1 +
r) * (1 + s) * (1 - t);
23 N[3] = 0.125 * (1 -
r) * (1 + s) * (1 - t);
27 template <
class T_X,
class T_N>
30 const double r = x[0];
31 const double s = x[1];
32 const double t = x[2];
34 dN[0] = -0.125 * (1.0 - s) * (1.0 - t);
36 dN[2] = 0.125 * (1.0 + s) * (1.0 - t);
40 dN[5] = -0.125 * (1.0 -
r) * (1.0 - t);
41 dN[6] = -0.125 * (1.0 +
r) * (1.0 - t);
46 dN[10] = -0.125 * (1.0 -
r) * (1.0 - s);
47 dN[11] = -0.125 * (1.0 +
r) * (1.0 - s);
48 dN[12] = -0.125 * (1.0 +
r) * (1.0 + s);
49 dN[13] = -0.125 * (1.0 -
r) * (1.0 + s);
static void computeGradShapeFunction(const T_X &x, T_N &dN)
static void computeShapeFunction(const T_X &x, T_N &N)