885 {
889
890 auto N_InLoop = getNinTheLoop();
891 auto sensee = getSkeletonSense();
892 auto nb_gauss_pts = getGaussPts().size2();
893 auto t_normal = getFTensor1NormalsAtGaussPts();
894
895 auto t_sigma_ptr =
897 if (N_InLoop == 0) {
900 }
902
903 for (int gg = 0; gg != nb_gauss_pts; gg++) {
904 t_traction(
i) = t_sigma_ptr(
i,
j) * sensee * (t_normal(
j) / t_normal.l2());
905 ++t_traction;
906 ++t_sigma_ptr;
907 ++t_normal;
908 }
909
911}
#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
FTensor::Tensor1< FTensor::PackPtr< T *, S >, Tensor_Dim > getFTensor1FromMat(ublas::matrix< T, L, A > &data)
Get tensor rank 1 (vector) form data matrix.
FTensor::Tensor2< FTensor::PackPtr< double *, 1 >, Tensor_Dim1, Tensor_Dim2 > getFTensor2FromMat(MatrixDouble &data)
Get tensor rank 2 (matrix) form data matrix.
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts