997 {
1000
1001 auto n_in_the_loop = getNinTheLoop();
1002 auto loop_size = getLoopSize();
1003 auto sense = getSkeletonSense();
1004 auto nb_gauss_pts = getGaussPts().size2();
1005 auto t_normal = getFTensor1NormalsAtGaussPts();
1006
1007 auto t_sigma =
dataAtPts->getFTensorApproxP(getGaussPts().size2());
1008 auto get_tracion =
1010 dataAtPts->tractionAtPts, nb_gauss_pts);
1011 if (!n_in_the_loop) {
1013 }
1014
1015 auto t_traction = get_tracion();
1016 for (int gg = 0; gg != nb_gauss_pts; gg++) {
1018 t_sigma(
i,
j) * sense * (t_normal(
j) / t_normal.l2()) / loop_size;
1019 ++t_traction;
1020 ++t_sigma;
1021 ++t_normal;
1022 }
1023
1025}
#define FTENSOR_INDEXES(DIM,...)
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
DataLayoutTraits< DataLayout::GaussByCoeffs > DL
decltype(GetFTensor1FromMatImpl< Tensor_Dim, S, DL, M >::get(std::declval< M & >(), 0, 0)) GetFTensor1FromMatType