756 const auto nb_gauss_pts = BoundaryEleOp::getGaussPts().size2();
763 sdf_vec.resize(nb_gauss_pts,
false);
764 grad_mat.resize(DIM, nb_gauss_pts,
false);
765 hess_mat.resize((DIM * (DIM + 1)) / 2, nb_gauss_pts,
false);
766 constraint_vec.resize(nb_gauss_pts,
false);
768 auto t_traction = getFTensor1FromMat<DIM>(contactTraction_mat);
771 auto t_grad_sdf = getFTensor1FromMat<DIM>(grad_mat);
772 auto t_hess_sdf = getFTensor2SymmetricFromMat<DIM>(hess_mat);
775 auto t_disp = getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp);
776 auto t_coords = BoundaryEleOp::getFTensor1CoordsAtGaussPts();
777 auto t_normal_at_pts = BoundaryEleOp::getFTensor1NormalsAtGaussPts();
782 auto ts_time = BoundaryEleOp::getTStime();
783 auto ts_time_step = BoundaryEleOp::getTStimeStep();
786 BoundaryEleOp::getFTensor1CoordsAtGaussPts(),
787 getFTensor1FromMat<DIM>(
commonDataPtr->contactDisp), nb_gauss_pts);
789 BoundaryEleOp::getFTensor1NormalsAtGaussPts(), nb_gauss_pts);
796 m_spatial_coords, m_normals_at_pts, block_id);
800 m_spatial_coords, m_normals_at_pts, block_id);
804 m_spatial_coords, m_normals_at_pts, block_id);
807 auto t_grad_sdf_v = getFTensor1FromMat<3>(m_grad_sdf);
808 auto t_hess_sdf_v = getFTensor2SymmetricFromMat<3>(m_hess_sdf);
822 for (
auto gg = 0; gg != nb_gauss_pts; ++gg) {
824 auto tn = -t_traction(
i) * t_grad_sdf_v(
i);
828 t_grad_sdf(
i) = t_grad_sdf_v(
i);
829 t_hess_sdf(
i,
j) = t_hess_sdf_v(
i,
j);