18 double node_xi[4][2] = {{-1.0, -1.0}, {1.0, -1.0}, {1.0, 1.0}, {-1.0, 1.0}};
19 return 0.25 * (1.0 + node_xi[node_id][0] * xi_2d[0]) * (1.0 + node_xi[node_id][1] * xi_2d[1]);
25 std::vector<double> area_xi(3, 0.0);
26 area_xi[0] = xi_2d[0];
27 area_xi[1] = xi_2d[1];
28 area_xi[2] = 1.0 - xi_2d[0] - xi_2d[1];
29 return area_xi[node_id];
35 double node_xi[8][3] = {{-1.0, -1.0, -1.0},
43 return 0.125 * (1.0 + node_xi[node_id][0] * xi_3d[0]) * (1.0 + node_xi[node_id][1] * xi_3d[1]) *
44 (1.0 + node_xi[node_id][2] * xi_3d[2]);
50 std::vector<double> vol_xi(4, 0.0);
51 for (
unsigned int i = 0; i < 3; ++i)
53 vol_xi[3] = 1.0 - xi_3d[0] - xi_3d[1] - xi_3d[2];
54 return vol_xi[node_id];
double linearTetShape3D(unsigned int node_id, std::vector< double > &xi_3d)
double linearHexShape3D(unsigned int node_id, std::vector< double > &xi_3d)
double linearQuadShape2D(unsigned int node_id, std::vector< double > &xi_2d)
double linearTriShape2D(unsigned int node_id, std::vector< double > &xi_2d)