705 {
709
710 auto N_InLoop = getNinTheLoop();
711 auto sensee = getSkeletonSense();
712 auto nb_gauss_pts = getGaussPts().size2();
713 auto t_normal = getFTensor1NormalsAtGaussPts();
714
715 auto t_sigma_ptr =
dataAtPts->getFTensorApproxP(getGaussPts().size2());
716 auto get_tracion =
719 if (N_InLoop == 0) {
721 }
722
723 auto t_traction = get_tracion();
724 for (int gg = 0; gg != nb_gauss_pts; gg++) {
725 t_traction(
i) = t_sigma_ptr(
i,
j) * sensee * (t_normal(
j) / t_normal.l2());
726 ++t_traction;
727 ++t_sigma_ptr;
728 ++t_normal;
729 }
730
732}
#define FTENSOR_INDEX(DIM, I)
#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