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