1097 const size_t nb_gauss_pts = AssemblyBoundaryEleOp::getGaussPts().size2();
1098 auto &locMat = AssemblyBoundaryEleOp::locMat;
1100 auto t_normal_at_pts = AssemblyBoundaryEleOp::getFTensor1NormalsAtGaussPts();
1101 auto t_traction = getFTensor1FromMat<DIM>(
commonDataPtr->contactTraction);
1102 auto t_coords = AssemblyBoundaryEleOp::getFTensor1CoordsAtGaussPts();
1104 auto t_w = AssemblyBoundaryEleOp::getFTensor0IntegrationWeight();
1105 auto t_row_base = row_data.getFTensor1N<3>();
1106 size_t nb_face_functions = row_data.getN().size2() / 3;
1109 BoundaryEleOp::getFTensor1CoordsAtGaussPts(),
1110 getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp), nb_gauss_pts);
1112 BoundaryEleOp::getFTensor1NormalsAtGaussPts(), nb_gauss_pts);
1114 auto t_normal = getFTensor1FromMat<3>(m_normals_at_pts);
1116 auto ts_time = AssemblyBoundaryEleOp::getTStime();
1117 auto ts_time_step = AssemblyBoundaryEleOp::getTStimeStep();
1124 m_spatial_coords, m_normals_at_pts, block_id);
1128 m_spatial_coords, m_normals_at_pts, block_id);
1131 auto t_grad_sdf = getFTensor1FromMat<3>(m_grad_sdf);
1133 for (
size_t gg = 0; gg != nb_gauss_pts; ++gg) {
1135 double jacobian = 1.;
1137 jacobian = 2. * M_PI * t_coords(0);
1139 const double alpha = t_w * jacobian * AssemblyBoundaryEleOp::getMeasure();
1141 auto tn = -t_traction(
i) * t_grad_sdf(
i);
1145 t_cP(
i,
j) = (
c * t_grad_sdf(
i)) * t_grad_sdf(
j);
1153 for (; rr != AssemblyBoundaryEleOp::nbRows / DIM; ++rr) {
1155 auto t_mat = getFTensor2FromArray<DIM, DIM, DIM>(locMat, DIM * rr);
1156 const double row_base = t_row_base(
i) * t_normal(
i);
1158 auto t_col_base = col_data.getFTensor1N<3>(gg, 0);
1159 for (
size_t cc = 0; cc != AssemblyBoundaryEleOp::nbCols / DIM; ++cc) {
1160 const double col_base = t_col_base(
i) * t_normal(
i);
1161 const double beta = alpha * row_base * col_base;
1163 t_mat(
i,
j) -= beta * t_res_dt(
i,
j);
1171 for (; rr < nb_face_functions; ++rr)