945 {
948
949 auto n_in_the_loop = getNinTheLoop();
950 auto loop_size = getLoopSize();
951 auto sense = getSkeletonSense();
952 auto nb_gauss_pts = getGaussPts().size2();
953 auto t_normal = getFTensor1NormalsAtGaussPts();
954
955 auto t_sigma =
dataAtPts->getFTensorApproxP(getGaussPts().size2());
956 auto get_tracion =
959 if (!n_in_the_loop) {
961 }
962
963 auto t_traction = get_tracion();
964 for (int gg = 0; gg != nb_gauss_pts; gg++) {
966 t_sigma(
i,
j) * sense * (t_normal(
j) / t_normal.l2()) / loop_size;
967 ++t_traction;
968 ++t_sigma;
969 ++t_normal;
970 }
971
973}
#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