167 int faceInRespectToTet = getFEMethod()->nInTheLoop;
168 int nb_face_gauss_pts = getGaussPts().size2();
171 if (
type == MBVERTEX) {
174 for (
int fgg = 0; fgg < nb_face_gauss_pts; fgg++) {
182 0.5 * getNormalsAtGaussPtss();
185 getCoordsAtGaussPts();
190 getCoordsAtGaussPts();
195 }
catch (
const std::exception &ex) {
196 std::ostringstream ss;
197 ss <<
"throw in method: " << ex.what() << std::endl;
198 SETERRQ(PETSC_COMM_SELF, 1, ss.str().c_str());
202 for (
int fgg = 0; fgg < nb_face_gauss_pts; fgg++) {
206 CommonData::MultiIndexData gauss_pt_data(gg, side,
type);
207 std::pair<CommonData::Container::iterator, bool>
p;
211 "data not inserted");
213 CommonData::MultiIndexData &p_data =
214 const_cast<CommonData::MultiIndexData &
>(*
p.first);
216 int nb_shape_fun = data.getN().size2();
217 shape_fun.resize(nb_shape_fun);
220 cblas_dcopy(nb_shape_fun, &data.getN()(fgg, 0), 1, &shape_fun[0],
223 p_data.iNdices = data.getIndices();
224 p_data.dofOrders.resize(data.getFieldDofs().size(),
false);
225 for (
unsigned int dd = 0;
dd < data.getFieldDofs().size();
dd++) {
226 p_data.dofOrders[
dd] = data.getFieldDofs()[
dd]->getDofOrder();
230 }
catch (
const std::exception &ex) {
231 std::ostringstream ss;
232 ss <<
"throw in method: " << ex.what() << std::endl;
233 SETERRQ(PETSC_COMM_SELF, 1, ss.str().c_str());
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
@ MOFEM_DATA_INCONSISTENCY
#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)
UBlasVector< double > VectorDouble
std::vector< std::vector< int > > inTetFaceGaussPtsNumber
std::vector< MatrixDouble > faceNormals
std::vector< VectorDouble > rAy
std::vector< MatrixDouble > coordsAtGaussPts
std::vector< VectorDouble > cOords
std::vector< MatrixDouble > hoCoordsAtGaussPts
std::vector< MatrixDouble > faceGaussPts