1022 {
1025
1026 auto n_in_the_loop = getNinTheLoop();
1027 auto loop_size = getLoopSize();
1028 auto sense = getSkeletonSense();
1029 auto nb_gauss_pts = getGaussPts().size2();
1030 auto t_normal = getFTensor1NormalsAtGaussPts();
1031
1032 auto t_sigma =
dataAtPts->getFTensorApproxP(getGaussPts().size2());
1033 auto get_tracion =
1035 dataAtPts->tractionAtPts, nb_gauss_pts);
1036 if (!n_in_the_loop) {
1038 }
1039
1040 auto t_traction = get_tracion();
1041 for (int gg = 0; gg != nb_gauss_pts; gg++) {
1043 t_sigma(
i,
j) * sense * (t_normal(
j) / t_normal.l2()) / loop_size;
1044 ++t_traction;
1045 ++t_sigma;
1046 ++t_normal;
1047 }
1048
1050}
#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