214 {
216
217 if (CN::Dimension(type) != 2)
219
220 int nb_gauss_pts = data.
getN().size1();
222
223 for (int gg = 0; gg < nb_gauss_pts; gg++) {
224 for (
int dd = 0;
dd < nb_dofs;
dd++) {
225
226 double w = getGaussPts()(2, gg);
227 const double n0 = getNormalsAtGaussPts(gg)[0];
228 const double n1 = getNormalsAtGaussPts(gg)[1];
229 const double n2 = getNormalsAtGaussPts(gg)[2];
230 if (getFEType() == MBTRI) {
232 }
233
237 w;
238 }
239 }
240
242}
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
const Tensor2_symmetric_Expr< const ddTensor0< T, Dim, i, j >, typename promote< T, double >::V, Dim, i, j > dd(const Tensor0< T * > &a, const Index< i, Dim > index1, const Index< j, Dim > index2, const Tensor1< int, Dim > &d_ijk, const Tensor1< double, Dim > &d_xyz)
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
const VectorDouble & getFieldData() const
get dofs values
const MatrixAdaptor getVectorN(const FieldApproximationBase base, const int gg)
get Hdiv of base functions at Gauss pts
