209 {
211
212 if (CN::Dimension(type) != 2)
214
215 int nb_gauss_pts = data.
getN().size1();
217
218 for (int gg = 0; gg < nb_gauss_pts; gg++) {
219 for (
int dd = 0;
dd < nb_dofs;
dd++) {
220
221 double w = getGaussPts()(2, gg);
222 const double n0 = getNormalsAtGaussPts(gg)[0];
223 const double n1 = getNormalsAtGaussPts(gg)[1];
224 const double n2 = getNormalsAtGaussPts(gg)[2];
225 if (getFEType() == MBTRI) {
227 }
228
232 w;
233 }
234 }
235
237}
#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
