1206 {
1208
1211 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
1213
1214 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
1215 const auto nb_row_base_functions = row_data.getN().size2();
1216
1217 auto t_c_dstrain =
1218 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->resCdStrain);
1219 auto t_c_dplastic_strain =
1220 getFTensor2SymmetricFromMat<SPACE_DIM>(
commonDataPtr->resCdPlasticStrain);
1221
1222 auto next = [&]() {
1223 ++t_c_dstrain;
1224 ++t_c_dplastic_strain;
1225 };
1226
1228
1229 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
1230 auto t_row_base = row_data.getFTensor0N();
1231 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
1232 const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
1233 ++t_w;
1234
1236 t_res_vec(L) =
1237 t_L(
i,
j, L) * (t_c_dplastic_strain(
i,
j) - t_c_dstrain(
i,
j));
1238 next();
1239
1242 size_t rr = 0;
1243 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
1244 const auto row_base = alpha * t_row_base;
1245 auto t_col_base = col_data.getFTensor0N(gg, 0);
1246 for (
size_t cc = 0; cc != AssemblyDomainEleOp::nbCols /
size_symm; cc++) {
1247 t_mat(L) += (row_base * t_col_base) * t_res_vec(L);
1248 ++t_mat;
1249 ++t_col_base;
1250 }
1251 ++t_row_base;
1252 }
1253 for (; rr != nb_row_base_functions; ++rr)
1254 ++t_row_base;
1255 }
1256
1258}
#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
auto get_mat_scalar_dtensor_sym(MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)